title: Quantum Chiral Higher Spin Gravity

abstract: An example of a higher spin gravity in four-dimensional flat space has recently been constructed in arXiv:1609.04655 [hep-th]. This theory is chiral and the action is written in the light-cone gauge. The theory has certain stringy features, e.g. admits Chan-Paton factors. We show that the theory is consistent, both at the classical and quantum level. Even though the interactions are non-trivial, due to the coupling conspiracy all tree level amplitudes vanish on-shell. The loop corrections also vanish. Therefore, the full quantum S-matrix is one and the theory is consistent with the numerous no-go theorems. This provides the first example of a (quantum) interacting higher spin gravity with an action. We argue that higher spin gravities in AdS space should display the same features.


title: An ice giant exoplanet interpretation of the anomaly in microlensing event OGLE-2011-BLG-0173

abstract: We analyze the microlensing event OGLE-2011-BLG-0173, which shows a small perturbation at the end of the microlensing event caused by the primary lens. We consider both binary lens and binary source models and we explore their degeneracies, some of which have not previously been recognized. There are two families of binary lens solutions, one with a mass ratio $q\approx4\times10^{-4}$ and a separation s~4.6 and the other with q~0.015 and s~0.22, i.e, both have companions in the planetary regime. We search for solutions by using Bayesian analysis that includes planet frequency as a prior and find that the s~4.6 family is the preferred one with ~4 M_Uranus mass planet on an orbit of ~10 AU. The degeneracies arise from a paucity of information on the anomaly, demonstrating that high-cadence observations are essential for characterizing wide-orbit microlensing planets. Hence, we predict that the planned WFIRST microlensing survey will be less prone to these degeneracies than the ongoing ground-based surveys. We discuss the known low-mass, wide-orbit companions and we notice that for the largest projected separations the mass ratios are either high (consistent with brown dwarf companions) or low (consistent with Uranus analogs), but intermediate mass ratios (Jupiter analogs on wide orbits) have not been detected to date, despite the fact that the sensitivity to such planets should be higher than that of Uranus analogs. This is therefore tentative evidence of the existence of a massive ice giant desert at wide separations. On the other hand, given their low intrinsic detection sensitivity, Uranus analogs may be ubiquitous.


title: Diagnosing the Structure of Massive Stars with Galactic Supernova Neutrinos

abstract: It has been suggested that whether a star explodes or not, and what kind of explosion properties it shows, is strongly dependent on the progenitor's core structure. We present the results from 101 axisymmetric core-collapse supernova simulations performed with progenitors spanning initial masses in the range from 10.8 to 75 solar masses, and focus on their connections to the compactness of the progenitor's core. Our simulations confirm a correlation between the neutrinos emitted during the accretion phase and the progenitor's compactness. We suggest that the ratio of observed neutrino events during the first hundreds of milliseconds can be used to infer the progenitor's inner mass density structure.


title: Critically Slow Learning in Flashcard Learning Models

abstract: Algorithmic education theory examines, among other things, the trade-off between reviewing old material and studying new material: time spent learning the new comes at the expense of reviewing and solidifying one's understanding of the old. This trade-off is captured in the `Slow Flashcard System' (SFS) -- a system that has been studied not only for its applications in educational software but also for its critical properties; it is a simple discrete deterministic system capable of remarkable complexity, with standing conjectures regarding its longterm behavior. Here we introduce a probabilistic model of SFS and further derive a continuous time, continuous space PDE model. These two models of SFS shed light on the longterm behavior of SFS and open new avenues of research.


title: Lagrangian Structure for Compressible Flow in the Half-space with the Navier Boundary Condition

abstract: We show the uniqueness of particle paths of a velocity field, which solves the compressible isentropic Navier-Stokes equations in the half-space $\mathbb{R}_+^3$ with the Navier boundary condition. In particular, by means of energy estimates and the assumption of small energy we prove that the velocity field satisfies the necessary regularity needed to prove the uniqueness of particle paths.


title: First-passage distributions for the one-dimensional Fokker-Planck equation

abstract: We present an analytical framework to study the first-passage (FP) and first-return (FR) distributions for the broad family of models described by the one-dimensional Fokker-Planck equation in finite domains, identifying general properties of these distributions for different classes of models. When in the Fokker-Planck equation the diffusion coefficient is positive (nonzero) and the drift term is bounded, as in the case of a Brownian walker, both distributions may exhibit a power-law decay with exponent -3/2 for intermediate times. We discuss how the influence of an absorbing state changes this exponent. The absorbing state is characterized by a vanishing diffusion coefficient and/or a diverging drift term. Remarkably, the exponent of the Brownian walker class of models is still found, as long as the departure and arrival regions are far enough from the absorbing state, but the range of times where the power law is observed narrows. Close enough to the absorbing point, though, a new exponent may appear. The particular value of the exponent depends on the behavior of the diffusion and the drift terms of the Fokker-Planck equation. We focus on the case of a diffusion term vanishing linearly at the absorbing point. In this case, the FP and FR distributions are similar to those of the voter model, characterized by a power law with exponent -2. As an illustration of the general theory, we compare it with exact analytical solutions and extensive numerical simulations of a two-parameter voter-like family models. We study the behavior of the FP and FR distributions by tuning the importance of the absorbing points throughout changes of the parameters. Finally, the possibility of inferring relevant information about the steady-sate probability distribution of a model from the FP and FR distributions is addressed.


title: Benchmarking the Capabilities and Limitations of SAT Solvers in Defeating Obfuscation Schemes

abstract: In this paper, we investigate the strength of six different SAT solvers in attacking various obfuscation schemes. Our investigation revealed that Glucose and Lingeling SAT solvers are generally suited for attacking small-to-midsize obfuscated circuits, while the MapleGlucose, if the system is not memory bound, is best suited for attacking mid-to-difficult obfuscation methods. Our experimental result indicates that when dealing with extremely large circuits and very difficult obfuscation problems, the SAT solver may be memory bound, and Lingeling, for having the most memory efficient implementation, is the best-suited solver for such problems. Additionally, our investigation revealed that SAT solver execution times may vary widely across different SAT solvers. Hence, when testing the hardness of an obfuscation method, although the increase in difficulty could be verified by one SAT solver, the pace of increase in difficulty is dependent on the choice of a SAT solver.


title: Efficient numerical simulations with Tensor Networks: Tensor Network Python (TeNPy)

abstract: Tensor product state (TPS) based methods are powerful tools to efficiently simulate quantum many-body systems in and out of equilibrium. In particular, the one-dimensional matrix-product (MPS) formalism is by now an established tool in condensed matter theory and quantum chemistry. In these lecture notes, we combine a compact review of basic TPS concepts with the introduction of a versatile tensor library for Python (TeNPy) [https://github.com/tenpy/tenpy]. As concrete examples, we consider the MPS based time-evolving block decimation and the density matrix renormalization group algorithm. Moreover, we provide a practical guide on how to implement abelian symmetries (e.g., a particle number conservation) to accelerate tensor operations.


title: Many-body delocalization in the presence of a quantum bath

abstract: Closed generic quantum many-body systems may fail to thermalize under certain conditions even after long times, a phenomenon called many-body localization (MBL). Numerous studies support the stability of the MBL phase in strongly disordered one-dimensional systems. However, the situation is much less clear when a small part of the system is ergodic, a scenario which also has important implications for the existence of many-body localization in higher dimensions. Here we address this question experimentally using a large-scale quantum simulator of ultracold bosons in a two-dimensional optical lattice. We prepare two-component mixtures of varying relative population and implement a disorder potential which is only experienced by one of the components. The second non-disordered ''clean'' component plays the role of a bath of adjustable size that is collisionally coupled to the ''dirty'' component. Our experiments show how the dynamics of the dirty component, which, when on its own, show strong evidence of localization, become affected by the coupling to the clean component. For a high clean population, the clean component appears to behave as an effective bath for the system which leads to its delocalization, while for a smaller clean population, the ability of the bath to destabilize the system becomes strongly reduced. Our results reveal how a finite-sized quantum system can bring another one towards thermalization, in a regime of complex interplay between disorder, tunneling and intercomponent interactions. They provide a new benchmark for effective theories aiming to capture the complex physics of MBL in the weakly localized regime.

categories: cond-mat.quant-gas cond-mat.dis-nn cond-mat.stat-mech cond-mat.str-el


title: Identifying Effects of Multivalued Treatments

abstract: Multivalued treatment models have typically been studied under restrictive assumptions: ordered choice, and more recently unordered monotonicity. We show how treatment effects can be identified in a more general class of models that allows for multidimensional unobserved heterogeneity. Our results rely on two main assumptions: treatment assignment must be a measurable function of threshold-crossing rules, and enough continuous instruments must be available. We illustrate our approach for several classes of models.


title: Superpositions of bright and dark solitons supporting the creation of balanced gain and loss optical potentials

abstract: Bright and dark solitons of the cubic nonlinear Schrodinger equation are used to construct complex-valued potentials with all-real spectrum. The real part of these potentials is equal to the intensity of a bright soliton while their imaginary part is defined by the product of such soliton with its concomitant, a dark soliton. Considering light propagation in Kerr media, the real part of the potential refers to the self-focusing of the signal and the imaginary one provides the system with balanced gain-and-loss rates.


title: Universality in Freezing Cellular Automata

abstract: Cellular Automata have been used since their introduction as a discrete tool of modelization. In many of the physical processes one may modelize thus (such as bootstrap percolation, forest fire or epidemic propagation models, life without death, etc), each local change is irreversible. The class of freezing Cellular Automata (FCA) captures this feature. In a freezing cellular automaton the states are ordered and the cells can only decrease their state according to this "freezing-order". We investigate the dynamics of such systems through the questions of simulation and universality in this class: is there a Freezing Cellular Automaton (FCA) that is able to simulate any Freezing Cellular Automata, i.e. an intrinsically universal FCA? We show that the answer to that question is sensitive to both the number of changes cells are allowed to make, and geometric features of the space. In dimension 1, there is no universal FCA. In dimension 2, if either the number of changes is at least 2, or the neighborhood is Moore, then there are universal FCA. On the other hand, there is no universal FCA with one change and Von Neumann neighborhood. We also show that monotonicity of the local rule with respect to the freezing-order (a common feature of bootstrap percolation) is also an obstacle to universality.


title: On improving the approximation ratio of the r-shortest common superstring problem

abstract: The Shortest Common Superstring problem (SCS) consists, for a set of strings S = {s_1,...,s_n}, in finding a minimum length string that contains all s_i, 1<= i <= n, as substrings. While a 2+11/30 approximation ratio algorithm has recently been published, the general objective is now to break the conceptual lower bound barrier of 2. This paper is a step ahead in this direction. Here we focus on a particular instance of the SCS problem, meaning the r-SCS problem, which requires all input strings to be of the same length, r. Golonev et al. proved an approximation ratio which is better than the general one for r<= 6. Here we extend their approach and improve their approximation ratio, which is now better than the general one for r<= 7, and less than or equal to 2 up to r = 6.


title: Machine Learning for Predictive On-Demand Deployment of UAVs for Wireless Communications

abstract: In this paper, a novel machine learning (ML) framework is proposed for enabling a predictive, efficient deployment of unmanned aerial vehicles (UAVs), acting as aerial base stations (BSs), to provide on-demand wireless service to cellular users. In order to have a comprehensive analysis of cellular traffic, an ML framework based on a Gaussian mixture model (GMM) and a weighted expectation maximization (WEM) algorithm is introduced to predict the potential network congestion. Then, the optimal deployment of UAVs is studied to minimize the transmit power needed to satisfy the communication demand of users in the downlink, while also minimizing the power needed for UAV mobility, based on the predicted cellular traffic. To this end, first, the optimal partition of service areas of each UAV is derived, based on a fairness principle. Next, the optimal location of each UAV that minimizes the total power consumption is derived. Simulation results show that the proposed ML approach can reduce the required downlink transmit power and improve the power efficiency by over 20%, compared with an optimal deployment of UAVs with no ML prediction.


title: Whispering gallery microsensors: a review

abstract: Optical whispering gallery mode (WGM) microresonators, confining resonant photons in a microscale resonator for long periods of time, could strongly enhance light-matter interaction, making it an ideal platform for all kinds of sensors. In this paper, an overview of optical sensors based on WGM microresonators is comprehensively summarized. First, the fundamental sensing mechanisms as well as several recently developed enhanced sensing techniques are introduced. Then, different types of WGM structures for sensing as well as microfluidics techniques are summarized. Furthermore, several important sensing parameters are discussed. Most importantly, a variety of WGM sensing applications are reviewed, including both traditional matter sensing and field sensing. Last, we give a brief summary and perspective of the WGM sensors and their role in future applications.


title: Adversarial Semantic Alignment for Improved Image Captions

abstract: In this paper we study image captioning as a conditional GAN training, proposing both a context-aware LSTM captioner and co-attentive discriminator, which enforces semantic alignment between images and captions. We empirically focus on the viability of two training methods: Self-critical Sequence Training (SCST) and Gumbel Straight-Through (ST) and demonstrate that SCST shows more stable gradient behavior and improved results over Gumbel ST, even without accessing discriminator gradients directly. We also address the problem of automatic evaluation for captioning models and introduce a new semantic score, and show its correlation to human judgement. As an evaluation paradigm, we argue that an important criterion for a captioner is the ability to generalize to compositions of objects that do not usually co-occur together. To this end, we introduce a small captioned Out of Context (OOC) test set. The OOC set, combined with our semantic score, are the proposed new diagnosis tools for the captioning community. When evaluated on OOC and MS-COCO benchmarks, we show that SCST-based training has a strong performance in both semantic score and human evaluation, promising to be a valuable new approach for efficient discrete GAN training.


title: Static friction boost in edge-driven incommensurate contacts

abstract: We present a numerical investigation of the size scaling of static friction in incommensurate two dimensional contacts performed for different lateral loading configurations. Results of model simulations show that both the absolute value of the force $F_s$ and the scaling exponent $\gamma$ strongly depend on the loading configuration adopted to drive the slider along the substrate. Under edge-loading a sharp increase of static friction is observed above a critical size corresponding to the appearance of a localized commensurate dislocation. Noticeably, the existence of sublinear scaling, which is a fingerprint of superlubricity, does not conflict with the possibility to observe shear-induced localized commensurate regions at the contact interface. Atomistic simulations of gold islands sliding over graphite corroborate these findings suggesting that similar elasticity effects should be at play in real frictional contacts.


title: A General Framework for Counterfactual Learning-to-Rank

abstract: Implicit feedback (e.g., click, dwell time) is an attractive source of training data for Learning-to-Rank, but its naive use leads to learning results that are distorted by presentation bias. For the special case of optimizing average rank for linear ranking functions, however, the recently developed SVM-PropRank method has shown that counterfactual inference techniques can be used to provably overcome the distorting effect of presentation bias. Going beyond this special case, this paper provides a general and theoretically rigorous framework for counterfactual learning-to-rank that enables unbiased training for a broad class of additive ranking metrics (e.g., Discounted Cumulative Gain (DCG)) as well as a broad class of models (e.g., deep networks). Specifically, we derive a relaxation for propensity-weighted rank-based metrics which is subdifferentiable and thus suitable for gradient-based optimization. We demonstrate the effectiveness of this general approach by instantiating two new learning methods. One is a new type of unbiased SVM that optimizes DCG -- called SVM PropDCG --, and we show how the resulting optimization problem can be solved via the Convex Concave Procedure (CCP). The other is Deep PropDCG, where the ranking function can be an arbitrary deep network. In addition to the theoretical support, we empirically find that SVM PropDCG significantly outperforms existing linear rankers in terms of DCG. Moreover, the ability to train non-linear ranking functions via Deep PropDCG further improves performance.


title: Modular Compact Modeling of Magnetic Tunnel Junction Devices

abstract: This paper describes a robust, modular, and physics- based circuit framework to model conventional and emerging Magnetic Tunnel Junction (MTJ) devices. Magnetization dynamics are described by the stochastic Landau-Lifshitz-Gilbert (sLLG) equation whose results are rigorously benchmarked with a Fokker-Planck Equation (FPE) description of magnet dynamics. We then show how sLLG is coupled to transport equations of MTJ-based devices in a unified circuit platform. Step by step, we illustrate how the physics-based MTJ model can be extended to include different spintronics phenomena, including spin-transfer-torque (STT), voltage-control of magnetic anisotropy (VCMA) and spin-orbit torque (SOT) phenomena by experimentally benchmarked examples. To demonstrate how our approach can be used in the exploration of novel MTJ-based devices, we also present a recently proposed MEMS resonator- driven spin-torque nano oscillator (STNO) that can reduce the phase noise of STNOs. We briefly elaborate on the use of our framework beyond conventional devices.


title: A General Framework for Relational Parametricity

abstract: Reynolds' original theory of relational parametricity was intended to capture the idea that polymorphically typed System F programs preserve all relations between inputs. But as Reynolds himself later showed, his theory can only be formalized in a meta-theory with an impredicative universe, such as the Calculus of Inductive Constructions. Abstracting from Reynolds' ideas, Dunphy and Reddy developed their well-known framework for parametricity that uses parametric limits in reflexive graph categories and aims to subsume a variety of parametric models. As we observe, however, their theory is not sufficiently general to subsume the very model that inspired parametricity, namely Reynolds' original model, expressed inside type theory. To correct this, we develop an abstract framework for relational parametricity that generalizes the notion of a reflexive graph categories and delivers Reynolds' model as a direct instance in a natural way. This framework is uniform with respect to a choice of meta-theory, which allows us to obtain the well-known PER model of Longo and Moggi as a direct instance in a natural way as well. In addition, we offer two novel relationally parametric models of System F: i) a categorical version of Reynolds' model, where types are functorial on isomorphisms and all polymorphic functions respect the functorial action, and ii) a proof-relevant categorical version of Reynolds' model (after Orsanigo), where, additionally, witnesses of relatedness are themselves suitably related. We show that, unlike previously existing frameworks for parametricity, ours recognizes both of these new models in a natural way. Our framework is thus descriptive, in that it accounts for well-known models, as well as prescriptive, in that it identifies abstract properties that good models of relational parametricity should satisfy and suggests new constructions of such models.


title: Exploiting Answer Set Programming with External Sources for Meta-Interpretive Learning

abstract: Meta-Interpretive Learning (MIL) learns logic programs from examples by instantiating meta-rules, which is implemented by the Metagol system based on Prolog. Viewing MIL-problems as combinatorial search problems, they can alternatively be solved by employing Answer Set Programming (ASP), which may result in performance gains as a result of efficient conflict propagation. However, a straightforward ASP-encoding of MIL results in a huge search space due to a lack of procedural bias and the need for grounding. To address these challenging issues, we encode MIL in the HEX-formalism, which is an extension of ASP that allows us to outsource the background knowledge, and we restrict the search space to compensate for a procedural bias in ASP. This way, the import of constants from the background knowledge can for a given type of meta-rules be limited to relevant ones. Moreover, by abstracting from term manipulations in the encoding and by exploiting the HEX interface mechanism, the import of such constants can be entirely avoided in order to mitigate the grounding bottleneck. An experimental evaluation shows promising results.


title: Ultralow-power chip-based soliton microcombs for photonic integration

abstract: The generation of dissipative Kerr solitons in optical microresonators has provided a route to compact frequency combs of high repetition rate, which have already been employed for optical frequency synthesizers, ultrafast ranging, coherent telecommunication and dual-comb spectroscopy. Silicon nitride (Si$_3$N$_4$) microresonators are promising for photonic integrated soliton microcombs. Yet to date, soliton formation in Si$_3$N$_4$ microresonators at electronically detectable repetition rates, typically less than 100 GHz, is hindered by the requirement of external power amplifiers, due to the low quality ($Q$) factors, as well as by thermal effects which necessitate the use of frequency agile lasers to access the soliton state. These requirements complicate future photonic integration, heterogeneous or hybrid, of soliton microcomb devices based on Si$_3$N$_4$ microresonators with other active or passive components. Here, using the photonic Damascene reflow process, we demonstrate ultralow-power single soliton formation in high-Q ($Q_0>15\times10^6$) Si$_3$N$_4$ microresonators with 9.8 mW input power (6.2 mW in the waveguide) for devices of electronically detectable, 99 GHz repetition rate. We show that solitons can be accessed via simple, slow laser piezo tuning, in many resonances in the same sample. These power levels are compatible with current silicon-photonics-based lasers for full photonic integration of soliton microcombs, at repetition rates suitable for applications such as ultrafast ranging and coherent communication. Our results show the technological readiness of Si$_3$N$_4$ optical waveguides for future all-on-chip soliton microcomb devices.


title: Probing heavy neutrino oscillations in rare W boson decays

abstract: In this work, we study the lepton number violating W boson and top quark decays via intermediate on-shell Majorana neutrinos Nj into three charged leptons and a light neutrino. We discuss the neutrino oscillation effects present in the decay due to the small mass gap between the heavy neutrino states. We focus on a scenario that contains at least two heavy Majorana neutrinos in the mass range 2-80 GeV. The results indicate that the modulation of the branching ratios as a function of the distance between the vertices may be detected in a future experiment such as High-Luminosity Large Hadron Collider. As a secondary result, the CP-violating phases could be explored.


title: Understanding Regularization to Visualize Convolutional Neural Networks

abstract: Variational methods for revealing visual concepts learned by convolutional neural networks have gained significant attention during the last years. Being based on noisy gradients obtained via back-propagation such methods require the application of regularization strategies. We present a mathematical framework unifying previously employed regularization methods. Within this framework, we propose a novel technique based on Sobolev gradients which can be implemented via convolutions and does not require specialized numerical treatment, such as total variation regularization. The experiments performed on feature inversion and activation maximization demonstrate the benefit of a unified approach to regularization, such as sharper reconstructions via the proposed Sobolev filters and a better control over reconstructed scales.


title: On p-adic Multidimensional Continued Fractions

abstract: Multidimensional continued fractions (MCFs) were introduced by Jacobi and Perron in order to generalize the classical continued fractions. In this paper, we propose an introductive fundamental study about MCFs in the field of the $p$--adic numbers $\mathbb Q_p$. First, we introduce them from a formal point of view, i.e., without considering a specific algorithm that produces the partial quotients of a MCF, and we perform a general study about their convergence in $\mathbb Q_p$. In particular, we derive some conditions about their convergence and we prove that convergent MCFs always strongly converge in $\mathbb Q_p$ contrarily to the real case where strong convergence is not ever guaranteed. Then, we focus on a specific algorithm that, starting from a $m$--tuple of numbers in $\mathbb Q_p$, produces the partial quotients of the corresponding MCF. We see that this algorithm is derived from a generalized $p$--adic Euclidean algorithm and we prove that it always terminates in a finite number of steps when it processes rational numbers.


title: Photoionization of Xe and Rn from the relativistic random-phase theory

abstract: Photoionization cross section $\sigma_{n\kappa}$, asymmetry parameter $\beta_{n\kappa}$, and polarization parameters $\xi_{n\kappa}$, $\eta_{n\kappa}$, $\zeta_{n\kappa}$ of Xe and Rn are calculated in the fully relativistic formalism. To deal with the relativistic and correlation effects, we adopt the relativistic random-phase theory with channel couplings among different subshells. Energy ranges for giant \emph{d}-resonance regions are especially considered.


title: Checking is Believing: Event-Aware Program Anomaly Detection in Cyber-Physical Systems

abstract: Securing cyber-physical systems (CPS) against malicious attacks is of paramount importance because these attacks may cause irreparable damages to physical systems. Recent studies have revealed that control programs running on CPS devices suffer from both control-oriented attacks (e.g., code-injection or code-reuse attacks) and data-oriented attacks (e.g., non-control data attacks). Unfortunately, existing detection mechanisms are insufficient to detect runtime data-oriented exploits, due to the lack of runtime execution semantics checking. In this work, we propose Orpheus, a new security methodology for defending against data-oriented attacks by enforcing cyber-physical execution semantics. We first present a general method for reasoning cyber-physical execution semantics of a control program (i.e., causal dependencies between the physical context and program control flows), including the event identification and dependence analysis. As an instantiation of Orpheus, we then present a new program behavior model, i.e., the event-aware finite-state automaton (eFSA). eFSA takes advantage of the event-driven nature of CPS control programs and incorporates event checking in anomaly detection. It detects data-oriented exploits if a specific physical event is missing along with the corresponding event dependent state transition. We evaluate our prototype's performance by conducting case studies under data-oriented attacks. Results show that eFSA can successfully detect different runtime attacks. Our prototype on Raspberry Pi incurs a low overhead, taking 0.0001s for each state transition integrity checking, and 0.063s~0.211s for the cyber-physical contextual consistency checking.


title: Greedy approximations by signed harmonic sums and the Thue--Morse sequence

abstract: Given a real number $\tau$, we study the approximation of $\tau$ by signed harmonic sums $\sigma_N(\tau) := \sum_{n \leq N}{s_n(\tau)}/n$, where the sequence of signs $(s_N(\tau))_{N \in\mathbb{N}}$ is defined "greedily" by setting $s_{N+1}(\tau) := +1$ if $\sigma_N(\tau) \leq \tau$, and $s_{N+1}(\tau) := -1$ otherwise. Precisely, we compute the limit points and the decay rate of the sequence $(\sigma_N(\tau)-\tau)_{N \in \mathbb{N}}$. Moreover, we give an accurate description of the behavior of the sequence of signs $(s_N(\tau))_{N\in\mathbb{N}}$, highlighting a surprising connection with the Thue--Morse sequence.


title: Several Topics in Experimental Mathematics

abstract: This thesis deals with applications of experimental mathematics to a number of problems. The first problem is related to random graph statistics. We consider a certain class of Galton-Watson random trees and look at the total height statistic. We provide an automated procedure for computing values of the moments of this statistic. Next, we investigate several problems related to lattice paths staying below a line of rational slope. These results are largely data-based. Using the generated data, we are able to find recurrences for the number of such paths for the cases of slopes 3/2 and 5/2. There is also investigation of a generalization of these problems to three dimensions. We also examine generalizations of Sister Celine's method and Gosper's algorithm for evaluating summations. For both, we greatly extend the classes of applicable functions and applications to proving, or reproving in an automated way, interesting combinatorial problems. For the generalization of Sister Celine's method, we allow summations of arbitrary products of hypergeometric terms and linear recurrent sequences with rational coefficients. We also show a partial result related to the bunk bed conjecture, a problem concerning random finite graphs. Let $G$ be a finite graph. Remove edges from $G\square K_2$ independently and with the same probability. In $G\square K_2$, there is an edge placed between all vertices of $G$ and the corresponding vertex in a copy of $G$. Then, label these vertices as either $(v,0)$ or $(v,1)$ for each $v\in V(G)$. The conjecture says that for any $x,y \in V(G)$, it is least as likely to have $(x,0)$ connected to $(y,0)$ as to have $(x,0)$ connected to $(y,1)$. We prove the conjecture in the case that only two of the edges going between the two copes of $G$ are retained.


title: Linear dynamics in reproducing kernel Hilbert spaces

abstract: Complementing earlier results on dynamics of unilateral weighted shifts, we obtain a sufficient (but not necessary, with supporting examples) condition for hypercyclicity, mixing and chaos for $M_z^*$, the adjoint of $M_z$, on vector-valued analytic reproducing kernel Hilbert spaces $\mathcal{H}$ in terms of the derivatives of kernel functions on the open unit disc $\mathbb{D}$ in $\mathbb{C}$. Here $M_z$ denotes the multiplication operator by the coordinate function $z$, that is \[ (M_z f) (w) = w f(w), \] for all $f \in \mathcal{H}$ and $w \in \mathbb{D}$. We analyze the special case of quasi-scalar reproducing kernel Hilbert spaces. We also present a complete characterization of hypercyclicity of $M_z^*$ on tridiagonal reproducing kernel Hilbert spaces and some special classes of vector-valued analytic reproducing kernel Hilbert spaces.


title: Properties of the entanglement Hamiltonian for finite free-fermion chains

abstract: We study the entanglement Hamiltonian for fermionic hopping models on rings and open chains and determine single-particle spectra, eigenfunctions and the form in real space. For the chain, we find a commuting operator as for the ring and compare with its properties in both cases. In particular, a scaling relation between the eigenvalues is found for large systems. We also show how the commutation property carries over to the critical transverse Ising model.


title: Coherence, entanglement and quantumness in closed and open systems with conserved charge, with an application to many-body localisation

abstract: While the scaling of entanglement in a quantum system can be used to distinguish many-body quantum phases, it is usually hard to quantify the amount of entanglement in mixed states of open quantum systems, while measuring entanglement experimentally, even for the closed systems, requires in general quantum state tomography. In this work we show how to remedy this situation in system with a fixed or conserved charge, e.g., density or magnetization, due to an emerging relation between quantum correlations and coherence. First, we show how, in these cases, the presence of multipartite entanglement or quantumness can be faithfully witnessed simply by detecting coherence in the quantum system, while bipartite entanglement or bipartite quantum discord are implied by asymmetry (block coherence) in the system. Second, we prove that the relation between quantum correlations and coherence is also quantitative. Namely, we establish upper and lower bounds on the amount of multipartite and bipartite entanglement in a many-body system with a fixed local charge, in terms of the amount of coherence and asymmetry present in the system. Importantly, both for pure and mixed quantum states, these bounds are expressed as closed formulas, and furthermore, for bipartite entanglement, are experimentally accessible by means of the multiple quantum coherence spectra. In particular, in one-dimensional systems, our bounds may detect breaking of the area law of entanglement entropy. We illustrate our results on the example of a many-body localized system, also in the presence of dephasing.

categories: cond-mat.stat-mech cond-mat.dis-nn cond-mat.quant-gas quant-ph


title: Measurement of Dispersion of a Clinical Proton Therapy Beam

abstract: Purpose: To measure dispersion of a clinical proton therapy beam. Methods and Materials: The proton center at our institution uses a Hitachi PROBEAT V synchrotron that delivers beams with energies ranging from 70 to 230 MeV to five treatment rooms, four of which have gantries and one fixed beam. For this study, the magnetic field strengths of all optics in the transport line were varied to measure the dispersion of 230 MeV beams in each of the four gantry rooms. Beam position was measured using the Spot Position Monitor (SPM), a multi-wire proportional counter (MWPC) in the nozzle approximately 0.5 meters upstream of the treatment isocenter. Measured beam positions as a function of momentum offset were fitted to a linear function to extract dispersion. In one gantry, dispersions for 70 and 140 MeV beams were also measured. Because spot position is continuously monitored during treatments and a position deviation will abort the beam delivery, the dispersion at the spot position monitor was used to set an upper limit of range uncertainty originating from momentum deviation. When dispersion was approximately zero, dose rate measurements were used to set an upper limit of momentum deviation that could reach the patient. Results: Dispersion was 0.45 m, 0.48 m and 0.30 m for 70 MeV, 140 MeV, and 230 MeV beams in a gantry, respectively, corresponding to 0.3 mm, 1.0 mm and 3.5 mm maximum range uncertainty during treatments. For 230 MeV beams, dispersion ranged from 0 to 0.73 m across the four gantries, corresponding to the upper limits of range uncertainty due to momentum deviation of 1.4 to 6.2 mm (0.4 to 1.9% in fractional range uncertainty). Conclusions: The measured dispersion was both energy- and gantry-dependent and ranged from 0 to 0.73 m.


title: Impacts of multiflagellarity on stability and speed of bacterial locomotion

abstract: Trajectories and conformations of uni- and multiflagellar bacteria are studied with a coarse-grained model of a cell comprised of elastic flagella connected to a cell body. The elasticities of both the hook protein (connecting cell body and flagellum) and flagella are varied. Flexibility plays contrasting roles for uni- and multiflagellar swimmers. For a uniflagellar swimmer, hook and/or flagellar buckling occurs above a critical flexibility relative to the torque exerted by the flagellar motor. Addition of a second flagellum greatly expands the parameter regime of stable locomotion, because flexible hooks that would lead to buckling instability in the uniflagellar case provide the flexibility required for flagellar bundling in the biflagellar case. Similar observations hold for tri- and quadriflagellar swimmers. Indeed the stability regimes for uni- and quadriflagellar swimming are virtually inverted -- to a first approximation what is stable in one case is unstable in the other. Swimming speed is also examined: it increases very weakly with number of flagella and a simple theory is developed that explains this observation.


title: Measuring uncertainty during respiratory rate estimation using pressure-sensitive mats

abstract: We develop and evaluate a respiratory rate estimation algorithm that utilizes data from pressure-sensitive mat (PSM) technology for continuous patient monitoring in neonatal intensive care units (NICU). An analysis of the random effect of drift and systematic effect of creep in the PSM data is presented, showing that these are essentially dependent on the applied load and contact surface. Uncertainty measurements are pivotal when estimating physiologic parameters. The standard uncertainty in the PSM data is here represented by the percent drift. Next, we evaluate the applicability of PSM technology to estimate RR in neonatal patient simulator trials under five mixed effects including internally and externally induced motion, mattress type, grunting, laying position, and different breathing rates. We analyze the limits of agreement on the mixed effects model to derive the uncertainty in the estimated RR obtained through two estimation techniques. In comparison with the gold standard RR values, we achieved a mean bias of 0.56 breaths per minute (bpm) with an error bounded by a 95% confidence interval of [-2.26, 3.37] bpm. These results meet the clinical accuracy requirements of RR within +/-5 bpm.


title: Comparing time and frequency domain estimation of neonatal respiratory rate using pressure-sensitive mats

abstract: Pressure-sensitive mats (PSM) have proved to be useful in the estimation of respiratory rates (RR) in adult patients. However, PSM technology has not been extensively applied to derive physiologic parameters in infant and neonatal patients. This research evaluates the applicability of the capacitive XSensor PSM technology to estimate a range of RR in neonatal patient simulator trials conducted under several experimental conditions. PSM data are analyzed in both the time and frequency domain and comparative results are presented. For the frequency-domain approach, in addition to estimating RR, a measure of confidence is also derived from the relative height of peaks in the periodogram. The study demonstrates that frequency domain analysis of mean-shifted PSM data achieves the best possible RR estimation, with zero percent error, as compared to the lowest achievable RMS error of 1.57 percent in the time domain. The frequency domain approach outperforms the time domain analysis whether examining raw data or those preprocessed by normalizing, detrending and median filtering.


title: Observation of a Topological Insulator Dirac Cone Reshaped by Non-magnetic Impurity Resonance

abstract: The massless Dirac electrons found at topological insulator surfaces are thought to be influenced very little by weak, non-magnetic disorder. However, a resonance effect of strongly perturbing non-magnetic impurities has been theoretically predicted to change the dispersion and physical nature of low-energy quasiparticles, resulting in unique particle-like states that lack microscopic translational symmetry. Here we report the direct observation of impurities reshaping the surface Dirac cone of the model 3D topological insulator Bi2Se3. For the first time, a pronounced kink-like dispersion feature is observed in disorder-enriched samples, and found to be closely associated with the anomaly caused by impurity resonance in the surface state density of states, as observed by dichroic angle resolved photoemission spectroscopy (ARPES). The experimental observation of these features, which closely resemble theoretical predictions, has significant implications for the properties of topological Dirac cones in applied scenarios that commonly feature point defect disorder at surfaces or interfaces.


title: PASTA: Python Algorithms for Searching Transition stAtes

abstract: Chemical reactions are often associated with an energy barrier along the reaction pathway which hinders the spontaneity of the reaction. Changing the energy barrier along the reaction pathway allows one to modulate the performance of a reaction. We present a module, Python Algorithms for Searching Transition stAtes (PASTA), to calculate the energy barrier and locate the transition state of a reaction efficiently. The module is written in python and can perform nudged elastic band, climbing image nudged elastic band and automated nudged elastic band calculations. These methods require the knowledge of the potential energy surface (and its gradient along some direction). This module is written such that it works in conjunction with density functional theory (DFT) codes to obtain this information. Presently it is interfaced with three well known DFT packages: Vienna Ab initio Simulation Package (VASP), Quantum Espresso and Spanish Initiative for Electronic Simulations with Thousands of Atoms (SIESTA). This module is easily extendable and can be interfaced with other DFT, force-field or empirical potential based codes. The uniqueness of the module lies in its user-friendliness. For users with limited computing resources, this module will be an effective tool as it allows to perform the calculations image by image. On the other hand, users with plentiful computing resources (such as users in a high performance computing environment) can perform the calculations for large number of images simultaneously. This module gives users complete flexibility, thereby enabling them to perform calculations on large systems making the best use of the available resources.


title: Implementation of Artifact Detection in Critical Care: A Methodological Review

abstract: Artifact Detection (AD) techniques minimize the impact of artifacts on physiologic data acquired in Critical Care Units (CCU) by assessing quality of data prior to Clinical Event Detection (CED) and Parameter Derivation (PD). This methodological review introduces unique taxonomies to synthesize over 80 AD algorithms based on these six themes: (1) CCU; (2) Physiologic Data Source; (3) Harvested data; (4) Data Analysis; (5) Clinical Evaluation; and (6) Clinical Implementation. Review results show that most published algorithms: (a) are designed for one specific type of CCU; (b) are validated on data harvested only from one Original Equipment Manufacturer (OEM) monitor; (c) generate Signal Quality Indicators (SQI) that are not yet formalised for useful integration in clinical workflows; (d) operate either in standalone mode or coupled with CED or PD applications; (e) are rarely evaluated in real-time; and (f) are not implemented in clinical practice. In conclusion, it is recommended that AD algorithms conform to generic input and output interfaces with commonly defined data: (1) type; (2) frequency; (3) length; and (4) SQIs. This shall promote (a) reusability of algorithms across different CCU domains; (b) evaluation on different OEM monitor data; (c) fair comparison through formalised SQIs; (d) meaningful integration with other AD, CED and PD algorithms; and (e) real-time implementation in clinical workflows.


title: Detailed study on the Fermi surfaces of the type-II Dirac semimetallic candidates PdTe2 and PtTe2

abstract: We present a detailed quantum oscillatory study on the Dirac type-II semimetallic candidates PdTe$_{2}$ and PtTe$_{2}$ \emph{via} the temperature and the angular dependence of the de Haas-van Alphen (dHvA) and Shubnikov-de Haas (SdH) effects. In high quality single crystals of both compounds, i.e. displaying carrier mobilities between $10^3$ and $10^4$ cm$^2$/Vs, we observed a large non-saturating magnetoresistivity (MR) which in PtTe$_2$ at a temperature $T = 1.3$ K, leads to an increase in the resistivity up to $5 \times 10^{4}$ % under a magnetic field $\mu_0 H = 62$ T. These high mobilities correlate with their light effective masses in the range of 0.04 to 1 bare electron mass according to our measurements. For PdTe$_{2}$ the experimentally determined Fermi surface cross-sectional areas show an excellent agreement with those resulting from band-structure calculations. Surprisingly, this is not the case for PtTe$_{2}$ whose agreement between calculations and experiments is relatively poor even when electronic correlations are included in the calculations. Therefore, our study provides a strong support for the existence of a Dirac type-II node in PdTe$_2$ and probably also for PtTe$_2$. Band structure calculations indicate that the topologically non-trivial bands of PtTe$_2$ do not cross the Fermi-level ($\varepsilon_F$). In contrast, for PdTe$_2$ the Dirac type-II cone does intersect $\varepsilon_F$, although our calculations also indicate that the associated cyclotron orbit on the Fermi surface is located in a distinct $k_z$ plane with respect to the one of the Dirac type-II node. Therefore it should yield a trivial Berry-phase.


title: Freeze-out temperature from net-Kaon fluctuations at RHIC

abstract: We compare the mean-over-variance ratio of the net-kaon distribution calculated within a state-of-the-art hadron resonance gas model to the latest experimental data from the Beam Energy Scan at RHIC by the STAR collaboration. Our analysis indicates that it is not possible to reproduce the experimental results using the freeze-out parameters from the existing combined fit of net-proton and net-electric charge mean-over-variance. The strange mesons need about 10-15 MeV higher temperatures than the light hadrons at the highest collision energies. In view of the future $\Lambda$ fluctuation measurements, we predict the $\Lambda$ variance-over-mean and skewness-times-variance at the light and strange chemical freeze-out parameters. We observe that the $\Lambda$ fluctuations are sensitive to the difference in the freeze-out temperatures established in this analysis. Our results have implications for other phenomenological models in the field of relativistic heavy ion collisions.


title: Concolic Testing for Deep Neural Networks

abstract: Concolic testing combines program execution and symbolic analysis to explore the execution paths of a software program. This paper presents the first concolic testing approach for Deep Neural Networks (DNNs). More specifically, we formalise coverage criteria for DNNs that have been studied in the literature, and then develop a coherent method for performing concolic testing to increase test coverage. Our experimental results show the effectiveness of the concolic testing approach in both achieving high coverage and finding adversarial examples.


title: Using Multi Expression Programming in Software Effort Estimation

abstract: Estimating the effort of software systems is an essential topic in software engineering, carrying out an estimation process reliably and accurately for a software forms a vital part of the software development phases. Many researchers have utilized different methods and techniques hopping to find solutions to this issue, such techniques include COCOMO, SEER-SEM,SLIM and others. Recently, Artificial Intelligent techniques are being utilized to solve such problems; different studies have been issued focusing on techniques such as Neural Networks NN, Genetic Algorithms GA, and Genetic Programming GP. This work uses one of the linear variations of GP, namely: Multi Expression Programming (MEP) aiming to find the equation that best estimates the effort of software. Benchmark datasets (based on previous projects) are used learning and testing. Results are compared with those obtained by GP using different fitness functions. Results show that MEP is far better in discovering effective functions for the estimation of about 6 datasets each comprising several projects.


title: On the probability distribution associated to commutator word map in finite groups \rom{2}

abstract: Let $P(G)$ denotes the set of sizes of fibers of non-trivial commutators of the commutator word map. Here, we prove that $|P(G)|=1$, for any finite group $G$ of nilpotency class $3$ with exactlly two conjugacy class sizes. We also show that for given $n\geq 1$, there exists a finite group $G$ of nilpotency class $2$ with exactlly two conjugacy class sizes such that $|P(G)|=n$.


title: New Methods of Studying Valley Fitness Landscapes

abstract: The word "valley" is a popular term used in intuitively describing fitness landscapes. What is a valley on a fitness landscape? How to identify the direction and location of a valley if it exists? However, such questions are seldom rigorously studied in evolutionary optimization especially when the search space is a high dimensional continuous space. This paper presents two methods of studying valleys on a fitness landscape. The first method is based on the topological homeomorphism. It establishes a rigorous definition of a valley. A valley is regarded as a one-dimensional manifold. The second method takes a different viewpoint from statistics. It provides an algorithm of identifying the valley direction and location using principle component analysis.


title: A new perspective on the powers of two descent for discrete logarithms in finite fields

abstract: A new proof is given for the correctness of the powers of two descent method for computing discrete logarithms. The result is slightly stronger than the original work, but more importantly we provide a unified geometric argument, eliminating the need to analyse all possible subgroups of $\mathrm{PGL}_2(\mathbb F_q)$. Our approach sheds new light on the role of $\mathrm{PGL}_2$, in the hope to eventually lead to a complete proof that discrete logarithms can be computed in quasi-polynomial time in finite fields of fixed characteristic.


title: Comparative Analysis and Framework Evaluating Web Single Sign-On Systems

abstract: We perform a comprehensive analysis and comparison of 14 web single sign-on (SSO) systems proposed and/or deployed over the last decade, including federated identity and credential/password management schemes. We identify common design properties and use them to develop a taxonomy for SSO schemes, highlighting the associated trade-offs in benefits (positive attributes) offered. We develop a framework to evaluate the schemes, in which we identify 14 security, usability, deployability, and privacy benefits. We also discuss how differences in priorities between users, service providers (SPs), and identity providers (IdPs) impact the design and deployment of SSO schemes.


title: Stationary Bernstein-Greene-Kruskal structures in a current carrying relativistic cold plasma

abstract: Nonlinear stationary structures formed in a cold plasma with immobile ions in the presence of a relativistic electron current beam have been investigated analytically in the collisionless limit. These are cold plasma version of the relativistic Berstein-Greene-Kruskal (BGK) waves. The structure profile is governed by the ratio of maximum electrostatic field energy density to the relativistic kinetic energy density of the electron beam, {\it i.e.}, $\kappa_{R} = E_{m}/ (8 \pi n_{0} (\gamma_{0}-1)m_{0} c^{2})^{1/2}$, where $E_{m}$ is the maximum electric field associated with the nonlinear structure and $\gamma_{0}$ is the Lorentz factor associated with the beam velocity. It is found that, in the linear limit, {\it i.e.}, $\kappa_{R} \ll 1/\sqrt{\gamma_{0}}$, the fluid variables, {\it viz}, density, electric field, and velocity vary harmonically in space. In the range $0 < \kappa_{R} \leq 1/\sqrt{\gamma_{0}}$, the fluid variables exhibit an-harmonic behavior. For values of $\kappa_{R} > 1/\sqrt{\gamma_{0}}$, the electric field shows finite discontinuities at specific spatial locations indicating the formation of negatively charged planes at these locations.


title: From thermal dissociation to condensation in the atmospheres of ultra hot Jupiters: WASP-121b in context

abstract: A new class of exoplanets has emerged: the ultra hot Jupiters, the hottest close-in gas giants. Most of them have weaker than expected spectral features in the $1.1-1.7\mu m$ bandpass probed by HST/WFC3 but stronger spectral features at longer wavelengths probed by Spitzer. This led previous authors to puzzling conclusions about the thermal structures and chemical abundances of these planets. Using the SPARC/MITgcm, we investigate how thermal dissociation, ionization, H$^-$ opacity and clouds shape the thermal structures and spectral properties of ultra hot Jupiters with a special focus on WASP-121b. We expand our findings to the whole population of ultra hot Jupiters through analytical quantification of the thermal dissociation and its influence on the strength of spectral features. We predict that most molecules are thermally dissociated and alkalies are ionized in the dayside photospheres of ultra hot Jupiters. This includes H$_{\rm 2}$O, TiO, VO, and H$_{\rm 2}$ but not CO, which has a stronger molecular bond. The vertical molecular gradient created by the dissociation significantly weakens the spectral features from water while the $4.5\mu m$ CO feature remain unchanged. The water band in the HST/WFC3 bandpass is further weakened by H$^-$ continuum opacity. Molecules are expected to recombine before reaching the limb, leading to order of magnitude variations of the chemical composition and cloud coverage between the limb and the dayside. Overall, molecular dissociation provides a qualitative understanding of the lack of strong spectral feature of water in the $1-2\mu m$ bandpass observed in most ultra hot Jupiters. Quantitatively, however, our model does not provide a satisfactory match to the WASP-121b emission spectrum. Together with WASP-33b and Kepler-33Ab, they seem the outliers among the population of ultra hot Jupiters in need of a more thorough understanding.


title: Syntactic Patterns Improve Information Extraction for Medical Search

abstract: Medical professionals search the published literature by specifying the type of patients, the medical intervention(s) and the outcome measure(s) of interest. In this paper we demonstrate how features encoding syntactic patterns improve the performance of state-of-the-art sequence tagging models (both linear and neural) for information extraction of these medically relevant categories. We present an analysis of the type of patterns exploited, and the semantic space induced for these, i.e., the distributed representations learned for identified multi-token patterns. We show that these learned representations differ substantially from those of the constituent unigrams, suggesting that the patterns capture contextual information that is otherwise lost.


title: Light-ray operators in conformal field theory

abstract: We argue that every CFT contains light-ray operators labeled by a continuous spin J. When J is a positive integer, light-ray operators become integrals of local operators over a null line. However for non-integer J, light-ray operators are genuinely nonlocal and give the analytic continuation of CFT data in spin described by Caron-Huot. A key role in our construction is played by a novel set of intrinsically Lorentzian integral transforms that generalize the shadow transform. Matrix elements of light-ray operators can be computed via the integral of a double-commutator against a conformal block. This gives a simple derivation of Caron-Huot's Lorentzian OPE inversion formula and lets us generalize it to arbitrary four-point functions. Furthermore, we show that light-ray operators enter the Regge limit of CFT correlators, and generalize conformal Regge theory to arbitrary four-point functions. The average null energy operator is an important example of a light-ray operator. Using our construction, we find a new proof of the average null energy condition (ANEC), and furthermore generalize the ANEC to continuous spin.


title: Reionization in Technicolor

abstract: We present the Technicolor Dawn simulations, a suite of cosmological radiation-hydrodynamic simulations of the first 1.2 billion years. By modeling a spatially-inhomogeneous UVB on-the-fly with 24 frequencies and resolving dark matter halos down to $10^8 M_\odot$ within 12 $h^{-1}$ Mpc volumes, our simulations unify observations of the intergalactic and circumgalactic media, galaxies, and reionization into a common framework. The only empirically-tuned parameter, the fraction $f_{\mathrm{esc,gal}}(z)$ of ionizing photons that escape the interstellar medium, is adjusted to match observations of the Lyman-$\alpha$ forest and the cosmic microwave background. With this single calibration, our simulations reproduce the history of reionization; the stellar mass-star formation rate relation of galaxies; the number density and metallicity of damped Lyman-$\alpha$ absorbers (DLAs) at $z\sim5$; the abundance of weak metal absorbers; the ultraviolet background (UVB) amplitude; and the Lyman-$\alpha$ flux power spectrum at $z=5.4$. The galaxy stellar mass and UV luminosity functions are underproduced by $\leq2\times$, suggesting an overly vigorous feedback model. The mean transmission in the Lyman-$\alpha$ forest is underproduced at $z<6$, indicating tension between measurements of the UVB amplitude and Lyman-$\alpha$ transmission. The observed SiIV column density distribution is reasonably well-reproduced ($\sim 1\sigma$ low). By contrast, CIV remains significantly underproduced despite being boosted by an intense $>4$ Ryd UVB. Solving this problem by increasing metal yields would overproduce both weak absorbers and DLA metallicities. Instead, the observed strength of high-ionization emission from high-redshift galaxies and absorption from their environments suggest that the ionizing flux from conventional stellar population models is too soft.


title: Control of Energy Storage in Home Energy Management Systems: Non-Simultaneous Charging and Discharging Guarantees

abstract: In this paper we provide non-simultaneous charging and discharging guarantees for a linear energy storage system (ESS) model for a model predictive control (MPC) based home energy management system (HEMS) algorithm. The HEMS optimally controls the residential load and residentially-owned power sources, such as photovoltaic (PV) power generation and energy storage, given residential customer preferences such as energy cost sensitivity and ESS lifetime. Under certain problem formulations with a linear ESS model, simultaneous charging and discharging can be observed as the optimal solution when there is high penetration of PV power. We present analysis for a proposed HEMS optimization formulation that ensures non-simultaneous ESS charging and discharging operation for a linear ESS model that captures both charging and discharging efficiency of the ESS. The energy storage system model behavior guarantees are shown for various electricity pricing schemes such as time of use (TOU) pricing and net metering. Simulation results demonstrating desirable ESS behavior are provided for each electricity pricing scheme.


title: Gravity effects on mixing with magnetic micro-convection in microfluidics

abstract: Mixing remains an important problem for development of successful microfluidic and lab-on-a-chip devices, where simple and predictable systems are particularly interesting. One is magnetic micro-convection, an instability happening on the interface of miscible magnetic and non-magnetic fluids in a Hele-Shaw cell under applied field. Previous work proved that Brinkman model quantitatively explains the experiments. However, a gravity caused convective motion complicated the tests. Here we first improve the experimental system to exclude the gravitational convective motion. Afterwards, we observe and quantify how gravity and laminar flow play an important role in stabilizing the perturbations that create the instability. Accordingly, we improve our theoretical model and perform linear analysis. Two dimensionless quantities explain the experimental observations of change in critical field needed for instability and characteristic size of the emerging pattern. Finally, we discuss the conditions at which gravity plays an important role in microfluidic systems.


title: Locality in the Fukaya category of a hyperk\"ahler manifold

abstract: Let $(M,I,J,K,g)$ be a hyperk\"ahler manifold. Then the complex manifold $(M,I)$ is holomorphic symplectic. We prove that for all real $x, y,$ with $x^2 + y^2 = 1$ except countably many, any finite energy $(xJ+yK)$-holomorphic curve with boundary in a collection of $I$-holomorphic Lagrangians must be constant. By an argument based on the Lojasiewicz inequality, this result holds no matter how the Lagrangians intersect each other. It follows that one can choose perturbations such that the holomorphic polygons of the associated Fukaya category lie in an arbitrarily small neighborhood of the Lagrangians. That is, the Fukaya category is local. We show that holomorphic Lagrangians are tautologically unobstructed in the sense of Fukaya-Oh-Ohta-Ono. Moreover, the Fukaya $A_\infty$ algebra of a holomorphic Lagrangian is formal. Our result also explains why the special Lagrangian condition holds without instanton corrections for holomorphic Lagrangians.


title: Momentum space treatment of inclusive neutrino scattering off the deuteron and trinucleons

abstract: The $\bar{\nu}_e + {^2{\rm H}} \rightarrow e^+ + n + n$, $\nu_e + {^2{\rm H}} \rightarrow e^- + p + p$, $\bar{\nu}_l + {^2{\rm H}} \rightarrow \bar{\nu}_l + {^2{\rm H}}$, $\nu_l + {^2{\rm H}} \rightarrow \nu_l + {^2{\rm H}}$, $\bar{\nu}_l + {^2{\rm H}} \rightarrow \bar{\nu}_l + p + n$, $\nu_l + {^2{\rm H}} \rightarrow \nu_l + p + n$, $\bar{\nu}_e + {^3{\rm He}} \rightarrow e^+ + {^3{\rm H}}$, $\bar{\nu}_l + {^3{\rm He}} \rightarrow \bar{\nu}_l + {^3{\rm He}}$, $\nu_l + {^3{\rm He}} \rightarrow \nu_l + {^3{\rm He}}$, $\bar{\nu}_l + {^3{\rm H}} \rightarrow \bar{\nu}_l + {^3{\rm H}}$, $\nu_l + {^3{\rm H}} \rightarrow \nu_l + {^3{\rm H}}$, $\bar{\nu}_e + {^3{\rm He}} \rightarrow e^+ + n + d$, $\bar{\nu}_e + {^3{\rm He}} \rightarrow e^+ + n + n + p$, $\bar{\nu}_l + {^3{\rm He}} \rightarrow \bar{\nu}_l + p + d$, $\bar{\nu}_l + {^3{\rm He}} \rightarrow \bar{\nu}_l + p + p +n$, $\nu_l + {^3{\rm H}} \rightarrow \nu_l + n + d$ and $\nu_l + {^3{\rm H}} \rightarrow \nu_l + n + n + p$ reactions ($l= e, \mu, \tau$) are studied consistently in momentum space for (anti)neutrino energies up to 300 MeV. For most of these processes we provide predictions for the total cross sections and in the case of the (anti)neutrino-$^3$He and (anti)neutrino-$^3$H inelastic scattering we compute examples of essential response functions, using the AV18 nucleon-nucleon potential and a single-nucleon weak current operator. For the reactions with the deuteron we study relativistic effects in the final state kinematics and compare two-nucleon scattering states obtained in momentum and coordinate spaces. Our results from momentum space are compared with the theoretical predictions by G.Shen et al., Phys. Rev. C 86, 035503 (2012). The observed disagreement can be attributed to the differences in kinematics and in the weak current operator.


title: Niebrzydowski Algebras and Trivalent Spatial Graphs

abstract: We introduce \textit{Niebrzydowski algebras}, algebraic structures with a ternary operation and a partially defined multiplication, with axioms motivated by the Reidemeister moves for $Y$-oriented trivalent spatial graphs and handlebody-links. As part of this definition, we identify generating sets of $Y$-oriented Reidemeister moves. We give some examples to demonstrate that the counting invariant can distinguish some $Y$-oriented trivalent spatial graphs and handlebody-links.


title: A Cyberinfrastructure for BigData Transportation Engineering

abstract: Big Data-driven transportation engineering has the potential to improve utilization of road infrastructure, decrease traffic fatalities, improve fuel consumption, decrease construction worker injuries, among others. Despite these benefits, research on Big Data-driven transportation engineering is difficult today due to the computational expertise required to get started. This work proposes BoaT, a transportation-specific programming language, and it's Big Data infrastructure that is aimed at decreasing this barrier to entry. Our evaluation that uses over two dozen research questions from six categories show that research is easier to realize as a BoaT computer program, an order of magnitude faster when this program is run, and exhibits 12-14x decrease in storage requirements.


title: FPGA Acceleration of Short Read Alignment

abstract: Aligning millions of short DNA or RNA reads, of 75 to 250 base pairs each, to a reference genome is a significant computation problem in bioinformatics. We present a flexible and fast FPGA-based short read alignment tool. Our aligner makes use of the processing power of FPGAs in conjunction with the greater host memory bandwidth and flexibility of software to improve performance and achieve a high level of configurability. This flexible design supports a variety of reference genome sizes without the performance degradation suffered by other software and FPGA-based aligners. It is also better able to support the features of new alignment algorithms, which frequently crop up in the rapidly evolving field of bioinformatics. We demonstrate these advantages in a case study where we align RNA-Seq data from a hypothesized mouse / human xenograft. In this case study, our aligner provides a speedup of 5.6x over BWA-SW with energy savings of 21%, while also reducing incorrect short read classification by 29%. To demonstrate the flexibility of our system we show that the speedup can be substantially improved while retaining most of the accuracy gains over BWA-SW. The speedup can be increased to 71.3x, while still enjoying a 28% incorrect classification improvement and 52% improvement in unaligned reads.


title: MV-YOLO: Motion Vector-aided Tracking by Semantic Object Detection

abstract: Object tracking is the cornerstone of many visual analytics systems. While considerable progress has been made in this area in recent years, robust, efficient, and accurate tracking in real-world video remains a challenge. In this paper, we present a hybrid tracker that leverages motion information from the compressed video stream and a general-purpose semantic object detector acting on decoded frames to construct a fast and efficient tracking engine. The proposed approach is compared with several well-known recent trackers on the OTB tracking dataset. The results indicate advantages of the proposed method in terms of speed and/or accuracy.Other desirable features of the proposed method are its simplicity and deployment efficiency, which stems from the fact that it reuses the resources and information that may already exist in the system for other reasons.


title: Conditional molecular design with deep generative models

abstract: Although machine learning has been successfully used to propose novel molecules that satisfy desired properties, it is still challenging to explore a large chemical space efficiently. In this paper, we present a conditional molecular design method that facilitates generating new molecules with desired properties. The proposed model, which simultaneously performs both property prediction and molecule generation, is built as a semi-supervised variational autoencoder trained on a set of existing molecules with only a partial annotation. We generate new molecules with desired properties by sampling from the generative distribution estimated by the model. We demonstrate the effectiveness of the proposed model by evaluating it on drug-like molecules. The model improves the performance of property prediction by exploiting unlabeled molecules, and efficiently generates novel molecules fulfilling various target conditions.


title: Quantum computational finance: Monte Carlo pricing of financial derivatives

abstract: Financial derivatives are contracts that can have a complex payoff dependent upon underlying benchmark assets. In this work, we present a quantum algorithm for the Monte Carlo pricing of financial derivatives. We show how the relevant probability distributions can be prepared in quantum superposition, the payoff functions can be implemented via quantum circuits, and the price of financial derivatives can be extracted via quantum measurements. We show how the amplitude estimation algorithm can be applied to achieve a quadratic quantum speedup in the number of steps required to obtain an estimate for the price with high confidence. This work provides a starting point for further research at the interface of quantum computing and finance.


title: A MultiMesh Finite Element Method for the Stokes Problem

abstract: The multimesh finite element method enables the solution of partial differential equations on a computational mesh composed by multiple arbitrarily overlapping meshes. The discretization is based on a continuous--discontinuous function space with interface conditions enforced by means of Nitsche's method. In this contribution, we consider the Stokes problem as a first step towards flow applications. The multimesh formulation leads to so called cut elements in the underlying meshes close to overlaps. These demand stabilization to ensure coercivity and stability of the stiffness matrix. We employ a consistent least-squares term on the overlap to ensure that the inf-sup condition holds. We here present the method for the Stokes problem, discuss the implementation, and verify that we have optimal convergence.


title: Perturbative traveling wave solution for a flux-limited reaction-diffusion morphogenesis equation

abstract: In this study, we investigate a porous medium-type flux limited reaction--diffusion equation that arises in morphogenesis modeling. This nonlinear partial differential equation is an extension of the generalized Fisher--Kolmogorov--Petrovsky--Piskunov (Fisher-KPP) equation in one-dimensional space. The approximate analytical traveling wave solution is found by using a perturbation method. We show that the morphogen concentration propagates as a sharp wave front where the wave speed has a saturated value. The numerical solutions of this equation are also provided to compare them with the analytical predictions. Finally, we qualitatively compare our theoretical results with those obtained in experimental studies.


title: Viability analysis of the first-order mean field games

abstract: The paper is concerned with the dependence of the solution of the deterministic mean field game on the initial distribution of players. The main object of study is the mapping which assigns to the initial time and the initial distribution of players the set of expected rewards of the representative player corresponding to solutions of mean field game. This mapping can be regarded as a value multifunction. We obtain the sufficient condition for a multifunction to be a value multifunction. It states that if a multifunction is viable with respect to the dynamics generated by the original mean field game, then it is a value multifunction. Furthermore, the infinitesimal variant of this condition is derived.


title: Identities from representation theory

abstract: We give a new Jacobi--Trudi-type formula for characters of finite-dimensional irreducible representations in type $C_n$ using characters of the fundamental representations and non-intersecting lattice paths. We give equivalent determinant formulas for the decomposition multiplicities for tensor powers of the spin representation in type $B_n$ and the exterior representation in type $C_n$. This gives a combinatorial proof of an identity of Katz and equates such a multiplicity with the dimension of an irreducible representation in type $C_n$. By taking certain specializations, we obtain identities for $q$-Catalan triangle numbers, the $q,t$-Catalan number of Stump, $q$-triangle versions of Motzkin and Riordan numbers, and generalizations of Touchard's identity. We use (spin) rigid tableaux and crystal base theory to show some formulas relating Catalan, Motzkin, and Riordan triangle numbers.


title: Algebraic dual polynomials for the equivalence of curl-curl problems

abstract: In this paper we will consider two curl-curl equation in two dimensions. One curl-curl problem for a scalar quantity $F$ and one problem for a vector field $\bf{E}$. For Dirichlet boundary conditions $\bf{n} \times \bf{E} =$ $ \hat{E}_{\dashv}$ on $\bf{E}$ and Neumann boundary conditions $\bf{n} \times \mbox{curl}$ $F=\hat{E}_{\dashv}$, we expect the solutions to satisfy $\bf{E}=\mbox{curl}$ $F$. When we use algebraic dual polynomial representations, these identities continue to hold at the discrete level. Equivalence will be proved and illustrated with a computational example.


title: Counting tropical rational curves with cross-ratio constraints

abstract: We enumerate rational curves in toric surfaces passing through points and satisfying cross-ratio constraints using tropical and combinatorial methods. Our starting point is arXiv:1509.07453, where a tropical-algebraic correspondence theorem was proved that relates counts of rational curves in toric varieties that satisfy point conditions and cross-ratio constraints to the analogous tropical counts. We proceed in two steps: based on tropical intersection theory we first study tropical cross-ratios and introduce degenerated cross-ratios. Second we provide a lattice path algorithm that produces all tropical curves satisfying such degenerated conditions explicitly. In a special case simpler combinatorial objects, so-called cross-ratio floor diagrams, are introduced which can be used to determine these enumerative numbers as well.


title: A Canonical Image Set for Examining and Comparing Image Processing Algorithms

abstract: The purpose of this paper is to introduce a set of four test images containing features and structures that can facilitate effective examination and comparison of image processing algorithms. More specifically, the images are designed to more explicitly expose the characteristic properties of algorithms for image compression, virtual resolution adjustment, and enhancement. This set was developed at the Naval Research Laboratory (NRL) in the late 1990s as a more rigorous alternative to Lena and other images that have come into common use for purely ad hoc reasons with little or no rigorous consideration of their suitability. The increasing number of test images appearing in the literature not only makes it more difficult to compare results from different papers, it also introduces the potential for cherry-picking to influence results. The key contribution of this paper is the proposal to establish {\em some} canonical set to ensure that published results can be analyzed and compared in a rigorous way from one paper to another, and consideration of the four NRL images is proposed for this purpose.


title: Fixation probabilities for the Moran process in evolutionary games with two strategies: graph shapes and large population asymptotics

abstract: This paper is based on the complete classification of evolutionary scenarios for the Moran process with two strategies given by Taylor et al. (B. Math. Biol. 66(6): 1621--1644, 2004). Their classification is based on whether each strategy is a Nash equilibrium and whether the fixation probability for a single individual of each strategy is larger or smaller than its value for neutral evolution. We improve on this analysis by showing that each evolutionary scenario is characterized by a definite graph shape for the fixation probability function. A second class of results deals with the behavior of the fixation probability when the population size tends to infinity. We develop asymptotic formulae that approximate the fixation probability in this limit and conclude that some of the evolutionary scenarios cannot exist when the population size is large.


title: Uncertainty Quantification in Non-Equilibrium Molecular Dynamics Simulations of Thermal Transport

abstract: Bulk thermal conductivity estimates based on predictions from non-equilibrium molecular dynamics (NEMD) using the so-called direct method are known to be severely under-predicted since finite simulation length-scales are unable to mimic bulk transport. Moreover, subjecting the system to a temperature gradient by means of thermostatting tends to impact phonon transport adversely. Additionally, NEMD predictions are tightly coupled with the choice of the inter-atomic potential and the underlying values associated with its parameters. In the case of silicon (Si), nominal estimates of the Stillinger-Weber (SW) potential parameters are largely based on a constrained regression approach aimed at agreement with experimental data while ensuring structural stability. However, this approach has its short-comings and it may not be ideal to use the same set of parameters to study a wide variety of Si-based systems subjected to different thermodynamic conditions. In this study, NEMD simulations are performed on a Si bar to investigate the impact of bar-length, and the applied thermal gradient on the discrepancy between predictions and the available measurement for bulk thermal conductivity at 300 K by constructing statistical response surfaces at different temperatures. The approach helps quantify the discrepancy, observed to be largely dependent on the system-size, with minimal computational effort. A computationally efficient approach based on derivative-based sensitivity measures to construct a reduced-order polynomial chaos surrogate for NEMD predictions is also presented. The surrogate is used to perform parametric sensitivity analysis, forward propagation of the uncertainty, and calibration of the important SW potential parameters in a Bayesian setting. It is found that only two (out of seven) parameters contribute significantly to the uncertainty in bulk thermal conductivity estimates for Si.


title: Risk-Averse Classification

abstract: We develop a new approach to solving classification problems, which is bases on the theory of coherent measures of risk and risk sharing ideas. The proposed approach aims at designing a risk-averse classifier. The new approach allows for associating distinct risk functional to each classes. The risk may be measured by different (non-linear in probability) measures, We analyze the structure of the new classifier design problem and establish its theoretical relation to known risk-neutral design problems. In particular, we show that the risk-sharing classification problem is equivalent to an implicitly defined optimization problem with unequal, implicitly defined but unknown, weights for each data point. We implement our methodology in a binary classification scenario on several different data sets and carry out numerical comparison with classifiers which are obtained using the Huber loss function and other loss functions known in the literature. We formulate specific risk-averse support vector machines in order to demonstrate the viability of our method.


title: Types for Information Flow Control: Labeling Granularity and Semantic Models

abstract: Language-based information flow control (IFC) tracks dependencies within a program using sensitivity labels and prohibits public outputs from depending on secret inputs. In particular, literature has proposed several type systems for tracking these dependencies. On one extreme, there are fine-grained type systems (like Flow Caml) that label all values individually and track dependence at the level of individual values. On the other extreme are coarse-grained type systems (like HLIO) that track dependence coarsely, by associating a single label with an entire computation context and not labeling all values individually. In this paper, we show that, despite their glaring differences, both these styles are, in fact, equally expressive. To do this, we show a semantics- and type-preserving translation from a coarse-grained type system to a fine-grained one and vice-versa. The forward translation isn't surprising, but the backward translation is: It requires a construct to arbitrarily limit the scope of a context label in the coarse-grained type system (e.g., HLIO's "toLabeled" construct). As a separate contribution, we show how to extend work on logical relation models of IFC types to higher-order state. We build such logical relations for both the fine-grained type system and the coarse-grained type system. We use these relations to prove the two type systems and our translations between them sound.


title: A Missing Information Loss function for implicit feedback datasets

abstract: Latent factor models for Recommender Systems with implicit feedback typically treat unobserved user-item interactions (i.e. missing information) as negative feedback. This is frequently done either through negative sampling (point--wise loss) or with a ranking loss function (pair-- or list--wise estimation). Since a zero preference recommendation is a valid solution for most common objective functions, regarding unknown values as actual zeros results in users having a zero preference recommendation for most of the available items. In this paper we propose a novel objective function, the \emph{Missing Information Loss} (MIL), that explicitly forbids treating unobserved user-item interactions as positive or negative feedback. We apply this loss to both traditional Matrix Factorization and user--based Denoising Autoencoder, and compare it with other established objective functions such as cross-entropy (both point- and pair-wise) or the recently proposed multinomial log-likelihood. MIL achieves competitive performance in ranking-aware metrics when applied to three datasets. Furthermore, we show that such a relevance in the recommendation is obtained while displaying popular items less frequently (up to a $20 \%$ decrease with respect to the best competing method). This debiasing from the recommendation of popular items favours the appearance of infrequent items (up to a $50 \%$ increase of long-tail recommendations), a valuable feature for Recommender Systems with a large catalogue of products.


title: Galactic Rotation Curves vs. Ultra-Light Dark Matter: Implications of the Soliton -- Host Halo Relation

abstract: Bosonic ultra-light dark matter (ULDM) would form cored density distributions at the center of galaxies. These cores, seen in numerical simulations, admit analytic description as the lowest energy bound state solution ("soliton") of the Schroedinger-Poisson equations. Numerical simulations of ULDM galactic halos found empirical scaling relations between the mass of the large-scale host halo and the mass of the central soliton. We discuss how the simulation results of different groups can be understood in terms of the basic properties of the soliton. Importantly, simulations imply that the energy per unit mass in the soliton and in the virialised host halo should be approximately equal. This relation lends itself to observational tests, because it predicts that the peak circular velocity, measured for the host halo in the outskirts of the galaxy, should approximately repeat itself in the central region. Contrasting this prediction to the measured rotation curves of well-resolved near-by galaxies, we show that ULDM in the mass range $m\sim (10^{-22}\div 10^{-21})$ eV, which has been invoked as a possible solution to the small-scale puzzles of $\Lambda$CDM, is in tension with the data. We suggest that a dedicated analysis of the Milky Way inner gravitational potential could probe ULDM up to $m\lesssim 10^{-19}$ eV.


title: CrowdHuman: A Benchmark for Detecting Human in a Crowd

abstract: Human detection has witnessed impressive progress in recent years. However, the occlusion issue of detecting human in highly crowded environments is far from solved. To make matters worse, crowd scenarios are still under-represented in current human detection benchmarks. In this paper, we introduce a new dataset, called CrowdHuman, to better evaluate detectors in crowd scenarios. The CrowdHuman dataset is large, rich-annotated and contains high diversity. There are a total of $470K$ human instances from the train and validation subsets, and $~22.6$ persons per image, with various kinds of occlusions in the dataset. Each human instance is annotated with a head bounding-box, human visible-region bounding-box and human full-body bounding-box. Baseline performance of state-of-the-art detection frameworks on CrowdHuman is presented. The cross-dataset generalization results of CrowdHuman dataset demonstrate state-of-the-art performance on previous dataset including Caltech-USA, CityPersons, and Brainwash without bells and whistles. We hope our dataset will serve as a solid baseline and help promote future research in human detection tasks.


title: Force Generation by Parallel Combinations of Fiber-Reinforced Fluid-Driven Actuators

abstract: The compliant structure of soft robotic systems enables a variety of novel capabilities in comparison to traditional rigid-bodied robots. A subclass of soft fluid-driven actuators known as fiber reinforced elastomeric enclosures (FREEs) is particularly well suited as actuators for these types of systems. FREEs are inherently soft and can impart spatial forces without imposing a rigid structure. Furthermore, they can be configured to produce a large variety of force and moment combinations. In this paper we explore the potential of combining multiple differently configured FREEs in parallel to achieve fully controllable multi-dimensional soft actuation. To this end, we propose a novel methodology to represent and calculate the generalized forces generated by soft actuators as a function of their internal pressure. This methodology relies on the notion of a state dependent fluid Jacobian that yields a linear expression for force. We employ this concept to construct the set of all possible forces that can be generated by a soft system in a given state. This force zonotope can be used to inform the design and control of parallel combinations of soft actuators. The approach is verified experimentally with the parallel combination of three carefully designed actuators constrained to a 2DOF testing rig. The force predictions matched measured values with a root-mean-square error of less than 1.5 N force and 8 x 10^(-3)Nm moment, demonstrating the utility of the presented methodology.


title: Spin Polarized Non-Relativistic Fermions in 1+1 Dimensions

abstract: We study by Monte Carlo methods the thermodynamics of a spin polarized gas of non-relativistic fermions in 1+1 dimensions. The main result of this work is that our action suffers no significant sign problem for any spin polarization in the region relevant for dilute degenerate fermi gases. This lack of sign problem allows us to study attractive spin polarized fermions non-perturbatively at spin polarizations not previously explored. For some parameters values we verify results previously obtained by methods which include an uncontrolled step like complex Langevin and/or analytical continuation from imaginary chemical potential. For others, larger values of the polarization, we deviate from these previous results.


title: Superconductors that do not expel magnetic flux

abstract: An ultrathin superconducting bilayer creates a coreless fractional vortex when only the second layer has a hole. The quantization is broken by the hole, and the normal core disappears. The magnetic flux is no longer confined near the normal core, and its density profile around the hole becomes similar to that of a cr\`eme caramel; the divergence of the magnetic flux density is truncated around the center. We propose basic design of a practical device to realize a coreless fractional vortex.


title: Hundredfold Enhancement of Light Emission via Defect Control in Monolayer Transition-Metal Dichalcogenides

abstract: Two dimensional (2D) transition-metal dichalcogenide (TMD) based semiconductors have generated intense recent interest due to their novel optical and electronic properties, and potential for applications. In this work, we characterize the atomic and electronic nature of intrinsic point defects found in single crystals of these materials synthesized by two different methods - chemical vapor transport and self-flux growth. Using a combination of scanning tunneling microscopy (STM) and scanning transmission electron microscopy (STEM), we show that the two major intrinsic defects in these materials are metal vacancies and chalcogen antisites. We show that by control of the synthetic conditions, we can reduce the defect concentration from above $10^{13} /cm^2$ to below $10^{11} /cm^2$. Because these point defects act as centers for non-radiative recombination of excitons, this improvement in material quality leads to a hundred-fold increase in the radiative recombination efficiency.


title: The effect of cosmic-ray acceleration on supernova blast wave dynamics

abstract: Non-relativistic shocks accelerate ions to highly relativistic energies provided that the orientation of the magnetic field is closely aligned with the shock normal (quasi-parallel shock configuration). In contrast, quasi-perpendicular shocks do not efficiently accelerate ions. We model this obliquity-dependent acceleration process in a spherically expanding blast wave setup with the moving-mesh code {\sc arepo} for different magnetic field morphologies, ranging from homogeneous to turbulent configurations. A Sedov-Taylor explosion in a homogeneous magnetic field generates an oblate ellipsoidal shock surface due to the slower propagating blast wave in the direction of the magnetic field. This is because of the efficient cosmic ray (CR) production in the quasi-parallel polar cap regions, which softens the equation of state and increases the compressibility of the post-shock gas. We find that the solution remains self-similar because the ellipticity of the propagating blast wave stays constant in time. This enables us to derive an effective ratio of specific heats for a composite of thermal gas and CRs as a function of the maximum acceleration efficiency. We finally discuss the behavior of supernova remnants expanding into a turbulent magnetic field with varying coherence lengths. For a maximum CR acceleration efficiency of about 15 per cent at quasi-parallel shocks (as suggested by kinetic plasma simulations), we find an average efficiency of about 5 per cent, independent of the assumed magnetic coherence length.


title: Bubble cloud dynamics in an ultrasound field

abstract: The dynamics of bubble clouds induced by high-intensity focused ultrasound are investigated in a regime where the cloud size is similar to the ultrasound wavelength. High-speed images show that the cloud is asymmetrical; the bubbles nearest the source grow to a larger radius than the distal ones. Similar structures of bubble clouds are observed in numerical simulations that mimic the laboratory experiment. To elucidate the structure, a parametric study is conducted for plane ultrasound waves with various amplitudes and diffuse clouds with different initial void fractions. Based on an analysis of the kinetic energy of liquid induced by bubble oscillations, a new scaling parameter is introduced to characterize the dynamics. The new parameter generalizes the cloud interaction parameter originally introduced by d'Agostino and Brennen (1989). The dynamic interaction parameter controls the energy localization and consequent anisoptropy of the cloud. Moreover, the amplitude of the far-field, bubble-scattered acoustics is likewise correlated with the proposed parameter. Findings of the present study not only shed light on the physics of cloud cavitation, but may also be of use to quantification of the effects of cavitation on outcomes of ultrasound therapies including HIFU-based lithotripsy.


title: Intrinsic timing jitter and latency in superconducting single photon nanowire detectors

abstract: We analyze the origin of the intrinsic timing jitter in superconducting nanowire single photon detectors (SNSPDs) in terms of fluctuations in the latency of the detector response, which is determined by the microscopic physics of the photon detection process. We demonstrate that fluctuations in the physical parameters which determine the latency give rise to the intrinsic timing jitter. We develop a general description of latency by introducing the explicit time dependence of the internal detection efficiency. By considering the dynamic Fano fluctuations together with static spatial inhomogeneities, we study the details of the connection between latency and timing jitter. We develop both a simple phenomenological model and a more general microscopic model of detector latency and timing jitter based on the solution of the generalized time-dependent Ginzburg-Landau equations for the 1D hotbelt geometry. While the analytical model is sufficient for qualitative interpretation of recent data, the general approach establishes the framework for a quantitative analysis of detector latency and the fundamental limits of intrinsic timing jitter. These theoretical advances can be used to interpret the results of recent experiments measuring the dependence of detection latency and timing jitter on photon energy to the few-picosecond level.


title: Wieferich Primes and a mod $p$ Leopoldt Conjecture

abstract: We consider questions in Galois cohomology which arise by considering mod $p$ Galois representations arising from automorphic forms. We consider a Galois cohomological analog for the standard heuristics about the distribution of Wieferich primes, i.e. prime $p$ such that $2^{p-1}$ is 1 mod $p^2$. Our analog relates to asking if in a compatible system of Galois representations, for almost all primes $p$, the residual mod $p$ representation arising from it has unobstructed deformation theory. This analog leads in particular to formulating a mod $p$ analog for almost all primes $p$ of the classical Leopoldt conjecture, which has been considered previously by G. Gras. Leopoldt conjectured that for a number field $F$, and a prime $p$, the $p$-adic regulator $R_{F,p}$ is non-zero. The mod $p$ analog is that for a fixed number field $F$, for almost all primes $p$, the $p$-adic regulator $R_{F,p}$ is a unit at $p$.