title: A new mechanism for symmetry breaking from nilmanifolds

abstract: We present a method to obtain a scalar potential at tree level from a pure gauge theory on nilmanifolds, a class of negatively-curved compact spaces, and discuss the spontaneous symmetry breaking mechanism induced in the residual Minkowski space after compactification at low energy. We show that the scalar potential is completely determined by the gauge symmetries and the geometry of the compact manifold. In order to allow for simple analytic calculations we consider three extra space dimensions as the minimal example of a nilmanifold, therefore considering a pure Yang-Mills theory in seven dimensions.


title: Local Group star formation in warm and self-interacting dark matter cosmologies

abstract: The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting dark matter (SIDM) models with the EAGLE galaxy formation code to study how galaxy formation times change with dark matter model. We are able to identify the same haloes in different simulations, since they share the same initial density field phases. We find that the stellar mass of galaxies depends systematically on resolution, and can differ by as much as a factor of two in haloes of a given dark matter mass. The evolution of the stellar populations in SIDM is largely identical to that of CDM, but in WDM early star formation is instead suppressed. The time at which LG haloes can begin to form stars through atomic cooling is delayed by $\sim$200~Myr in WDM models compared to CDM. It will be necessary to measure stellar ages of old populations to a precision of better than 100~Myr, and to address degeneracies with the redshift of reionization -- and potentially other baryonic processes -- in order to use these observables to distinguish between dark matter models.


title: Magneto-gravity wave packet dynamics in strongly magnetised cores of evolved stars

abstract: Magnetic fields are believed to be generated in the cores of massive main sequence stars, and these may survive on to later stages of evolution. Observations of depressed dipole modes in red giant stars have been touted as evidence for these fields, but the predictions of existing magnetic theories have difficulty accommodating several aspects, including the need to return a fraction of wave energy from the core to the envelope, and the persistent gravity-like character of affected modes. In this work we perform a Hamiltonian ray tracing study investigating the dynamics of magneto-gravity waves in full spherical geometry, using realistic stellar models and magnetic field configurations. This technique applies in the limit where wavelengths are much shorter than scales of background variation. We conduct a comprehensive exploration of parameter space, examining the roles of wave frequency, spherical harmonic degree, wavevector polarisation, incoming latitude, field strength, field radius, and evolutionary state. We demonstrate that even in the presence of a strong field, there exist trajectories where waves remain predominantly gravity-like in character, and these are able to undergo reflection out of the core much like pure gravity waves. The remaining trajectories are ones where waves acquire significant Alfven character, becoming trapped and eventually dissipated. Orientation effects, i.e. wavevector polarisation and incoming latitude, are found to be crucial factors in determining the outcome (trapped versus reflected) of individual wave packets. The allowance for partial energy return from the core offers a solution to the conundrum faced by the magnetic hypothesis.


title: Obliquity Constraints on an Extrasolar Planetary-Mass Companion

abstract: We place the first constraints on the obliquity of a planetary-mass companion (PMC) outside of the Solar System. Our target is the directly imaged system 2MASS J01225093-2439505 (2M0122), which consists of a 120 Myr 0.4 M_sun star hosting a 12-27 M_J companion at 50 AU. We constrain all three of the system's angular momentum vectors: how the companion spin axis, the stellar spin axis, and the orbit normal are inclined relative to our line of sight. To accomplish this, we measure projected rotation rates (vsini) for both the star and the companion using new near-infrared high-resolution spectra with NIRSPEC at Keck Observatory. We combine these with a new stellar photometric rotation period from TESS and a published companion rotation period from HST to obtain spin axis inclinations for both objects. We also fitted multiple epochs of astrometry, including a new observation with NIRC2/Keck, to measure 2M0122b's orbital inclination. The three line-of-sight inclinations place limits on the true de-projected companion obliquity and stellar obliquity. We find that while the stellar obliquity marginally prefers alignment, the companion obliquity tentatively favors misalignment. We evaluate possible origin scenarios. While collisions, secular spin-orbit resonances, and Kozai-Lidov oscillations are unlikely, formation by gravitational instability in a gravito-turbulent disk - the scenario favored for brown dwarf companions to stars - appears promising.


title: Electron Heating in Perpendicular Low-Beta Shocks

abstract: Collisionless shocks heat electrons in the solar wind, interstellar blast waves, and hot gas permeating galaxy clusters. How much shock heating goes to electrons instead of ions, and what plasma physics controls electron heating? We simulate 2-D perpendicular shocks with a fully kinetic particle-in-cell code. For magnetosonic Mach number $\mathcal{M}_\mathrm{ms} \sim 1$-$10$ and plasma beta $\beta_\mathrm{p} \lesssim 4$, the post-shock electron/ion temperature ratio $T_\mathrm{e}/T_\mathrm{i}$ decreases from $1$ to $0.1$ with increasing $\mathcal{M}_\mathrm{ms}$. In a representative $\mathcal{M}_\mathrm{ms}=3.1$, $\beta_\mathrm{p}=0.25$ shock, electrons heat above adiabatic compression in two steps: ion-scale $E_\parallel = \vec{E} \cdot \hat{b}$ accelerates electrons into streams along $\vec{B}$, which then relax via two-stream-like instability. The $\vec{B}$-parallel heating is mostly induced by waves; $\vec{B}$-perpendicular heating is mostly adiabatic compression by quasi-static fields.


title: Finite-Temperature Quarkyonic Matter with an Excluded Volume Model

abstract: We introduce a theoretical framework to analyze the thermal properties of a recently proposed "excluded volume" model of quarkyonic matter. This entails proposing finite temperature distribution functions and entropy functionals constrained by internal theoretical consistency. We use the formalism developed here to analyze the effect of temperature on the quark onset as a function of baryon density. We find that the baryon density at which the quarkyonic phase emerges decreases only slightly with temperature, an effect produced by large $N_c$. We furthermore find that a finite temperature increases the density of quarks produced. The framework presented here can be used to compute the equation of state of hot and dense quarkyonic matter which is of relevance to neutron star mergers.


title: White dwarf deflagrations for Type Iax supernovae: Polarisation signatures from the explosion and companion interaction

abstract: Growing evidence suggests that Type Iax supernovae might be the result of thermonuclear deflagrations of Chandrasekhar-mass white dwarfs in binary systems. We carry out Monte Carlo radiative transfer simulations and predict spectropolarimetric features originating from the supernova explosion and subsequent ejecta interaction with the companion star. Specifically, we calculate viewing-angle dependent flux and polarisation spectra for a 3D model simulating the deflagration of a Chandrasekhar-mass white dwarf and, for a second model, simulating the ejecta interaction with a main-sequence star. We find that the intrinsic signal is weakly polarised and only mildly viewing-angle dependent, owing to the overall spherical symmetry of the explosion and the depolarising contribution of iron-group elements dominating the ejecta composition. The interaction with the companion star carves out a cavity in the ejecta and produces a detectable, but modest signal that is significant only at relatively blue wavelengths ($\lesssim$ 5000 $\unicode{x212B}$). In particular, increasingly fainter and redder spectra are predicted for observer orientations further from the cavity, while a modest polarisation signal $P\sim0.2$ per cent is found at blue wavelengths for orientations 30$^\circ$ and 45$^\circ$ away from the cavity. We find a reasonable agreement between the interaction model viewed from these orientations and spectropolarimetric data of SN 2005hk and interpret the maximum-light polarisation signal seen at blue wavelengths for this event as a possible signature of the ejecta-companion interaction. We encourage further polarimetric observations of SNe Iax to test whether our results can be extended and generalised to the whole SN Iax class.


title: Fingerprints of giant planets in the composition of solar twins

abstract: The Sun shows a $\sim 10$% depletion in refractory elements relative to nearby solar twins. It has been suggested that this depletion is a signpost of planet formation. The exoplanet statistics are now good enough to show that the origin of this depletion does not arise from the sequestration of refractory material inside the planets themselves. This conclusion arises because most sun-like stars host close-in planetary systems that are on average more massive than the Sun's. Using evolutionary models for the protoplanetary discs that surrounded the young Sun and solar twins we demonstrate that the origin of the depletion likely arises due to the trapping of dust exterior to the orbit of a forming giant planet. In this scenario a forming giant planet opens a gap in the gas disc, creating a pressure trap. If the planet forms early enough, while the disc is still massive, the planet can trap $\gtrsim 100 $M$_\oplus$ of dust exterior to its orbit, preventing the dust from accreting onto the star in contrast to the gas. Forming giant planets can create refractory depletions of $\sim 5-15$%, with the larger values occurring for initial conditions that favour giant planet formation (e.g. more massive discs, that live longer). The incidence of solar-twins that show refractory depletion matches both the occurrence of giant planets discovered in exoplanet surveys and "transition" discs that show similar depletion patterns in the material that is accreting onto the star.


title: Fluid & Gyrokinetic turbulence in open field-line, helical plasmas

abstract: Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently has it become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas (T. N. Bernard, et. al., Phys. of Plasmas 26, 042301 (2019)). Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs which can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, as well as improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.


title: Dynamic Incentive-aware Learning: Robust Pricing in Contextual Auctions

abstract: Motivated by pricing in ad exchange markets, we consider the problem of robust learning of reserve prices against strategic buyers in repeated contextual second-price auctions. Buyers' valuations for an item depend on the context that describes the item. However, the seller is not aware of the relationship between the context and buyers' valuations, i.e., buyers' preferences. The seller's goal is to design a learning policy to set reserve prices via observing the past sales data, and her objective is to minimize her regret for revenue, where the regret is computed against a clairvoyant policy that knows buyers' heterogeneous preferences. Given the seller's goal, utility-maximizing buyers have the incentive to bid untruthfully in order to manipulate the seller's learning policy. We propose learning policies that are robust to such strategic behavior. These policies use the outcomes of the auctions, rather than the submitted bids, to estimate the preferences while controlling the long-term effect of the outcome of each auction on the future reserve prices. When the market noise distribution is known to the seller, we propose a policy called Contextual Robust Pricing (CORP) that achieves a T-period regret of $O(d\log(Td) \log (T))$, where $d$ is the dimension of {the} contextual information. When the market noise distribution is unknown to the seller, we propose two policies whose regrets are sublinear in $T$.


title: Revealing the tidal scars of the Small Magellanic Cloud

abstract: Due to their close proximity, the Large and Small Magellanic Clouds (SMC/LMC) provide natural laboratories for understanding how galaxies form and evolve. With the goal of determining the structure and dynamical state of the SMC, we present new spectroscopic data for $\sim$ 3000 SMC red giant branch stars observed using the AAOmega spectrograph at the Anglo-Australian Telescope. We complement our data with further spectroscopic measurements from previous studies that used the same instrumental configuration and proper motions from the \textit{Gaia} Data Release 2 catalogue. Analysing the photometric and stellar kinematic data, we find that the SMC centre of mass presents a conspicuous offset from the velocity centre of its associated $\mbox{H\,{\sc i}}$ gas, suggesting that the SMC gas is likely to be far from dynamical equilibrium. Furthermore, we find evidence that the SMC is currently undergoing tidal disruption by the LMC within 2\,kpc of the centre of the SMC, and possibly all the way in to the very core. This is evidenced by a net outward motion of stars from the SMC centre along the direction towards the LMC and apparent tangential anisotropy at all radii. The latter is expected if the SMC is undergoing significiant tidal stripping, as we demonstrate using a suite of $N$-body simulations of the SMC/LMC system disrupting around the Milky Way. These results suggest that dynamical models for the SMC that assume a steady state will need to be revisited.


title: Optical Spectroscopy of nearby type1-LINERs

abstract: We present the highlights from our recent study of 22 local (z$<$0.025) type-1 LINERs from the Palomar Survey, on the basis of optical long-slit spectroscopic observations taken with TWIN/CAHA, ALFOSC/NOT and HST/STIS (Cazzoli et al. 2018). Our goals were threefold: (a) explore the AGN- nature of these LINERs by studying the broad (BLR-originated) H$\alpha$ component; (b) derive a reliable interpretation for the multiple narrow components of emission lines by studying their kinematics and ionisation mechanism (via standard BPTs); (c) probe the neutral gas in the nuclei of these LINERs for the first time. Hence, kinematics and fluxes of a set of emission lines, from H$\beta$ to [SII], and the NaD doublet in absorption have been modelled and measured, after the subtraction of the underlying light from the stellar component.


title: Dynamical evolution of fractal structures in star-forming regions

abstract: The Q-parameter is used extensively to quantify the spatial distributions of stars and gas in star-forming regions as well as older clusters and associations. It quantifies the amount of structure using the ratio of the average length of a minimum spanning tree, mbar, to the average length within the complete graph, sbar. The interpretation of the Q-parameter often relies on comparing observed values of Q, mbar and sbar to idealised synthetic geometries, where there is little or no match between the observed star-forming regions and the synthetic regions. We measure Q, mbar, and sbar over 10 Myr in N-body simulations which are compared to IC 348, NGC 1333, and the ONC. For each star-forming region we set up simulations that approximate their initial conditions for a combination of different virial rations and fractal dimensions. We find that dynamical evolution of idealised fractal geometries can account for the observed Q, mbar, and sbar values in nearby star-forming regions. In general, an initially fractal star-forming region will tend to evolve to become more smooth and centrally concentrated. However, we show that initial conditions, as well as where the edge of the region is defined, can cause significant differences in the path that a star-forming region takes across the mbar-sbar plot as it evolves. We caution that the observed Q-parameter should not be directly compared to idealised geometries. Instead, it should be used to determine the degree to which a star-forming region is either spatially substructured or smooth and centrally concentrated.


title: Optics-free imaging of complex, non-sparse QR-codes with Deep Neural Networks

abstract: We demonstrate optics-free imaging of complex QR-codes using a bare image sensor and a trained artificial neural network (ANN). The ANN is trained to interpret the raw sensor data for human visualization. The image sensor is placed at a specified gap from the QR code. We studied the robustness of our approach by experimentally testing the output of the ANNs with system perturbations of this gap, and the translational and rotational alignments of the QR code to the image sensor. Our demonstration opens us the possibility of using completely optics-free cameras for application-specific imaging of complex, non-sparse objects.


title: Full 5D characterisation of the Sagittarius stream with Gaia DR2 RR Lyrae

abstract: The Sagittarius stream is one of the best tools that we currently have to estimate the mass and shape of our Galaxy. However, assigning membership and obtaining the phase-space distribution of the stars that form the tails is quite challenging. Our goal is to produce a catalogue of RR Lyrae stars of Sagittarius and obtain an empiric measurement of the trends along the stream in sky position, distance and tangential velocities. We generate two initial samples from the Gaia DR2 RR Lyrae catalogue: one, selecting only the stars within \pm20deg of the orbital plane of Sagittarius (Strip) and the other, the result of applying the Pole Count Map (nGC3) algorithm. We then use the model-independent, deterministic method developed in this work to remove most of the contamination by detecting and isolating the stream in distance and proper motions. The output is two empiric catalogues: the Strip sample (higher-completeness, lower-purity) which contains 11 677 stars, and the nGC3 sample (higher-purity, lower-completeness) with 6 608 stars. We characterise the changes along the stream in all the available dimensions, the 5 astrometric ones plus the metallicity, covering more than 2pi rad in the sky and obtain new estimates for the apocentres and the mean [Fe/H] of the RR Lyrae population. Also, we show the first map of the two components of the tangential velocity, thanks to the combination of distances and proper motions. Finally, we detect the bifurcation in the leading arm and report no significant difference between the two branches, either in metallicity, kinematics or distance. We provide the largest sample of RR Lyrae candidates of Sagittarius, which can be used as an input for a spectroscopic follow-up or as a reference for the new generation of models of the stream through the interpolators in distance and velocity that we have constructed.


title: End-to-End Entity Linking and Disambiguation leveraging Word and Knowledge Graph Embeddings

abstract: Entity linking - connecting entity mentions in a natural language utterance to knowledge graph (KG) entities is a crucial step for question answering over KGs. It is often based on measuring the string similarity between the entity label and its mention in the question. The relation referred to in the question can help to disambiguate between entities with the same label. This can be misleading if an incorrect relation has been identified in the relation linking step. However, an incorrect relation may still be semantically similar to the relation in which the correct entity forms a triple within the KG; which could be captured by the similarity of their KG embeddings. Based on this idea, we propose the first end-to-end neural network approach that employs KG as well as word embeddings to perform joint relation and entity classification of simple questions while implicitly performing entity disambiguation with the help of a novel gating mechanism. An empirical evaluation shows that the proposed approach achieves a performance comparable to state-of-the-art entity linking while requiring less post-processing.


title: Integrable deformation of $\mathbb{CP}^n$ and generalised Kaehler geometry

abstract: We build on the results of arXiv:1912.11036 for generalised frame fields on generalised quotient spaces and study integrable deformations for $\mathbb{CP}^n$. In particular we show how, when the target space of the Principal Chiral Model is a complex projective space, a two-parameter deformation can be introduced in principle. The second parameter can however be removed via a diffeomorphism, which we construct explicitly, in accordance with the results stemming from a thorough integrability analysis we carry out. We also elucidate how the deformed target space can be seen as an instance of generalised Kaehler, or equivalently bi-Hermitian, geometry. In this respect, we find the generic form of the pure spinors for $\mathbb{CP}^n$ and the explicit expression for the generalised Kaehler potential for $n=1,2$.


title: Symmetry protection and giant Fermi arcs from multifold fermions in binary, ternary, and quaternary compounds

abstract: Higher-fold chiral fermions that go beyond two-fold Weyl fermions have recently been reported in crystalline systems. Here, we focus on such excitations in several binary, ternary and quaternary alloys/compounds with CoGe, BiSbPt and KMgBO3 as the representative examples that belong to the crystal space group (SG) 198. We found distinct three-fold, four-fold and six-fold chiral fermions in the bulk via Density Functional computations. We provide general symmetry arguments for the protection of these degeneracies with special emphasis on the four-fold fermions. Our surface spectra simulations show that the size of Fermi arcs resulting from these chiral fermions are large, robust and untouched from the bulk states due to the near absence of trivial bulk Fermi pockets. All these features make these systems -- especially CoGe and KMgBO3 -- promising topological semimetal candidates to realize higher-fold fermions in future photo-emission and transport experiments.


title: Quantum Algorithms for Simulating the Lattice Schwinger Model

abstract: The Schwinger model (quantum electrodynamics in 1+1 dimensions) is a testbed for the study of quantum gauge field theories. We give scalable, explicit digital quantum algorithms to simulate the lattice Schwinger model in both NISQ and fault-tolerant settings. In particular, we perform a tight analysis of low-order Trotter formula simulations of the Schwinger model, using recently derived commutator bounds, and give upper bounds on the resources needed for simulations in both scenarios. In lattice units, we find a Schwinger model on $N/2$ physical sites with coupling constant $x^{-1/2}$ and electric field cutoff $x^{-1/2}\Lambda$ can be simulated on a quantum computer for time $2xT$ using a number of $T$-gates or CNOTs in $\widetilde{O}( N^{3/2} T^{3/2} \sqrt{x} \Lambda )$ for fixed operator error. This scaling with the truncation $\Lambda$ is better than that expected from algorithms such as qubitization or QDRIFT. Furthermore, we give scalable measurement schemes and algorithms to estimate observables which we cost in both the NISQ and fault-tolerant settings by assuming a simple target observable---the mean pair density. Finally, we bound the root-mean-square error in estimating this observable via simulation as a function of the diamond distance between the ideal and actual CNOT channels. This work provides a rigorous analysis of simulating the Schwinger model, while also providing benchmarks against which subsequent simulation algorithms can be tested.


title: Analysis of lower bounds for quantum control times and their relation to the quantum speed limit

abstract: Limitations to the speed of evolution of quantum systems, typically referred to as quantum speed limits (QSLs), have important consequences for quantum control problems. However, in its standard formulation, is not straightforward to obtain meaningful QSL bounds for time-dependent Hamiltonians with unknown control parameters. In this paper we present a short introductory overview of quantum speed limit for unitary dynamics and its connection to quantum control. We then analyze potential methods for obtaining new bounds on control times inspired by the QSL. We finally extend the work in [Poggi, Lombardo and Wisniacki EPL 104 40005 (2013)] by studying the properties and limitations of these new bounds in the context of a driven two-level quantum system


title: Nonreciprocal Optomechanical Entanglement against Backscattering Losses

abstract: We propose how to achieve nonreciprocal quantum entanglement of light and motion and reveal its counterintuitive robustness against random losses. We find that by splitting the counterpropagating lights of a spinning resonator via the Sagnac effect, photons and phonons can be entangled strongly in a chosen direction but fully uncorrelated in the other. This makes it possible both to realize quantum nonreciprocity even in the absence of any classical nonreciprocity and also to achieve significant entanglement revival against backscattering losses in practical devices. Our work provides a way to protect and engineer quantum resources by utilizing diverse nonreciprocal devices, for building noise-tolerant quantum processors, realizing chiral networks, and backaction-immune quantum sensors.


title: Fluid critical behavior at liquid-gas phase transition: Analytic method for microscopic description

abstract: The behavior of fluids in the vicinity of the liquid-gas critical point is studied within the cell fluid model framework. The analytic method for deriving the equation of state of a cell fluid model in the low-temperature region (T<Tc) is developed using the renormalization group transformation within the collective variables approach. Mathematical description within the grand canonical ensemble is illustrated by an example of the Morse interaction potential possessing the Fourier transform. A specific feature of the proposed method lies in the possibility to use exclusively microscopic characteristics of a fluid (parameters of the interaction potential) for calculating macroscopic quantities (pressure and other thermodynamic quantities) without involving the hard-spheres reference system. The grand partition function, thermodynamic potential, and equation of state of the model near the critical point are derived taking into account the non-Gaussian (quartic) distribution of order parameter fluctuations. A nonlinear equation, which links the density and the chemical potential, is presented and solved. Graphs of the dependence of the density on the chemical potential are plotted for various values of the relative temperature. The numerical estimates of the critical-point parameters for potassium, obtained in addition to the estimates for sodium, are given. The calculated critical-point parameters for liquid alkali metals (sodium and potassium) are in accord with experimental data. The coexistence curve for sodium is plotted and compared with other authors' data in the immediate vicinity of Tc, where theoretical and experimental researches are difficult to carry out. The differences between the obtained results and the earlier published results for T>Tc are discussed.


title: Flavour anomalies from a split dark sector

abstract: We investigate solutions to the flavour anomalies in $B$ decays based on loop diagrams of a "split" dark sector characterised by the simultaneous presence of heavy particles at the TeV scale and light particles around and below the $B$-meson mass scale. We show that viable parameter space exists for solutions based on penguin diagrams with a vector mediator, while minimal constructions relying on box diagrams are in strong tension with the constraints from the LHC, LEP, and the anomalous magnetic moment of the muon. In particular, we highlight a regime where the mediator lies close to the $B$-meson mass, naturally realising a resonance structure and a $q^2$-dependent effective coupling. We perform a full fit to the relevant flavour observables and analyse the constraints from intensity frontier experiments. Besides new measurements of the anomalous magnetic moment of the muon, we find that decays of the $B$ meson, $B_s$-mixing, missing energy searches at Belle-II, and LHC searches for top/bottom partners can robustly test these scenarios in the near future.


title: TxSim:Modeling Training of Deep Neural Networks on Resistive Crossbar Systems

abstract: Resistive crossbars have attracted significant interest in the design of Deep Neural Network (DNN) accelerators due to their ability to natively execute massively parallel vector-matrix multiplications within dense memory arrays. However, crossbar-based computations face a major challenge due to a variety of device and circuit-level non-idealities, which manifest as errors in the vector-matrix multiplications and eventually degrade DNN accuracy. To address this challenge, there is a need for tools that can model the functional impact of non-idealities on DNN training and inference. Existing efforts towards this goal are either limited to inference, or are too slow to be used for large-scale DNN training. We propose TxSim, a fast and customizable modeling framework to functionally evaluate DNN training on crossbar-based hardware considering the impact of non-idealities. The key features of TxSim that differentiate it from prior efforts are: (i) It comprehensively models non-idealities during all training operations (forward propagation, backward propagation, and weight update) and (ii) it achieves computational efficiency by mapping crossbar evaluations to well-optimized BLAS routines and incorporates speedup techniques to further reduce simulation time with minimal impact on accuracy. TxSim achieves orders-of-magnitude improvement in simulation speed over prior works, and thereby makes it feasible to evaluate training of large-scale DNNs on crossbars. Our experiments using TxSim reveal that the accuracy degradation in DNN training due to non-idealities can be substantial (3%-10%) for large-scale DNNs, underscoring the need for further research in mitigation techniques. We also analyze the impact of various device and circuit-level parameters and the associated non-idealities to provide key insights that can guide the design of crossbar-based DNN training accelerators.


title: Fundamental Issues Regarding Uncertainties in Artificial Neural Networks

abstract: Artificial Neural Networks (ANNs) implement a specific form of multi-variate extrapolation and will generate an output for any input pattern, even when there is no similar training pattern. Extrapolations are not necessarily to be trusted, and in order to support safety critical systems, we require such systems to give an indication of the training sample related uncertainty associated with their output. Some readers may think that this is a well known issue which is already covered by the basic principles of pattern recognition. We will explain below how this is not the case and how the conventional (Likelihood estimate of) conditional probability of classification does not correctly assess this uncertainty. We provide a discussion of the standard interpretations of this problem and show how a quantitative approach based upon long standing methods can be practically applied. The methods are illustrated on the task of early diagnosis of dementing diseases using Magnetic Resonance Imaging.


title: Stochastic Makespan Minimization in Structured Set Systems

abstract: We study stochastic combinatorial optimization problems where the objective is to minimize the expected maximum load (a.k.a.\ the makespan). In this framework, we have a set of $n$ tasks and $m$ resources, where each task $j$ uses some subset of the resources. Tasks have random sizes $X_j$, and our goal is to non-adaptively select $t$ tasks to minimize the expected maximum load over all resources, where the load on any resource $i$ is the total size of all selected tasks that use $i$. For example, when resources are points and tasks are intervals in a line, we obtain an $O(\log\log m)$-approximation algorithm. Our technique is also applicable to other problems with some geometric structure in the relation between tasks and resources; e.g., packing paths, rectangles, and "fat" objects. Our approach uses a strong LP relaxation using the cumulant generating functions of the random variables. We also show that this LP has an $\Omega(\log^* m)$ integrality gap, even for the problem of selecting intervals on a line; here $\log^* m$ is the iterated logarithm function.


title: A metric version of Poincar\'e's theorem concerning biholomorphic inequivalence of domains

abstract: We show that if $Y_j\subset \mathbb{C}^{n_j}$ is a bounded strongly convex domain with $C^3$-boundary for $j=1,\dots,q$, and $X_j\subset \mathbb{C}^{m_j}$ is a bounded convex domain for $j=1,\ldots,p$, then the product domain $\prod_{j=1}^p X_j\subset \mathbb{C}^m$ cannot be isometrically embedded into $\prod_{j=1}^q Y_j\subset \mathbb{C}^n$ under the Kobayashi distance, if $p>q$. This result generalises Poincar\'e's theorem which says that there is no biholomorphic map from the polydisc onto the Euclidean ball in $\mathbb{C}^n$ for $n\geq 2$. The method of proof only relies on the metric geometry of the spaces and will be derived from a result for products of proper geodesic metric spaces with the sup-metric. In fact, the main goal of the paper is to establish a general criterion, in terms of certain asymptotic geometric properties of the individual metric spaces, that yields an obstruction for the existence of an isometric embedding between product metric spaces.


title: Experimental Study on the Aerodynamic Sealing of Air Curtains

abstract: Controlling the air quality is of the utmost importance in today buildings. Vertical air curtains are often used to separate two different climatic zones with a view to reduce heat transfer. In fact, this research work proposes an air curtain aimed to ensure a proper separation between two zones, a clean one and a contaminated one. The methodology of this research includes: (i) small-scale tests on water models to ensure that the contamination does not pass through the air curtain, and (ii) an analytical development integrating the main physical characteristics of plane jets. In the solution developed, the airflow is extracted from the contaminated compartment to reduce the curtain airflow rejected to the exterior of the compartment. In this research work, it was possible to determine the minimum exhaust flow necessary to ensure the aerodynamic sealing of the air curtain. This article addresses the methodology used to perform the small-scale water tests and the corresponding results.


title: In-Medium Jet Modification Measured by PHENIX Via two-particle Correlations and High $p_T$ Hadrons in $A+A$ Collisions

abstract: The first evidence of jet quenching was observed at RHIC via suppression of single high $p_T$ hadron $R_{AA}$ and the disappearance of the away-side jet peak in two-particle correlations. Since then, hadron $R_{AA}$ and two-particle correlations continue to be useful probes of the QGP in heavy ion collisions, since the particles involved are fragments of the jets produced in the initial hard scattering. PHENIX recently improved the width measurements extracted from $\pi^0$-hadron correlations after removing the higher order flow terms in the underlying event subtraction. Measurements of the away-side jet correlated with high $p_T$ neutral pions show an increase in low momentum particle production at wide angles consistent with theoretical expectations for energy loss. The system size dependence of energy loss is further investigated at PHENIX by measuring the absolute yield and $R_{AA}$ for various hadron species at high $p_T$ in several collision systems including $U+U$ and $Cu+Au$. These proceedings will present the newest PHENIX $R_{AA}$ and two-particle correlation measurements and their role in our understanding of jet quenching and medium response in heavy ion collisions


title: 2-Dimensional Palindromes with $k$ Mismatches

abstract: This paper extends the problem of 2-dimensional palindrome search into the area of approximate matching. Using the Hamming distance as the measure, we search for 2D palindromes that allow up to $k$ mismatches. We consider two different definitions of 2D palindromes and describe efficient algorithms for both of them. The first definition implies a square, while the second definition (also known as a \emph{centrosymmetric factor}), can be any rectangular shape. Given a text of size $n \times m$, the time complexity of the first algorithm is $O(nm (\log m + \log n + k))$ and for the second algorithm it is $O(nm(\log m + k) + occ)$ where $occ$ is the size of the output.


title: SHIFT: A Highly Realistic Financial Market Simulation Platform

abstract: This paper presents a new financial market simulator that may be used as a tool in both industry and academia for research in market microstructure. It allows multiple automated traders and/or researchers to simultaneously connect to an exchange-like environment, where they are able to asynchronously trade several financial assets at the same time. In its current iteration, this order-driven market implements the basic rules of U.S. equity markets, supporting both market and limit orders, and executing them in a first-in-first-out fashion. We overview the system architecture and we present possible use cases. We demonstrate how a set of automated agents is capable of producing a price process with characteristics similar to the statistics of real price from financial markets. Finally, we detail a market stress scenario and we draw, what we believe to be, interesting conclusions about crash events.


title: Smoothing Graphons for Modelling Exchangeable Relational Data

abstract: Modelling exchangeable relational data can be described by \textit{graphon theory}. Most Bayesian methods for modelling exchangeable relational data can be attributed to this framework by exploiting different forms of graphons. However, the graphons adopted by existing Bayesian methods are either piecewise-constant functions, which are insufficiently flexible for accurate modelling of the relational data, or are complicated continuous functions, which incur heavy computational costs for inference. In this work, we introduce a smoothing procedure to piecewise-constant graphons to form {\em smoothing graphons}, which permit continuous intensity values for describing relations, but without impractically increasing computational costs. In particular, we focus on the Bayesian Stochastic Block Model (SBM) and demonstrate how to adapt the piecewise-constant SBM graphon to the smoothed version. We initially propose the Integrated Smoothing Graphon (ISG) which introduces one smoothing parameter to the SBM graphon to generate continuous relational intensity values. We then develop the Latent Feature Smoothing Graphon (LFSG), which improves on the ISG by introducing auxiliary hidden labels to decompose the calculation of the ISG intensity and enable efficient inference. Experimental results on real-world data sets validate the advantages of applying smoothing strategies to the Stochastic Block Model, demonstrating that smoothing graphons can greatly improve AUC and precision for link prediction without increasing computational complexity.


title: Uncoupled Majorana fermions in open quantum systems: On the efficient simulation of non-equilibrium stationary states of quadratic Fermi models

abstract: A decomposition of the non-equilibrium stationary state of a quadratic Fermi system influenced by linear baths is obtained and used to establish a simulation protocol in terms of tensor states. The scheme is then applied to examine the occurrence of uncoupled Majorana fermions in Kitaev chains subject to baths on the ends. The resulting phase diagram is compared against the topological characterization of the equilibrium chain and the protocol efficiency is studied with respect to this model


title: Low-mass planet migration in three dimensional wind-driven inviscid discs: a negative corotation torque

abstract: We present simulations of low-mass planet--disc interactions in inviscid three-dimensional discs. We show that a wind-driven laminar accretion flow through the surface layers of the disc does not significantly modify the migration torque experienced by embedded planets. More importantly, we find that 3D effects lead to a dramatic change in the behaviour of the dynamical corotation torque compared to earlier 2D theory and simulations. Although it was previously shown that the dynamical corotation torque could act to slow and essentially stall the inward migration of a low-mass planet, our results in 3D show that the dynamical corotation torque has the complete opposite effect and speeds up inward migration. Our numerical experiments implicate buoyancy resonances as the cause. These have two effects: (i) they exert a direct torque on the planet, whose magnitude relative to the Lindblad torque is measured in our simulations to be small; (ii) they torque the gas librating on horseshoe orbits in the corotation region and drive evolution of its vortensity, leading to the negative dynamical corotation torque. This indicates that at low turbulent viscosity, the detailed vertical thermal structure of the protoplanetary disc plays an important role in determining the migration behaviour of embedded planets. If this result holds up under a more refined treatment of disc thermal evolution, then it has important implications for understanding the formation and early evolution of planetary systems.


title: On the Information Leakage of Camera Fingerprint Estimates

abstract: Camera fingerprints based on sensor PhotoResponse Non-Uniformity (PRNU) have gained broad popularity in forensic applications due to their ability to univocally identify the camera that captured a certain image. This fingerprint of a given sensor is extracted through some estimation method that requires a few images known to be taken with such sensor. In this paper, we show that the fingerprints extracted in this way leak a considerable amount of information from those images used in the estimation, thus constituting a potential threat to privacy. We propose to quantify the leakage via two measures: one based on the Mutual Information, and another based on the output of a membership inference test. Experiments with practical fingerprint estimators on a real-world image dataset confirm the validity of our measures and highlight the seriousness of the leakage and the importance of implementing techniques to mitigate it. Some of these techniques are presented and briefly discussed.


title: sBSNN: Stochastic-Bits Enabled Binary Spiking Neural Network with On-Chip Learning for Energy Efficient Neuromorphic Computing at the Edge

abstract: In this work, we propose stochastic Binary Spiking Neural Network (sBSNN) composed of stochastic spiking neurons and binary synapses (stochastic only during training) that computes probabilistically with one-bit precision for power-efficient and memory-compressed neuromorphic computing. We present an energy-efficient implementation of the proposed sBSNN using 'stochastic bit' as the core computational primitive to realize the stochastic neurons and synapses, which are fabricated in 90nm CMOS process, to achieve efficient on-chip training and inference for image recognition tasks. The measured data shows that the 'stochastic bit' can be programmed to mimic spiking neurons, and stochastic Spike Timing Dependent Plasticity (or sSTDP) rule for training the binary synaptic weights without expensive random number generators. Our results indicate that the proposed sBSNN realization offers possibility of up to 32x neuronal and synaptic memory compression compared to full precision (32-bit) SNN and energy efficiency of 89.49 TOPS/Watt for two-layer fully-connected SNN.


title: Topologically sensitive metaheuristics

abstract: This paper proposes topologically sensitive metaheuristics, and describes conceptual design of topologically sensitive Variable Neighborhood Search method (TVNS) and topologically sensitive Electromagnetism Metaheuristic (TEM).


title: Voronoi-based similarity distances between arbitrary crystal lattices

abstract: This paper develops a new continuous approach to a similarity between periodic lattices of ideal crystals. Quantifying a similarity between crystal structures is needed to substantially speed up the Crystal Structure Prediction, because the prediction of many target properties of crystal structures is computationally slow and is essentially repeated for many nearly identical simulated structures. The proposed distances between arbitrary periodic lattices of crystal structures are invariant under all rigid motions, satisfy the metric axioms and continuity under atomic perturbations. The above properties make these distances ideal tools for clustering and visualizing large datasets of crystal structures. All the conclusions are rigorously proved and justified by experiments on real and simulated crystal structures reported in the Nature 2017 paper "Functional materials discovery using energy-structure-function maps".


title: E2 distribution and statistical regularity in polygonal planar tessellations

abstract: From solar supergranulation to salt flat in Bolivia, from veins on leaves to cells on Drosophila wing discs, polygon-based networks exhibit great complexities, yet similarities persist and statistical distributions can be remarkably consistent. Based on analysis of 99 polygonal tessellations of a wide variety of physical origins, this work demonstrates the ubiquity of an exponential distribution in the squared norm of the deformation tensor, $E^{2}$, which directly leads to the ubiquitous presence of Gamma distributions in polygon aspect ratio. The $E^{2}$ distribution in turn arises as a $\chi^{2}$-distribution, and an analytical framework is developed to compute its statistics. $E^{2}$ is closely related to many energy forms, and its Boltzmann-like feature allows the definition of a pseudo-temperature. Together with normality in other key variables such as vertex displacement, this work reveals regularities universally present in all systems alike


title: Data-driven super-parameterization using deep learning: Experimentation with multi-scale Lorenz 96 systems and transfer-learning

abstract: To make weather/climate modeling computationally affordable, small-scale processes are usually represented in terms of the large-scale, explicitly-resolved processes using physics-based or semi-empirical parameterization schemes. Another approach, computationally more demanding but often more accurate, is super-parameterization (SP), which involves integrating the equations of small-scale processes on high-resolution grids embedded within the low-resolution grids of large-scale processes. Recently, studies have used machine learning (ML) to develop data-driven parameterization (DD-P) schemes. Here, we propose a new approach, data-driven SP (DD-SP), in which the equations of the small-scale processes are integrated data-drivenly using ML methods such as recurrent neural networks. Employing multi-scale Lorenz 96 systems as testbed, we compare the cost and accuracy (in terms of both short-term prediction and long-term statistics) of parameterized low-resolution (LR), SP, DD-P, and DD-SP models. We show that with the same computational cost, DD-SP substantially outperforms LR, and is better than DD-P, particularly when scale separation is lacking. DD-SP is much cheaper than SP, yet its accuracy is the same in reproducing long-term statistics and often comparable in short-term forecasting. We also investigate generalization, finding that when models trained on data from one system are applied to a system with different forcing (e.g., more chaotic), the models often do not generalize, particularly when the short-term prediction accuracy is examined. But we show that transfer-learning, which involves re-training the data-driven model with a small amount of data from the new system, significantly improves generalization. Potential applications of DD-SP and transfer-learning in climate/weather modeling and the expected challenges are discussed.

categories: physics.ao-ph nlin.CD physics.comp-ph physics.flu-dyn physics.geo-ph stat.ML


title: Perturbative Quantum Field Theory and Homotopy Algebras

abstract: We review the homotopy algebraic perspective on perturbative quantum field theory: classical field theories correspond to homotopy algebras such as $A_\infty$- and $L_\infty$-algebras. Furthermore, their scattering amplitudes are encoded in minimal models of these homotopy algebras at tree level and their quantum relatives at loop level. The translation between Lagrangian field theories and homotopy algebras is provided by the Batalin--Vilkovisky formalism. The minimal models are computed recursively using the homological perturbation lemma, which induces useful recursion relations for the computation of scattering amplitudes. After explaining how the homological perturbation lemma produces the usual Feynman diagram expansion, we use our techniques to verify an identity for the Berends--Giele currents which implies the Kleiss--Kuijf relations.


title: Unsupervised Discovery, Control, and Disentanglement of Semantic Attributes with Applications to Anomaly Detection

abstract: Our work focuses on unsupervised and generative methods that address the following goals: (a) learning unsupervised generative representations that discover latent factors controlling image semantic attributes, (b) studying how this ability to control attributes formally relates to the issue of latent factor disentanglement, clarifying related but dissimilar concepts that had been confounded in the past, and (c) developing anomaly detection methods that leverage representations learned in (a). For (a), we propose a network architecture that exploits the combination of multiscale generative models with mutual information (MI) maximization. For (b), we derive an analytical result (Lemma 1) that brings clarity to two related but distinct concepts: the ability of generative networks to control semantic attributes of images they generate, resulting from MI maximization, and the ability to disentangle latent space representations, obtained via total correlation minimization. More specifically, we demonstrate that maximizing semantic attribute control encourages disentanglement of latent factors. Using Lemma 1 and adopting MI in our loss function, we then show empirically that, for image generation tasks, the proposed approach exhibits superior performance as measured in the quality and disentanglement trade space, when compared to other state of the art methods, with quality assessed via the Frechet Inception Distance (FID), and disentanglement via mutual information gap. For (c), we design several systems for anomaly detection exploiting representations learned in (a), and demonstrate their performance benefits when compared to state-of-the-art generative and discriminative algorithms. The above contributions in representation learning have potential applications in addressing other important problems in computer vision, such as bias and privacy in AI.


title: Entanglement and the measurement problem

abstract: The entangled "measurement state" (MS), predicted by von Neumann to arise during quantum measurement, seems to display paradoxical properties such as multiple macroscopic outcomes. But analysis of interferometry experiments using entangled photon pairs shows that entangled states differ surprisingly from simple superposition states. Based on standard quantum theory, this paper shows that (i) the MS does not represent multiple detector readings but instead represents nonparadoxical multiple statistical correlations between system states and detector readings, (ii) exactly one outcome actually occurs, and (iii) when one outcome occurs, the other possible outcomes simultaneously collapse nonlocally. Point (iii) resolves an issue first raised in 1927 by Einstein who demonstrated that quantum theory requires instantaneous state collapse. This conundrum's resolution requires nonlocal correlations, which from today's perspective implies the MS must be an entangled state. Thus, contrary to previous presumed proofs of the measurement problem's insolubility, we find the MS to be the collapsed state and just what we expect upon measurement.


title: Dynamic Set Cover: Improved Amortized and Worst-Case Update Time

abstract: In the dynamic minimum set cover problem, a challenge is to minimize the update time while guaranteeing close to the optimal $\min(O(\log n), f)$ approximation factor. (Throughout, $m$, $n$, $f$, and $C$ are parameters denoting the maximum number of sets, number of elements, frequency, and the cost range.) In the high-frequency range, when $f=\Omega(\log n)$, this was achieved by a deterministic $O(\log n)$-approximation algorithm with $O(f \log n)$ amortized update time [Gupta et al. STOC'17]. In the low-frequency range, the line of work by Gupta et al. [STOC'17], Abboud et al. [STOC'19], and Bhattacharya et al. [ICALP'15, IPCO'17, FOCS'19] led to a deterministic $(1+\epsilon)f$-approximation algorithm with $O(f \log (Cn)/\epsilon^2)$ amortized update time. In this paper we improve the latter update time and provide the first bounds that subsume (and sometimes improve) the state-of-the-art dynamic vertex cover algorithms. We obtain: 1. $(1+\epsilon)f$-approximation ratio in $O(f\log^2 (Cn)/\epsilon^3)$ worst-case update time: No non-trivial worst-case update time was previously known for dynamic set cover. Our bound subsumes and improves by a logarithmic factor the $O(\log^3 n/\text{poly}(\epsilon))$ worst-case update time for unweighted dynamic vertex cover (i.e., when $f=2$ and $C=1$) by Bhattacharya et al. [SODA'17]. 2. $(1+\epsilon)f$-approximation ratio in $O\left((f^2/\epsilon^3)+(f/\epsilon^2) \log C\right)$ amortized update time: This result improves the previous $O(f \log (Cn)/\epsilon^2)$ update time bound for most values of $f$ in the low-frequency range, i.e. whenever $f=o(\log n)$. It is the first that is independent of $m$ and $n$. It subsumes the constant amortized update time of Bhattacharya and Kulkarni [SODA'19] for unweighted dynamic vertex cover (i.e., when $f = 2$ and $C = 1$).


title: A Sample Complexity Separation between Non-Convex and Convex Meta-Learning

abstract: One popular trend in meta-learning is to learn from many training tasks a common initialization for a gradient-based method that can be used to solve a new task with few samples. The theory of meta-learning is still in its early stages, with several recent learning-theoretic analyses of methods such as Reptile [Nichol et al., 2018] being for convex models. This work shows that convex-case analysis might be insufficient to understand the success of meta-learning, and that even for non-convex models it is important to look inside the optimization black-box, specifically at properties of the optimization trajectory. We construct a simple meta-learning instance that captures the problem of one-dimensional subspace learning. For the convex formulation of linear regression on this instance, we show that the new task sample complexity of any initialization-based meta-learning algorithm is $\Omega(d)$, where $d$ is the input dimension. In contrast, for the non-convex formulation of a two layer linear network on the same instance, we show that both Reptile and multi-task representation learning can have new task sample complexity of $\mathcal{O}(1)$, demonstrating a separation from convex meta-learning. Crucially, analyses of the training dynamics of these methods reveal that they can meta-learn the correct subspace onto which the data should be projected.


title: Quantum simulations of materials on near-term quantum computers

abstract: Quantum computers hold promise to enable efficient simulations of the properties of molecules and materials; however, at present they only permit ab initio calculations of a few atoms, due to a limited number of qubits. In order to harness the power of near-term quantum computers for simulations of larger systems, it is desirable to develop hybrid quantum-classical methods where the quantum computation is restricted to a small portion of the system. This is of particular relevance for molecules and solids where an active region requires a higher level of theoretical accuracy than its environment. Here we present a quantum embedding theory for the calculation of strongly-correlated electronic states of active regions, with the rest of the system described within density functional theory. We demonstrate the accuracy and effectiveness of the approach by investigating several defect quantum bits in semiconductors that are of great interest for quantum information technologies. We perform calculations on quantum computers and show that they yield results in agreement with those obtained with exact diagonalization on classical architectures, paving the way to simulations of realistic materials on near-term quantum computers.


title: TanksWorld: A Multi-Agent Environment for AI Safety Research

abstract: The ability to create artificial intelligence (AI) capable of performing complex tasks is rapidly outpacing our ability to ensure the safe and assured operation of AI-enabled systems. Fortunately, a landscape of AI safety research is emerging in response to this asymmetry and yet there is a long way to go. In particular, recent simulation environments created to illustrate AI safety risks are relatively simple or narrowly-focused on a particular issue. Hence, we see a critical need for AI safety research environments that abstract essential aspects of complex real-world applications. In this work, we introduce the AI safety TanksWorld as an environment for AI safety research with three essential aspects: competing performance objectives, human-machine teaming, and multi-agent competition. The AI safety TanksWorld aims to accelerate the advancement of safe multi-agent decision-making algorithms by providing a software framework to support competitions with both system performance and safety objectives. As a work in progress, this paper introduces our research objectives and learning environment with reference code and baseline performance metrics to follow in a future work.


title: Small-mass naked singularities censored by the Higgs field

abstract: We investigate a vacuum decay around an over-spinning naked singularity by using the Israel junction condition. We found that if the Higgs field develops the second minimum at higher energy scale, a spinning small-mass naked singularity could cause the vacuum decay around it within the cosmic age. An event horizon may form around the singularity due to the angular momentum transport from the singularity to a vacuum bubble wall. The newly formed event horizon leads to the increase of Bekenstein-Hawking entropy, which contributes to the enhancement of the vacuum decay rate. We conclude that small-mass naked singularities may be hidden by the event horizon within the cosmological time.


title: An elementary proof of 1D LSM theorems

abstract: The Lieb-Schultz-Mattis (LSM) theorem and its generalizations forbids the existence of a unique gapped ground state in the presence of certain lattice and internal symmetries and thus imposes powerful constraints on the low energy properties of quantum many-body systems. We provide an elementary proof of a class of generalized LSM theorems in 1D using matrix product state representations and the representation theory of groups.


title: Relating relative R\'enyi entropies and Wigner-Yanase-Dyson skew information to generalized multiple quantum coherences

abstract: Quantum coherence is a crucial resource for quantum information processing. By employing the language of coherence orders largely applied in NMR systems, quantum coherence has been currently addressed in terms of multiple quantum coherences (MQCs). Here we investigate the $\alpha$-MQCs, a novel class of multiple quantum coherences which is based on $\alpha$-relative purity, an information-theoretic quantifier analogous to quantum fidelity and closely related to R\'{e}nyi relative entropy of order $\alpha$. Our framework enables linking $\alpha$-MQCs to Wigner-Yanase-Dyson skew information (WYDSI), an asymmetry monotone finding applications in quantum thermodynamics and quantum metrology. Furthermore, we derive a family of bounds on $\alpha$-MQCs, particularly showing that $\alpha$-MQC define a lower bound to quantum Fisher information (QFI). We illustrate these ideas for quantum systems described by single-qubit states, two-qubit Bell-diagonal states, and a wide class of multiparticle mixed states. Finally, we investigate the time evolution of the $\alpha$-MQC spectrum and the overall signal of relative purity, by simulating the time reversal dynamics of a many-body all-to-all Ising Hamiltonian and comment on applications to physical platforms such as NMR systems, trapped ions, and ultracold atoms.


title: Limits on Absorption from a 332-MHz survey for Fast Radio Bursts

abstract: Fast Radio Bursts (FRBs) are bright, extragalactic radio pulses whose origins are still unknown. Until recently, most FRBs have been detected at frequencies greater than 1 GHz with a few exceptions at 800 MHz. The recent discoveries of FRBs at 400 MHz from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope has opened up possibilities for new insights about the progenitors while many other low frequency surveys in the past have failed to find any FRBs. Here, we present results from a FRB survey recently conducted at the Jodrell Bank Observatory at 332 MHz with the 76-m Lovell telescope for a total of 58 days. We did not detect any FRBs in the survey and report a 90$\%$ upper limit of 5500 FRBs per day per sky for a Euclidean Universe above a fluence threshold of 46 Jy ms. We discuss the possibility of absorption as the main cause of non-detections in low frequency (< 800 MHz) searches and invoke different absorption models to explain the same. We find that Induced Compton Scattering alone cannot account for absorption of radio emission and that our simulations favour a combination of Induced Compton Scattering and Free-Free Absorption to explain the non-detections. For a free-free absorption scenario, our constraints on the electron density are consistent with those expected in the post-shock region of the ionized ejecta in Super-Luminous SuperNovae (SLSNe).


title: On the Bertini regularity theorem for arithmetic varieties

abstract: Let $\mathcal{X}$ be a regular projective arithmetic variety equipped with an ample hermitian line bundle $\overline{\mathcal{L}}$. We prove that the proportion of global sections $\sigma$ with $\left\lVert \sigma \right\rVert_{\infty}<1$ of $\overline{\mathcal{L}}^{\otimes d}$ whose divisor does not have a singular point on the fiber $\mathcal{X}_p$ over any prime $p<e^{\varepsilon d}$ tends to $\zeta_{\mathcal{X}}(1+\dim \mathcal{X})^{-1}$ as $d\rightarrow \infty$.


title: Big Quantum cohomology of orbifold spheres

abstract: We construct a Kodaira-Spencer map from the big quantum cohomology of a sphere with three orbifold points to the Jacobian ring of the mirror Landau-Ginzburg potential function. This is constructed via the Lagrangian Floer theory of the Seidel Lagrangian and we show that Kodaira-Spencer map is a ring isomorphism.


title: Firing rate of the leaky integrate-and-fire neuron with stochastic conductance-based synaptic inputs with short decay times

abstract: We compute the firing rate of a leaky integrate-and-fire (LIF) neuron with stochastic conductance-based inputs in the limit when synaptic decay times are much shorter than the membrane time constant. A comparison of our analytical results to numeric simulations is presented for a range of biophysically-realistic parameters.


title: Information Directed Sampling for Linear Partial Monitoring

abstract: Partial monitoring is a rich framework for sequential decision making under uncertainty that generalizes many well known bandit models, including linear, combinatorial and dueling bandits. We introduce information directed sampling (IDS) for stochastic partial monitoring with a linear reward and observation structure. IDS achieves adaptive worst-case regret rates that depend on precise observability conditions of the game. Moreover, we prove lower bounds that classify the minimax regret of all finite games into four possible regimes. IDS achieves the optimal rate in all cases up to logarithmic factors, without tuning any hyper-parameters. We further extend our results to the contextual and the kernelized setting, which significantly increases the range of possible applications.


title: Arithmetic statistics on cubic surfaces

abstract: In this paper we compute the distributions of various markings on smooth cubic surfaces defined over the finite field $\mathbb{F}_q$, for example the distribution of pairs of points, `tritangents' or `double sixes'. We also compute the (rational) cohomology of certain associated bundles and covers over complex numbers.


title: The Sampled Moran Genealogy Process

abstract: We define the Sampled Moran Genealogy Process, a continuous-time Markov process on the space of genealogies with the demography of the classical Moran process, sampled through time. To do so, we begin by defining the Moran Genealogy Process using a novel representation. We then extend this process to include sampling through time. We derive exact conditional and marginal probability distributions for the sampled process under a stationarity assumption, and an exact expression for the likelihood of any sequence of genealogies it generates. This leads to some interesting observations pertinent to existing phylodynamic methods in the literature. Throughout, our proofs are original and make use of strictly forward-in-time calculations and are exact for all population sizes and sampling processes.


title: NOMA Design with Power-Outage Tradeoff for Two-User Systems

abstract: This letter proposes a modified non-orthogonal multiple-access (NOMA) scheme for systems with a multi-antenna base station (BS) and two single-antenna users, where NOMA transmissions are conducted only when the absolute correlation coefficient (CC) between the user channels exceeds a threshold and the BS uses matched-filter (MF) precoding along the user with the stronger average channel gain. We derive the average minimal transmit power to guarantee the signal-to-interference-plus-noise-ratio (SINR) levels of both users. Our results show that the average minimal power grows logarithmically in the reciprocal of the CC threshold and a non-zero threshold is necessary for the modified NOMA scheme to have finite average minimal transmit power. Further, for the massive MIMO scenario, we derive the scaling laws of the average transmit power and outage probability with respect to the antenna numbers, as well as their tradeoff law. Simulation results are shown to validate our theoretical results.


title: Comment on "Quantum Time Crystals and Interacting Gauge Theories in Atomic Bose-Einstein Condensates"

abstract: In a recent letter [Phys. Rev. Lett. 123, 250402], \"Ohberg and Wright describe a Bose-Einstein condensate trapped on a ring in the presence of the density-dependent gauge potential. It is claimed that the ground state of the system corresponds to a rotating chiral bright soliton and consequently it forms a genuine time crystal which minimizes its energy by performing periodic motion. We show that the energy of the chiral soliton in the laboratory frame is not correctly calculated in the letter. The correct energy becomes minimal if the soliton does not move.

categories: cond-mat.quant-gas cond-mat.mes-hall physics.atom-ph quant-ph


title: Analysis of Discriminator in RKHS Function Space for Kullback-Leibler Divergence Estimation

abstract: Several scalable sample-based methods to compute the Kullback Leibler (KL) divergence between two distributions have been proposed and applied in large-scale machine learning models. While they have been found to be unstable, the theoretical root cause of the problem is not clear. In this paper, we study a generative adversarial network based approach that uses a neural network discriminator to estimate KL divergence. We argue that, in such case, high fluctuations in the estimates are a consequence of not controlling the complexity of the discriminator function space. We provide a theoretical underpinning and remedy for this problem by first constructing a discriminator in the Reproducing Kernel Hilbert Space (RKHS). This enables us to leverage sample complexity and mean embedding to theoretically relate the error probability bound of the KL estimates to the complexity of the discriminator in RKHS. Based on this theory, we then present a scalable way to control the complexity of the discriminator for a reliable estimation of KL divergence. We support both our proposed theory and method to control the complexity of the RKHS discriminator through controlled experiments.


title: IoT Based Real Time Noise Mapping System for Urban Sound Pollution Study

abstract: This paper describes the development of a system that enables real time data visualization via a webapp regarding sound intensity using multiple node devices connected through internet. The prototypes were realized using ATmega328 (Arduino Nano) and ESP8266 hardware modules, NodeMCU Arduino wrapper library, Google maps and firebase API along with JavaScript webapp. System architecture is such that multiple node devices will be installed in different locations of the target area. On each node device, an Arduino Nano interfaced with a Sound Sensor measures the ambient sound intensity and ESP8266 Wi-Fi module transmits the data to a database via web API. On the webapp, it plots all the real-time data from the devices over Google maps according to the locations of the node devices. The logged data that is collected can then be used to carry out researches regarding sound pollution in targeted areas.


title: Low-temperature formation of platinum silicides on polycrystalline silicon

abstract: Phase formation and transitions in platinum silicide films on polycrystalline Si in the process of Pt magnetron sputtering and heat treatments at different temperatures are studied using X-ray photoelectron spectroscopy, high-resolution transmission electron microscopy and X-ray diffractometry. Phases of amorphous and polycrystalline Pt$_3$Si and Pt$_2$Si are shown to form during the room-temperature Pt deposition on poly-Si beneath the Pt layer. The relaxation of the interfacial film and partial transformation of Pt$_3$Si into Pt$_2$Si occurs as a result of the thermal treatment for 30 minutes in the temperature range from 125 to 300$^{\circ}$C. The Pt$_3$Si and Pt$_2$Si phases crystallize to PtSi due to annealing at the temperatures from 320 to 480$^{\circ}$C for the same time period.


title: A Reliable Millimeter-Wave Quadrature Interferometer

abstract: A simple, low-cost millimeter-wave (70 GHz) interferometer with a phase accuracy better than plus or minur 2 degrees, and a response time of 10 ns is described. The simplicity of this interferometer makes it ideal for measurement of electron densities in laboratory or semiconductor processing plasmas. The relative high accuracy and low cost are attained by a homodyne system with a quadrature intermediate frequency (IF) mixer, now commercially available in millimeter-wave frequency bands. The design and construction of magnetic shielding for the system isolator, which may be required around magnetically-confined plasmas, is also described.


title: Depletion Forces in Athermally Sheared Mixtures of Frictionless Disks and Rods in Two Dimensions

abstract: We carry out numerical simulations to study the behavior of an athermal mixture of frictionless circular disks and elongated rods in two dimensions, under three different types of global linear deformation at a finite strain rate: (i) simple shearing, (ii) pure shearing, and (iii) isotropic compression. We find that the fluctuations induced by such deformations lead to depletion forces that cause rods to group in parallel oriented clusters for the cases of simple and pure shear, but not for isotropic compression. For simple shearing, we find that as the fraction of rods increases, this clustering increases, leading to an increase in the average rate of rotation of the rods, and a decrease in the magnitude of their nematic ordering.


title: End-to-End Models for the Analysis of System 1 and System 2 Interactions based on Eye-Tracking Data

abstract: While theories postulating a dual cognitive system take hold, quantitative confirmations are still needed to understand and identify interactions between the two systems or conflict events. Eye movements are among the most direct markers of the individual attentive load and may serve as an important proxy of information. In this work we propose a computational method, within a modified visual version of the well-known Stroop test, for the identification of different tasks and potential conflicts events between the two systems through the collection and processing of data related to eye movements. A statistical analysis shows that the selected variables can characterize the variation of attentive load within different scenarios. Moreover, we show that Machine Learning techniques allow to distinguish between different tasks with a good classification accuracy and to investigate more in depth the gaze dynamics.


title: Computing the Relative Value of Spatio-Temporal Data in Wholesale and Retail Data Marketplaces

abstract: Spatio-temporal information is used for driving a plethora of intelligent transportation, smart-city, and crowd-sensing applications. Since data is now considered a valuable production factor, data marketplaces have appeared to help individuals and enterprises bring it to market to satisfy the ever-growing demand. In such marketplaces, several sources may need to combine their data in order to meet the requirements of different applications. In this paper we study the problem of estimating the relative value of different spatio-temporal datasets combined in wholesale and retail marketplaces for the purpose of predicting demand in metropolitan areas. Using as case studies large datasets of taxi rides from Chicago and New York, we ask questions such as "When does it make sense for different taxi companies to combine their data?", and "How should different companies be compensated for the data that they share?". We then turn our attention to the even harder problem of establishing the relative value of the data brought to retail marketplaces by individual drivers. Overall, we show that simplistic but popular approaches for estimating the relative value of data, such as using volume, or the ``leave-one-out'' heuristic, are inaccurate. Instead, more complex notions of value from economics and game-theory, such as the Shapley value need to be employed if one wishes to capture the complex effects of mixing different datasets on the accuracy of forecasting algorithms. Applying the Shapley value to large datasets from many sources is, of course, computationally challenging. We resort to structured sampling and manage to compute accurately the importance of thousands of data sources. We show that the relative value of the data held by different taxi companies and drivers may differ substantially, and that its relative ranking may change from district to district within a metropolitan area.


title: On the stability of uniform motion

abstract: We show that the uniform motion of a homogeneous distribution of electric charge can be stable or unstable depending on its geometry. When the electrodynamic body is perturbed from a state of rest, it starts to perform fast oscillations, irrespective of the frequency of the perturbation. This nonlinear oscillation is the result of the feedback interaction between Coulombian and radiative fields. The resulting spontaneous symmetry breaking of the Lorentz group implies that the principle of inertia only holds on average and suggests that the default state of matter is not necessarily uniform motion, but self-oscillation as well. We propose that the excitability of electrodynamic bodies under external perturbations, which leads to limit cycle oscillations, is at the basis of the wave particle duality and its related quantum effects.


title: Quantum Cognitive Triad. Semantic geometry of context representation

abstract: The paper describes an algorithm for semantic representation of behavioral contexts relative to a dichotomic decision alternative. The contexts are represented as quantum qubit states in two-dimensional Hilbert space visualized as points on the Bloch sphere. The azimuthal coordinate of this sphere functions as a one-dimensional semantic space in which the contexts are accommodated according to their subjective relevance to the considered uncertainty. The contexts are processed in triples defined by knowledge of a subject about a binary situational factor. The obtained triads of context representations function as stable cognitive structure at the same time allowing a subject to model probabilistically-variative behavior. The developed algorithm illustrates an approach for quantitative subjectively-semantic modeling of behavior based on conceptual and mathematical apparatus of quantum theory.


title: Three-dimensional phonon population anisotropy in silicon nanomembranes

abstract: Nanoscale single-crystals possess modified phonon dispersions due to the truncation of the crystal. The introduction of surfaces alters the population of phonons relative to the bulk and introduces anisotropy arising from the breaking of translational symmetry. Such modifications exist throughout the Brillouin zone, even in structures with dimensions of several nanometers, posing a challenge to the characterization of vibrational properties and leading to uncertainty in predicting the thermal, optical, and electronic properties of nanomaterials. Synchrotron x-ray thermal diffuse scattering studies find that freestanding Si nanomembranes with thicknesses as large as 21 nm exhibit a higher scattering intensity per unit thickness than bulk silicon. In addition, the anisotropy arising from the finite thickness of these membranes produces particularly intense scattering along reciprocal-space directions normal to the membrane surface compared to corresponding in-plane directions. These results reveal the dimensions at which calculations of materials properties and device characteristics based on bulk phonon dispersions require consideration of the nanoscale size of the crystal.


title: Behavior Cloning in OpenAI using Case Based Reasoning

abstract: Learning from Observation (LfO), also known as Behavioral Cloning, is an approach for building software agents by recording the behavior of an expert (human or artificial) and using the recorded data to generate the required behavior. jLOAF is a platform that uses Case-Based Reasoning to achieve LfO. In this paper we interface jLOAF with the popular OpenAI Gym environment. Our experimental results show how our approach can be used to provide a baseline for comparison in this domain, as well as identify the strengths and weaknesses when dealing with environmental complexity.


title: Photoemission spectrum in paramagnetic FeO under pressure: towards an \textit{ab-initio} description

abstract: In this work we provide an exhaustive study of the photemission spectrum of paramagnetic FeO under pressure using a refined version of our recently derived many-body effective energy theory (MEET). We show that, within a nonmagnetic description of the paramagnetic phase, the MEET gives an overall good description of the photoemission spectrum at ambient pressure as well as the changes it undergoes by increasing pressure. In particular at ambient pressure the band gap opens between the mixed Fe $t_{2g}$ and O $2p$ states and the Fe 4s states and, moreover, a $d$-$d$ gap opens, which is compatible with a high-spin configuration (hence nonzero local magnetic moments as observed in experiment), whereas decreasing pressure the band gap tends to close, $t_{2g}$ states tend to become fully occupied and $e_{g}$ fully unoccupied, which is compatible with a low-spin configuration (hence a collapse of the magnetic moments as observed in experiment). This is a remarkable result, since, within a nonmagnetic description of the paramagnetic phase, the MEET is capable to correctly describe the photoemission spectrum and the spin configuration at ambient as well as high pressure. For comparison we report the band gap values obtained using density-functional theory with a hybrid functional containing screened exchange (HSE06) and a variant of the $GW$ method (self-consistent COHSEX), which are reliable for the description of the antiferromagnetic phase. Both methods open a gap at ambient pressure, although, by construction, they give a low-spin configuration; increasing pressure they correctly describes the band gap closing. We also report the photoemission spectrum of the metallic phase obtained with one-shot fully-dynamical $GW$ on top of LDA, which gives a spectrum very similar to DMFT results from literature.


title: Expansivity and unique shadowing

abstract: Let $f\colon X\to X$ be a continuous function on a compact metric space. We show that shadowing is equivalent to backwards shadowing and two-sided shadowing when the map $f$ is onto. Using this we go on to show that, for expansive surjective maps the properties shadowing, two-sided shadowing, s-limit shadowing and two-sided s-limit shadowing are equivalent. We show that $f$ is positively expansive and has shadowing if and only if it has unique shadowing (i.e.\ each pseudo-orbit is shadowed by a unique point), extending a result implicit in Walter's proof that positively expansive maps with shadowing are topologically stable. We use the aforementioned result on two-sided shadowing to find an equivalent characterisation of shadowing and expansivity and extend these results to the notion of $n$-expansivity due to Morales.


title: Motor-free Contractility in Active Gels

abstract: Animal cells form contractile structures to promote various functions, from cell motility to cell division. Force generation in these structures is often due to molecular motors such as myosin that require polar substrates for their function. Here, we propose a motor-free mechanism that can generate contraction in biopolymer networks without the need for polarity. This mechanism is based on active binding/unbinding of crosslinkers that breaks the principle of detailed balance, together with the asymmetric force-extension response of semiflexible biopolymers. We find that these two ingredients can generate steady state contraction via a non-thermal, ratchet-like process. We calculate the resulting force-velocity relation using both coarse-grained and microscopic models.


title: Geometric Fusion via Joint Delay Embeddings

abstract: We introduce geometric and topological methods to develop a new framework for fusing multi-sensor time series. This framework consists of two steps: (1) a joint delay embedding, which reconstructs a high-dimensional state space in which our sensors correspond to observation functions, and (2) a simple orthogonalization scheme, which accounts for tangencies between such observation functions, and produces a more diversified geometry on the embedding space. We conclude with some synthetic and real-world experiments demonstrating that our framework outperforms traditional metric fusion methods.


title: Tetragonal mixed system $Cs_2CuCl_{4-x}Br_x$ complemented by the tetragonal phase realisation of $Cs_2CuCl_4$

abstract: Realisation of the tetragonal phase of $Cs_2CuCl_4$ is possible using specific crystal growth conditions at a temperature below $281K$. This work deals with the comparison of the magnetic susceptibility and the magnetization of this new tetragonal compound with the magnetic behaviour of tetragonal $Cs_2CuCl_{2.9}Br_{1.1}$, $Cs_2CuCl_{2.5}Br_{1.5}$, $Cs_2CuCl_{2.2}Br_{1.8}$ and presents consistent results for such quasi $2-D$ antiferromagnets. Structural investigation at low temperature for $Cs_2CuCl_{2.2}Br_{1.8}$ shows no phase transition. The structure remains in the tetragonal symmetry $I4/mmm$. Furthermore, several magnetic reflections corresponding to the propagation vector $k = (0, 0, 0)$ are observed for this tetragonal compound through neutron diffraction experiments below the magnetic phase transition at $T_N = 11.3K$ confirming its antiferromagnetic nature.


title: Interactive Summarizing -- Automatic Slide Localization Technology as Generative Learning Tool

abstract: Making a summary is a common learning strategy in lecture learning. It is an effective way for learners to engage in both traditional and video lectures. Video summarization is an effective technology applied to enhance learners' summarizing experience in a video lecture. In this article, we propose to apply cutting-edge automatic slide localization technology to lecture video learning experience. An interactive summarizing model is designed to explain how learners are engaged in the video lecture learning process supported by convolutional neural network and the possibility of related learning analytics.


title: Classical and Bayesian Analyses of a Mixture of Exponential and Lomax Distributions

abstract: The exponential and the Lomax distributions are widely used in life testing experiments in mixture models. A mixture model of exponential distribution and Lomax distribution is proposed. Parameters of the proposed model are estimated using classical and Bayesian procedures under type-I right censoring. Expressions for Bayes estimators are derived assuming noninformative (uniform and Jeffreys) priors under symmetric and asymmetric loss functions. Posterior predictive distributions of a future observation are derived and predictive estimates are obtained. Extensive Monte Carlo simulations are carried out to investigate performance of the estimators in terms of sample sizes, censoring times and mixing proportions. The analysis of mixture model is carried out using a data set of lifetime of transmitter receivers. Interesting properties of estimators are observed and discussed.


title: The firefighter problem on polynomial and intermediate growth groups

abstract: We prove that any Cayley graph $G$ with degree $d$ polynomial growth does not satisfy $\{f(n)\}$-containment for any $f=o(n^{d-2})$. This settles the asymptotic behaviour of the firefighter problem on such graphs as it was known that $Cn^{d-2}$ firefighters are enough, answering and strengthening a conjecture of Develin and Hartke. We also prove that intermediate growth Cayley graphs do not satisfy polynomial containment, and give explicit lower bounds depending on the growth rate of the group. These bounds can be further improved when more geometric information is available, such as for Grigorchuk's group.


title: Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics

abstract: Hierarchical code coupling strategies make it possible to combine the results of individual numerical solvers into a self-consistent symplectic solution. We explore the possibility of allowing such a coupling strategy to be non-intrusive. In that case, the underlying numerical implementation is not affected by the coupling itself, but its functionality is carried over in the interface. This method is efficient for solving the equations of motion for a self-gravitating system over a wide range of scales. We adopt a dedicated integrator for solving each particular part of the problem and combine the results to a self-consistent solution. In particular, we explore the possibilities of combining the evolution of one or more microscopic systems that are embedded in a macroscopic system. The here presented generalizations of Bridge include higher-order coupling strategies (from the classic 2nd order up to 10th-order), but we also demonstrate how multiple bridges can be nested and how additional processes can be introduced at the bridge time-step to enrich the physics, for example by incorporating dissipative processes. Such augmentation allows for including additional processes in a classic Newtonian N-body integrator without alterations to the underlying code. These additional processes include for example the Yarkovsky effect, dynamical friction or relativistic dynamics. Some of these processes operate on all particles whereas others apply only to a subset. The presented method is non-intrusive in the sense that the underlying methods remain operational without changes to the code (apart from adding the get- and set-functions to enable the bridge operator). As a result, the fundamental integrators continue to operate with their internal time step and preserve their local optimizations and parallelism. ... abridged ...


title: Measurement of Identical Particle Entanglement and the Influence of Antisymmetrisation

abstract: We explore the relationship between symmetrisation and entanglement through measurements on few-particle systems in a multi-well potential. In particular, considering two or three trapped atoms, we measure and distinguish correlations arising from two different physical origins: antisymmetrisation of the fermionic wavefunction and interaction between particles. We quantify this through the entanglement negativity of states, and the introduction of an antisymmetric negativity, which allows us to understand the role that symmetrisation plays in the measured entanglement properties. We apply this concept both to pure theoretical states and to experimentally reconstructed density matrices of two or three mobile particles in an array of optical tweezers.


title: Automatic Segmentation, Feature Extraction and Comparison of Healthy and Stroke Cerebral Vasculature

abstract: Accurate segmentation of cerebral vasculature and a quantitative assessment of cerebrovascular morphology is critical to various diagnostic and therapeutic purposes and is pertinent to studying brain health and disease. However, this is still a challenging task due to the complexity of the vascular imaging data. We propose an automated method for cerebral vascular segmentation without the need of any manual intervention as well as a method to skeletonize the binary volume to extract vascular geometric features which can characterize vessel structure. We combine a probabilistic vessel-enhancing filtering with an active-contour technique to segment magnetic resonance and computed tomography angiograms (MRA and CTA) and subsequently extract the vessel centerlines and diameters to calculate the geometrical properties of the vasculature. Our method was validated using a 3D phantom of the Circle-of-Willis region with 84% mean Dice Similarity and 85% mean Pearson Correlation with minimal modified Hausdorff distance error. We applied this method to a dataset of healthy subjects and stroke patients and present a quantitative comparison between them. We found significant differences in the geometric features including total length (2.88 +/- 0.38 m for healthy and 2.20 +/- 0.67 m for stroke), volume (40.18 +/- 25.55 ml for healthy and 34.43 +/- 21.83 ml for stroke), tortuosity (3.24 +/- 0.88 rad/cm for healthy and 5.80 +/- 0.92 rad/cm for stroke) and fractality (box dimension 1.36 +/- 0.28 for healthy vs. 1.69 +/- 0.20 for stroke). This technique can be applied on any imaging modality and can be used in the future to automatically obtain the 3D segmented vasculature for diagnosis and treatment planning of Stroke and other cerebrovascular diseases (CVD) in the clinic and also to study the morphological changes caused by various CVD.


title: Zero-inertia Systems: Sufficient Conditions for Phasor Modeling

abstract: Time-domain simulations are a critical tool for power system operators. Depending on the instability mechanism under consideration and the system characteristics, such as the time constants of controllers, either phasor or Electro-Magnetic Transient (EMT) models should be employed. On the one hand, EMT models provide a detailed-modeling of the system dynamics, thus increase the reliability of stability analysis; on the other end, using these models increase the computational times of simulations, slowing down the security assessment process. To decrease computational time, system operators could resort to phasor-mode simulations for a (hopefully large) subset of disturbances. This paper investigates the appropriateness of phasor-approximation models on simulating events related to power supply and balance stability in zero-inertia systems. First, the stability boundaries, which each model is able to identify, are analyzed; then sufficient conditions for control parameters are derived, which allow using phasor-approximation models to monitor power sharing among grid-forming converter-based resources. Time-domain simulations are performed in PowerFactory DigSilent to verify the results.


title: Mobility and Blockage-aware Communications in Millimeter-Wave Vehicular Networks

abstract: Mobility may degrade the performance of next-generation vehicular networks operating at the millimeter-wave spectrum: frequent mis-alignment and blockages require repeated beam training and handover, and incur enormous overhead. Nevertheless, mobility induces temporal correlations in the communication beams and in blockage events. In this paper, an adaptive design of beam training, data transmission and handover is proposed, that learns and exploits these temporal correlations to reduce the beam training overhead and optimally trade-off throughput and power consumption. At each time-slot, the serving base station (BS) decides to perform either beam training, data communication, or handover when blockage is detected, under uncertainty in the system state. The decision problem is cast as a partially observable Markov decision process, and the goal is to maximize the throughput delivered to the UE, under an average power constraint. To address the high dimensional optimization, an approximate constrained point-based value iteration (C-PBVI) method is developed, which simultaneously optimizes the primal and dual functions to meet the power constraint. Numerical results demonstrate a good match between the analysis and a simulation based on 2D mobility and 3D analog beamforming via uniform planar arrays at both BSs and UE, and reveal that C-PBVI performs near-optimally, and outperforms a baseline scheme with periodic beam training by 38% in spectral efficiency. Motivated by the structure of the C-PBVI policy, two heuristics are proposed, that trade complexity with sub-optimality, and achieve only 4% and 15% loss in spectral efficiency.


title: Hours Worked and the U.S. Distribution of Real Annual Earnings 1976-2019

abstract: We examine the impact of annual hours worked on annual earnings by decomposing changes in the real annual earnings distribution into composition, structural and hours effects. We do so via a nonseparable simultaneous model of hours, wages and earnings. Using the Current Population Survey for the survey years 1976--2019, we find that changes in the female distribution of annual hours of work are important in explaining movements in inequality in female annual earnings. This captures the substantial changes in their employment behavior over this period. Movements in the male hours distribution only affect the lower part of their earnings distribution and reflect the sensitivity of these workers' annual hours of work to cyclical factors.


title: Computationally efficient transfinite patches with fullness control

abstract: Transfinite patches provide a simple and elegant solution to the problem of representing non-four-sided continuous surfaces, which are useful in a variety of applications, such as curve network based design. Real-time responsiveness is essential in this context, and thus reducing the computation cost is an important concern. The Midpoint Coons (MC) patch presented in this paper is a fusion of two previous transfinite schemes, combining the speed of one with the superior control mechanism of the other. This is achieved using a new constrained parameterization based on generalized barycentric coordinates and transfinite blending functions.