title: Global Defensive Alliances in the Lexicographic Product of Paths and Cycles

abstract: A set $S$ of vertices of graph $G$ is a \textit{defensive alliance} of $G$ if for every $v \in S$, it holds $|N[v] \cap S| \geq |N[v]-S|$. An alliance $S$ is called $global$ if it is also a dominating set. In this paper, we determine the exact values of the global defensive alliance number of lexicographic products of path and cycles.


title: Precision QCD at the LHC: from the structure of the proton to all-order resummations

abstract: Experiments at the LHC are collecting a wealth of data with an unprecedented level of precision. As a consequence, the theoretical error is now starting to lag behind the experimental one, and a ceaseless effort is required to reduce the theory uncertainty to match the precision of the data. At hadron colliders, QCD predictions are obtained by convoluting perturbative parton-level results with non-perturbative PDFs, whose precise determination is crucial to reach percent-accurate theoretical predictions. At the parton level, cross sections are obtained through a perturbative expansion in the strong coupling. In some cases, large terms appear at all orders and spoil the convergence of the series. The perturbative description is rescued by resumming the series to all orders, thereby making theory calculations accurate in regions where a fixed-order treatment is not sufficient. In this thesis I first present two global PDF sets where fixed-order calculations are supplemented by threshold and high-energy resummation, respectively. In the first case, it is found that including resummation into PDFs can compensate for the enhancement in the partonic cross sections, with implications for high-mass resonance searches. In the second case, resummation quantitatively improves the description of the HERA structure functions, thus providing evidence of the onset of a new dynamical regime of QCD in the HERA data. I then focus on Higgs production in gluon fusion. The effect of threshold resummation on the total cross section is found to significantly improve the convergence of the perturbative series and to provide a robust method for estimating missing higher order uncertainty. Finally, I present predictions for the Higgs transverse-momentum spectrum both in the inclusive case and within fiducial cuts, exploiting a novel approach where transverse-momentum resummation is performed in direct space.


title: Science with an ngVLA: Imaging Cold Gas to 1 kpc Scales in High-Redshift Galaxies with the ngVLA

abstract: The next generation Very Large Array (ngVLA) will revolutionize our understanding of the distant Universe via the detection of cold molecular gas in the first galaxies. Its impact on studies of galaxy characterization via detailed gas dynamics will provide crucial insight on dominant physical drivers for star-formation in high redshift galaxies, including the exchange of gas from scales of the circumgalactic medium down to resolved clouds on mass scales of $\sim10^{5}\,M_\odot$. In this study, we employ a series of high-resolution, cosmological, hydrodynamic zoom simulations from the MUFASA simulation suite and a CASA simulator to generate mock ngVLA observations of a $z\sim4.5$ gas rich star-forming galaxy. Using the DESPOTIC radiative transfer code that encompasses simultaneous thermal, chemical, and statistical equilibrium in calculating the molecular and atomic level transitions of CO from ALMA for comparison. We find that observations of CO(1-0) are especially important for tracing the systemic redshift of the galaxy and the total mass of the well-shielded molecular gas reservoir, while even CO(2-1) can predominantly trace denser gas regions distinct from CO(1-0). The factor of 100 times improvement in mapping speed for the ngVLA beyond the Jansky VLA and the proposed ALMA Band 1 will make these detailed, high-resolution imaging and kinematic studies of CO(1-0) routine at $z\sim2-5$.


title: A systematic investigation of classical causal inference strategies under mis-specification due to network interference

abstract: We systematically investigate issues due to mis-specification that arise in estimating causal effects when (treatment) interference is informed by a network available pre-intervention, i.e., in situations where the outcome of a unit may depend on the treatment assigned to other units. We develop theory for several forms of interference through the concept of exposure neighborhood, and develop the corresponding semi-parametric representation for potential outcomes as a function of the exposure neighborhood. Using this representation, we extend the definition of two popular classes of causal estimands, marginal and average causal effects, to the case of network interference. We characterize the bias and variance one incurs when combining classical randomization strategies (namely, Bernoulli, Completely Randomized, and Cluster Randomized designs) and estimators (namely, difference-in-means and Horvitz-Thompson) used to estimate average treatment effect and on the total treatment effect, under misspecification due to interference. We illustrate how difference-in-means estimators can have arbitrarily large bias when estimating average causal effects, depending on the form and strength of interference, which is unknown at design stage. Horvitz-Thompson (HT) estimators are unbiased when the correct weights are specified. Here, we derive the HT weights for unbiased estimation of different estimands, and illustrate how they depend on the design, the form of interference, which is unknown at design stage, and the estimand. More importantly, we show that HT estimators are in-admissible for a large class of randomization strategies, in the presence of interference. We develop new model-assisted and model-dependent strategies to improve HT estimators, and we develop new randomization strategies for estimating the average treatment effect and total treatment effect.


title: Merge: An Architecture for Interconnected Testbed Ecosystems

abstract: In the cybersecurity research community, there is no one-size-fits-all solution for merging large numbers of heterogeneous resources and experimentation capabilities from disparate specialized testbeds into integrated experiments. The current landscape for cyber-experimentation is diverse, encompassing many fields including critical infrastructure, enterprise IT, cyber-physical systems, cellular networks, automotive platforms, IoT and industrial control systems. Existing federated testbeds are constricted in design to predefined domains of applicability, lacking the systematic ability to integrate the burgeoning number of heterogeneous devices or tools that enable their effective use for experimentation. We have developed the Merge architecture to dynamically integrate disparate testbeds in a logically centralized way that allows researchers to effectively discover, and use the resources and capabilities provided the by evolving ecosystem of distributed testbeds for the development of rigorous and high-fidelity cybersecurity experiments.


title: Simulating the Collapse of a Thick Accretion Disk due to a Type I X-ray Burst from a Neutron Star

abstract: We use two-dimensional, general relativistic, viscous, radiation hydrodynamic simulations to study the impact of a Type I X-ray burst on a hot and geometrically thick accretion disk surrounding an unmagnetized, non-rotating neutron star. The disk is initially consistent with a system in its low/hard spectral state, and is subject to a burst which rises to a peak luminosity of $10^{38}$ erg s$^{-1}$ in $2.05$ s. At the peak of the burst, the temperature of the disk has dropped by more than three orders of magnitude and its scale height has gone down by more than one order of magnitude. The simulations show that these effects predominantly happen due to Compton cooling of the hot plasma, and clearly illustrate the potential cooling effects of bursts on accretion disk coronae. In addition, we demonstrate the presence of Poynting-Robertson drag, though it only enhances the mass accretion rate onto the neutron star by a factor of $\sim 3$-$4$ compared to a simulation with no burst. Simulations such as these are important for building a general understanding of the response of an accretion disk to an intense X-ray impulse, which, in turn, will be crucial for deciphering burst spectra. Detailed analysis of such spectra offers the potential to measure neutron star radii, and hence constrain the neutron star equation of state, but only if the contributions coming from the impacted disk and its associated corona can be understood.


title: A novel chiral spin texture: Antiferromagnetic Skyrmionium

abstract: Exotic spin textures viz. chiral domain wall, vortices, skyrmion, skyrmionium, etc. have recently emerged as active field of research because of their potential applications in high density data storage technology and logic gate computing. Magnetic skyrmionium is a skyrmion like soliton, which carries zero topological quantum number. Skyrmioniums are superior to conventional skyrmions in ferromagnets due to their negligible skyrmion hall effect and higher velocity. The physical properties of both skyrmion and skyrmionium have been investigated rigorously in ferromagnetic systems. Recent observations hint that such chiral spin structures in antiferromagnetic (AFM) systems are more promising in comparison to the ferromagnetic ones because of their robustness towards external perturbation, absence of Skyrmion hall effect, etc. However skyrmionium in AFM materials are not reported in literature so far. In this work, we demonstrate that skyrmionium can be created and stabilized in AFM materials by application of spin polarized current in an experimentally feasible geometry. We have further studied the dynamics of AFM skyrmionium by applying spin polarized current.


title: Too Many Hats

abstract: A puzzle about prisoners trying to identify the color of a hat on their head leads to a version where there are k more hats than prisoners. This generalized puzzle is related to the independence number of the arrangement graph A(m, n) and to Steiner systems and other designs. A natural conjecture is that perfect hat-guessing strategies exist in all cases, where "perfect" means that the success probability is 1/(k+1). This is true when k = 1, but we show that it is false when k = 2. Further, we present a strategy with success rate at least 1/O(k log k), independent of the number of prisoners.


title: Quantile Regression Under Memory Constraint

abstract: This paper studies the inference problem in quantile regression (QR) for a large sample size $n$ but under a limited memory constraint, where the memory can only store a small batch of data of size $m$. A natural method is the na\"ive divide-and-conquer approach, which splits data into batches of size $m$, computes the local QR estimator for each batch, and then aggregates the estimators via averaging. However, this method only works when $n=o(m^2)$ and is computationally expensive. This paper proposes a computationally efficient method, which only requires an initial QR estimator on a small batch of data and then successively refines the estimator via multiple rounds of aggregations. Theoretically, as long as $n$ grows polynomially in $m$, we establish the asymptotic normality for the obtained estimator and show that our estimator with only a few rounds of aggregations achieves the same efficiency as the QR estimator computed on all the data. Moreover, our result allows the case that the dimensionality $p$ goes to infinity. The proposed method can also be applied to address the QR problem under distributed computing environment (e.g., in a large-scale sensor network) or for real-time streaming data.


title: Toric generalized K$\ddot{a}$hler structures. I

abstract: This is a sequel of \cite{Wang}, which provides a general formalism for this paper. We mainly investigate thoroughly a subclass of toric generalized K$\ddot{a}$hler manifolds of symplectic type introduced by Boulanger in \cite{Bou}. We find torus actions on such manifolds are all \emph{strong Hamiltonian} in the sense of \cite{Wang}. For each such a manifold, we prove that besides the ordinary two complex structures $J_\pm$ associated to the biHermitian description, there is a \emph{third} canonical complex structure $J_0$ underlying the geometry, which makes the manifold toric K$\ddot{a}$hler. We find the other generalized complex structure besides the symplectic one is always a B-transform of a generalized complex structure induced from a $J_0$-holomorphic Poisson structure $\beta$ characterized by an anti-symmetric constant matrix. Stimulated by the above results, we introduce a \emph{generalized Delzant construction} which starts from a Delzant polytope with $d$ faces of codimension 1, the standard K$\ddot{a}$hler structure of $\mathbb{C}^d$ and an anti-symmetric $d\times d$ matrix. This construction is used to produce non-abelian examples of strong Hamiltonian actions.


title: A Robust Feature-aware Sparse Mesh Representation

abstract: The sparse representation of signals defined on Euclidean domains has been successfully applied in signal processing. Bringing the power of sparse representations to non-regular domains is still a challenge, but promising approaches have started emerging recently. In this paper, we investigate the problem of sparsely representing discrete surfaces and propose a new representation that is capable of providing tools for solving different geometry processing problems. The sparse discrete surface representation is obtained by combining innovative approaches into an integrated method. First, to deal with irregular mesh domains, we devised a new way to subdivide discrete meshes into a set of patches using a feature-aware seed sampling. Second, we achieve good surface approximation with over-fitting control by combining the power of a continuous global dictionary representation with a modified Orthogonal Marching Pursuit. The discrete surface approximation results produced were able to preserve the shape features while being robust to over-fitting. Our results show that the method is quite promising for applications like surface re-sampling and mesh compression.


title: Connectivity-Preserving Swarm Teleoperation With A Tree Network

abstract: During swarm teleoperation, the human operator may threaten the distance-dependent inter-robot communications and, with them, the connectivity of the slave swarm. To prevent the harmful component of the human command from disconnecting the swarm network, this paper develops a constructive strategy to dynamically modulate the interconnections of, and the locally injected damping at, all slave robots. By Lyapunov-based set invariance analysis, the explicit law for updating that control gains has been rigorously proven to synchronize the slave swarm while preserving all interaction links in the tree network. By properly limiting the impact of the user command rather than rejecting it entirely, the proposed control law enables the human operator to guide the motion of the slave swarm to the extent to which it does not endanger the connectivity of the swarm network. Experiment results demonstrate that the proposed strategy can maintain the connectivity of the tree network during swarm teleoperation.


title: Magnetohydrodynamical equilibria with current singularities and continuous rotational transform

abstract: We revisit the Hahm-Kulsrud-Taylor (HKT) problem, a classic prototype problem for studying resonant magnetic perturbations and 3D magnetohydrodynamical equilibria. We employ the boundary-layer techniques developed by Rosenbluth, Dagazian, and Rutherford (RDR) for the internal $m=1$ kink instability, while addressing the subtle difference in the matching procedure for the HKT problem. Pedagogically, the essence of RDR's approach becomes more transparent in the simplified slab geometry of the HKT problem. We then compare the boundary-layer solution, which yields a "DC" current singularity at the resonant surface, to the numerical solution obtained using a flux-preserving Grad-Shafranov solver. The remarkable agreement between the solutions demonstrates the validity and universality of RDR's approach. In addition, we show that RDR's approach consistently preserves the rotational transform, which hence stays continuous, contrary to a recent claim that RDR's solution contains a discontinuity in the rotational transform.


title: The NANOGrav 11-year Data Set: Pulse Profile Variability

abstract: Access to 50 years of data has led to the discovery of pulsar emission and rotation variability on timescales of months and years. Most of this long-term variability has been seen in long-period pulsars, with relatively little focus on recycled millisecond pulsars. We have analyzed a 38-pulsar sub-set of the 45 millisecond pulsars in the NANOGrav 11-year data set, in order to review their pulse profile stability. The most variability, on any timescale, is seen in PSRs J1713+0747, B1937+21 and J2145-0750. The strongest evidence for long-timescale pulse profile changes is seen in PSRs B1937+21 and J1643-1224. We have focused our analyses on these four pulsars in an attempt to elucidate the causes of their profile variability. Effects of scintillation seem to be responsible for the profile modifications of PSR J2145-0750. We see evidence that imperfect polarization calibration contributes to the profile variability of PSRs J1713+0747 and B1937+21, along with radio frequency interference around 2 GHz, but find that propagation effects also have an influence. The changes seen in PSR J1643-1224 have been reported previously, yet elude explanation beyond their astrophysical nature. Regardless of cause, unmodeled pulse profile changes are detrimental to the accuracy of pulsar timing and must be incorporated into the timing models where possible.


title: Lower bounds for fluctuations in first-passage percolation for general distributions

abstract: In first-passage percolation (FPP), one assigns i.i.d.~weights to the edges of the cubic lattice $\mathbb{Z}^d$ and analyzes the induced weighted graph metric. If $T(x,y)$ is the distance between vertices $x$ and $y$, then a primary question in the model is: what is the order of the fluctuations of $T(0,x)$? It is expected that the variance of $T(0,x)$ grows like the norm of $x$ to a power strictly less than 1, but the best lower bounds available are (only in two dimensions) of order $\log \|x\|$. This result was found in the '90s and there has not been any improvement since. In this paper, we address the problem of getting stronger fluctuation bounds: to show that $T(0,x)$ is with high probability not contained in an interval of size $o(\log \|x\|)^{1/2}$, and similar statements for FPP in thin cylinders. Such statements have been proved for special edge-weight distributions, and here we obtain such bounds for general edge-weight distributions. The methods involve inducing a fluctuation in the number of edges in a box whose weights are of "hi-mode" (large).


title: Improved general-purpose five-point model for water: TIP5P/2018

abstract: A new five point potential for liquid water, TIP5P/2018, is presented along with the techniques used to derive its charges from ab initio per-molecule electrostatic potentials in the liquid phase using the split charge equilibration (SQE) of Nistor et al. [J. Chem. Phys. 125, 094108 (2006)]. By taking the density and diffusion dependence on temperature as target properties, significant improvements to the behavior of isothermal compressibility were achieved along with improvements to other thermodynamic and rotational properties. While exhibiting a dipole moment close to ab initio values, TIP5P/2018 suffers from a too small quadrupole moment due to the charge assignment procedure and results in an overestimation of the dielectric constant.


title: BabyAI: A Platform to Study the Sample Efficiency of Grounded Language Learning

abstract: Allowing humans to interactively train artificial agents to understand language instructions is desirable for both practical and scientific reasons, but given the poor data efficiency of the current learning methods, this goal may require substantial research efforts. Here, we introduce the BabyAI research platform to support investigations towards including humans in the loop for grounded language learning. The BabyAI platform comprises an extensible suite of 19 levels of increasing difficulty. The levels gradually lead the agent towards acquiring a combinatorially rich synthetic language which is a proper subset of English. The platform also provides a heuristic expert agent for the purpose of simulating a human teacher. We report baseline results and estimate the amount of human involvement that would be required to train a neural network-based agent on some of the BabyAI levels. We put forward strong evidence that current deep learning methods are not yet sufficiently sample efficient when it comes to learning a language with compositional properties.


title: Spectra of ultrabroadband squeezed pulses and the finite-time Unruh-Davies effect

abstract: We study spectral properties of quantum radiation of ultimately short duration. In particular, we introduce a continuous multimode squeezing operator for the description of subcycle pulses of entangled photons generated by a coherent-field driving in a thin nonlinear crystal with second order susceptibility. We find the ultrabroadband spectra of the emitted quantum radiation perturbatively in the strength of the driving field. These spectra can be related to the spectra expected in an Unruh-Davies experiment with a finite time of acceleration. In the time domain, we describe the corresponding behavior of the normally ordered electric field variance.


title: 3C 17: The BCG of a galaxy cluster at z=0.22

abstract: Gemini Multi Object Spectrograph medium-resolution spectra and photometric data of 39 objects in the field of the radio galaxy 3C 17 are presented. Based on the new data, a previously uncataloged cluster of galaxies is identified at a mean redshift of z=0.220 +/-0.003, a projected virial radius of 0.37 Mpc, and a velocity dispersion of sigma_v = 821 +/-171 km/s. The brightest member of this cluster is 3C 17 with M_r = -22.45 mag. The surface brightness profile of 3C 17 is best fit with two components (Exponential + Sersic) characteristic of brightest cluster galaxies. The spectrum of 3C 17 is dominated by broad emission lines Halpha + N[II] and Hbeta + [OIII]. Analysis of Chandra data shows extended emission around the cluster core that supports the existence of hot gas cospatial with 3C 17. The discovery of a cluster of galaxies around 3C 17 better explains the sharply bent morphology of the radio jet given that it propagates through a dense intracluster medium.


title: Fractional supersymmetry algebra and lacunary Hermite polynomials

abstract: We consider a realization of fractional supersymmetric of quantum mechanics of order $r$, where the Hamiltonian and supercharges involve reflection operators. It is shown that the Hamiltonian has $r$-fold degenerate spectrum and the eigenvalues of hermitian supercharges are zeros of the associated Hermite polynomials of Askey and Wimp. Also it is shown that the associated eigenfunctions involve lacunary Hermite polynomials.


title: FPT algorithms to recognize well covered graphs

abstract: Given a graph $G$, let $vc(G)$ and $vc^+(G)$ be the sizes of a minimum and a maximum minimal vertex covers of $G$, respectively. We say that $G$ is well covered if $vc(G)=vc^+(G)$ (that is, all minimal vertex covers have the same size). Determining if a graph is well covered is a coNP-complete problem. In this paper, we obtain $O^*(2^{vc})$-time and $O^*(1.4656^{vc^+})$-time algorithms to decide well coveredness, improving results of Boria et. al. (2015). Moreover, using crown decomposition, we show that such problems admit kernels having linear number of vertices. In 2018, Alves et. al. (2018) proved that recognizing well covered graphs is coW[2]-hard when the independence number $\alpha(G)=n-vc(G)$ is the parameter. Contrasting with such coW[2]-hardness, we present an FPT algorithm to decide well coveredness when $\alpha(G)$ and the degeneracy of the input graph $G$ are aggregate parameters. Finally, we use the primeval decomposition technique to obtain a linear time algorithm for extended $P_4$-laden graphs and $(q,q-4)$-graphs, which is FPT parameterized by $q$, improving results of Klein et al (2013).


title: Quantum computing from a mathematical perspective: a description of the quantum circuit model

abstract: This paper is an essentially self-contained and rigorous description of the fundamental principles of quantum computing from a mathematical perspective. It is intended to help mathematicians who want to get a grasp of this quickly growing discipline and find themselves taken aback by the language gap between mathematics and the pioneering fields on the matter: computer science and quantum physics.


title: Interpolating between Optimal Transport and MMD using Sinkhorn Divergences

abstract: Comparing probability distributions is a fundamental problem in data sciences. Simple norms and divergences such as the total variation and the relative entropy only compare densities in a point-wise manner and fail to capture the geometric nature of the problem. In sharp contrast, Maximum Mean Discrepancies (MMD) and Optimal Transport distances (OT) are two classes of distances between measures that take into account the geometry of the underlying space and metrize the convergence in law. This paper studies the Sinkhorn divergences, a family of geometric divergences that interpolates between MMD and OT. Relying on a new notion of geometric entropy, we provide theoretical guarantees for these divergences: positivity, convexity and metrization of the convergence in law. On the practical side, we detail a numerical scheme that enables the large scale application of these divergences for machine learning: on the GPU, gradients of the Sinkhorn loss can be computed for batches of a million samples.


title: A Comparison of the Trojan Y Chromosome Strategy to Harvesting Models for Eradication of Non-Native Species

abstract: The Trojan Y Chromosome Strategy (TYC) is a promising eradication method for biological control of non-native species. The strategy works by manipulating the sex ratio of a population through the introduction of \textit{supermales} that guarantee male offspring. In the current manuscript, we compare the TYC method with a pure harvesting strategy. We also analyze a hybrid harvesting model that mirrors the TYC strategy. The dynamic analysis leads to results on stability, global boundedness of solutions and bifurcations of the model. Several conclusions about the different strategies are established via optimal control methods. In particular, the results affirm that either a pure harvesting or hybrid strategy may work better than the TYC method at controlling an invasive species population.


title: Exploring Adversarial Examples in Malware Detection

abstract: The convolutional neural network (CNN) architecture is increasingly being applied to new domains, such as malware detection, where it is able to learn malicious behavior from raw bytes extracted from executables. These architectures reach impressive performance with no feature engineering effort involved, but their robustness against active attackers is yet to be understood. Such malware detectors could face a new attack vector in the form of adversarial interference with the classification model. Existing evasion attacks intended to cause misclassification on test-time instances, which have been extensively studied for image classifiers, are not applicable because of the input semantics that prevents arbitrary changes to the binaries. This paper explores the area of adversarial examples for malware detection. By training an existing model on a production-scale dataset, we show that some previous attacks are less effective than initially reported, while simultaneously highlighting architectural weaknesses that facilitate new attack strategies for malware classification. Finally, we explore how generalizable different attack strategies are, the trade-offs when aiming to increase their effectiveness, and the transferability of single-step attacks.


title: Eigenvalue comparisons in Steklov eigenvalue problem and some other eigenvalue estimates

abstract: In this paper, two interesting eigenvalue comparison theorems for the first non-zero Steklov eigenvalue of the Laplacian have been established for manifolds with radial sectional curvature bounded from above. Besides, sharper bounds for the first non-zero eigenvalue of the Wentzell eigenvalue problem of the weighted Laplacian, which can be seen as a natural generalization of the classical Steklov eigenvalue problem, have been obtained.


title: Steepest-Entropy-Ascent Quantum Thermodynamics Models in Materials Science

abstract: Steepest-entropy-ascent quantum thermodynamics, or SEAQT, is a unified approach of quantum mechanics and thermodynamics that avoids many of the inconsistencies that can arise between the two theories. Given a set of energy levels, i.e., energy eigenstructure, accessible to a given physical system, SEAQT predicts the unique kinetic path from any initial non-equilibrium state to stable equilibrium by solving a master equation that directs the system along the path of steepest entropy ascent. There are no intrinsic limitations on the length and time scales the method can treat so it is well-suited for calculations where the dynamics over multiple spacial scales need to be taken into account within a single framework. In this paper, the theoretical framework and its advantages are described, and several applications are presented to illustrate the use of the SEAQT equation of motion and the construction of a simplified, reduced-order, energy eigenstructure.


title: Optimal Mesh Size for Inverse Medium Scattering Problems

abstract: An optimal mesh size of the sampling region can help to reduce computational burden in practical applications. In this work, we investigate optimal choices of mesh sizes for the identifications of medium obstacles from either the far-field or near-field data in two and three dimensions. The results would have applications in the reconstruction process of inverse scattering problems.


title: Social learning of prescribing behavior can promote population optimum of antibiotic use

abstract: The rise and spread of antibiotic resistance causes worsening medical cost and mortality especially for life-threatening bacteria infections, thereby posing a major threat to global health. Prescribing behavior of physicians is one of the important factors impacting the underlying dynamics of resistance evolution. It remains unclear when individual prescribing decisions can lead to the overuse of antibiotics on the population level, and whether population optimum of antibiotic use can be reached through an adaptive social learning process that governs the evolution of prescribing norm. Here we study a behavior-disease interaction model, specifically incorporating a feedback loop between prescription behavior and resistance evolution. We identify the conditions under which antibiotic resistance can evolve as a result of the tragedy of the commons in antibiotic overuse. Furthermore, we show that fast social learning that adjusts prescribing behavior in prompt response to resistance evolution can steer out cyclic oscillations of antibiotic usage quickly towards the stable population optimum of prescribing. Our work demonstrates that provision of prompt feedback to prescribing behavior with the collective consequences of treatment decisions and costs that are associated with resistance helps curb the overuse of antibiotics.


title: Log-symmetric regression models for correlated errors with an application to mortality data

abstract: Log-symmetric regression models are particularly useful when the response variable is continuous, strictly positive and asymmetric. In this paper, we proposed a class of log-symmetric regression models in the context of correlated errors. The proposed models provide a novel alternative to the existing log-symmetric regression models due to its flexibility in accommodating correlation. We discuss some properties, parameter estimation by the conditional maximum likelihood method and goodness of fit of the proposed model. We also provide expressions for the observed Fisher information matrix. A Monte Carlo simulation study is presented to evaluate the performance of the conditional maximum likelihood estimators. Finally, a full analysis of a real-world mortality data set is presented to illustrate the proposed approach.


title: A note on positive $\mathcal{AN}$ operators

abstract: We show that positive absolutely norm attaining operators can be characterized by a simple property of their spectra. This result clarifies and simplifies a result of Ramesh. As an application we characterize weighted shift operators which are absolutely norm attaining.


title: Performance Improvement in Noisy Linear Consensus Networks with Time-Delay

abstract: We analyze performance of a class of time-delay first-order consensus networks from a graph topological perspective and present methods to improve it. The performance is measured by network's square of H-2 norm and it is shown that it is a convex function of Laplacian eigenvalues and the coupling weights of the underlying graph of the network. First, we propose a tight convex, but simple, approximation of the performance measure in order to achieve lower complexity in our design problems by eliminating the need for eigen-decomposition. The effect of time-delay reincarnates itself in the form of non-monotonicity, which results in nonintuitive behaviors of the performance as a function of graph topology. Next, we present three methods to improve the performance by growing, re-weighting, or sparsifying the underlying graph of the network. It is shown that our suggested algorithms provide near-optimal solutions with lower complexity with respect to existing methods in literature.


title: Isomorphisms of $AC(\sigma)$ spaces for linear graphs

abstract: We show that among compact subsets of the plane which are drawings of linear graphs, two sets $\sigma$ and $\tau$ are homeomorphic if and only if the corresponding spaces of absolutely continuous functions (in the sense of Ashton and Doust) are isomorphic as Banach algebras. This gives an analogue for this class of sets to the well-known result of Gelfand and Kolmogorov for $C(\Omega)$ spaces.


title: Sample-Free Learning of Input Grammars for Comprehensive Software Fuzzing

abstract: Generating valid test inputs for a program is much easier if one knows the input language. We present first successes for a technique that, given a program P without any input samples or models, learns an input grammar that represents the syntactically valid inputs for P -- a grammar which can then be used for highly effective test generation for P . To this end, we introduce a test generator targeted at input parsers that systematically explores parsing alternatives based on dynamic tracking of constraints; the resulting inputs go into a grammar learner producing a grammar that can then be used for fuzzing. In our evaluation on subjects such as JSON, URL, or Mathexpr, our PYGMALION prototype took only a few minutes to infer grammars and generate thousands of valid high-quality inputs.


title: Deep Learning vs. Human Graders for Classifying Severity Levels of Diabetic Retinopathy in a Real-World Nationwide Screening Program

abstract: Deep learning algorithms have been used to detect diabetic retinopathy (DR) with specialist-level accuracy. This study aims to validate one such algorithm on a large-scale clinical population, and compare the algorithm performance with that of human graders. 25,326 gradable retinal images of patients with diabetes from the community-based, nation-wide screening program of DR in Thailand were analyzed for DR severity and referable diabetic macular edema (DME). Grades adjudicated by a panel of international retinal specialists served as the reference standard. Across different severity levels of DR for determining referable disease, deep learning significantly reduced the false negative rate (by 23%) at the cost of slightly higher false positive rates (2%). Deep learning algorithms may serve as a valuable tool for DR screening.


title: Qubit Allocation for Noisy Intermediate-Scale Quantum Computers

abstract: In the era of noisy-intermediate-scale quantum computers, we expect to see quantum devices with increasing numbers of qubits emerge in the foreseeable future. To practically run quantum programs, logical qubits have to be mapped to the physical qubits by a qubit allocation algorithm. However, on present day devices, qubits differ by their error rate and connectivity. Here, we establish and demonstrate on current experimental devices a new allocation algorithm that combines the simulated annealing method with local search of the solution space using Dijkstra's algorithm. Our algorithm takes into account the weighted connectivity constraints of both the quantum hardware and the quantum program being compiled. New quantum programs will enable unprecedented developments in physics, chemistry, and materials science and our work offers an important new pathway toward optimizing compilers for quantum programs.


title: A similarity measure for second order properties of non-stationary functional time series with applications to clustering and testing

abstract: Due to the surge of data storage techniques, the need for the development of appropriate techniques to identify patterns and to extract knowledge from the resulting enormous data sets, which can be viewed as collections of dependent functional data, is of increasing interest in many scientific areas. We develop a similarity measure for spectral density operators of a collection of functional time series, which is based on the aggregation of Hilbert-Schmidt differences of the individual time-varying spectral density operators. Under fairly general conditions, the asymptotic properties of the corresponding estimator are derived and asymptotic normality is established. The introduced statistic lends itself naturally to quantify (dis)-similarity between functional time series, which we subsequently exploit in order to build a spectral clustering algorithm. Our algorithm is the first of its kind in the analysis of non-stationary (functional) time series and enables to discover particular patterns by grouping together `similar' series into clusters, thereby reducing the complexity of the analysis considerably. The algorithm is simple to implement and computationally feasible. As a further application we provide a simple test for the hypothesis that the second order properties of two non-stationary functional time series coincide.


title: CURE-OR: Challenging Unreal and Real Environments for Object Recognition

abstract: In this paper, we introduce a large-scale, controlled, and multi-platform object recognition dataset denoted as Challenging Unreal and Real Environments for Object Recognition (CURE-OR). In this dataset, there are 1,000,000 images of 100 objects with varying size, color, and texture that are positioned in five different orientations and captured using five devices including a webcam, a DSLR, and three smartphone cameras in real-world (real) and studio (unreal) environments. The controlled challenging conditions include underexposure, overexposure, blur, contrast, dirty lens, image noise, resizing, and loss of color information. We utilize CURE-OR dataset to test recognition APIs-Amazon Rekognition and Microsoft Azure Computer Vision- and show that their performance significantly degrades under challenging conditions. Moreover, we investigate the relationship between object recognition and image quality and show that objective quality algorithms can estimate recognition performance under certain photometric challenging conditions. The dataset is publicly available at https://ghassanalregib.com/cure-or/.


title: Linearized Analysis of Barotropic Perturbations around Spherically Symmetric Gaseous Stars Governed by the Euler-Poisson Equations

abstract: The hydrodynamic evolution of self-gravitating gaseous stars is governed by the Euler-Poisson equations. We study the structure of the linear approximation of barotropic perturbations around spherically symmetric equilibria based on functional analytic tools. In contrast to folklore, we show that the spectrum of the linearized operator for general perturbations is not of the Sturm-Liouville type unless the perturbations are restricted. In particular, we prove that it is of the Sturm-Liouville type for irrotational perturbations.


title: Approximate functional equation and upper bounds for the Barnes double zeta-function

abstract: As one of the asymptotic formulas of the zeta-function, Hardy and Littlewood gave asymptotic formulas called the approximate functional equation. In this paper, we prove an approximate functional equation of the Barnes double zeta-function $ \zeta_2 (s, \alpha ; v, w ) = \sum_{m=0}^\infty \sum_{n=0}^\infty (\alpha+vm+wn)^{-s} $. Also, applying this approximate functional equation and the van der Corput method, we obtain upper bounds for $ \zeta_2(1/2 + it, \alpha ; v, w) $ and $ \zeta_2(3/2 + it, \alpha ; v, w) $ with respect to $ t $ as $ t \rightarrow \infty $.


title: Strong entanglement criterion involving momentum weak values

abstract: In recent years weak values have been used to explore interesting quantum features in novel ways. In particular, the real part of the weak value of the momentum operator has been widely studied, mainly in connection with (nonlocal) Bohmian trajectories. Here we focus on the imaginary part and its role in relation with the entanglement of a bipartite system. We establish an entanglement criterion based on weak momentum correlations, that allows to discern whether the entanglement is encoded in the amplitude and/or in the phase of the wave function. Our results throw light on the physical role of the real and imaginary parts of the weak values, and stress the relevance of the latter in the multi-particle scenario


title: Dynamic Automatic Differentiation of GPU Broadcast Kernels

abstract: We show how forward-mode automatic differentiation (AD) can be employed within larger reverse-mode computations to dynamically differentiate broadcast operations in a GPU-friendly manner. Our technique fully exploits the broadcast Jacobian's inherent sparsity structure, and unlike a pure reverse-mode approach, this "mixed-mode" approach does not require a backwards pass over the broadcasted operation's subgraph, obviating the need for several reverse-mode-specific programmability restrictions on user-authored broadcast operations. Most notably, this approach allows broadcast fusion in primal code despite the presence of data-dependent control flow. We discuss an experiment in which a Julia implementation of our technique outperformed pure reverse-mode TensorFlow and Julia implementations for differentiating through broadcast operations within an HM-LSTM cell update calculation.


title: Stochastic Primal-Dual Q-Learning

abstract: In this work, we present a new model-free and off-policy reinforcement learning (RL) algorithm, that is capable of finding a near-optimal policy with state-action observations from arbitrary behavior policies. Our algorithm, called the stochastic primal-dual Q-learning (SPD Q-learning), hinges upon a new linear programming formulation and a dual perspective of the standard Q-learning. In contrast to previous primal-dual RL algorithms, the SPD Q-learning includes a Q-function estimation step, thus allowing to recover an approximate policy from the primal solution as well as the dual solution. We prove a first-of-its-kind result that the SPD Q-learning guarantees a certain convergence rate, even when the state-action distribution is time-varying but sub-linearly converges to a stationary distribution. Numerical experiments are provided to demonstrate the off-policy learning abilities of the proposed algorithm in comparison to the standard Q-learning.


title: Singular link concordance implies link homotopy in codimension $\ge 3$

abstract: We prove the analogue of the Concordance Implies Isotopy in Codimension $\ge 3$ Theorem for link maps, together with some other its singular analogues. In the case of spherical link maps, a stronger result was independently obtained by P. Teichner (by different methods).


title: Versatile laser-free trapped-ion entangling gates

abstract: We present a general theory for laser-free entangling gates with trapped-ion hyperfine qubits, using either static or oscillating magnetic-field gradients combined with a pair of uniform microwave fields symmetrically detuned about the qubit frequency. By transforming into a `bichromatic' interaction picture, we show that either ${\hat{\sigma}_{\phi}\otimes\hat{\sigma}_{\phi}}$ or ${\hat{\sigma}_{z}\otimes\hat{\sigma}_{z}}$ geometric phase gates can be performed. The gate basis is determined by selecting the microwave detuning. The driving parameters can be tuned to provide intrinsic dynamical decoupling from qubit frequency fluctuations. The ${\hat{\sigma}_{z}\otimes\hat{\sigma}_{z}}$ gates can be implemented in a novel manner which eases experimental constraints. We present numerical simulations of gate fidelities assuming realistic parameters.


title: SOFIA FIFI-LS Observations of Sgr B1: Ionization Structure and Sources of Excitation

abstract: The current paradigm of Galactic Center (GC) gas motions and star formation envisions sequential star formation in streams of gas as they pass near the supermassive black hole, Sgr A*. This is based on the relative positions of dense molecular clouds, the very young star-forming region Sgr B2, the much older region Sgr C, and the several Myr old Arches and Quintuplet Clusters. Because Sgr B1 is found with Sgr B2 in a common envelope of molecular gas and far-infrared emission, the two sources are thought to be physically related, even though there are indicators of a significantly greater age for Sgr B1. To clarify the status of Sgr B1, we have mapped it with the FIFI-LS spectrometer on SOFIA in the far-infrared lines of [O III] 52 and 88 micron. From the ratios of these lines and lines measured with the Spitzer Infrared Spectrograph, we find that there are at least eight separate sub-regions that must contain the stars that excite the gas. We infer spectral energy distributions (SEDs) of the ionizing sources from models and find they are in agreement only with SEDs of late O stars augmented at the highest frequencies with interstellar X-rays from fast shocks. We suggest that although the gas, from its velocity structure, must be part of the very young Sgr B2 complex, the stars that are ionizing the gas were not formed there but are the remnants of a previous generation of star formation in the GC.


title: Ferromagnetic spin fluctuations in the filled skutterudite SrFe4As12 revealed by 75As NMR-NQR measurements

abstract: $^{75}$As nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) measurements have been carried out to investigate the magnetic and electronic properties of the filled skutterudite metallic compound SrFe$_4$As$_{12}$. The temperature dependence of Knight shift $K$ determined by the NQR spectrum under a small magnetic field ($\le$ 0.5 T) shows the similar temperature dependence of the magnetic susceptibility $\chi$ which exhibits a broad maximum at $T^\ast$ $\sim$ 50 K. The nuclear spin-lattice relaxation rate divided by temperature, 1/$T_1T$, increases with decreasing temperature and exhibits a broad maximum at $T$ $\sim$ 70 K, similar to the case of $\chi$. The temperature dependence of $K$ and $1/T_1T$ is reasonably explained by a simple model where we assume a concave-shaped band structure near the Fermi energy. Based on a Korringa ratio analysis using the $T_1$ and $K$ data, ferromagnetic spin fluctuations are found to exist in SrFe$_4$As$_{12}$. These results indicate that SrFe$_4$As$_{12}$ can be characterized to be a metal with ferromagnetic correlations and also the peculiar band structure responsible for the suppression of $1/T_1T$ and $K$ at low temperatures.


title: Compositional Verification for Autonomous Systems with Deep Learning Components

abstract: As autonomy becomes prevalent in many applications, ranging from recommendation systems to fully autonomous vehicles, there is an increased need to provide safety guarantees for such systems. The problem is difficult, as these are large, complex systems which operate in uncertain environments, requiring data-driven machine-learning components. However, learning techniques such as Deep Neural Networks, widely used today, are inherently unpredictable and lack the theoretical foundations to provide strong assurance guarantees. We present a compositional approach for the scalable, formal verification of autonomous systems that contain Deep Neural Network components. The approach uses assume-guarantee reasoning whereby {\em contracts}, encoding the input-output behavior of individual components, allow the designer to model and incorporate the behavior of the learning-enabled components working side-by-side with the other components. We illustrate the approach on an example taken from the autonomous vehicles domain.


title: Anisotropic liquid drop models

abstract: We introduce and study certain variants of Gamow's liquid drop model in which an anisotropic surface energy replaces the perimeter. After existence and nonexistence results are established, the shape of minimizers is analyzed. Under suitable regularity and ellipticity assumptions on the surface tension, Wulff shapes are minimizers in this problem if and only if the surface energy is isotropic. In sharp contrast, Wulff shapes are the unique minimizers for certain crystalline surface tensions. We also introduce and study several related liquid drop models with anisotropic repulsion for which the Wulff shape is the minimizer in the small mass regime.


title: Open Vocabulary Learning on Source Code with a Graph-Structured Cache

abstract: Machine learning models that take computer program source code as input typically use Natural Language Processing (NLP) techniques. However, a major challenge is that code is written using an open, rapidly changing vocabulary due to, e.g., the coinage of new variable and method names. Reasoning over such a vocabulary is not something for which most NLP methods are designed. We introduce a Graph-Structured Cache to address this problem; this cache contains a node for each new word the model encounters with edges connecting each word to its occurrences in the code. We find that combining this graph-structured cache strategy with recent Graph-Neural-Network-based models for supervised learning on code improves the models' performance on a code completion task and a variable naming task --- with over $100\%$ relative improvement on the latter --- at the cost of a moderate increase in computation time.


title: PYTHIA and the preoccupied proton

abstract: The PYTHIA Monte Carlo (PMC) has been applied broadly to simulations of high-energy $p$-$p$ and $p$-$\bar p$ collisions. The PMC is based on several assumptions, such as that most hadrons result from jet production (multiple parton interactions or MPIs), that $p$-$p$ centrality is relevant and that color reconnection (CR) strongly influences fragmentation to jets. An alternative description is provided by the two-component (soft + hard) model (TCM) of hadron production. TCM analysis of $p$-Pb ensemble-mean-$p_t$ data reveals centrality trends quite different from those estimated via a geometric Glauber model based on the eikonal approximation. Glauber estimates of binary-collision number are three times TCM estimates. Detailed study of $p$-Pb data conflicts with a basic Glauber assumption -- that a projectile proton may interact simultaneously with multiple target nucleons. Instead, in both $p$-$p$ and $p$-A collisions, a $p$-N collision once initiated is exclusive of other possible interactions (during that collision), and within the collision any pair of participant partons may interact -- a $p$-N collision is thus "all or nothing." In this presentation the PMC is challenged by an assortment of contradictory data, and evidence for $p$-N exclusivity is reviewed to make a case for the "preoccupied proton" of the title.


title: Reduction of Parameter Redundancy in Biaffine Classifiers with Symmetric and Circulant Weight Matrices

abstract: Currently, the biaffine classifier has been attracting attention as a method to introduce an attention mechanism into the modeling of binary relations. For instance, in the field of dependency parsing, the Deep Biaffine Parser by Dozat and Manning has achieved state-of-the-art performance as a graph-based dependency parser on the English Penn Treebank and CoNLL 2017 shared task. On the other hand, it is reported that parameter redundancy in the weight matrix in biaffine classifiers, which has O(n^2) parameters, results in overfitting (n is the number of dimensions). In this paper, we attempted to reduce the parameter redundancy by assuming either symmetry or circularity of weight matrices. In our experiments on the CoNLL 2017 shared task dataset, our model achieved better or comparable accuracy on most of the treebanks with more than 16% parameter reduction.


title: Dyson Brownian Motion for General $\beta$ and Potential at the Edge

abstract: In this paper, we compare the solutions of Dyson Brownian motion with general $\beta$ and potential $V$ and the associated McKean-Vlasov equation near the edge. Under suitable conditions on the initial data and potential $V$, we obtain the optimal rigidity estimates of particle locations near the edge for short time $t=\text{o}(1)$. Our argument uses the method of characteristics along with a careful estimate involving an equation of the edge. With the rigidity estimates as an input, we prove a central limit theorem for mesoscopic statistics near the edge which, as far as we know, have been done for the first time in this paper. Additionally, combining with \cite{LandonEdge}, our rigidity estimates are used to give a proof of the local ergodicity of Dyson Brownian motion for general $\beta$ and potential at the edge, i.e. the distribution of extreme particles converges to Tracy-Widom $\beta$ distribution in short time.


title: Unsupervised Anomalous Data Space Specification

abstract: Computer algorithms are written with the intent that when run they perform a useful function. Typically any information obtained is unknown until the algorithm is run. However, if the behavior of an algorithm can be fully described by precomputing just once how this algorithm will respond when executed on any input, this precomputed result provides a complete specification for all solutions in the problem domain. We apply this idea to a previous anomaly detection algorithm, and in doing so transform it from one that merely detects individual anomalies when asked to discover potentially anomalous values, into an algorithm also capable of generating a complete specification for those values it would deem to be anomalous. This specification is derived by examining no more than a small training data, can be obtained in very small constant time, and is inherently far more useful than results obtained by repeated execution of this tool. For example, armed with such a specification one can ask how close an anomaly is to being deemed normal, and can validate this answer not by exhaustively testing the algorithm but by examining if the specification so generated is indeed correct. This powerful idea can be applied to any algorithm whose runtime behavior can be recovered from its construction and so has wide applicability.


title: An Efficient Data Retrieval Parallel Reeb Graph Algorithm

abstract: The Reeb graph of a scalar function defined on a domain gives a topologically meaningful summary of that domain. Reeb graphs have been shown in the past decade to be of great importance in geometric processing, image processing, computer graphics, and computational topology. The demand for analyzing large data sets has increased in the last decade. Hence the parallelization of topological computations needs to be more fully considered. We propose a parallel augmented Reeb graph algorithm on triangulated meshes with and without a boundary. That is, in addition to our parallel algorithm for computing a Reeb graph, we describe a method for extracting the original manifold data from the Reeb graph structure. We demonstrate the running time of our algorithm on standard datasets. As an application, we show how our algorithm can be utilized in mesh segmentation algorithms.


title: Zero-Forcing Per-Group Precoding (ZF-PGP) for Robust Optimized Downlink Massive MIMO Performance

abstract: In this paper, we propose a new, combined Zero-Forcing Per-Group Precoding (ZF-PGP) method that achieves very high gains in comparison to ZF Precoding techniques, while it simultaneously offers individually separate streams to reach individual User Equipment (UE), i.e., it obliterates the need for coordinated, joint decoding by the group's UEs. Although ZF-PGP in general experiences a performance loss in comparison to the near-optimal Per-Group Precoding within groups (PGP-WG), it is much simpler to implement than PGP-WG. The method works extremely well in conjunction with Combined Frequency and Space Division and Multiplexing (CFSDM) in correlated channels or in cases in which the group's total number of antennas is greater than the group's available Virtual Channel Model Beams (VCMB) number. As we show by multiple examples, for a Uniform Linear Array (ULA), by combining the proposed technique with massive Multiple-Input Multiple-Output (MIMO) and Joint Spatial Division and Multiplexing for Finite Alphabets (JSDM-FA) type hybrid precoding or CFSDM a very efficient, near-optimal design with reduced complexity at the UEs is achieved, while each UE decodes its received data independently, without interference from data destined for other UEs. Due to JSDM-FA Statistical Channel State Information (SCSI) inner beamforming, the groups are easily identified. We show that the combined ZF-PGP scheme achieves very high gains in comparison to ZF and Regularized ZF (RZF) precoding and other modified sub-optimal techniques investigated, while its complexity is still relatively low due to PGP-WG. We also demonstrate the robustness of the proposed precoding techniques to channel estimation errors, showing minimal performance loss.


title: DeepXS: Fast approximation of MSSM electroweak cross sections at NLO

abstract: We present a deep learning solution to the prediction of particle production cross sections over a complicated, high-dimensional parameter space. We demonstrate the applicability by providing state-of-the-art predictions for the production of charginos and neutralinos at the Large Hadron Collider (LHC) at the next-to-leading order in the phenomenological MSSM-19 and explicitly demonstrate the performance for $pp\to\tilde{\chi}^+_1\tilde{\chi}^-_1,$ $\tilde{\chi}^0_2\tilde{\chi}^0_2$ and $\tilde{\chi}^0_2\tilde{\chi}^\pm_1$ as a proof of concept which will be extended to all SUSY electroweak pairs. We obtain errors that are lower than the uncertainty from scale and parton distribution functions with mean absolute percentage errors of well below $0.5\,\%$ allowing a safe inference at the next-to-leading order with inference times that improve the Monte Carlo integration procedures that have been available so far by a factor of $\mathcal{O}(10^7)$ from $\mathcal{O}(\rm{min})$ to $\mathcal{O}(\mu\rm{s})$ per evaluation.


title: Adaptive Communication Strategies to Achieve the Best Error-Runtime Trade-off in Local-Update SGD

abstract: Large-scale machine learning training, in particular distributed stochastic gradient descent, needs to be robust to inherent system variability such as node straggling and random communication delays. This work considers a distributed training framework where each worker node is allowed to perform local model updates and the resulting models are averaged periodically. We analyze the true speed of error convergence with respect to wall-clock time (instead of the number of iterations), and analyze how it is affected by the frequency of averaging. The main contribution is the design of AdaComm, an adaptive communication strategy that starts with infrequent averaging to save communication delay and improve convergence speed, and then increases the communication frequency in order to achieve a low error floor. Rigorous experiments on training deep neural networks show that AdaComm can take $3 \times$ less time than fully synchronous SGD, and still reach the same final training loss.


title: Minor-closed graph classes with bounded layered pathwidth

abstract: We prove that a minor-closed class of graphs has bounded layered pathwidth if and only if some apex-forest is not in the class. This generalises a theorem of Robertson and Seymour, which says that a minor-closed class of graphs has bounded pathwidth if and only if some forest is not in the class.


title: A Comparative Analysis of Registration Tools: Traditional vs Deep Learning Approach on High Resolution Tissue Cleared Data

abstract: Image registration plays an important role in comparing images. It is particularly important in analyzing medical images like CT, MRI, PET, etc. to quantify different biological samples, to monitor disease progression and to fuse different modalities to support better diagnosis. The recent emergence of tissue clearing protocols enable us to take images at cellular level resolution. Image registration tools developed for other modalities are currently unable to manage images of such extreme high resolution. The recent popularity of deep learning based methods in the computer vision community justifies a rigorous investigation of deep-learning based methods on tissue cleared images along with their traditional counterparts. In this paper, we investigate and compare the performance of a deep learning based registration method with traditional optimization based methods on samples from tissue-clearing methods. From the comparative results it is found that a deep-learning based method outperforms all traditional registration tools in terms of registration time and has achieved promising registration accuracy.


title: Heteroskedastic PCA: Algorithm, Optimality, and Applications

abstract: A general framework for principal component analysis (PCA) in the presence of heteroskedastic noise is introduced. We propose an algorithm called HeteroPCA, which involves iteratively imputing the diagonal entries of the sample covariance matrix to remove estimation bias due to heteroskedasticity. This procedure is computationally efficient and provably optimal under the generalized spiked covariance model. A key technical step is a deterministic robust perturbation analysis on singular subspaces, which can be of independent interest. The effectiveness of the proposed algorithm is demonstrated in a suite of problems in high-dimensional statistics, including singular value decomposition (SVD) under heteroskedastic noise, Poisson PCA, and SVD for heteroskedastic and incomplete data.


title: Enabling Grasp Action: Generalized Evaluation of Grasp Stability via Contact Stiffness from Contact Mechanics Insight

abstract: Performing a grasp is a pivotal capability for a robotic gripper. We propose a new evaluation approach of grasping stability via constructing a model of grasping stiffness based on the theory of contact mechanics. First, the mathematical models are built to explore soft contact and the general grasp stiffness between a finger and an object. Next, the grasping stiffness matrix is constructed to reflect the normal, tangential and torsion stiffness coefficients. Finally, we design two grasping cases to verify the proposed measurement criterion of grasping stability by comparing different grasping configurations. Specifically, a standard grasping index is used and compared with the minimum eigenvalue index of the constructed grasping stiffness we built. The comparison result reveals a similar tendency between them for measuring the grasping stability and thus, validates the proposed approach.


title: Newly observed $\Omega(2012)$ state and strong decays of the low-lying $\Omega$ excitations

abstract: Stimulated by the newly discovered $\Omega(2012)$ resonance at Belle II, in this work we have studied the OZI allowed strong decays of the low-lying $1P$- and $1D$-wave $\Omega$ baryons within the $^3P_0$ model. It is found that $\Omega(2012)$ is most likely to be a $1P$-wave $\Omega$ state with $J^P=3/2^-$. We also find that the $\Omega(2250)$ state could be assigned as a $1D$-wave state with $J^P=5/2^+$. The other missing $1P$- and $1D$-wave $\Omega$ baryons may have large potentials to be observed in their main decay channels.


title: Power and Efficiency of a Thermal Engine with a Coherent Bath

abstract: We consider a quantum engine driven by repeated weak interactions with a heat bath of identical three-level atoms. This model was first introduced by Scully et al. [Science, 2003], who showed that coherence between the energy-degenerate ground states serves as a thermodynamic resource that allows operation of a thermal cycle with a coherence-dependent thermalisation temperature. We consider a similar engine out of the quasistatic limit and find that the ground-state coherence also determines the rate of thermalisation, therefore increasing the output power and the engine efficiency only when the thermalisation temperature is reduced; revealing a more nuanced perspective of coherence as a resource. This allows us to optimise the output power by adjusting the coherence and relative stroke durations.


title: Evidence for quasi-one-dimensional charge density wave in CuTe by angle-resolved photoemission spectroscopy

abstract: We report the electronic structure of CuTe with a high charge density wave (CDW) transition temperature Tc = 335 K by angle-resolved photoemission spectroscopy (ARPES). An anisotropic charge density wave gap with a maximum value of 190 meV is observed in the quasi-one-dimensional band formed by Te px orbitals. The CDW gap can be filled by increasing temperature or electron doping through in situ potassium deposition. Combining the experimental results with calculated electron scattering susceptibility and phonon dispersion, we suggest that both Fermi surface nesting and electron-phonon coupling play important roles in the emergence of the CDW.


title: Effective $n$-type Doping of Mg$_3$Sb$_2$ with Group-3 Elements

abstract: The recent discovery of high thermoelectric performance in Mg$_3$Sb$_2$ has been critically enabled by the success in $n$-type doping of this material, which is achieved under Mg-rich growth conditions, typically with chalcogens (Se, Te) as extrinsic dopants. Using first-principles defect calculations, we previously predicted that higher electron concentrations ($\sim10^{20}$ cm$^{-3}$) can be achieved in Mg$_3$Sb$_2$ by doping with La instead of Se or Te. Subsequent experiments showed that free electron concentration in La-doped Mg$_3$Sb$_{2-x}$Bi$_x$ indeed exceeds those in the Te-doped material. Herein, we further investigate $n$-type doping of Mg$_3$Sb$_2$ and predict that, in addition to La, other group-3 elements (Sc, Y) are also effective as $n$-type dopants; Y is as good as La while Sc slightly less. Overall, we find that doping with any group-3 elements should lead to higher free electron concentrations than doping with chalcogens.


title: Sequenced-Replacement Sampling for Deep Learning

abstract: We propose sequenced-replacement sampling (SRS) for training deep neural networks. The basic idea is to assign a fixed sequence index to each sample in the dataset. Once a mini-batch is randomly drawn in each training iteration, we refill the original dataset by successively adding samples according to their sequence index. Thus we carry out replacement sampling but in a batched and sequenced way. In a sense, SRS could be viewed as a way of performing "mini-batch augmentation". It is particularly useful for a task where we have a relatively small images-per-class such as CIFAR-100. Together with a longer period of initial large learning rate, it significantly improves the classification accuracy in CIFAR-100 over the current state-of-the-art results. Our experiments indicate that training deeper networks with SRS is less prone to over-fitting. In the best case, we achieve an error rate as low as 10.10%.


title: Learning Multi-Layer Transform Models

abstract: Learned data models based on sparsity are widely used in signal processing and imaging applications. A variety of methods for learning synthesis dictionaries, sparsifying transforms, etc., have been proposed in recent years, often imposing useful structures or properties on the models. In this work, we focus on sparsifying transform learning, which enjoys a number of advantages. We consider multi-layer or nested extensions of the transform model, and propose efficient learning algorithms. Numerical experiments with image data illustrate the behavior of the multi-layer transform learning algorithm and its usefulness for image denoising. Multi-layer models provide better denoising quality than single layer schemes.


title: On the structure of $\lambda$-Cantor set with overlaps

abstract: Given $\lambda\in(0, 1)$, let $E_\lambda$ be the self-similar set generated by the iterated function system $\{x/3,(x+\lambda)/3,(x+2)/3\}$. Then $E_\lambda$ is a self-similar set with overlaps. We obtain the necessary and sufficient condition for $E_\lambda$ to be totally self-similar, which is a concept first introduced by Broomhead, Montaldi, and Sidorov in 2004. When $E_\lambda$ is totally self-similar, all its generating IFSs are investigated, and the size of the set of points having finite triadic codings is determined. Besides, we give some properties of the spectrum of $E_\lambda$ and show that the spectrum of $E_\lambda$ vanishes if and only if $\lambda$ is irrational.


title: Ordered Size Ramsey Number of Paths

abstract: An ordered graph is a simple graph with an ordering on its vertices. Define the ordered path $P_n$ to be the monotone increasing path with $n$ edges. The ordered size Ramsey number $\tilde{r}(P_r,P_s)$ is the minimum number $m$ for which there exists an ordered graph $H$ with $m$ edges such that every two-coloring of the edges of $H$ contains a red copy of $P_r$ or a blue copy of $P_s$. For $2\leq r\leq s$, we show $\frac{1}{8}r^2s\leq \tilde{r}(P_r,P_s)\leq Cr^2s(\log s)^3$, where $C>0$ is an absolute constant. This problem is motivated by the recent results of Buci\'c-Letzter-Sudakov and Letzter-Sudakov for oriented graphs.


title: Domain-Invariant Projection Learning for Zero-Shot Recognition

abstract: Zero-shot learning (ZSL) aims to recognize unseen object classes without any training samples, which can be regarded as a form of transfer learning from seen classes to unseen ones. This is made possible by learning a projection between a feature space and a semantic space (e.g. attribute space). Key to ZSL is thus to learn a projection function that is robust against the often large domain gap between the seen and unseen classes. In this paper, we propose a novel ZSL model termed domain-invariant projection learning (DIPL). Our model has two novel components: (1) A domain-invariant feature self-reconstruction task is introduced to the seen/unseen class data, resulting in a simple linear formulation that casts ZSL into a min-min optimization problem. Solving the problem is non-trivial, and a novel iterative algorithm is formulated as the solver, with rigorous theoretic algorithm analysis provided. (2) To further align the two domains via the learned projection, shared semantic structure among seen and unseen classes is explored via forming superclasses in the semantic space. Extensive experiments show that our model outperforms the state-of-the-art alternatives by significant margins.


title: Decoupling carrier concentration and electron-phonon coupling in oxide heterostructures observed with resonant inelastic x-ray scattering

abstract: We report the observation of multiple phonon satellite features in ultra thin superlattices of form $n$SrIrO$_3$/$m$SrTiO$_3$ using resonant inelastic x-ray scattering. As the values of $n$ and $m$ vary the energy loss spectra show a systematic evolution in the relative intensity of the phonon satellites. Using a closed-form solution for the cross section, we extract the variation in the electron-phonon coupling strength as a function of $n$ and $m$. Combined with the negligible carrier doping into the SrTiO$_3$ layers, these results indicate that tuning of the electron-phonon coupling can be effectively decoupled from doping. This work showcases both a feasible method to extract the electron-phonon coupling in superlattices and unveils a potential route for tuning this coupling which is often associated with superconductivity in SrTiO$_3$-based systems.


title: Finite Groups with a Prescribed Number of Cyclic Subgroups II

abstract: T\u{a}rn\u{a}uceanu described the finite groups $G$ having exactly $|G|-1$ cyclic subgroups. In "Finite Groups with a Prescribed Number of Cyclic Subgroups,", we used elementary methods to completely characterize those finite groups $G$ having exactly $|G|-\Delta$ cyclic subgroups for $\Delta=2, 3, 4$ and $5$. In this paper, we prove that for any $\Delta >0$ if $G$ has exactly $|G|-\Delta$ cyclic subgroups, then $|G|\le 8\Delta$ and therefore the number of such $G$ is finite. We then use the computer program GAP to find all $G$ with exactly $|G|-\Delta$ cyclic subgroups for $\Delta=1,\ldots,32$.


title: Transferrable Feature and Projection Learning with Class Hierarchy for Zero-Shot Learning

abstract: Zero-shot learning (ZSL) aims to transfer knowledge from seen classes to unseen ones so that the latter can be recognised without any training samples. This is made possible by learning a projection function between a feature space and a semantic space (e.g. attribute space). Considering the seen and unseen classes as two domains, a big domain gap often exists which challenges ZSL. Inspired by the fact that an unseen class is not exactly `unseen' if it belongs to the same superclass as a seen class, we propose a novel inductive ZSL model that leverages superclasses as the bridge between seen and unseen classes to narrow the domain gap. Specifically, we first build a class hierarchy of multiple superclass layers and a single class layer, where the superclasses are automatically generated by data-driven clustering over the semantic representations of all seen and unseen class names. We then exploit the superclasses from the class hierarchy to tackle the domain gap challenge in two aspects: deep feature learning and projection function learning. First, to narrow the domain gap in the feature space, we integrate a recurrent neural network (RNN) defined with the superclasses into a convolutional neural network (CNN), in order to enforce the superclass hierarchy. Second, to further learn a transferrable projection function for ZSL, a novel projection function learning method is proposed by exploiting the superclasses to align the two domains. Importantly, our transferrable feature and projection learning methods can be easily extended to a closely related task -- few-shot learning (FSL). Extensive experiments show that the proposed model significantly outperforms the state-of-the-art alternatives in both ZSL and FSL tasks.


title: A common trajectory recapitulated by urban economies

abstract: Is there a general economic pathway recapitulated by individual cities over and over? Identifying such evolution structure, if any, would inform models for the assessment, maintenance, and forecasting of urban sustainability and economic success as a quantitative baseline. This premise seems to contradict the existing body of empirical evidences for path-dependent growth shaping the unique history of individual cities. And yet, recent empirical evidences and theoretical models have amounted to the universal patterns, mostly size-dependent, thereby expressing many of urban quantities as a set of simple scaling laws. Here, we provide a mathematical framework to integrate repeated cross-sectional data, each of which freezes in time dimension, into a frame of reference for longitudinal evolution of individual cities in time. Using data of over 100 millions employment in thousand business categories between 1998 and 2013, we decompose each city's evolution into a pre-factor and relative changes to eliminate national and global effects. In this way, we show the longitudinal dynamics of individual cities recapitulate the observed cross-sectional regularity. Larger cities are not only scaled-up versions of their smaller peers but also of their past. In addition, our model shows that both specialization and diversification are attributed to the distribution of industry's scaling exponents, resulting a critical population of 1.2 million at which a city makes an industrial transition into innovative economies.


title: The Novel Symmetry Constraint and Binary Nonlinearization of the Super Generalized Broer-Kaup Hierarchy with Self-consistent Sources and Conservation Laws

abstract: The super generalized Broer-Kaup(gBK) hierarchy and its super Hamiltonian structure are established based on a loop super Lie algebra and super-trace identity. Then the self-consistent sources, the conservation laws, the novel symmetry constraint and the binary nonlinearization of the super gBK hierarchy are generated, respectively. In addition, the integrals of motion required for Liouville integrability are explicitly given.


title: Zero and Few Shot Learning with Semantic Feature Synthesis and Competitive Learning

abstract: Zero-shot learning (ZSL) is made possible by learning a projection function between a feature space and a semantic space (e.g.,~an attribute space). Key to ZSL is thus to learn a projection that is robust against the often large domain gap between the seen and unseen class domains. In this work, this is achieved by unseen class data synthesis and robust projection function learning. Specifically, a novel semantic data synthesis strategy is proposed, by which semantic class prototypes (e.g., attribute vectors) are used to simply perturb seen class data for generating unseen class ones. As in any data synthesis/hallucination approach, there are ambiguities and uncertainties on how well the synthesised data can capture the targeted unseen class data distribution. To cope with this, the second contribution of this work is a novel projection learning model termed competitive bidirectional projection learning (BPL) designed to best utilise the ambiguous synthesised data. Specifically, we assume that each synthesised data point can belong to any unseen class; and the most likely two class candidates are exploited to learn a robust projection function in a competitive fashion. As a third contribution, we show that the proposed ZSL model can be easily extended to few-shot learning (FSL) by again exploiting semantic (class prototype guided) feature synthesis and competitive BPL. Extensive experiments show that our model achieves the state-of-the-art results on both problems.


title: Angular Distributions of Thomson Scattering in Combined Laser and Magnetic fields

abstract: Angular distributions of Thomson scattering are researched in the combined fields with a circularly polarization of laser field and a strong uniform magnetic field. The trajectories of the electron and the dependence of it on the initial phase are also given. It is found that the angular distributions with respect to the azimuthal angle show twofold symmetry whatever the laser intensity, the order of harmonics, the resonance parameter, and the initial axial momentum are. On the other hand, the radiation with respect to the polar angle is mainly distributed in two regions which are roughly symmetric of the laser propagation. In addition, the larger the laser intensity, the resonance parameters, the initial axial momentum are, the closer radiation is to the laser propagation direction. Besides, a new possibility of X-ray production is indicated. That is to say, with appropriate choice of laser and electron parameters, the high frequency part of the Thomson scattering radiation can reach the frequency range of X-ray (10^17 Hz - 10^18 Hz).


title: Regime-Switching Jump Diffusions with Non-Lipschitz Coefficients and Countably Many Switching States: Existence and Uniqueness, Feller, and Strong Feller Properties

abstract: This work focuses on a class of regime-switching jump diffusion processes, which is a two component Markov processes $(X(t),\Lambda(t))$, where $\Lambda(t)$ is a component representing discrete events taking values in a countably infinite set. Considering the corresponding stochastic differential equations, our main focus is on treating those with non-Lipschitz coefficients. We first show that there exists a unique strong solution to the corresponding stochastic differential equation. Then Feller and strong Feller properties are investigated.


title: Multilinear SVD for Millimeter Wave Channel Parameter Estimation

abstract: Fifth generation (5G) cellular standards are set to utilize millimeter wave (mmWave) frequencies, which enable data speeds greater than 10 Gbps and sub-centimeter localization accuracy. These capabilities rely on accurate estimates of the channel parameters, which we define as the angle of arrival, angle of departure, and path distance for each path between the transmitter and receiver. Estimating the channel parameters in a computationally efficient manner poses a challenge because it requires estimation of parameters from a high-dimensional measurement -- particularly for multi-carrier systems since each subcarrier must be estimated separately. Additionally, channel parameter estimation must be able to handle hybrid beamforming, which uses a combination of digital and analog beamforming to reduce the number of required analog to digital converters. This paper introduces a channel parameter estimation technique based on the multilinear singular value decomposition (MSVD), a tensor analog of the singular value decomposition, for massive multiple input multiple output (MIMO) multi-carrier systems with hybrid beamforming. The MSVD tensor estimation approach provides a computationally efficient method and is shown to closely match the Cramer-Rao bound (CRB) of parameter estimates through simulations. Limitations of channel parameter estimation and communication waveform effects are also studied.


title: A note on spanning trees of connected $K_{1,t}$-free graphs whose stems have a few leaves

abstract: Let $T$ be a tree, a vertex of degree one is called a leaf. The set of leaves of $T$ is denoted by $Leaf(T)$. The subtree $T-Leaf(T)$ of $T$ is called the stem of $T$ and denoted by $Stem(T).$ In this note, we give a sharp sufficient condition to show that a $K_{1,t}-$free graph has a spanning tree whose stem has a few leaves. By applying the main result, we give improvements of previous related results.


title: Optimal hedging under fast-varying stochastic volatility

abstract: In a market with a rough or Markovian mean-reverting stochastic volatility there is no perfect hedge. Here it is shown how various delta-type hedging strategies perform and can be evaluated in such markets in the case of European options. A precise characterization of the hedging cost, the replication cost caused by the volatility fluctuations, is presented in an asymptotic regime of rapid mean reversion for the volatility fluctuations. The optimal dynamic asset based hedging strategy in the considered regime is identified as the so-called `practitioners' delta hedging scheme. It is moreover shown that the performances of the delta-type hedging schemes are essentially independent of the regularity of the volatility paths in the considered regime and that the hedging costs are related to a vega risk martingale whose magnitude is proportional to a new market risk parameter. It is also shown via numerical simulations that the proposed hedging schemes which derive from option price approximations in the regime of rapid mean reversion, are robust: the `practitioners' delta hedging scheme that is identified as being optimal by our asymptotic analysis when the mean reversion time is small seems to be optimal with arbitrary mean reversion times.


title: Multi-Domain Pose Network for Multi-Person Pose Estimation and Tracking

abstract: Multi-person human pose estimation and tracking in the wild is important and challenging. For training a powerful model, large-scale training data are crucial. While there are several datasets for human pose estimation, the best practice for training on multi-dataset has not been investigated. In this paper, we present a simple network called Multi-Domain Pose Network (MDPN) to address this problem. By treating the task as multi-domain learning, our methods can learn a better representation for pose prediction. Together with prediction heads fine-tuning and multi-branch combination, it shows significant improvement over baselines and achieves the best performance on PoseTrack ECCV 2018 Challenge without additional datasets other than MPII and COCO.


title: Quality Assessment for Tone-Mapped HDR Images Using Multi-Scale and Multi-Layer Information

abstract: Tone mapping operators and multi-exposure fusion methods allow us to enjoy the informative contents of high dynamic range (HDR) images with standard dynamic range devices, but also introduce distortions into HDR contents. Therefore methods are needed to evaluate tone-mapped image quality. Due to the complexity of possible distortions in a tone-mapped image, information from different scales and different levels should be considered when predicting tone-mapped image quality. So we propose a new no-reference method of tone-mapped image quality assessment based on multi-scale and multi-layer features that are extracted from a pre-trained deep convolutional neural network model. After being aggregated, the extracted features are mapped to quality predictions by regression. The proposed method is tested on the largest public database for TMIQA and compared to existing no-reference methods. The experimental results show that the proposed method achieves better performance.


title: The Morphology of H{\alpha} emission in CALIFA galaxies

abstract: We have determined the H{\alpha} emission line radial profiles for a sample of 86 face-on galaxies observed in the CALIFA survey and analyzed with the Pipe3D pipeline. From a visual analysis we propose a two step classification of these profiles. Initially, they were divided in two classes with respect to the maximum of the H{\alpha} emission: C (for central) profiles have the maximum line emission at the galaxy center, whereas EX (for extended) profiles have the maximum H{\alpha} emission outside the galaxy center. After, we divided the C galaxies in two classes, CE and CL (where E and L stands here for 'early' and 'late'), through the value of $c_r$, the concentration index in the r-band. We analyzed the profile class dependence of several galaxy parameters, as well as examined the nature of line emission through the BPT diagram. We notice that almost 75% of the sample is in the C class. Elliptical and S0 galaxies dominate the CE class, with spiral galaxies being found mostly in the CL and EX classes. We also notice that spirals in each of these classes have different properties, with CL objects seeming less evolved than those in the EX class.


title: Quasi-free-standing monolayer hexagonal boron nitride on Ni

abstract: The electronic structure of monolayer hexagonal boron nitride grown on Ni by the diffusion and precipitation method was studied by x-ray absorption spectroscopy, emission spectroscopy, x-ray photoelectron spectroscopy and micro-ultraviolet photoemission spectroscopy. No indication of hybridization between h-BN pi and Ni 3d orbitals was observed. That is, the monolayer h-BN was found to be in the quasi-free-standing state. These results are in striking contrast to those of previous studies in which h-BN was strongly bound to the Ni surface by the orbital hybridization. The absence of hybridization is attributed to absence of a Ni(111) surface in this study. The lattice-matched Ni(111) surface is considered to be essential to orbital hybridization between h-BN and Ni.


title: Flow Network Models for Online Scheduling Real-time Tasks on Multiprocessors

abstract: We consider the flow network model to solve the multiprocessor real-time task scheduling problems. Using the flow network model or its generic form, linear programming (LP) formulation, for the problems is not new. However, the previous works have limitations, for example, that they are classified as offline scheduling techniques since they establish a flow network model or an LP problem considering a very long time interval. In this study, we propose how to construct the flow network model for online scheduling periodic real-time tasks on multiprocessors. Our key idea is to construct the flow network only for the active instances of tasks at the current scheduling time, while guaranteeing the existence of an optimal schedule for the future instances of the tasks. The optimal scheduling is here defined to ensure that all real-time tasks meet their deadlines when the total utilization demand of the given tasks does not exceed the total processing capacity. We then propose the flow network model-based polynomial-time scheduling algorithms. Advantageously, the flow network model allows the task workload to be collected unfairly within a certain time interval without losing the optimality. It thus leads us to designing three unfair-but-optimal scheduling algorithms on both continuous and discrete-time models. Especially, our unfair-but-optimal scheduling algorithm on a discrete-time model is, to the best of our knowledge, the first in the problem domain. We experimentally demonstrate that it significantly alleviates the scheduling overheads, i.e., the reduced number of preemptions with the comparable number of task migrations across processors.


title: Reply to "Comment on 'Can accretion disk properties observationally distinguish black holes from naked singularities?'"

abstract: In the Comment on "Can accretion disk properties observationally distinguish black holes from naked singularities?", by Bertrand Chauvineau, Phys. Rev. D {\bf 98}, 088501 (2018), the author did show that the metric used in Z. Kov\'{a}cs and T. Harko, Phys. Rev. D {\bf 82}, 124047 (2010), and initially introduced in K. D. Krori and D. R. Bhattacharjee, J. Math. Phys. \textbf{23}, 637 (1982) and K. K. Nandi, P. M. Alsing, J. C. Evans, and T. B. Nayak, Phys. Rev. D \textbf{63}, 084027 (2001), does not satisfy the Einstein gravitational field equations with a minimally coupled scalar field. In our reply we would like to point out that this result is actually not new, but it was already published in the literature. Moreover, a rotating solution that generalizes the Kerr metric for a nonminimally coupled scalar field does exist. We briefly discuss the nature of the singularities for the generalized metric, and point out that it can be used as a testing ground to differentiate black holes from naked singularities. We also mention the existence of some other typing or technical errors existing in the literature.


title: Time-dependent density functional theory for interaction of ultrashort light pulse with thin materials

abstract: We present a comprehensive theoretical description for an irradiation of an ultrashort light pulse normally on thin materials based on first-principles time-dependent density functional theory. As the most elaborate scheme, we develop a microscopic description solving Maxwell equations for light electromagnetic fields and the time-dependent Kohn-Sham equation for electron dynamics simultaneously in the time domain using a common spatial grid. We call it the microscopic Maxwell-TDDFT scheme. We test this scheme for silicon thin films of various thickness, from a few atomic layers to a few tens of nm. We show that the microscopic Maxwell-TDDFT scheme provides a satisfactory description incorporating the electronic structure of thin films in the first-principles level, multiple reflections of the electromagnetic fields at the surfaces, and nonlinear light-matter interaction when the incident light pulse is strong. However, the calculation becomes expensive as the thickness increases. We then consider two limiting cases of extremely thin and sufficiently thick films and develop approximate schemes. For the extremely thin case including two-dimensional atomic-layered materials, a two-dimensional macroscopic electromagnetism is developed: a two-dimensional susceptibility is introduced for a weak field, while time evolution equation is derived for an intense field. For sufficiently thick films, the microscopic Maxwell-TDDFT scheme is expected to coincide with a description utilizing ordinary macroscopic electromagnetism. We numerically confirm it comparing the calculated results: For a weak field, a comparison is made with a description using the bulk dielectric susceptibility. For a strong field, a comparison is made with a multiscale Maxwell-TDDFT scheme which the authors' group developed previously.


title: A Matrix Chernoff Bound for Strongly Rayleigh Distributions and Spectral Sparsifiers from a few Random Spanning Trees

abstract: Strongly Rayleigh distributions are a class of negatively dependent distributions of binary-valued random variables [Borcea, Branden, Liggett JAMS 09]. Recently, these distributions have played a crucial role in the analysis of algorithms for fundamental graph problems, e.g. Traveling Salesman Problem [Gharan, Saberi, Singh FOCS 11]. We prove a new matrix Chernoff bound for Strongly Rayleigh distributions. As an immediate application, we show that adding together the Laplacians of $\epsilon^{-2} \log^2 n$ random spanning trees gives an $(1\pm \epsilon)$ spectral sparsifiers of graph Laplacians with high probability. Thus, we positively answer an open question posed in [Baston, Spielman, Srivastava, Teng JACM 13]. Our number of spanning trees for spectral sparsifier matches the number of spanning trees required to obtain a cut sparsifier in [Fung, Hariharan, Harvey, Panigraphi STOC 11]. The previous best result was by naively applying a classical matrix Chernoff bound which requires $\epsilon^{-2} n \log n$ spanning trees. For the tree averaging procedure to agree with the original graph Laplacian in expectation, each edge of the tree should be reweighted by the inverse of the edge leverage score in the original graph. We also show that when using this reweighting of the edges, the Laplacian of single random tree is bounded above in the PSD order by the original graph Laplacian times a factor $\log n$ with high probability, i.e. $L_T \preceq O(\log n) L_G$. We show a lower bound that almost matches our last result, namely that in some graphs, with high probability, the random spanning tree is $\it{not}$ bounded above in the spectral order by $\frac{\log n}{\log\log n}$ times the original graph Laplacian. We also show a lower bound that in $\epsilon^{-2} \log n$ spanning trees are necessary to get a $(1\pm \epsilon)$ spectral sparsifier.


title: On the lower bounds of Davenport constant

abstract: Let $G = C_{n_1} \oplus \cdots \oplus C_{n_r}$ with $1 < n_1 | \cdots | n_r$ be a finite abelian group. The Davenport constant $\mathsf D(G)$ is the smallest integer $t$ such that every sequence $S$ over $G$ of length $|S|\ge t$ has a non-empty zero-sum subsequence. It is a starting point of zero-sum theory but only has a trivial lower bound $\mathsf D^*(G) = n_1 + \cdots + n_r - r + 1$, which equals $\mathsf D(G)$ over $p$-groups. We investigate the non-dispersive sequences over group $C_n^r$, thereby revealing the growth of $\mathsf D(G)-\mathsf D^*(G)$ over non-$p$-groups $G = C_n^r \oplus C_{kn}$ with $n,k \ne 1$. We give a general lower bound of $\mathsf D(G)$ over non-$p$-groups and show that, let $G$ be abelian groups with $\exp(G)=m$ and rank $r$, fix $m>0$ a non-prime-power, then for each $N>0$ there exists an $\varepsilon>0$ such that if $|G|/m^r<\varepsilon $, then $\mathsf D(G)-\mathsf D^*(G)>N$.


title: Excitations of a vortex line in an elongated dipolar condensate

abstract: We characterise the properties of a vortex line in an elongated dipolar Bose-Einstein condensate. Increasing the strength of the dipole-dipole interactions (DDIs) relative to the short ranged contact interactions we find that the system crosses over to a self-bound vortex droplet stabilized from collapse by quantum fluctuations. We calculate the quasiparticle excitation spectrum of the vortex state, which is important in characterizing the vortex response, and assessing its stability. When the DDIs are sufficiently strong we find that the vortex is dynamically unstable to quadrupolar modes.


title: Harmonic maps in connection of phase transitions with higher dimensional potential wells

abstract: This is in the sequel of authors' paper \cite{LPW} in which we had set up a program to verify rigorously some formal statements associated with the multiple component phase transitions with higher dimensional wells. The main goal here is to establish a regularity theory for minimizing maps with a rather non-standard boundary condition at the sharp interface of the transition. We also present a proof, under simplified geometric assumptions, of existence of local smooth gradient flows under such constraints on interfaces which are in the motion by the mean-curvature. In a forthcoming paper, a general theory for such gradient flows and its relation to Keller-Rubinstein-Sternberg's work \cite{KRS1, KRS2} on the fast reaction, slow diffusion and motion by the mean curvature would be addressed.