title: BMI: A Behavior Measurement Indicator for Fuel Poverty Using Aggregated Load Readings from Smart Meters

abstract: Fuel poverty affects between 50 and 125 million households in Europe and is a significant issue for both developed and developing countries globally. This means that fuel poor residents are unable to adequately warm their home and run the necessary energy services needed for lighting, cooking, hot water, and electrical appliances. The problem is complex but is typically caused by three factors; low income, high energy costs, and energy inefficient homes. In the United Kingdom (UK), 4 million families are currently living in fuel poverty. Those in series financial difficulty are either forced to self-disconnect or have their services terminated by energy providers. Fuel poverty contributed to 10,000 reported deaths in England in the winter of 2016-2107 due to homes being cold. While it is recognized by governments as a social, public health and environmental policy issue, the European Union (EU) has failed to provide a common definition of fuel poverty or a conventional set of indicators to measure it. This chapter discusses current fuel poverty strategies across the EU and proposes a new and foundational behavior measurement indicator designed to directly assess and monitor fuel poverty risks in households using smart meters, Consumer Access Device (CAD) data and machine learning. By detecting Activities of Daily Living (ADLS) through household appliance usage, it is possible to spot the early signs of financial difficulty and identify when support packages are required.


title: MajorityNets: BNNs Utilising Approximate Popcount for Improved Efficiency

abstract: Binarized neural networks (BNNs) have shown exciting potential for utilising neural networks in embedded implementations where area, energy and latency constraints are paramount. With BNNs, multiply-accumulate (MAC) operations can be simplified to XnorPopcount operations, leading to massive reductions in both memory and computation resources. Furthermore, multiple efficient implementations of BNNs have been reported on field-programmable gate array (FPGA) implementations. This paper proposes a smaller, faster, more energy-efficient approximate replacement for the XnorPopcountoperation, called XNorMaj, inspired by state-of-the-art FPGAlook-up table schemes which benefit FPGA implementations. Weshow that XNorMaj is up to 2x more resource-efficient than the XnorPopcount operation. While the XNorMaj operation has a minor detrimental impact on accuracy, the resource savings enable us to use larger networks to recover the loss.


title: Haupt's theorem for strata of abelian differentials

abstract: Let S be a closed topological surface. Haupt's theorem provides necessary and sufficient conditions for a complex-valued character of the first integer homology group of S to be realized by integration against a complex-valued 1-form that is holomorphic with respect to some complex structure on S. We prove a refinement of this theorem that takes into account the divisor data of the 1-form.


title: Chemistry on quantum computers with virtual quantum subspace expansion

abstract: Simulating chemical systems on quantum computers has been limited to a few electrons in a minimal basis. We demonstrate experimentally that the virtual quantum subspace expansion [Phys. Rev. X 10, 011004 (2020)] can achieve full basis accuracy for hydrogen and lithium dimers, comparable to simulations requiring twenty or more qubits. We developed an approach to minimize the impact of experimental noise on the stability of the generalized eigenvalue problem, a crucial component of the quantum algorithm. In addition, we were able to obtain an accurate potential energy curve for the nitrogen dimer in a quantum simulation on a classical computer.


title: The estimation error of general first order methods

abstract: Modern large-scale statistical models require to estimate thousands to millions of parameters. This is often accomplished by iterative algorithms such as gradient descent, projected gradient descent or their accelerated versions. What are the fundamental limits to these approaches? This question is well understood from an optimization viewpoint when the underlying objective is convex. Work in this area characterizes the gap to global optimality as a function of the number of iterations. However, these results have only indirect implications in terms of the gap to statistical optimality. Here we consider two families of high-dimensional estimation problems: high-dimensional regression and low-rank matrix estimation, and introduce a class of `general first order methods' that aim at efficiently estimating the underlying parameters. This class of algorithms is broad enough to include classical first order optimization (for convex and non-convex objectives), but also other types of algorithms. Under a random design assumption, we derive lower bounds on the estimation error that hold in the high-dimensional asymptotics in which both the number of observations and the number of parameters diverge. These lower bounds are optimal in the sense that there exist algorithms whose estimation error matches the lower bounds up to asymptotically negligible terms. We illustrate our general results through applications to sparse phase retrieval and sparse principal component analysis.


title: Fermi liquid theory and divergences of the two-particle irreducible vertex in the periodic Anderson lattice

abstract: Here we analyze the divergences of the irreducible vertex function in dynamical mean field theory, which may indicate either a non-physical breakdown of the perturbation theory or a response to some physical phenomenon. To investigate this question, we construct a quasiparticle vertex from the diverging irreducible vertex functions. This vertex describes the scattering between quasiparticles and quasiholes in a Fermi liquid. We show that the quasparticle vertex \textit{does not} diverge in the charge channel, wherein the irreducible vertex \textit{does} diverge; and we show that the quasiparticle vertex does diverge in the spin channel, wherein the irreducible vertex does not diverge. This divergence occurs at the Mott transition wherein the Fermi liquid theory breaks down. Both Hubbard and Anderson lattices are investigated. In general, our results support that the divergences of the irreducible vertex function do not indicate a non-physical failure of the perturbation theory. Instead, the divergences are the mathematical consequence of inverting a matrix (the local charge susceptibility) which accumulates increasingly negative diagonal elements as the Hubbard interaction suppresses charge fluctuations. Indeed, we find that the first symmetric divergences of the irreducible vertex in both Hubbard and Anderson lattices occurs near the maximum magnitude of the (negative) vertex-connected part of the charge susceptibility, which suppresses charge fluctuations.


title: Heterogeneous surface charge confining an electrolyte solution

abstract: The structure of dilute electrolyte solutions close to a surface carrying a spatially inhomogeneous surface charge distribution is investigated by means of classical density functional theory (DFT) within the approach of fundamental measure theory (FMT). For electrolyte solutions the influence of these inhomogeneities is particularly strong because the corresponding characteristic length scale is the Debye length, which is large compared to molecular sizes. Here a fully three-dimensional investigation is performed, which accounts explicitly for the solvent particles, and thus provides insight into effects caused by ion-solvent coupling. The present study introduces a versatile framework to analyze a broad range of types of surface charge heterogeneities even beyond the linear response regime. This reveals a sensitive dependence of the number density profiles of the fluid components and of the electrostatic potential on the magnitude of the charge as well as on details of the surface charge patterns at small scales.


title: BlueFlood: Concurrent Transmissions for Multi-Hop Bluetooth 5 -- Modeling and Evaluation

abstract: Bluetooth is an omnipresent technology, available on billions of devices today. While it has been traditionally limited to peer-to-peer communication and star networks, the recent Bluetooth Mesh standard extends it to multi-hop networking. In addition, the Bluetooth 5 standard introduces new modes to allow for increased reliability. In this paper, we evaluate the feasibility of concurrent transmissions (CT) in Bluetooth via modeling and controlled experiments and then devise an efficient network-wide data dissemination protocol, BlueFlood, based on CT for multi-hop Bluetooth networks. First, we model and analyze how CT distorts the received waveform and characterize the Bit Error Rate of a Frequency-Shift Keying receiver to show that CT is feasible over Bluetooth. Second, we verify our analytic results with a controlled experimental study of CT over Bluetooth PHY. Third, we present BlueFlood, a fast and efficient network-wide data dissemination in multi-hop Bluetooth networks. In our experimental evaluation, in two testbeds deployed in university buildings, we show that BlueFlood achieves 99.9% end-to-end delivery ratio with a duty-cycle of 0.4% for periodic dissemination of advertising packets of 38 bytes with 200 milliseconds intervals at 2 Mbps. Moreover, we show that BlueFlood can be received by off-the-shelf devices such as smartphones, paving a seamless integration with existing technologies.


title: Two types of alternating spin-$\frac12$ chains and their field-induced transitions in $\varepsilon$-LiVOPO$_4$

abstract: Thermodynamic properties, $^{31}$P nuclear magnetic resonance (NMR) measurements, and density-functional band-structure calculations for $\varepsilon$-LiVOPO$_4$ are reported. This quantum magnet features a singlet ground state and comprises two types of alternating spin-$\frac12$ chains that manifest themselves by the double maxima in the susceptibility and magnetic specific heat, and by the two-step magnetization process with an intermediate $\frac12$-plateau. From thermodynamic data and band-structure calculations, we estimate the leading couplings of $J_1\simeq 20$ K and $J_2\simeq 60$ K and the alternation ratios of $\alpha_1=J_1'/J_1\simeq 0.6$ and $\alpha_2=J_2'/J_2\simeq 0.3$ within the two chains, respectively. The zero-field spin gap $\Delta_0/k_{\rm B}\simeq 7.3$ K probed by thermodynamic and NMR measurements is caused by the $J_1$-$J_1'$ spin chains and can be closed in the applied field of $\mu_{0}H_{\rm c1}\simeq 5.6$ T, giving rise to a field-induced long-range order. The NMR data reveal predominant three-dimensional spin-spin correlations at low temperatures. Field-induced magnetic ordering transition observed above $H_{c1}$ is attributed to the Bose-Einstein condensation of triplons in the sublattice formed by the $J_1$-$J_1'$ chains with weaker exchange couplings.


title: Nonlinear Characterization of Tissue Viscoelasticity with Acoustoelastic Attenuation of Shear-Waves

abstract: Shear-wave elastography (SWE) measures shear-wave speed (SWS), which is related to the underlying shear modulus of soft tissue. SWE methods generally assume that soft tissue viscoelasticity is independent of mechanical loading, however, soft tissues are known to have viscoelasticity changing nonlinearly with pre-compression. Hence, characterization by SWS alone is insufficient, where nonlinear properties can be seen as confounders but may also be utilized as additional bio-markers. The viscoelastic nature of a medium is fully characterized by its storage and loss moduli, which are related to SWS and shear-wave attenuation (SWA). In this work, we study SWA characteristics as a function of applied strain to measure nonlinear viscoelastic parameters in soft tissues. For this purpose, we apply incremental quasi-static compression on the samples while measuring SWS and SWA, from which we derive storage and loss moduli to estimate nonlinear viscoelastic parameters as a function of applied strain using acoustoelasticity (AE) theory. Results from tissue-mimicking phantoms with varying oil percentages and ex-vivo porcine liver experiments demonstrate the feasibility of the proposed approach. In both these experiments, SWA was observed to decrease with applied strain. For 10 % compression in ex-vivo livers, shear-wave attenuation decreased on average by 28 % (93 Np/m), while SWS increased on average by 20 % (0.26 m/s).


title: Deep Reinforcement Learning for FlipIt Security Game

abstract: Reinforcement learning has shown much success in games such as chess, backgammon and Go. However, in most of these games, agents have full knowledge of the environment at all times. In this paper, we describe a deep learning model in which agents successfully adapt to different classes of opponents and learn the optimal counter-strategy using reinforcement learning in a game under partial observability. We apply our model to FlipIt, a two-player security game in which both players, the attacker and the defender, compete for ownership of a shared resource and only receive information on the current state of the game upon making a move. Our model is a deep neural network combined with Q-learning and is trained to maximize the defender's time of ownership of the resource. Despite the noisy information, our model successfully learns a cost-effective counter-strategy outperforming its opponent's strategies and shows the advantages of the use of deep reinforcement learning in game theoretic scenarios. We also extend FlipIt to a larger action-spaced game with the introduction of a new lower-cost move and generalize the model to $n$-player FlipIt.


title: Addressing $R_{D^{(*)}}$, $R_{K^{(*)}}$, muon $g-2$ and ANITA anomalies in a minimal $R$-parity violating supersymmetric framework

abstract: We analyze the recent hints of lepton flavor universality violation in both charged-current and neutral-current rare decays of $B$-mesons in an $R$-parity violating supersymmetric scenario. Motivated by simplicity and minimality, we had earlier postulated the third-generation superpartners to be the lightest (calling the scenario "RPV3") and explicitly showed that it preserves gauge coupling unification and of course has the usual attribute of naturally addressing the Higgs radiative stability. Here we show that both $R_{D^{(*)}}$ and $R_{K^{(*)}}$ flavor anomalies can be addressed in this RPV3 framework. Interestingly, this scenario may also be able to accommodate two other seemingly disparate anomalies, namely, the longstanding discrepancy in the muon $(g-2)$, as well as the recent anomalous upgoing ultra-high energy ANITA events. Based on symmetry arguments, we consider three different benchmark points for the relevant RPV3 couplings and carve out the regions of parameter space where all (or some) of these anomalies can be simultaneously explained. We find it remarkable that such overlap regions exist, given the plethora of precision low-energy and high-energy experimental constraints on the minimal model parameter space. The third-generation superpartners needed in this theoretical construction are all in the 1-10 TeV range, accessible at the LHC and/or next-generation hadron collider. We also discuss some testable predictions for the lepton-flavor-violating decays of the tau-lepton and $B$-mesons for the current and future $B$-physics experiments, such as LHCb and Belle II. Complementary tests of the flavor anomalies in the high-$p_T$ regime in collider experiments such as the LHC are also discussed.


title: First Order Methods take Exponential Time to Converge to Global Minimizers of Non-Convex Functions

abstract: Machine learning algorithms typically perform optimization over a class of non-convex functions. In this work, we provide bounds on the fundamental hardness of identifying the global minimizer of a non convex function. Specifically, we design a family of parametrized non-convex functions and employ statistical lower bounds for parameter estimation. We show that the parameter estimation problem is equivalent to the problem of function identification in the given family. We then claim that non convex optimization is at least as hard as function identification. Jointly, we prove that any first order method can take exponential time to converge to a global minimizer.


title: Quantum to classical crossover of Floquet engineering in correlated quantum systems

abstract: Light-matter coupling involving classical and quantum light offers a wide range of possibilities to tune the electronic properties of correlated quantum materials. Two paradigmatic results are the dynamical localization of electrons and the ultrafast control of spin dynamics, which have been discussed within classical Floquet engineering and in the deep quantum regime where vacuum fluctuations modify the properties of materials. Here we discuss how these two extreme limits are interpolated by a cavity which is driven to the excited states. In particular, this is achieved by formulating a Schrieffer-Wolff transformation for the cavity-coupled system, which is mathematically analogous to its Floquet counterpart. Some of the extraordinary results of Floquet-engineering, such as the sign reversal of the exchange interaction or electronic tunneling, which are not obtained by coupling to a dark cavity, can already be realized with a single-photon state (no coherent states are needed). The analytic results are verified and extended with numerical simulations on a two-site Hubbard model coupled to a driven cavity mode. Our results generalize the well-established Floquet-engineering of correlated electrons to the regime of quantum light. It opens up a new pathway of controlling properties of quantum materials with high tunability and low energy dissipation.


title: Applying Tensor Decomposition to image for Robustness against Adversarial Attack

abstract: Nowadays the deep learning technology is growing faster and shows dramatic performance in computer vision areas. However, it turns out a deep learning based model is highly vulnerable to some small perturbation called an adversarial attack. It can easily fool the deep learning model by adding small perturbations. On the other hand, tensor decomposition method widely uses for compressing the tensor data, including data matrix, image, etc. In this paper, we suggest combining tensor decomposition for defending the model against adversarial example. We verify this idea is simple and effective to resist adversarial attack. In addition, this method rarely degrades the original performance of clean data. We experiment on MNIST, CIFAR10 and ImageNet data and show our method robust on state-of-the-art attack methods.


title: Social Welfare and Price of Anarchy in Preemptive Priority Queues

abstract: Consider an unobservable $M|G|1$ queue with preemptive-resume scheduling and two priority classes. Customers are strategic and may join the premium class for a fee. We analyze the resulting equilibrium outcomes, equilibrium stability, and social welfare. We find that for service distributions with coefficient of variation greater than 1, there exists a unique and stable mixed equilibrium at low loads. We also establish a tight bound on the price of anarchy, which is $4/3$.


title: The Implicit and Explicit Regularization Effects of Dropout

abstract: Dropout is a widely-used regularization technique, often required to obtain state-of-the-art for a number of architectures. This work demonstrates that dropout introduces two distinct but entangled regularization effects: an explicit effect (also studied in prior work) which occurs since dropout modifies the expected training objective, and, perhaps surprisingly, an additional implicit effect from the stochasticity in the dropout training update. This implicit regularization effect is analogous to the effect of stochasticity in small mini-batch stochastic gradient descent. We disentangle these two effects through controlled experiments. We then derive analytic simplifications which characterize each effect in terms of the derivatives of the model and the loss, for deep neural networks. We demonstrate these simplified, analytic regularizers accurately capture the important aspects of dropout, showing they faithfully replace dropout in practice.


title: Spatial Kasner solution and an infinite slab with constant energy density

abstract: We study the solutions of the Einstein equations in the presence of a thick infinite slab with constant energy density. When there is an isotropy in the plane of the slab, we find an explicit exact solution that matches with the Rindler and Weyl-Levi-Civita spacetimes outside the slab. We also show that there are solutions that can be matched with general anisotropic Kasner spacetime outside the slab. In any case, it is impossible to avoid the presence of the Kasner type singularities in contrast to the well-known case of spherical symmetry, where by matching the internal Schwarzschild solution with the external one, the singularity in the center of coordinates can be eliminated.


title: Multivariate Haar systems in Besov function spaces

abstract: We determine all cases for which the $d$-dimensional Haar wavelet system $H^d$ on the unit cube $I^d$ is a conditional or unconditional Schauder basis in the classical isotropic Besov function spaces ${B}_{p,q,1}^s(I^d)$, $0<p,q<\infty$, $0\le s < 1/p$, defined in terms of first-order $L_p$ moduli of smoothness. We obtain similar results for the tensor-product Haar system $\tilde{H}^d$, and characterize the parameter range for which the dual of ${B}_{p,q,1}^s(I^d)$ is trivial for $0<p<1$.


title: Homotopical rigidity of the pre-Lie operad

abstract: We show that the celebrated operad of pre-Lie algebras is very rigid: it has no "non-obvious" degrees of freedom from either of the three points of view: deformations of maps to and from the "three graces of operad theory", homotopy automorphisms, and operadic twisting. Examining the latter, it is possible to answer two questions of Markl from 2005, including a Lie-theoretic version of the Deligne conjecture.


title: MMC-Based Distributed Maximum Power Point Tracking for Photovoltaic Systems

abstract: This paper proposes a novel topology for grid connected photovoltaic (PV) system based on modular multilevel converter (MMC). In this topology, a PV array is connected to capacitors of each submodule (SM) of the MMC through a DC-DC boost converter with maximum power point tracking (MPPT) control. This topology will maximize the efficiency of the system in the case of partial shading conditions, as it can regulate the SM capacitor voltages independently from each other to realize distributed MPPT. A model predictive control is used to track the AC output current, balance the SMs capacitor voltages, and to mitigate the circulating current. The proposed PV generation topology with 7 level MMC system validity has been verified by simulations via MATLAB/Simulink toolbox under normal operation, partial shading and PV array failure.


title: Automatic Perturbation Analysis for Scalable Certified Robustness and Beyond

abstract: Linear relaxation based perturbation analysis (LiRPA) for neural networks, which computes provable linear bounds of output neurons given a certain amount of input perturbation, has become a core component in robustness verification and certified defense. The majority of LiRPA-based methods focus on simple feed-forward networks and need particular manual derivations and implementations when extended to other architectures. In this paper, we develop an automatic framework to enable perturbation analysis on any neural network structures, by generalizing existing LiRPA algorithms such as CROWN to operate on general computational graphs. The flexibility, differentiability and ease of use of our framework allow us to obtain state-of-the-art results on LiRPA based certified defense on fairly complicated networks like DenseNet, ResNeXt and Transformer that are not supported by prior works. Our framework also enables loss fusion, a technique that significantly reduces the computational complexity of LiRPA for certified defense. For the first time, we demonstrate LiRPA based certified defense on Tiny ImageNet and Downscaled ImageNet where previous approaches cannot scale to due to the relatively large number of classes. Our work also yields an open-source library for the community to apply LiRPA to areas beyond certified defense without much LiRPA expertise, e.g., we create a neural network with a probably flat optimization landscape by applying LiRPA to network parameters. Our opensource library is available at https://github.com/KaidiXu/auto_LiRPA.


title: VegasFlow: accelerating Monte Carlo simulation across multiple hardware platforms

abstract: We present VegasFlow, a new software for fast evaluation of high dimensional integrals based on Monte Carlo integration techniques designed for platforms with hardware accelerators. The growing complexity of calculations and simulations in many areas of science have been accompanied by advances in the computational tools which have helped their developments. VegasFlow enables developers to delegate all complicated aspects of hardware or platform implementation to the library so they can focus on the problem at hand. This software is inspired on the Vegas algorithm, ubiquitous in the particle physics community as the driver of cross section integration, and based on Google's powerful TensorFlow library. We benchmark the performance of this library on many different consumer and professional grade GPUs and CPUs.


title: Building the largest spectroscopic sample of ultra-compact massive galaxies with the Kilo Degree Survey

abstract: Ultra-compact massive galaxies UCMGs, i.e. galaxies with stellar masses $M_{*} > 8 \times 10^{10} M_{\odot}$ and effective radii $R_{e} < 1.5$ kpc, are very rare systems, in particular at low and intermediate redshifts. Their origin as well as their number density across cosmic time are still under scrutiny, especially because of the paucity of spectroscopically confirmed samples. We have started a systematic census of UCMG candidates within the ESO Kilo Degree Survey, together with a large spectroscopic follow-up campaign to build the largest possible sample of confirmed UCMGs. This is the third paper of the series and the second based on the spectroscopic follow-up program. Here, we present photometrical and structural parameters of 33 new candidates at redshifts $0.15 \lesssim z \lesssim 0.5$ and confirm 19 of them as UCMGs, based on their nominal spectroscopically inferred $M_{*}$ and $R_{e}$. This corresponds to a success rate of $\sim 58\%$, nicely consistent with our previous findings. The addition of these 19 newly confirmed objects, allows us to fully assess the systematics on the system selection, and finally reduce the number density uncertainties. Moreover, putting together the results from our current and past observational campaigns and some literature data, we build the largest sample of UCMGs ever collected, comprising 92 spectroscopically confirmed objects at $0.1 \lesssim z \lesssim 0.5$. This number raises to 116, allowing for a $3\sigma$ tolerance on the $M_{*}$ and $R_{e}$ thresholds for the UCMG definition. For all these galaxies we have estimated the velocity dispersion values at the effective radii which have been used to derive a preliminary mass-velocity dispersion correlation.


title: Stabilization of uni-directional water-wave trains over an uneven bottom

abstract: We study the evolution of nonlinear surface gravity water-wave packets developing from modulational instability over an uneven bottom. A nonlinear Schr\"odinger equation (NLSE) with coefficients varying in space along propagation is used as a reference model. Based on a low-dimensional approximation obtained by considering only three complex harmonic modes, we discuss how to stabilize a one-dimensional pattern in the form of train of large peaks sitting on a background and propagating over a significant distance. Our approach is based on a gradual depth variation, while its conceptual framework is the theory of autoresonance in nonlinear systems and leads to a quasi-frozen state. Three main stages are identified: amplification from small sideband amplitudes, separatrix crossing, and adiabatic conversion to orbits oscillating around an elliptic fixed point. Analytical estimates on the three stages are obtained from the low-dimensional approximation and validated by NLSE simulations. Our result will contribute to understand dynamical stabilization of nonlinear wave packets and the persistence of large undulatory events in hydrodynamics and other nonlinear dispersive media.


title: Porous media equations with multiplicative space-time white noise

abstract: The existence of martingale solutions for stochastic porous media equations driven by nonlinear multiplicative space-time white noise is established in spatial dimension one. The Stroock-Varopoulos inequality is identified as a key tool in the derivation of the corresponding estimates.


title: Predicting excited states from ground state wavefunction by supervised quantum machine learning

abstract: Excited states of molecules lie in the heart of photochemistry and chemical reactions. The recent development in quantum computational chemistry leads to inventions of a variety of algorithms that calculate the excited states of molecules on near-term quantum computers, but they require more computational burdens than the algorithms for calculating the ground states. In this study, we propose a scheme of supervised quantum machine learning which predicts the excited-state properties of molecules only from their ground state wavefunction resulting in reducing the computational cost for calculating the excited states. Our model is comprised of a quantum reservoir and a classical machine learning unit which processes the measurement results of single-qubit Pauli operators with the output state from the reservoir. The quantum reservoir effectively transforms the single-qubit operators into complicated multi-qubit ones which contain essential information of the system, so that the classical machine learning unit may decode them appropriately. The number of runs for quantum computers is saved by training only the classical machine learning unit, and the whole model requires modest resources of quantum hardware that may be implemented in current experiments. We illustrate the predictive ability of our model by numerical simulations for small molecules with and without noise inevitable in near-term quantum computers. The results show that our scheme well reproduces the first and second excitation energies as well as the transition dipole moment between the ground states and excited states only from the ground state as an input. We expect our contribution will enhance the applications of quantum computers in the study of quantum chemistry and quantum materials.


title: Cities as they could be: Artificial Life and Urban Systems

abstract: The metaphor of cities as organisms has a long history in urban planning, and a few urban modeling approaches have explicitly been linked to Artificial Life. We propose in that paper to explore the extent of Artificial Life and Artificial Intelligence application to urban issues, by constructing and exploring a citation network of around 225,000 papers. It shows that most of the literature is indeed application of methodologies and a rather strong modularity of approaches. We finally develop ALife concepts which have a strong potential for the development of new urban theories.


title: Machine Learning ${\cal N}=8, D=5$ Gauged Supergravity

abstract: Type IIB string theory on a 5-sphere gives rise to ${\cal N}=8, SO(6)$ gauged supergravity in five dimensions. Motivated by the fact that this is the context of the most widely studied example of the AdS/CFT correspondence, we undertake an investigation of its critical points. The scalar manifold is an $E_{6(6)}/USp(8)$ coset, and the challenge is that it is 42-dimensional. We take a Machine Learning approach to the problem using TensorFlow, and this results in a substantial increase in the number of known critical points. Our list of 32 critical points contains all five of the previously known ones, including an ${\cal N}=2$ supersymmetric point identified by Khavaev, Pilch and Warner.


title: A Self-Tuning Actor-Critic Algorithm

abstract: Reinforcement learning algorithms are highly sensitive to the choice of hyperparameters, typically requiring significant manual effort to identify hyperparameters that perform well on a new domain. In this paper, we take a step towards addressing this issue by using metagradients to automatically adapt hyperparameters online by meta-gradient descent (Xu et al., 2018). We apply our algorithm, Self-Tuning Actor-Critic (STAC), to self-tune all the differentiable hyperparameters of an actor-critic loss function, to discover auxiliary tasks, and to improve off-policy learning using a novel leaky V-trace operator. STAC is simple to use, sample efficient and does not require a significant increase in compute. Ablative studies show that the overall performance of STAC improved as we adapt more hyperparameters. When applied to the Arcade Learning Environment (Bellemare et al. 2012), STAC improved the median human normalized score in 200M steps from 243% to 364%. When applied to the DM Control suite (Tassa et al., 2018), STAC improved the mean score in 30M steps from 217 to 389 when learning with features, from 108 to 202 when learning from pixels, and from 195 to 295 in the Real-World Reinforcement Learning Challenge (Dulac-Arnold et al., 2020).


title: Coincidence of dimensions in closed ordered differential fields

abstract: Let $\mathcal K=\langle\mathcal R, \delta\rangle$ be a closed ordered differential field, in the sense of M. Singer, and $C$ its field of constants. In this note, we prove that, for sets definable in the pair $\mathcal M=\langle \mathcal R, C\rangle$, the $\delta$-dimension and the large dimension coincide. As an application, we characterize the definable sets in $\mathcal K$ that are internal to $C$ as those sets that are definable in $\mathcal M$ and have $\delta$-dimension $0$. We further show that, for sets definable in $\mathcal K$, having $\delta$-dimension $0$ does not generally imply co-analyzability in $C$ (in contrast to the case of transseries). We also point out that the coincidence of dimensions also holds in the context of differentially closed fields and in the context of transseries.


title: Optical Magnetism and Huygens' Surfaces in Arrays of Atoms Induced by Cooperative Responses

abstract: By utilizing strong optical resonant interactions in arrays of atoms with electric dipole transitions, we show how to synthesize collective optical responses that correspond to those formed by arrays of magnetic dipoles and other multipoles. Optically active magnetism with the strength comparable with that of electric dipole transitions is achieved in collective excitation eigenmodes of the array. By controlling the atomic level shifts, an array of spectrally overlapping, crossed electric and magnetic dipoles can be excited, providing a physical realization of a nearly-reflectionless quantum Huygens' surface with the full $2\pi$ phase control of the transmitted light that allows for extreme wavefront engineering even at a single photon level. We illustrate this by transforming a plane wave into a vortex beam.

categories: physics.atom-ph cond-mat.quant-gas physics.optics quant-ph


title: Compensated Isocurvature Perturbations in the Galaxy Power Spectrum

abstract: We investigate the potential of the galaxy power spectrum to constrain compensated isocurvature perturbations (CIPs), primordial fluctuations in the baryon density that are compensated by fluctuations in CDM density to ensure an unperturbed total matter density. We show that CIPs contribute to the galaxy overdensity at linear order, and if they are close to scale-invariant, their effects are nearly perfectly degenerate with the local PNG parameter $f_{\rm nl}$ if they correlate with the adiabatic perturbations. This degeneracy can however be broken by analyzing multiple galaxy samples with different bias parameters, or by taking CMB priors on $f_{\rm nl}$ into account. Parametrizing the amplitude of the CIP power spectrum as $P_{\sigma\sigma} = A^2P_{\mathcal{R}\mathcal{R}}$ (where $P_{\mathcal{R}\mathcal{R}}$ is the adiabatic power spectrum) we find, for a number of fiducial galaxy samples in a simplified forecast setup, that constraints on $A$, relative to those on $f_{\rm nl}$, of order $\sigma_{A}/\sigma_{f_{\rm nl}} \approx 1-2$ are achievable for CIPs correlated with adiabatic perturbations, and $\sigma_{A}/\sigma_{f_{\rm nl}} \approx 5$ for the uncorrelated case. These values are independent of survey volume, and suggest that current galaxy data are already able to improve significantly on the tightest existing constraints on CIPs from the CMB. Future galaxy surveys that aim to achieve $\sigma_{f_{\rm nl}} \sim 1$ have the potential to place even stronger bounds on CIPs.


title: Non-Liquid Cellular States

abstract: The existence of quantum non-liquid states and fracton orders, both gapped and gapless states, challenges our understanding of phases of entangled matter. We generalize the cellular topological states to liquid or non-liquid cellular states. We propose a mechanism to construct more general non-abelian states by gluing gauge-symmetry-breaking vs gauge-symmetry-extension interfaces as extended defects in a cellular network, including defects of higher-symmetries, in any dimension. Our approach also naturally incorporates the anyonic particle/string condensations and composite string (related to particle-string or p-string)/membrane condensations. This approach shows gluing the familiar extended topological quantum field theory or conformal field theory data via topology, geometry, and renormalization consistency criteria (via certain modified group cohomology or cobordism theory data) in a tensor network can still guide us to analyze the non-liquid states. (Part of the abelian construction can be understood from the K-matrix Chern-Simons theory approach and the coupled-layer-by-junction constructions.) This approach may also lead us toward a unifying framework for quantum systems of both higher-symmetries and sub-system/sub-dimensional symmetries.

categories: cond-mat.str-el hep-th math-ph math.MP quant-ph


title: The mathematics of Bitcoin

abstract: We survey recent results on the mathematical stability of Bitcoin protocol. Profitability and probability of a double spend are estimated in closed form with classical special functions. The stability of Bitcoin mining rules is analyzed and several theorems are proved using martingale and combinatorics techniques. In particular, the empirical observation of the stability of the Bitcoin protocol is proved. This survey article on the mathematics of Bitcoin is published by the Newsletter of the European Mathematical Society, vol.115, 2020, p.31-37. Continuation of arXiv:1601.05254 (EMS Newsletter, 100, 2016 p.32).


title: A Nonlinear Extension of Korovkin's Theorem

abstract: In this paper we extend the classical Korovkin theorems to the framework of comonotone additive, sublinear and monotone operators. Based on the theory of Choquet capacities, several concrete examples illustrating our results are also discussed.


title: Toward Active and Passive Confidentiality Attacks On Cryptocurrency Off-Chain Networks

abstract: Cryptocurrency off-chain networks such as Lightning (e.g., Bitcoin) or Raiden (e.g., Ethereum) aim to increase the scalability of traditional on-chain transactions. To support nodes in learning about possible paths to route their transactions, these networks need to provide gossip and probing mechanisms. This paper explores whether these mechanisms may be exploited to infer sensitive information about the flow of transactions, and eventually harm privacy. In particular, we identify two threats, related to an active and a passive adversary. The first is a probing attack: here the adversary aims to detect the maximum amount which is transferable in a given direction over a target channel by actively probing it and differentiating the response messages it receives. The second is a timing attack: the adversary discovers how close the destination of a routed payment actually is, by acting as a passive man-in-the middle and analyzing the time deltas between sent messages and their corresponding responses. We then analyze the limitations of these attacks and propose remediations for scenarios in which they are able to produce accurate results.


title: Volterra-Choquet nonlinear operators

abstract: In this paper we study to what extend some properties of the classical linear Volterra operators could be transferred to the nonlinear Volterra-Choquet operators, obtained by replacing the classical linear integral with respect to the Lebesgue measure, by the nonlinear Choquet integral with respect to a nonadditive set function. Compactness, Lipschitz and cyclicity properties are studied.


title: Elliptic Solutions for Higher Order KdV Equations

abstract: We study higher order KdV equations from the GL(2,$\mathbb{R}$) $\cong$ SO(2,1) Lie group point of view. We find elliptic solutions of higher order KdV equations up to the ninth order. We argue that the main structure of the trigonometric/hyperbolic/elliptic $N$-soliton solutions for higher order KdV equations is the same as that of the original KdV equation. Pointing out that the difference is only the time dependence, we find $N$-soliton solutions of higher order KdV equations can be constructed from those of the original KdV equation by properly replacing the time-dependence. We discuss that there always exist elliptic solutions for all higher order KdV equations.


title: Dynamic contrast in scanning microscopic OCT

abstract: While optical coherence tomography (OCT) provides a resolution down to 1 micrometer it has difficulties to visualize cellular structures due to a lack of scattering contrast. By evaluating signal fluctuations, a significant contrast enhancement was demonstrated using time-domain full-field OCT (FF-OCT), which makes cellular and subcellular structures visible. The putative cause of the dynamic OCT signal is ATP-dependent motion of cellular structures in a sub-micrometer range, which provides histology-like contrast. Here we demonstrate dynamic contrast with a scanning frequency-domain OCT (FD-OCT). Given the inherent sectional imaging geometry, scanning FD-OCT provides depth-resolved images across tissue layers, a perspective known from histopathology, much faster and more efficiently than FF-OCT. Both, shorter acquisition times and tomographic depth-sectioning reduce the sensitivity of dynamic contrast for bulk tissue motion artifacts and simplify their correction in post-processing. The implementation of dynamic contrast makes microscopic FD-OCT a promising tool for histological analysis of unstained tissues.


title: Generating EEG features from Acoustic features

abstract: In this paper we demonstrate predicting electroencephalograpgy (EEG) features from acoustic features using recurrent neural network (RNN) based regression model and generative adversarial network (GAN). We predict various types of EEG features from acoustic features. We compare our results with the previously studied problem on speech synthesis using EEG and our results demonstrate that EEG features can be generated from acoustic features with lower root mean square error (RMSE), normalized RMSE values compared to generating acoustic features from EEG features (ie: speech synthesis using EEG) when tested using the same data sets.


title: Reduction theory for connections over the formal punctured disc

abstract: We give a purely algebraic treatment of reduction theory for connections over the formal punctured disc. Our proofs apply to arbitrary connected linear algebraic groups over an algebraically closed field of characteristic 0. We also state and prove some new quantitative results.


title: Cell Mechanics Based Computational Classification of Red Blood Cells Via Machine Intelligence Applied to Morpho-Rheological Markers

abstract: Despite fluorescent cell-labelling being widely employed in biomedical studies, some of its drawbacks are inevitable, with unsuitable fluorescent probes or probes inducing a functional change being the main limitations. Consequently, the demand for and development of label-free methodologies to classify cells is strong and its impact on precision medicine is relevant. Towards this end, high-throughput techniques for cell mechanical phenotyping have been proposed to get a multidimensional biophysical characterization of single cells. With this motivation, our goal here is to investigate the extent to which an unsupervised machine learning methodology, which is applied exclusively on morpho-rheological markers obtained by real-time deformability and fluorescence cytometry (RT-FDC), can address the difficult task of providing label-free discrimination of reticulocytes from mature red blood cells. We focused on this problem, since the characterization of reticulocytes (their percentage and cellular features) in the blood is vital in multiple human disease conditions, especially bone-marrow disorders such as anemia and leukemia. Our approach reports promising label-free results in the classification of reticulocytes from mature red blood cells, and it represents a step forward in the development of high-throughput morpho-rheological-based methodologies for the computational categorization of single cells. Besides, our methodology can be an alternative but also a complementary method to integrate with existing cell-labelling techniques.


title: Distributed Momentum for Byzantine-resilient Learning

abstract: Momentum is a variant of gradient descent that has been proposed for its benefits on convergence. In a distributed setting, momentum can be implemented either at the server or the worker side. When the aggregation rule used by the server is linear, commutativity with addition makes both deployments equivalent. Robustness and privacy are however among motivations to abandon linear aggregation rules. In this work, we demonstrate the benefits on robustness of using momentum at the worker side. We first prove that computing momentum at the workers reduces the variance-norm ratio of the gradient estimation at the server, strengthening Byzantine resilient aggregation rules. We then provide an extensive experimental demonstration of the robustness effect of worker-side momentum on distributed SGD.


title: Controlled Online Optimization Learning (COOL): Finding the ground state of spin Hamiltonians with reinforcement learning

abstract: Reinforcement learning (RL) has become a proven method for optimizing a procedure for which success has been defined, but the specific actions needed to achieve it have not. We apply the so-called "black box" method of RL to what has been referred as the "black art" of simulated annealing (SA), demonstrating that an RL agent based on proximal policy optimization can, through experience alone, arrive at a temperature schedule that surpasses the performance of standard heuristic temperature schedules for two classes of Hamiltonians. When the system is initialized at a cool temperature, the RL agent learns to heat the system to "melt" it, and then slowly cool it in an effort to anneal to the ground state; if the system is initialized at a high temperature, the algorithm immediately cools the system. We investigate the performance of our RL-driven SA agent in generalizing to all Hamiltonians of a specific class; when trained on random Hamiltonians of nearest-neighbour spin glasses, the RL agent is able to control the SA process for other Hamiltonians, reaching the ground state with a higher probability than a simple linear annealing schedule. Furthermore, the scaling performance (with respect to system size) of the RL approach is far more favourable, achieving a performance improvement of one order of magnitude on L=14x14 systems. We demonstrate the robustness of the RL approach when the system operates in a "destructive observation" mode, an allusion to a quantum system where measurements destroy the state of the system. The success of the RL agent could have far-reaching impact, from classical optimization, to quantum annealing, to the simulation of physical systems.


title: High-throughput calculations of magnetic topological materials

abstract: The discoveries of intrinsically magnetic topological materials, including semimetals with a large anomalous Hall effect and axion insulators, have directed fundamental research in solid-state materials. Topological quantum chemistry has enabled the understanding of and the search for paramagnetic topological materials. Using magnetic topological indices obtained from magnetic topological quantum chemistry (MTQC), here we perform a high-throughput search for magnetic topological materials based on first-principles calculations. We use as our starting point the Magnetic Materials Database on the Bilbao Crystallographic Server, which contains more than 549 magnetic compounds with magnetic structures deduced from neutron-scattering experiments, and identify 130 enforced semimetals (for which the band crossings are implied by symmetry eigenvalues), and topological insulators. For each compound, we perform complete electronic structure calculations, which include complete topological phase diagrams using different values of the Hubbard potential. Using a custom code to find the magnetic co-representations of all bands in all magnetic space groups, we generate data to be fed into the algorithm of MTQC to determine the topology of each magnetic material. Several of these materials display previously unknown topological phases, including symmetry-indicated magnetic semimetals, three-dimensional anomalous Hall insulators and higher-order magnetic semimetals. We analyse topological trends in the materials under varying interactions: 60 per cent of the 130 topological materials have topologies sensitive to interactions, and the others have stable topologies under varying interactions. We provide a materials database for future experimental studies and open-source code for diagnosing topologies of magnetic materials.


title: The ALPINE-ALMA [CII] Survey: Size of Individual Star-Forming Galaxies at z=4-6 and their Extended Halo Structure

abstract: We present the physical extent of [CII] 158um line-emitting gas from 46 star-forming galaxies at z=4-6 from the ALMA Large Program to INvestigate CII at Early Times (ALPINE). Using exponential profile fits, we measure the effective radius of the [CII] line (r_e,[CII]) for individual galaxies and compare them with the rest-frame ultra-violet (UV) continuum (r_e,UV) from Hubble Space Telescope images. The effective radius r_e,[CII] exceeds r_e,UV by factors of ~2-3 and the ratio of r_e,[CII]/r_e,UV increases as a function of M_star. We do not find strong evidence that [CII] line, the rest-frame UV, and FIR continuum are always displaced over ~ 1-kpc scale from each other. We identify 30% of isolated ALPINE sources as having an extended [CII] component over 10-kpc scales detected at 4.1$\sigma$-10.9$\sigma$ beyond the size of rest-frame UV and far-infrared (FIR) continuum. One object has tentative rotating features up to ~10-kpc, where the 3D model fit shows the rotating [CII]-gas disk spread over 4 times larger than the rest-frame UV-emitting region. Galaxies with the extended [CII] line structure have high star-formation rate (SFR), stellar mass (M_star), low Lya equivalent-width, and more blue-shifted (red-shifted) rest-frame UV metal absorption (Lya line), as compared to galaxies without such extended [CII] structures. Although we cannot rule out the possibility that a selection bias towards luminous objects may be responsible for such trends, the star-formation driven outflow also explains all these trends. Deeper observations are essential to test whether the extended [CII] line structures are ubiquitous to high-z star-forming galaxies.


title: New Limits on Coloured Three Jet Resonances

abstract: We consider experimental limits on colour triplet fermions that decay dominantly to three jets via a scalar mediator that can be on- or off-shell. These fermions arise in top-partner models that can solve the hierarchy problem, and limits on this scenario are weaker than those on traditional top-partner models because of the messy all-hadronic final state with significant backgrounds. We do find, however, that while there are no dedicated searches for this scenario, especially in case of an on-shell mediator, the suite of LHC all-hadronic searches still constrains a significant portion of the parameter space. In particular, we find that searches for pair production of di-jet and tri-jet resonances are complementary, covering different regions of parameter space. We also find that if the final state is rich in b-jets, current limits do not change significantly relative to the scenario with all light jets, and we describe how modifications of current search strategies can improve limits in that case.


title: Double X/Peanut Structures in Barred Galaxies -- Insights from an $N$--body Simulation

abstract: Boxy, peanut- or X-shaped "bulges" are observed in a large fraction of barred galaxies viewed in, or close to, edge-on projection, as well as in the Milky Way. They are the product of dynamical instabilities occurring in stellar bars, which cause the latter to buckle and thicken vertically. Recent studies have found nearby galaxies that harbour two such features arising at different radial scales, in a nested configuration. In this paper we explore the formation of such double peanuts, using a collisionless N-body simulation of a pure disc evolving in isolation within a live dark matter halo, which we analyse in a completely analogous way to observations of real galaxies. In the simulation we find a stable double configuration consisting of two X/peanut structures associated to the same galactic bar - rotating with the same pattern speed - but with different morphology, formation time, and evolution. The inner, conventional peanut-shaped structure forms early via the buckling of the bar, and experiences little evolution once it stabilises. This feature is consistent in terms of size, strength and morphology, with peanut structures observed in nearby galaxies. The outer structure, however, displays a strong X, or "bow-tie", morphology. It forms just after the inner peanut, and gradually extends in time (within 1 to 1.5 Gyr) to almost the end of the bar, a radial scale where ansae occur. We conclude that, although both structures form, and are dynamically coupled to, the same bar, they are supported by inherently different mechanisms.


title: Three Dusty Star Forming Galaxies at $z\sim1.5$: Mergers and Disks on the Main Sequence

abstract: The main sequence of galaxies, a correlation between the star formation rates and stellar masses of galaxies, has been observed out to $z\sim4$. Galaxies within the scatter of the correlation are typically interpreted to be secularly evolving while galaxies with star formation rates elevated above the main sequence are interpreted to be undergoing interactions or to be Toomre-unstable disks with starbursting clumps. In this paper we investigate the recent merger histories of three dusty star forming galaxies, identified by their bright submillimeter emission at $z\sim1.5$. We analyze rest-frame optical and UV imaging, rest-frame optical emission line kinematics using slit spectra obtained with MOSFIRE on Keck I, and calculate Gini and M$_{20}$ statistics for each galaxy and conclude two are merger-driven while the third is an isolated disk galaxy. The disk galaxy lies $\sim$4$\times$ above the main sequence, one merger lies within the scatter of the main sequence, and one merger lies $\sim$4$\times$ below the main sequence. This hints that the location of a galaxy with respect to the main sequence may not be a useful discriminator of the recent star formation history of high-M$_{\star}$ galaxies at $z\sim1$.


title: From parallel to perpendicular -- On the orientation of magnetic fields in molecular clouds

abstract: We present synthetic dust polarization maps of simulated molecular clouds (MCs) with the goal to systematically explore the origin of the relative orientation of the magnetic field ($\bf{B}$) with respect to the MC sub-structures identified in density ($n$; 3D) and column density ($N$; 2D). The polarization maps are generated with the radiative transfer code POLARIS, including self-consistently calculated efficiencies for radiative torque alignment. The MCs are formed in two sets of 3D MHD simulations: in (i) colliding flows (CF), and (ii) the SILCC-Zoom simulations. In 3D, for the CF simulations with an initial field strength below $\sim$5 $\mu$G, $\bf{B}$ is oriented parallel or randomly with respect to the $n$-structures. For CF runs with stronger initial fields and all SILCC-Zoom simulations, which have an initial field strength of 3 $\mu$G, a flip from parallel to perpendicular orientation occurs at high densities of $n_\text{trans}$ $\simeq$ 10$^2$ - 10$^3$ cm$^{-3}$. We suggest that this flip happens if the MC's mass-to-flux ratio, $\mu$, is close to or below the critical value of 1. This corresponds to a field strength around 3 - 5 $\mu$G. In 2D, we use the Projected Rayleigh Statistics (PRS) to study the orientation of $\bf{B}$. If present, the flip in orientation occurs at $N_\text{trans}$ $\simeq$ 10$^{21 - 21.5}$ cm$^{-2}$, similar to the observed transition value from sub- to supercritical magnetic fields in the ISM. However, projection effects can reduce the power of the PRS method: Depending on the MC or LOS, the projected maps of the SILCC-Zoom simulations do not always show the flip, although expected from the 3D morphology. Such projection effects can explain the variety of recently observed field configurations, in particular within a single MC. Finally, we do not find a correlation between the observed orientation of $\bf{B}$ and the $N$-PDF.


title: Carrier transport theory for twisted bilayer graphene in the metallic regime

abstract: Understanding the normal-metal state transport in twisted bilayer graphene near magic angle is of fundamental importance as it provides insights into the mechanisms responsible for the observed strongly correlated insulating and superconducting phases. Here we provide a rigorous theory for phonon-dominated transport in twisted bilayer graphene describing its unusual signatures in the resistivity (including the variation with electron density, temperature, and twist angle) showing good quantitative agreement with recent experiments. We contrast this with the alternative Planckian dissipation mechanism that we show is incompatible with available experimental data. An accurate treatment of the electron-phonon scattering requires us to go well beyond the usual treatment, including both interband and intraband processes, considering the finite-temperature dynamical screening of the electron-phonon matrix element, and going beyond the linear Dirac dispersion. In addition to explaining the observations in currently available experimental data, we make concrete predictions that can be tested in ongoing experiments.


title: Venturing beyond the ISCO: Detecting X-ray emission from the plunging regions around black holes

abstract: We explore how X-ray reverberation around black holes may reveal the presence of the innermost stable circular orbit (ISCO), predicted by General Relativity, and probe the dynamics of the plunging region between the ISCO and the event horizon. Being able to directly detect the presence of the ISCO and probe the dynamics of material plunging through the event horizon represents a unique test of general relativity in the strong field regime. X-ray reverberation off of the accretion disc and material in the plunging region is modelled using general relativistic ray tracing simulations. X-ray reverberation from the plunging region has a minimal effect on the time-averaged X-ray spectrum and the overall lag-energy spectrum, but is manifested in the lag in the highest frequency Fourier components, above 0.01 c^3 (GM)^-1 (scaled for the mass of the black hole) in the 2-4keV energy band for a non-spinning black hole or the 1-2keV energy band for a maximally spinning black hole. The plunging region is distinguished from disc emission not just by the energy shifts characteristic of plunging orbits, but by the rapid increase in ionisation of material through the plunging region. Detection requires measurement of time lags to an accuracy of 20 per cent at these frequencies. Improving accuracy to 12 per cent will enable constraints to be placed on the dynamics of material in the plunging region and distinguish plunging orbits from material remaining on stable circular orbits, confirming the existence of the ISCO, a prime discovery space for future X-ray missions.


title: The white dwarf planet WD J0914+1914 b: barricading potential rocky pollutants?

abstract: An ice giant planet was recently reported orbiting white dwarf WD J0914+1914 at an approximate distance of 0.07 au. The striking non-detection of rocky pollutants in this white dwarf's photosphere contrasts with the observations of nearly every other known white dwarf planetary system. Here, I analyze the prospects for exterior extant rocky asteroids, boulders, cobbles and pebbles to radiatively drift inward past the planet due to the relatively high luminosity (0.1 L_Sun) of this particularly young (13 Myr) white dwarf. Pebbles and cobbles drift too slowly from Poynting-Robertson drag to bypass the planet, but boulders and asteroids are subject to the much stronger Yarkovsky effect. In this paper, I (i) place lower limits on the timescales for these objects to reach the planet's orbit, (ii) establish 3 m as the approximate limiting radius above which a boulder drifts too slowly to avoid colliding with the planet, and (iii) compute bounds on the fraction of boulders which succeed in traversing mean motion resonances and the planet's Hill sphere to eventually pollute the star. Overall, I find that the planet acts as a barrier against rather than a facilitator for radiatively-driven rocky pollution, suggesting that future rocky pollutants would most likely originate from distant scattering events.


title: Systematizing the Effective Theory of Self-Interacting Dark Matter

abstract: If dark matter has strong self-interactions, future astrophysical and cosmological observations, together with a clearer understanding of baryonic feedback effects, might be used to extract the velocity dependence of the dark matter scattering rate. To interpret such data, we should understand what predictions for this quantity are made by various models of the underlying particle nature of dark matter. In this paper, we systematically compute this function for fermionic dark matter with light bosonic mediators of vector, scalar, axial vector, and pseudoscalar type. We do this by matching to the nonrelativistic effective theory of self-interacting dark matter and then computing the spin-averaged viscosity cross section nonperturbatively by solving the Schrodinger equation, thus accounting for any possible Sommerfeld enhancement of the low-velocity cross section. In the pseudoscalar case, this requires a coupled-channel analysis of different angular momentum modes. We find, contrary to some earlier analyses, that nonrelativistic effects only provide a significant enhancement for the cases of light scalar and vector mediators. Scattering from light pseudoscalar and axial vector mediators is well described by tree-level quantum field theory.


title: Noncommutative generalized Gibbs ensemble in isolated integrable quantum systems

abstract: The generalized Gibbs ensemble (GGE), which involves multiple conserved quantities other than the Hamiltonian, has served as the statistical-mechanical description of the long-time behavior for several isolated integrable quantum systems. The GGE may involve a noncommutative set of conserved quantities in view of the maximum entropy principle, and show that the GGE thus generalized (noncommutative GGE, NCGGE) gives a more qualitatively accurate description of the long-time behaviors than that of the conventional GGE. Providing a clear understanding of why the (NC)GGE well describes the long-time behaviors, we construct, for noninteracting models, the exact NCGGE that describes the long-time behaviors without an error even at finite system size. It is noteworthy that the NCGGE involves nonlocal conserved quantities, which can be necessary for describing long-time behaviors of local observables. We also give some extensions of the NCGGE and demonstrate how accurately they describe the long-time behaviors of few-body observables.


title: Low-Noise GaAs Quantum Dots for Quantum Photonics

abstract: Quantum dots are both excellent single-photon sources and hosts for single spins. This combination enables the deterministic generation of Raman-photons -- bandwidth-matched to an atomic quantum-memory -- and the generation of photon cluster states, a resource in quantum communication and measurement-based quantum computing. GaAs quantum dots in AlGaAs can be matched in frequency to a rubidium-based photon memory, and have potentially improved electron spin coherence compared to the widely used InGaAs quantum dots. However, their charge stability and optical linewidths are typically much worse than for their InGaAs counterparts. Here, we embed GaAs quantum dots into an $n$-$i$-$p$-diode specially designed for low-temperature operation. We demonstrate ultra-low noise behaviour: charge control via Coulomb blockade, close-to lifetime-limited linewidths, and no blinking. We observe high-fidelity optical electron-spin initialisation and long electron-spin lifetimes for these quantum dots. Our work establishes a materials platform for low-noise quantum photonics close to the red part of the spectrum.


title: New material platform for superconducting transmon qubits with coherence times exceeding 0.3 milliseconds

abstract: The superconducting transmon qubit is a leading platform for quantum computing and quantum science. Building large, useful quantum systems based on transmon qubits will require significant improvements in qubit relaxation and coherence times, which are orders of magnitude shorter than limits imposed by bulk properties of the constituent materials. This indicates that relaxation likely originates from uncontrolled surfaces, interfaces, and contaminants. Previous efforts to improve qubit lifetimes have focused primarily on designs that minimize contributions from surfaces. However, significant improvements in the lifetime of two-dimensional transmon qubits have remained elusive for several years. Here, we fabricate two-dimensional transmon qubits that have both lifetimes and coherence times with dynamical decoupling exceeding 0.3 milliseconds by replacing niobium with tantalum in the device. We have observed increased lifetimes for seventeen devices, indicating that these material improvements are robust, paving the way for higher gate fidelities in multi-qubit processors.


title: Effective momentum-momentum coupling in a correlated electronic system

abstract: We present a way of partly reincorporate the effects of the localized bonding electrons on the dynamics of their itinerant counterparts in Hubbard-like Hamiltonians. This is done by relaxing the constraint that the former should be entirely frozen in the chemical bonds between the underlying lattice sites through the employment of a Born-Oppenheimer-like ansatz for the wavefunction of the whole electronic system. Accordingly, the latter includes itinerant as well as bonding electron coordinates. Going beyond the adiabatic approximation, we show that the net effect of virtual transitions of bonding electrons between their ground and excited states is to furnish the itinerant electrons with an effective inter-electronic momentum-momentum interaction. Although we have applied these ideas to the specific case of rings, our assumptions can be generalized to higher dimensional systems sharing the required properties of which we have made use herein.


title: Convective dynamics with mixed temperature boundary conditions: why thermal relaxation matters and how to accelerate it

abstract: Astrophysical simulations of convection frequently impose different thermal boundary conditions at the top and the bottom of the domain in an effort to more accurately model natural systems. In this work, we study Rayleigh-Benard convection (RBC) under the Boussinesq approximation. We examine simulations with mixed temperature boundary conditions in which the flux is fixed at the bottom boundary and the temperature is fixed at the top ("FT"). We aim to understand how FT boundaries change the nature of the convective solution compared to the traditional choice of thermal boundaries, in which the temperature is fixed at the top and bottom of the domain ("TT"). We demonstrate that the timescale of thermal relaxation for FT simulations is dependent upon the initial conditions. "Classic" initial conditions which employ a hydrostatically -- and thermally -- balanced linear temperature profile exhibit a long thermal relaxation. This long relaxation is not seen in FT simulations which use a TT simulation's nonlinear state as initial conditions ("TT-to-FT"). In the thermally relaxed, statistically stationary state, the mean behavior of an FT simulation corresponds to an equivalent simulation with TT boundaries, and time- and volume-averaged flow statistics like the Nusselt number and the Peclet number are indistinguishable between FT and TT simulations. FT boundaries are fundamentally asymmetric, and we examine the asymmetries that these boundaries produce in the flow. We find that the fixed-flux boundary produces more extreme temperature events than the fixed-temperature boundary. However, these near-boundary asymmetries do not measurably break the symmetry in the convective interior. We briefly explore rotating RBC to demonstrate that our findings with respect to thermal relaxation carry over to this more complex case, and to show the power of TT-to-FT initial conditions.


title: Flat Topological Bands and Eigenstate Criticality in a Quasiperiodic Insulator

abstract: The effects of downfolding a Brillouin zone can open gaps and quench the kinetic energy by flattening bands. Quasiperiodic systems are extreme examples of this process, which leads to new phases and critical eigenstates. We analytically and numerically investigate these effects in a two dimensional topological insulator with a quasiperiodic potential and discover a complex phase diagram. We study the nature of the resulting eigenstate quantum phase transitions; a quasiperiodic potential can make a trivial insulator topological and induce topological insulator-to-metal phase transitions through a unique universality class distinct from random systems. This wealth of critical behavior occurs concomitantly with the quenching of the kinetic energy, resulting in flat topological bands that could serve as a platform to realize the fractional quantum Hall effect without a magnetic field.


title: An ultra-massive white dwarf with a mixed hydrogen-carbon atmosphere as a likely merger remnant

abstract: White dwarfs are dense, cooling stellar embers consisting mostly of carbon and oxygen, or oxygen and neon (with a few percent carbon) at higher initial stellar masses. These stellar cores are enveloped by a shell of helium which in turn is usually surrounded by a layer of hydrogen, generally prohibiting direct observation of the interior composition. However, carbon is observed at the surface of a sizeable fraction of white dwarfs, sometimes with traces of oxygen, and it is thought to be dredged-up from the core by a deep helium convection zone. In these objects only traces of hydrogen are found as large masses of hydrogen are predicted to inhibit hydrogen/helium convective mixing within the envelope. We report the identification of WDJ055134.612+413531.09, an ultra-massive (1.14 $M_\odot$) white dwarf with a unique hydrogen/carbon mixed atmosphere (C/H=0.15 in number ratio). Our analysis of the envelope and interior indicates that the total hydrogen and helium mass fractions must be several orders of magnitude lower than predictions of single star evolution: less than $10^{-9.5}$ and $10^{-7.0}$, respectively. Due to the fast kinematics ($129\pm5$ km/s relative to the local standard of rest), large mass, and peculiar envelope composition, we argue that WDJ0551+4135 is consistent with formation from the merger of two white dwarfs in a tight binary system.


title: Estimating the impact of treatment compliance over time on smoking cessation using data from ecological momentary assessments (EMA)

abstract: The Wisconsin Smoker's Health Study (WSHS2) was a longitudinal trial conducted to compare the effectiveness of two commonly used smoking cessation treatments, varenicline and combination nicotine replacement therapy (cNRT) with the less intense standard of care, nicotine patch. The main outcome of the WSHS2 study was that all three treatments had equivalent treatment effects. However, in-depth analysis of the compliance data collected via ecological momentary assessment (EMA) were not analyzed. Compliance to the treatment regimens may represent a confounder as varenicline and cNRT are more intense treatments and would likely have larger treatment effects if all subjects complied. In order to estimate the causal compliance effect, we view the counterfactual, the outcome that would have been observed if the subject was allocated to the treatment counter to the fact, as a missing data problem and proceed to impute the counterfactual. Our contribution to the methodological literature lies in the extension of this idea to a more general analytic approach that includes mediators and confounders of the mediator-outcome relationship. Simulation results suggest that our method works well and application to the WSHS2 data suggest that the treatment effect of nicotine patch, varenicline, and cNRT are equivalent after accounting for differences in treatment compliance.


title: Policy-Aware Model Learning for Policy Gradient Methods

abstract: This paper considers the problem of learning a model in model-based reinforcement learning (MBRL). We examine how the planning module of an MBRL algorithm uses the model, and propose that the model learning module should incorporate the way the planner is going to use the model. This is in contrast to conventional model learning approaches, such as those based on maximum likelihood estimate, that learn a predictive model of the environment without explicitly considering the interaction of the model and the planner. We focus on policy gradient type of planning algorithms and derive new loss functions for model learning that incorporate how the planner uses the model. We call this approach Policy-Aware Model Learning (PAML). We theoretically analyze a generic model-based policy gradient algorithm and provide a convergence guarantee for the optimized policy. We also empirically evaluate PAML on some benchmark problems, showing promising results.


title: Experimental Demonstration of Dynamic Thermal Regulation using Vanadium Dioxide Thin Films

abstract: We present an experimental demonstration of passive, dynamic thermal regulation in a solid-state system with temperature-dependent thermal emissivity switching. We achieve this effect using a multilayered device, comprised of a vanadium dioxide (VO2) thin film on a silicon substrate with a gold back reflector. We experimentally characterize the optical properties of the VO2 film and use the results to optimize device design. Using a calibrated, transient calorimetry experiment we directly measure the temperature fluctuations arising from a time-varying heat load. Under laboratory conditions, we find that the device regulates temperature better than a constant emissivity sample. We use the experimental results to validate our thermal model, which can be used to predict device performance under the conditions of outer space. In this limit, thermal fluctuations are halved with reference to a constant-emissivity sample.


title: Declarative Stream Runtime Verification (hLola)

abstract: Stream Runtime Verification is a formal dynamic analysis technique that generalizes runtime verification algorithms from temporal logics like LTL to stream monitoring, allowing to compute richer verdicts than Booleans (including quantitative and arbitrary data). In this paper we study the problem of implementing an SRV engine that is truly extensible to arbitrary data theories, and we propose a solution as a Haskell embedded domain specific language. In spite of the theoretical clean separation in SRV between temporal dependencies and data computations, previous engines include ad-hoc implementations of a few data types, requiring complex changes to incorporate new data theories. We propose here an SRV language called hLola that borrows general Haskell types and embeds them transparently into an eDSL. This novel technique, which we call lift deep embedding, allows for example, the use of higher-order functions for static stream parameterization. We describe the Haskell implementation of hLola and illustrate simple extensions implemented using libraries, which require long and error-prone additions in other ad-hoc SRV formalisms.


title: Dynamic Beveridge Curve Accounting

abstract: We develop a dynamic decomposition of the empirical Beveridge curve, i.e., the level of vacancies conditional on unemployment. Using a standard model, we show that three factors can shift the Beveridge curve: reduced-form matching efficiency, changes in the job separation rate, and out-of-steady-state dynamics. We find that the shift in the Beveridge curve during and after the Great Recession was due to all three factors, and each factor taken separately had a large effect. Comparing the pre-2010 period to the post-2010 period, a fall in matching efficiency and out-of-steady-state dynamics both pushed the curve upward, while the changes in the separation rate pushed the curve downward. The net effect was the observed upward shift in vacancies given unemployment. In previous recessions changes in matching efficiency were relatively unimportant, while dynamics and the separation rate had more impact. Thus, the unusual feature of the Great Recession was the deterioration in matching efficiency, while separations and dynamics have played significant, partially offsetting roles in most downturns. The importance of these latter two margins contrasts with much of the literature, which abstracts from one or both of them. We show that these factors affect the slope of the empirical Beveridge curve, an important quantity in recent welfare analyses estimating the natural rate of unemployment.


title: Vibrational dressing in Kinetically Constrained Rydberg Spin Systems

abstract: Quantum spin systems with kinetic constraints have become paradigmatic for exploring collective dynamical behaviour in many-body systems. Here we discuss a facilitated spin system which is inspired by recent progress in the realization of Rydberg quantum simulators. This platform allows to control and investigate the interplay between facilitation dynamics and the coupling of spin degrees of freedom to lattice vibrations. Developing a minimal model, we show that this leads to the formation of polaronic quasiparticle excitations which are formed by many-body spin states dressed by phonons. We investigate in detail the properties of these quasiparticles, such as their dispersion relation, effective mass and the quasiparticle weight. Rydberg lattice quantum simulators are particularly suited for studying this phonon-dressed kinetically constrained dynamics as their exaggerated length scales permit the site-resolved monitoring of spin and phonon degrees of freedom.

categories: cond-mat.quant-gas cond-mat.stat-mech physics.atom-ph quant-ph


title: Avoiding monochromatic sub-paths in uniform hypergraph paths and cycles

abstract: We present a recursive formula for the number of ways to color $j$ vertices blue in an r-uniform hyperpath of size $n$ while avoiding a blue monochromatic sub-hyperpath of length k. We use this result to solve the corresponding problem for $(r-1)$-tight r-uniform paths and loose r-uniform cycles. This generalizes some well known results from reliability engineering and analysis.


title: Theoretical study of the pressure-induced structure, phase transition, mechanical and electronic properties in the V-N system

abstract: Stable compounds in the V-N system are systematically searched and four new high-pressure phases are found, including C2/m-V$_9$N, Pbam-V$_5$N$_2$, Pnma-V$_2$N and I4/mcm-VN$_2$. V$_2$N undergoes a phase transition from $\varepsilon$-Fe$_2$N-type V$_2$N (P$\bar{3}$1m) to $\zeta$-Fe$_2$N-type V$_2$N (Pbcn) at 10 GPa and to Fe$_2$C-type V$_2$N (Pnnm) at 59 GPa, then to Pnma-V$_2$N at 96 GPa. Low-temperature tetragonal VN is theoretically proved to belong to space group P$\bar{4}$2m. The estimated Vickers hardnesses and fracture toughness of WC-type VN are around 37 GPa and 4.3-6.1 MPa m$^{1/2}$, respectively. Al$_2$Cu-type VN$_2$ (I4/mcm) with a Vickers hardness of 25-27 GPa and fracture toughness of 3.6-6.6 MPa m$^{1/2}$ also shows excellent mechanical properties. Elastic properties of WC-type mononitrides of transition metals from IVB group (Ti, Zr and Hf), VB group (V, Nb and Ta) and VIB (Cr, Mo and W) are calculated and compared. Both the bond strength and structural configuration determine the mechanical properties of a material.


title: Tits type alternative for groups acting on toric affine varieties

abstract: Given a toric affine algebraic variety $X$ and a collection of one-parameter unipotent subgroups $U_1,\ldots,U_s$ of $\mathop{\rm Aut}(X)$ which are normalized by the torus acting on $X$, we show that the group $G$ generated by $U_1,\ldots,U_s$ verifies the following alternative of Tits' type: either $G$ is a unipotent algebraic group, or it contains a non-abelian free subgroup. We deduce that if $G$ is $2$-transitive on a $G$-orbit in $X$, then $G$ contains a non-abelian free subgroup, and so, is of exponential growth.


title: Experimental determination of the orientation of tilted magnetic anisotropy by measuring Hall voltage

abstract: The orientation of the tilted magnetic anisotropy has crucial importance in many spintronic devices. However, it is very challenging to determine it especially in very small structures produced by lithography. Here, we propose a new experimental method to directly and accurately measure both the polar and azimuthal angles of a tilted magnetic anisotropy. By using the proposed experimental method, we have successfully determined the out-of-plane and in-plane angles of the tilted magnetic anisotropy in a micro-structured multilayer thin film. The orientation of the tilted magnetic anisotropy in the sample has also been confirmed by the theoretical simulations proving the accuracy of the method.


title: Unsupervised Machine Learning of Quenched Gauge Symmetries: A Proof-of-Concept Demonstration

abstract: In condensed matter physics, one of the goals of machine learning is the classification of phases of matter. The consideration of a system's symmetries can significantly assist the machine in this goal. We demonstrate the ability of an unsupervised machine learning protocol, the Principal Component Analysis method, to detect hidden quenched gauge symmetries introduced via the so-called Mattis gauge transformation. Our work reveals that unsupervised machine learning can identify hidden properties of a model and may therefore provide new insights into the models themselves.


title: Regularisation Can Mitigate Poisoning Attacks: A Novel Analysis Based on Multiobjective Bilevel Optimisation

abstract: Machine Learning (ML) algorithms are vulnerable to poisoning attacks, where a fraction of the training data is manipulated to deliberately degrade the algorithms' performance. Optimal poisoning attacks, which can be formulated as bilevel optimisation problems, help to assess the robustness of learning algorithms in worst-case scenarios. However, current attacks against algorithms with hyperparameters typically assume that these hyperparameters remain constant ignoring the effect the attack has on them. We show that this approach leads to an overly pessimistic view of the robustness of the algorithms. We propose a novel optimal attack formulation that considers the effect of the attack on the hyperparameters by modelling the attack as a multiobjective bilevel optimisation problem. We apply this novel attack formulation to ML classifiers using $L_2$ regularisation and show that, in contrast to results previously reported, $L_2$ regularisation enhances the stability of the learning algorithms and helps to mitigate the attacks. Our empirical evaluation on different datasets confirms the limitations of previous strategies, evidences the benefits of using $L_2$ regularisation to dampen the effect of poisoning attacks and shows how the regularisation hyperparameter increases with the fraction of poisoning points.


title: InfiniWolf: Energy Efficient Smart Bracelet for Edge Computing with Dual Source Energy Harvesting

abstract: This work presents InfiniWolf, a novel multi-sensor smartwatch that can achieve self-sustainability exploiting thermal and solar energy harvesting, performing computationally high demanding tasks. The smartwatch embeds both a System-on-Chip (SoC) with an ARM Cortex-M processor and Bluetooth Low Energy (BLE) and Mr. Wolf, an open-hardware RISC-V based parallel ultra-low-power processor that boosts the processing capabilities on board by more than one order of magnitude, while also increasing energy efficiency. We demonstrate its functionality based on a sample application scenario performing stress detection with multi-layer artificial neural networks on a wearable multi-sensor bracelet. Experimental results show the benefits in terms of energy efficiency and latency of Mr. Wolf over an ARM Cortex-M4F micro-controllers and the possibility, under specific assumptions, to be self-sustainable using thermal and solar energy harvesting while performing up to 24 stress classifications per minute in indoor conditions.


title: Purcell enhancement of a single silicon carbide color center with coherent spin control

abstract: Silicon carbide has recently been developed as a platform for optically addressable spin defects. In particular, the neutral divacancy in the 4H polytype displays an optically addressable spin-1 ground state and near-infrared optical emission. Here, we present the Purcell enhancement of a single neutral divacancy coupled to a photonic crystal cavity. We utilize a combination of nanolithographic techniques and a dopant-selective photoelectrochemical etch to produce suspended cavities with quality factors exceeding 5,000. Subsequent coupling to a single divacancy leads to a Purcell factor of ~50, which manifests as increased photoluminescence into the zero-phonon line and a shortened excited-state lifetime. Additionally, we measure coherent control of the divacancy ground state spin inside the cavity nanostructure and demonstrate extended coherence through dynamical decoupling. This spin-cavity system represents an advance towards scalable long-distance entanglement protocols using silicon carbide that require the interference of indistinguishable photons from spatially separated single qubits.


title: Finding archetypal patterns for binary questionnaires

abstract: Archetypal analysis is an exploratory tool that explains a set of observations as mixtures of pure (extreme) patterns. If the patterns are actual observations of the sample, we refer to them as archetypoids. For the first time, we propose to use archetypoid analysis for binary observations. This tool can contribute to the understanding of a binary data set, as in the multivariate case. We illustrate the advantages of the proposed methodology in a simulation study and two applications, one exploring objects (rows) and the other exploring items (columns). One is related to determining student skill set profiles and the other to describing item response functions.


title: Asymptotically safe gravity

abstract: This is an introduction to asymptotically safe quantum gravity, explaining the main idea of asymptotic safety and how it could solve the problem of predictivity in quantum gravity. In the first part, the concept of an asymptotically safe fixed point is discussed within the functional Renormalization Group framework for gravity, which is also briefly reviewed. A concise overview of key results on asymptotically safe gravity is followed by a short discussion of important open questions. The second part highlights how the interplay with matter provides observational consistency tests for all quantum-gravity models, followed by an overview of the state of results on asymptotic safety and its implications in gravity-matter models. Finally, effective asymptotic safety is briefly discussed as a scenario in which asymptotically safe gravity could be connected to other approaches to quantum gravity.


title: Library Adoption Dynamics in Software Teams

abstract: When a group of people strives to understand new information, struggle ensues as various ideas compete for attention. Steep learning curves are surmounted as teams learn together. To understand how these team dynamics play out in software development, we explore Git logs, which provide a complete change history of software repositories. In these repositories, we observe code additions, which represent successfully implemented ideas, and code deletions, which represent ideas that have failed or been superseded. By examining the patterns between these commit types, we can begin to understand how teams adopt new information. We specifically study what happens after a software library is adopted by a project, i.e. when a library is used for the first time in the project. We find that a variety of factors, including team size, library popularity, and prevalence on Stack Overflow are associated with how quickly teams learn and successfully adopt new software libraries.


title: The monoatomic FPU system as a limit of a diatomic FPU system

abstract: We consider a diatomic infinite Fermi-Pasta-Ulam (FPU) system with light and heavy particles. For a small mass ratio, we prove error estimates for the approximation of the dynamics of this system by the dynamics of the monoatomic FPU system. The light particles are squeezed by the heavy particles at the mean value of their displacements. The error estimates are derived by means of the energy method and hold for sufficiently long times, for which the dynamics of the monoatomic FPU system is observed. The approximation result is restricted to sufficiently small displacements of the heavy particles relatively to each other.

categories: nlin.SI math-ph math.AP math.DS math.MP nlin.PS


title: Overcoming the tradeoff between confinement and focal distance using virtual ultrasonic optical waveguides

abstract: Conventional optical lenses have been used to focus light from outside without disturbing the medium. The focused spot size is proportional to the focal distance in a conventional lens, resulting in a tradeoff between depth of penetration in the target medium and spatial resolution. We have shown that virtual ultrasonically sculpted gradient-index (GRIN) optical waveguides can be formed in the target medium to guide and steer light without disturbing the medium. Here, we demonstrate that such virtual waveguides can relay an externally focused beam of light through the medium beyond the focal distance of an external physical lens to extend the penetration depth without compromising the spot size. Moreover, the spot size can be tuned by reconfiguring the virtual waveguide. We show that these virtual GRIN waveguides can be formed in transparent as well as turbid media to enhance the confinement and contrast ratio of the focused beam of light at the target location. This method can be extended to realize complex optical systems of external physical lenses and in situ virtual waveguides to extend the reach and flexibility of optical methods.


title: On the formation and recurrent shedding of ligaments in droplet aerobreakup

abstract: The breakup of water droplets when exposed to high-speed gas flows is investigated using both high-magnification shadowgraphy experiments as well as fully three-dimensional numerical simulations, which account for viscous as well as capillary effects. After thorough validation of the simulations with respect to the experiments, we elucidate the ligament formation process and the effect of surface tension. By Fourier decomposition of the flow field, we observe the development of specific azimuthal modes, which destabilize the liquid sheet surrounding the droplet. Eventually, the liquid sheet is ruptured, which leads to the formation of ligaments. We further observe the ligament formation and shedding to be a recurrent process. While the first ligament shedding weakly depends on the Weber number, subsequent shedding processes seem to be driven primarily by inertia and the vortex shedding in the wake of the deformed droplet.


title: Analysis of Skew Quadrupole Compensation in RF-Photoinjectors

abstract: In this paper, we present a detailed analysis of the effects of quadrupole components in radiofrequency (RF) photoinjectors, which could occur either in the RF cavity or in the gun solenoid, and their compensation using a skew quadrupole corrector. The analysis is generalized to the case in which the correction cannot be applied at the same location as the quadrupole. A practical guide to setting the corrector to eliminate correlation in the x-y space is discussed. Particle tracking simulations and measurements of emittance and beam matrix elements obtained at the UCLA Pegasus photoinjector are presented to validate the analysis.


title: Dynamic Nuclear Polarization by optical Stark effect in periodically-pumped gallium arsenide

abstract: Optical pump-probe time- and magnetic-field-resolved Kerr rotation measurements provide a window into the mechanisms that generate dynamic nuclear polarization in bulk gallium arsenide. Previously, we have reported an unexpected dependence of the direction of the nuclear polarization on the sweep direction of the applied external magnetic field. In this paper, we present numerical calculations based on a model for this nuclear polarization due to the optical orientation and optical Stark effect produced by a train of ultrafast optical pulses. We demonstrate the correspondence of the model to our experimental measurements for different laser wavelengths and magnetic field sweep rates. Finally, we show that the model reproduces the sweep direction dependence and provides an explanation for this behavior.


title: Dynamic Skyline Queries on Encrypted Data Using Result Materialization

abstract: Skyline computation is an increasingly popular query, with broad applicability in domains such as healthcare, travel and finance. Given the recent trend to outsource databases and query evaluation, and due to the proprietary and sometimes highly sensitivity nature of the data (e.g., in healthcare), it is essential to evaluate skylines on encrypted datasets. Several research efforts acknowledged the importance of secure skyline computation, but existing solutions suffer from at least one of the following shortcomings: (i) they only provide ad-hoc security; (ii) they are prohibitively expensive; or (iii) they rely on unrealistic assumptions, such as the presence of multiple non-colluding parties in the protocol. Inspired from solutions for secure nearest-neighbors (NN) computation, we conjecture that the most secure and efficient way to compute skylines is through result materialization. However, this approach is significantly more challenging for skylines than for NN queries. We exhaustively study and provide algorithms for pre-computation of skyline results, and we perform an in-depth theoretical analysis of this process. We show that pre-computing results while minimizing storage overhead is NP-hard, and we provide dynamic programming and greedy heuristics that solve the problem more efficiently, while maintaining storage at reasonable levels. Our algorithms are novel and applicable to plain-text skyline computation, but we focus on the encrypted setting where materialization reduces the cost of skyline computation from hours to seconds. Extensive experiments show that we clearly outperform existing work in terms of performance, and our security analysis proves that we obtain a smaller (and quantifiable) data leakage than competitors.


title: Learning Nonparametric Human Mesh Reconstruction from a Single Image without Ground Truth Meshes

abstract: Nonparametric approaches have shown promising results on reconstructing 3D human mesh from a single monocular image. Unlike previous approaches that use a parametric human model like skinned multi-person linear model (SMPL), and attempt to regress the model parameters, nonparametric approaches relax the heavy reliance on the parametric space. However, existing nonparametric methods require ground truth meshes as their regression target for each vertex, and obtaining ground truth mesh labels is very expensive. In this paper, we propose a novel approach to learn human mesh reconstruction without any ground truth meshes. This is made possible by introducing two new terms into the loss function of a graph convolutional neural network (Graph CNN). The first term is the Laplacian prior that acts as a regularizer on the reconstructed mesh. The second term is the part segmentation loss that forces the projected region of the reconstructed mesh to match the part segmentation. Experimental results on multiple public datasets show that without using 3D ground truth meshes, the proposed approach outperforms the previous state-of-the-art approaches that require ground truth meshes for training.


title: Dispersing Fermi-Ulam Models

abstract: We study a natural class of Fermi-Ulam Models that features good hyperbolicity properties and that we call dispersing Fermi-Ulam models. Using tools inspired by the theory of hyperbolic billiards we prove, under very mild complexity assumptions, a Growth Lemma for our systems. This allows us to obtain ergodicity of dispersing Fermi-Ulam Models. It follows that almost every orbit of such systems is oscillatory.