title: Active Brownian particles: mapping to equilibrium polymers and exact computation of moments

abstract: It is well known that path probabilities of Brownian motion correspond to the equilibrium configurational probabilities of flexible Gaussian polymers, while those of active Brownian motion correspond to in-extensible semiflexible polymers. Here we investigate the properties of the equilibrium polymer that corresponds to the trajectories of particles acted on simultaneously by both Brownian as well as active noise. Through this mapping we can see interesting crossovers in mechanical properties of the polymer with changing contour length. The polymer end-to-end distribution exhibits Gaussian behaviour for short lengths, which changes to the form of semiflexible filaments at intermediate lengths, to finally go back to a Gaussian form for long contour lengths. By performing a Laplace transform of the governing Fokker-Planck equation of the active Brownian particle, we discuss a direct method to derive exact expressions for all the moments of the relevant dynamical variables, in arbitrary dimensions. These are verified via numerical simulations and used to describe interesting qualitative features such as, for example, dynamical crossovers. Finally we discuss the kurtosis of the ABP's position which we compute exactly and show that it can be used to differentiate between active Brownian particles and active Ornstein-Uhlenbeck process.


title: High quality factor photonic cavity for dark matter axion searches

abstract: Searches for dark matter axion involve the use of microwave resonant cavities operating in a strong magnetic field. Detector sensitivity is directly related to the cavity quality factor, which is limited, however, by the presence of the external magnetic field. In this paper we present a cavity of novel design whose quality factor is not affected by a magnetic field. It is based on a photonic structure by the use of sapphire rods. The quality factor at cryogenic temperature is in excess of $5 \times 10^5$ for a selected mode.


title: Restricted Hilbert transform for non-Hermitian management of fields

abstract: Non-Hermitian systems exploiting the synergy between gain and loss have recently become the focus of interest to discover novel physical phenomena. The spatial symmetry breaking in such systems allows tailoring the wave propagation at will. Inspired by such property, we propose a feasible approach based on local Hilbert transform to control the field flows in two- or higher dimensional non-Hermitian systems, restricting the complex refractive index within practical limits. We propose an iterative procedure to reduce the dimensionality of complex refractive index parameter space to two, one or zero dimensions. The proposed method provides a flexible way to systematically design locally PT-symmetric systems realizable with a limited collection of realistic materials.


title: Minimality and identifiability of discrete-time SARX systems

abstract: The paper studies the problem of minimality and identifiability for Switched Auto-Regressive eXogenous (abbreviated by SARX) systems. We propose formal definitions of the concepts of identifiability and minimality for SARX models. Based on these formalizations, we derive conditions for minimality and identifiability of SARX systems. In particular, we show that polynomially parametrized SARX systems are generically identifiable.


title: Study Zitterbewegung effect in a Quasi One-dimensional Relativistic Quantum Plasma by DHW formalization

abstract: Using Dirac equation together with the Wigner distribution function,the trembling motion,known as Zitterbewegung effect,of moving electrons in quasi-one-dimensional relativistic quantum plasma is theoretically investigated.The relativistic Wigner matrix is used to calculate the mean values of the position and velocity operators for a Dirac gas of electrons.It is found that the oscillatory behavior of measurable quantities could be associated with the Zitterbewegung effect which manifests itself as a damping interference pattern stemming from mixing the positive and negative dispersion modes of Dirac particles.

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


title: Metric preheating and radiative decay in single-field inflation

abstract: At the end of inflation, the inflaton oscillates at the bottom of its potential and these oscillations trigger a parametric instability for scalar fluctuations with wavelength $\lambda$ comprised in the instability band $(3H m)^{-1/2} <\lambda < H^{-1}$, where $H$ is the Hubble parameter and $m$ the curvature of the potential at its minimum. This "metric preheating" instability, which proceeds in the narrow resonance regime, leads to various interesting phenomena such as early structure formation, production of gravitational waves and formation of primordial black holes. In this work we study its fate in the presence of interactions with additional degrees of freedom, in the form of perturbative decay of the inflaton into a perfect fluid. Indeed, in order to ensure a complete transition from inflation to the radiation-dominated era, metric preheating must be considered together with perturbative reheating. We find that the decay of the inflaton does not alter the instability structure until the fluid dominates the universe content. As an application, we discuss the impact of the inflaton decay on the production of primordial black holes from the instability. We stress the difference between scalar field and perfect fluid fluctuations and explain why usual results concerning the formation of primordial black holes from perfect fluid inhomogeneities cannot be used, clarifying some recent statements made in the literature.


title: Forecasts for detecting the gravitational-wave memory effect with Advanced LIGO and Virgo

abstract: The detection of gravitational waves (GWs) from binary black holes (BBHs) has allowed the theory of general relativity to be tested in a previously unstudied regime: that of strong curvature and high GW luminosities. One distinctive and measurable effect associated with this aspect of the theory is the nonlinear GW memory effect. The GW memory effect is characterized by its effect on freely falling observers: the proper distance between their locations differs before and after a burst of GWs passes by their locations. Gravitational-wave interferometers, like the LIGO and Virgo detectors, can measure features of this effect from a single BBH merger, but previous work has shown that it will require an event that is significantly more massive and closer than any previously detected GW event. Finding evidence for the GW memory effect within the entire population of BBH mergers detected by LIGO and Virgo is more likely to occur sooner. A prior study has shown that the GW memory effect could be detected in a population of BBHs consisting of binaries like the first GW150914 event after roughly one-hundred events. In this paper, we compute forecasts of the time it will take the advanced LIGO and Virgo detectors (when the detectors are operating at their design sensitivities) to find evidence for the GW memory effect in a population of BBHs that is consistent with the measured population of events in the first two observing runs of the LIGO detectors. We find that after five years of data collected by the advanced LIGO and Virgo detectors the signal-to-noise ratio for the nonlinear GW memory effect in the population will be about three (near a previously used threshold for detection). We point out that the different approximation methods used to compute the GW memory effect can lead to notably different signal-to-noise ratios.


title: Comparing clusterings and numbers of clusters by aggregation of calibrated clustering validity indexes

abstract: A key issue in cluster analysis is the choice of an appropriate clustering method and the determination of the best number of clusters. Different clusterings are optimal on the same data set according to different criteria, and the choice of such criteria depends on the context and aim of clustering. Therefore, researchers need to consider what data analytic characteristics the clusters they are aiming at are supposed to have, among others within-cluster homogeneity, between-clusters separation, and stability. Here, a set of internal clustering validity indexes measuring different aspects of clustering quality is proposed, including some indexes from the literature. Users can choose the indexes that are relevant in the application at hand. In order to measure the overall quality of a clustering (for comparing clusterings from different methods and/or different numbers of clusters), the index values are calibrated for aggregation. Calibration is relative to a set of random clusterings on the same data. Two specific aggregated indexes are proposed and compared with existing indexes on simulated and real data.


title: A Robust Sensorless Controller-Observer Strategy for PMSMs with Unknown Resistance and Mechanical Model

abstract: In this work, we present a mixed sensorless strategy for Permanent Magnets Synchronous Machines, combining a torque/current controller and an observer for position, speed, flux and stator resistance. The proposed co-design is motivated by the need of an appropriate signal injection technique, in order to guarantee full state observability. Neither the typical constant or slowly-varying speed assumptions, nor a priori mechanical model information are used in the observer design. Instead, the rotor speed is modeled as an unknown input disturbance with constant (unknown) sign and uniformly non-zero magnitude. With the proposed architecture, it is shown that the torque tracking and signal injection tasks can be achieved and asymptotically decoupled. Because of these features, we refer to this strategy as a sensorless controller-observer with no mechanical model. Employing a gradient descent resistance/back-EMF estimation, combined with the unit circle formalism to describe the rotor position, we rigorously prove regional practical asymptotic stability of the overall structure, with a domain of attraction that can be made arbitrarily large, not including a lower dimensional manifold. The effectiveness of this design is further validated with numerical simulations, related to a challenging application of UAV propellers control.


title: Discontinuous Constituent Parsing with Pointer Networks

abstract: One of the most complex syntactic representations used in computational linguistics and NLP are discontinuous constituent trees, crucial for representing all grammatical phenomena of languages such as German. Recent advances in dependency parsing have shown that Pointer Networks excel in efficiently parsing syntactic relations between words in a sentence. This kind of sequence-to-sequence models achieve outstanding accuracies in building non-projective dependency trees, but its potential has not been proved yet on a more difficult task. We propose a novel neural network architecture that, by means of Pointer Networks, is able to generate the most accurate discontinuous constituent representations to date, even without the need of Part-of-Speech tagging information. To do so, we internally model discontinuous constituent structures as augmented non-projective dependency structures. The proposed approach achieves state-of-the-art results on the two widely-used NEGRA and TIGER benchmarks, outperforming previous work by a wide margin.


title: Impact of Atmospheric and Model Physics Perturbations On a High-Resolution Ensemble Data Assimilation System of the Red Sea

abstract: The Ensemble Adjustment Kalman Filter (EAKF) of the Data Assimilation Research Testbed (DART) is implemented to assimilate observations of satellite sea surface temperature, altimeter sea surface height and in situ ocean temperature and salinity profiles into an eddy-resolving 4 km Massachusetts Institute of Technology general circulation model (MITgcm) of the Red Sea. We investigate the impact of three different ensemble generation strategies (1) Iexp, uses ensemble of ocean states to initialize the model on 1st January, 2011 and inflates filter error covariance by 10 percent, (2) IAexp, adds ensemble of atmospheric forcing to Iexp, and (3) IAPexp, adds perturbed model physics to IAexp. The assimilation experiments are run for one year, starting from the same initial ensemble and assimilating data every three days. Results demonstrate that the Iexp mainly improved the model outputs with respect to assimilation free MITgcm run in the first few months, before showing signs of dynamical imbalances in the ocean estimates, particularly in the data-sparse subsurface layers. The IAexp yielded substantial improvements throughout the assimilation period with almost no signs of imbalances, including the subsurface layers. It further well preserved the model mesoscale features resulting in an improved forecasts for eddies, both in terms of intensity and location. Perturbing model physics in IAPexp slightly improved the forecast statistics and also the placement of basin scale eddies. Increasing hydrographic coverage further improved the results of IAPexp compared to IAexp in the subsurface layers. Switching off multiplicative inflation in IAexp and IAPexp leads to further improvements, especially in the subsurface layers.


title: Long-time asymptotics of the one-dimensional damped nonlinear Klein-Gordon equation

abstract: For the one-dimensional nonlinear damped Klein-Gordon equation \[ \partial_{t}^{2}u+2\alpha\partial_{t}u-\partial_{x}^{2}u+u-|u|^{p-1}u=0 \quad \mbox{on $\mathbb{R}\times\mathbb{R}$,}\] with $\alpha>0$ and $p>2$, we prove that any global finite energy solution either converges to $0$ or behaves asymptotically as $t\to \infty$ as the sum of $K\geq 1$ decoupled solitary waves. In the multi-soliton case $K\geq 2$, the solitary waves have alternate signs and their distances are of order $\log t$.


title: Analyzing the Dependency of ConvNets on Spatial Information

abstract: Intuitively, image classification should profit from using spatial information. Recent work, however, suggests that this might be overrated in standard CNNs. In this paper, we are pushing the envelope and aim to further investigate the reliance on spatial information. We propose spatial shuffling and GAP+FC to destroy spatial information during both training and testing phases. Interestingly, we observe that spatial information can be deleted from later layers with small performance drops, which indicates spatial information at later layers is not necessary for good performance. For example, test accuracy of VGG-16 only drops by 0.03% and 2.66% with spatial information completely removed from the last 30% and 53% layers on CIFAR100, respectively. Evaluation on several object recognition datasets (CIFAR100, Small-ImageNet, ImageNet) with a wide range of CNN architectures (VGG16, ResNet50, ResNet152) shows an overall consistent pattern.


title: Lipschitz stability at the boundary for time-harmonic diffuse optical tomography

abstract: We study the inverse problem in Optical Tomography of determining the optical properties of a medium $\Omega\subset\mathbb{R}^n$, with $n\geq 3$, under the so-called diffusion approximation. We consider the time-harmonic case where $\Omega$ is probed with an input field that is modulated with a fixed harmonic frequency $\omega=\frac{k}{c}$, where $c$ is the speed of light and $k$ is the wave number. We prove a result of Lipschitz stability of the absorption coefficient $\mu_a$ at the boundary $\partial\Omega$ in terms of the measurements in the case when the scattering coefficient $\mu_s$ is assumed to be known and $k$ belongs to certain intervals depending on some a-priori bounds on $\mu_a$, $\mu_s$.


title: Light-by-Light Scattering and Spacetime Noncommutativity

abstract: The aim of this article is to explore a potential usability of a photon-photon self-interaction from the noncommutative quantum electrodynamics (NCQED) in the case of light-by-light scattering ($\gamma\gamma\to\gamma\gamma$) in ultraperipheral Pb+Pb collisions, a reaction measured recently by ATLAS and planned for future experiments in hadron-hadron colliders. We compute the total cross section from both, the full one-loop standard model (SM) and the tree-level NCQED amplitudes, in the equivalent photon approximation with impact parameters, for various noncommutative scales, $\rm\Lambda_{NC}$, and incoming nuclear spin-energy combinations. We find that NCQED contribution to the cross section has considerable increase at diphoton invariant mass range higher than $\rm\Lambda_{NC}$, while the SM contribution is strongly suppressed in such region. Our results show that the current ATLAS $\rm\sqrt{s_{NN}} = 5.02$ TeV experiment can only probe $\rm\Lambda_{NC}<100$ GeV region. On the other hand, future hadron-hadron collider proposals could have the potential to extend to $\rm\Lambda_{NC}\lesssim 300$ GeV region, making the performance of $\rm Pb+Pb(\gamma\gamma)\to Pb+Pb\gamma\gamma$ scattering on testing space-time noncommutativity close to that of the previously proposed photon-photon mode of linear electron-positron collider.


title: Divergence-preserving reconstructions on polygons and a really pressure-robust virtual element method for the Stokes problem

abstract: Non divergence-free discretisations for the incompressible Stokes problem may suffer from a lack of pressure-robustness characterised by large discretisations errors due to irrotational forces in the momentum balance. This paper argues that also divergence-free virtual element methods (VEM) on polygonal meshes are not really pressure-robust as long as the right-hand side is not discretised in a careful manner. To be able to evaluate the right-hand side for the test functions, some explicit interpolation of the virtual test functions is needed that can be evaluated pointwise everywhere. The standard discretisation via an $L^2$-best approximation does not preserve the divergence and so destroys the orthogonality between divergence-free test functions and possibly eminent gradient forces in the right-hand side. To repair this orthogonality and restore pressure-robustness another divergence-preserving reconstruction is suggested based on Raviart-Thomas approximations on local subtriangulations of the polygons. All findings are proven theoretically and are demonstrated numerically in two dimensions. The construction is also interesting for hybrid high-order methods on polygonal or polyhedral meshes.


title: Vertical bonding distances and interfacial band structure of PTCDA on a Sn-Ag surface alloy

abstract: Molecular materials enable a vast variety of functionalities for novel electronic and spintronic devices. The unique possibility to alter or substitute organic molecules or metallic substrates offers the opportunity to modify and optimize interfacial properties for almost any desired field of application. For this reason, we extend the successful approach to control molecular interfaces by surface alloying. We present a comprehensive characterization of the structural and electronic properties of the interface formed between the prototypical molecule PTCDA and a Sn-Ag surface alloy grown on an Ag(111) single crystal surface. We monitor the changes of adsorption height of the surface alloy atoms and electronic valence band structure upon adsorption of one layer of PTCDA using the normal incidence x-ray standing wave technique in combination with momentum-resolved photoelectron spectroscopy. We find that the vertical buckling and the surface band structure of the SnAg$_2$ surface alloy is not altered by the adsorption of one layer of PTCDA, in contrast to our recent study of PTCDA on a PbAg$_2$ surface alloy [Phys. Rev. Lett. 117, 096805 (2016)] . In addition, the vertical adsorption geometry of PTCDA and the interfacial energy level alignment indicate the absence of any chemical interaction between the molecule and the surface alloy. We attribute the different interactions at these PTCDA/surface alloy interfaces to the presence or absence of local $\sigma$-bonds between the PTCDA oxygen atoms and the surface atoms. Combining our findings with results from literature, we are able to propose an empiric rule for engineering the surface band structure of alloys by adsorption of organic molecules.


title: Resonance libration and width at arbitrary inclination

abstract: We apply the analytical disturbing function for arbitrary inclination derived in our previous work to characterize resonant width and libration of mean motion resonances at arbitrary inclination obtained from direct numerical simulations of the three-body problem. We examine the 2:1 and 3:1 inner Jupiter and 1:2 and 1:3 outer Neptune resonances and their possible asymmetric librations using a new analytical pendulum model of resonance that includes the simultaneous libration of multiple arguments and their second harmonics. The numerically-derived resonance separatrices are obtained using the Mean Exponential Growth factor of Nearby Orbits (MEGNO chaos indicator). We find that the analytical and numerical estimates are in agreement and that resonance width is determined by the first few fundamental resonance modes that librate simultaneously on the resonant timescale. Our results demonstrate that the new pendulum model may be used to ascertain resonance width analytically, and more generally, that the disturbing function for arbitrary inclination is a powerful analytical tool that describes resonance dynamics of low as well as high inclination asteroids in the solar system.


title: Observation of the Eckhaus Instability in Whispering-Gallery Mode Resonators

abstract: The Eckhaus instability is a secondary instability of nonlinear spatiotemporal patterns in which high-wavenumber periodic solutions become unstable against small-wavenumber perturbations. We show in this letter that this instability can take place in Kerr combs generated with ultra-high $Q$ whispering-gallery mode resonators. In our experiment, sub-critical Turing patterns (rolls) undergo Eckhaus instabilities upon changes in the laser detuning leading to cracking patterns with long-lived transients. In the spectral domain, this results in a metastable Kerr comb dynamics with a timescale that can be larger than one minute. This ultra-slow timescale is at least seven orders of magnitude larger than the intracavity photon lifetime, and is in sharp contrast with all the transient behaviors reported so far in cavity nonlinear optics, that are typically only few photon lifetimes long (i.e., in the ps range). We show that this phenomenology is well explained by the Lugiato-Lefever model, as the result of an Eckhaus instability. Our theoretical analysis is found to be in excellent agreement with the experimental measurements.


title: Exo-MerCat: a merged exoplanet catalog with Virtual Observatory connection

abstract: The heterogeneity of papers dealing with the discovery and characterization of exoplanets makes every attempt to maintain a uniform exoplanet catalog almost impossible. Four sources currently available online (NASA Exoplanet Archive, Exoplanet Orbit Database, Exoplanet Encyclopaedia, and Open Exoplanet Catalogue) are commonly used by the community, but they can hardly be compared, due to discrepancies in notations and selection criteria. Exo-MerCat is a Python code that collects and selects the most precise measurement for all interesting planetary and orbital parameters contained in the four databases, accounting for the presence of multiple aliases for the same target. It can download information about the host star as well by the use of Virtual Observatory ConeSearch connections to the major archives such as SIMBAD and those available in VizieR. A Graphical User Interface is provided to filter data based on the user's constraints and generate automatic plots that are commonly used in the exoplanetary community. With Exo-MerCat, we retrieved a unique catalog that merges information from the four main databases, standardizing the output and handling notation differences issues. Exo-MerCat can correct as many issues that prevent a direct correspondence between multiple items in the four databases as possible, with the available data. The catalog is available as a VO resource for everyone to use and it is periodically updated, according to the update rates of the source catalogs.


title: Class numbers and representations by ternary quadratic forms with congruence conditions

abstract: In this paper, we are interested in the interplay between integral ternary quadratic forms and class numbers. This is partially motivated by a question of Petersson.


title: Conformal quantum mechanics of causal diamonds

abstract: It is shown that a general radial conformal Killing vector in Minkowski space-time can be associated to a generator of time evolution in conformal quantum mechanics. Among these conformal Killing vectors one finds a class which maps causal diamonds in Minkowski space-time into themselves. The flow of such Killing vectors describes worldlines of accelerated observers with a finite lifetime within the causal diamond. Time evolution of static diamond observers is equivalent to time evolution in conformal quantum mechanics governed by a hyperbolic Hamiltonian and covering only a segment of the time axis. This indicates that the Unruh temperature perceived by static diamond observers in the vacuum state of inertial observers in Minkowski space can be obtained from the behaviour of the two-point functions of conformal quantum mechanics.


title: Higher order cumulants of net baryon-number distributions at non-zero $\mu_B$

abstract: Using recent results on higher order cumulants of conserved charge fluctuations from lattice QCD, we construct mean, variance, skewness, kurtosis, hyper-skewness and hyper-kurtosis of net-baryon number distributions for small baryon chemical potentials $\mu_B$. For the strangeness neutral case ($\mu_S=0$) at fixed ratio of electric charge to baryon number density ($\frac{n_Q}{n_B}=0.4$), which is appropriate for a comparison with heavy ion collisions, we present results for $\kappa_B \sigma_B^2$, $S_B \sigma_B^3/M_B$, $\kappa^{H}_{B}\sigma_{B}^4$ and $S^{H}_{B}\sigma^5_{B}/M_{B}$ on the crossover line for the chiral transition, $T_{pc}(\mu_B)$. Continuum extrapolations for this pseudo-critical transition line have recently been reported by HotQCD up to baryon chemical potentials $\mu_B\simeq 300$ MeV [arXiv:1812.08235]. These cumulant ratios are of direct relevance for comparisons with corresponding ratios measured by STAR in the BES-I and II runs at beam energies $\sqrt{s_{NN}}\ge 20$ GeV. In particular, we point out that recent high statistics results on skewness and kurtosis of net-baryon number distributions obtained by STAR at $\sqrt{s_{NN}} = 54.4$ GeV put strong constraints on freeze-out parameters and are consistent with predictions from thermal QCD.


title: Non-interacting many-particle quantum transport between finite reservoirs

abstract: We present a formalism to study many-particle quantum transport across a lattice locally connected to two finite, non-stationary (bosonic or fermionic) reservoirs, both of which are in a thermal state. We show that, for conserved total particle number, a system of nonlinear quantum-classical master equations describes the concurrent many-particle time evolution on the lattice and in the reservoirs. The finiteness of the reservoirs makes a macroscopic current emerge, which decreases exponentially in time, and asymptotically drives the many-particle configuration into an equilibrium state where the particle flow ceases. We analytically derive the time scale of this equilibration process, and, furthermore, investigate the imprint of many-particle interferences on the transport process.


title: How accurately can the Extended Thomas-Fermi method describe the inner crust of a neutron star?

abstract: We perform a systematic comparison between the results obtained by solving fully self-consistently the Hartree-Fock-Bogoliubov equations, and those obtained using the semi-classical Extended Thomas-Fermi method, for various Wigner-Seitz cells within the inner crust of a neutron star. The lack of pairing correlations in the semi-classical approach leads to a large discrepancy between the two approaches. This discrepancy is well beyond the error of the quantum-mechanical calculation, and is related to spurious shell effects in the neutron gas.


title: What prevents Finnish women from applying to software engineering roles? A preliminary analysis of survey data

abstract: Finland is considered a country with a good track record in gender equality. Whilst statistics support the notion that Finland is performing well compared to many other countries in terms of workplace equality, there are still many areas for improvement. This paper focuses on the problems that some women face in obtaining software engineering roles. We report a preliminary analysis of survey data from 252 respondents. These are mainly women who have shown an interest in gaining programming roles by joining the Mimmit koodaa initiative, which aims to increase equality and diversity within the software industry. The survey sought to understand what early experiences may influence later career choices and feelings of efficacy and confidence needed to pursue technology-related careers. These initial findings reveal that women's feelings of computing self-efficacy and attitudes towards software engineering are shaped by early experiences. More negative experiences decrease the likelihood of working in software engineering roles in the future, despite expressing an interest in the field.


title: The Casimir force, causality and the Gurzhi model

abstract: An extended Drude model, termed as the Gurzhi model, which takes into account the electron-phonon and electron-electron interactions, is applied to calculate the Casimir force between two metallic plates. It is shown that although the dielectric permittivity of the Gurzhi model has a first order pole in the upper half-plane of complex frequencies and, thus, violates the causality principle, it can be used in a restricted frequency interval in combination with the experimental permittivity determined by the optical data for the complex index of refraction. The imaginary part of the Gurzhi dielectric permittivity of Au at low frequencies demonstrates better agreement with the permittivity given by the optical data than the simple Drude model. The Casimir pressure between two Au plates is computed using the Gurzhi, Drude and plasma model approaches, taking into account the optical data, as well as with the simple Drude and plasma models. The contribution of the electron-electron scattering to the Casimir pressure is estimated. Although a comparison with the measurement data of two precise experiments show that the Gurzhi model does not resolve the Casimir puzzle, the obtained results suggest further clarification of this fundamental problem.


title: Patterns for Name Analysis and Type Analysis with JastAdd

abstract: In the last two decades, tools have been implemented to more formally specify the semantic analysis phase of a compiler instead of relying on handwritten code. In this paper, we introduce patterns and a method to translate a formal definition of a language's type system into a specification for JastAdd, which is one of the aforementioned tools based on Reference Attribute Grammars. This methodological approach will help language designers and compiler engineers to more systematically use such tools for semantic analysis. As an example, we use a simple, yet complete imperative language and provide an outlook on how the method can be extended to cover further language constructs or even type inference.


title: Constructing Multipartite Bell inequalities from stabilizers

abstract: Bell inequality with self-testing property has played an important role in quantum information field with both fundamental and practical applications. However, it is generally challenging to find Bell inequalities with self-testing property for multipartite states and actually there are not many known candidates. In this work, we propose a systematical framework to construct Bell inequalities from stabilizers which are maximally violated by general stabilizer states, with two observables for each local party. We show that the constructed Bell inequalities can self-test any stabilizer state which is essentially device-independent, if and only if these stabilizers can uniquely determine the state in a device-dependent manner. This bridges the gap between device-independent and device-dependent verification methods. Our framework can provide plenty of Bell inequalities for self-testing stabilizer states. Among them, we give two families of Bell inequalities with different advantages: (1) a family of Bell inequalities with a constant ratio of quantum and classical bounds using 2N correlations, (2) Single pair inequalities improving on all previous robustness self-testing bounds using N+1 correlations, which are both efficient and suitable for realizations in multipartite systems. Our framework can not only inspire more fruitful multipartite Bell inequalities from conventional verification methods, but also pave the way for their practical applications.


title: Hindered surface diffusion of bonded molecular clusters mediated by surface defects

abstract: The design of low dimensional materials through surface assisted self-assembly requires a better understanding of the factors that limit and control surface diffusion. We reveal how substrate surface defects hinder the mobility of sub-monolayer organic adsorbates on a metal surface with the model CuPc/Cu(111) system. Post-deposition annealing bonds CuPc molecules into dendritelike clusters that are often mobile at room temperature. Surface defects on Cu(111) create energetic barriers that prevent CuPc cluster motion on the metal surface. This phenomenon was unveiled by the motion of small clusters that show rigid-body diffusion solely in the available space in between defects. When clusters are sufficiently surrounded by defects, they become completely pinned in place and become immobilized.


title: Quantum transport between finite reservoirs

abstract: When driven by a potential bias between two finite reservoirs, the particle current across a quantum system evolves from an initial loading through a coherent, followed by a metastable phase, and ultimately fades away upon equilibration. We formulate a theory which fully accounts for the associated, distinct time scales, and identifies the parameter dependence of the decay rate which ultimately controls the convergence towards equilibrium. Our formalism guarantees total particle number conservation and fundamental consistency between macroscopic and internal currents flowing in the system. We furthermore establish a clear imprint of the fermionic or bosonic particle character on the resulting conductance.


title: Geosocial Location Classification: Associating Type to Places Based on Geotagged Social-Media Posts

abstract: Associating type to locations can be used to enrich maps and can serve a plethora of geospatial applications. An automatic method to do so could make the process less expensive in terms of human labor, and faster to react to changes. In this paper we study the problem of Geosocial Location Classification, where the type of a site, e.g., a building, is discovered based on social-media posts. Our goal is to correctly associate a set of messages posted in a small radius around a given location with the corresponding location type, e.g., school, church, restaurant or museum. We explore two approaches to the problem: (a) a pipeline approach, where each message is first classified, and then the location associated with the message set is inferred from the individual message labels; and (b) a joint approach where the individual messages are simultaneously processed to yield the desired location type. We tested the two approaches over a dataset of geotagged tweets. Our results demonstrate the superiority of the joint approach. Moreover, we show that due to the unique structure of the problem, where weakly-related messages are jointly processed to yield a single final label, linear classifiers outperform deep neural network alternatives.


title: Zendoo: a zk-SNARK Verifiable Cross-Chain Transfer Protocol Enabling Decoupled and Decentralized Sidechains

abstract: Sidechains are an appealing innovation devised to enable blockchain scalability and extensibility. The basic idea is simple yet powerful: construct a parallel chain -- sidechain -- with desired features, and provide a way to transfer coins between the mainchain and the sidechain. In this paper, we introduce Zendoo, a construction for Bitcoin-like blockchain systems that allows the creation and communication with sidechains of different types without knowing their internal structure. We consider a parent-child relationship between the mainchain and sidechains, where sidechain nodes directly observe the mainchain while mainchain nodes only observe cryptographically authenticated certificates from sidechain maintainers. We use zk-SNARKs to construct a universal verifiable transfer mechanism that is used by sidechains. Moreover, we propose a specific sidechain construction, named Latus, that can be built on top of this infrastructure, and realizes a decentralized verifiable blockchain system for payments. We leverage the use of recursive composition of zk-SNARKs to generate succinct proofs of sidechain state progression that are used to generate certificates' validity proofs. This allows the mainchain to efficiently verify all operations performed in the sidechain without knowing any details about those operations.


title: A^1-Euler classes: six functors formalisms, dualities, integrality and linear subspaces of complete intersections

abstract: We equate various Euler classes of algebraic vector bundles, including those of [BM, KW, DJK], and one suggested by M.J. Hopkins, A. Raksit, and J.-P. Serre. We establish integrality results for this Euler class, and give formulas for local indices at isolated zeros, both in terms of 6-functor formalism of coherent sheaves and as an explicit recipe in commutative algebra of Scheja and Storch. As an application, we compute the Euler classes associated to arithmetic counts of d-planes on complete intersections in P^n in terms of topological Euler numbers over R and C.


title: Rank $2r$ iterative least squares: efficient recovery of ill-conditioned low rank matrices from few entries

abstract: We present a new, simple and computationally efficient iterative method for low rank matrix completion. Our method is inspired by the class of factorization-type iterative algorithms, but substantially differs from them in the way the problem is cast. Precisely, given a target rank $r$, instead of optimizing on the manifold of rank $r$ matrices, we allow our interim estimated matrix to have a specific over-parametrized rank $2r$ structure. Our algorithm, denoted R2RILS for rank $2r$ iterative least squares, has low memory requirements, and at each iteration it solves a computationally cheap sparse least-squares problem. We motivate our algorithm by its theoretical analysis for the simplified case of a rank-1 matrix. Empirically, R2RILS is able to recover ill conditioned low rank matrices from very few observations -- near the information limit, and it is stable to additive noise.


title: Role of IR-Improvement in Precision LHC/FCC Physics and in Quantum Gravity

abstract: IR-improvement based on amplitude-level resummation allows one to control unintegrable results in quantum field theory with arbitrary precision in principle. We illustrate such resummation in specific examples in precision LHC and FCC physics and in quantum gravity.


title: A generalization of the Freidlin-Wentcell theorem on averaging of Hamiltonian systems

abstract: In this paper, we generalize the classical Freidlin-Wentzell's theorem for random perturbations of Hamiltonian systems. In stead of the two-dimensional standard Brownian motion, the coefficient for the noise term is no longer the identity matrix but a state-dependent matrix plus a state-dependent matrix that converges uniformly to 0 on any compact sets as $\epsilon$ tends to 0. We also take the drift term into consideration where the drfit term also contains two parts, the state-dependent mapping and a state-dependent mapping that converges uniformly to 0 on any compact sets as $\epsilon$ tends to 0. In the proof, we use the result of generalized differential operator. We also adapt a new way to prove the weak convergence inside the edge by constructing an auxiliary process and apply Girsanov's theorem in the proof of gluing condition.


title: TPPO: A Novel Trajectory Predictor with Pseudo Oracle

abstract: Forecasting pedestrian trajectories in dynamic scenes remains a critical problem in various applications, such as autonomous driving and socially aware robots. Such forecasting is challenging due to human-human and human-object interactions and future uncertainties caused by human randomness. Generative model-based methods handle future uncertainties by sampling a latent variable. However, few studies explored the generation of the latent variable. In this work, we propose the Trajectory Predictor with Pseudo Oracle (TPPO), which is a generative model-based trajectory predictor. The first pseudo oracle is pedestrians' moving directions, and the second one is the latent variable estimated from ground truth trajectories. A social attention module is used to aggregate neighbors' interactions based on the correlation between pedestrians' moving directions and future trajectories. This correlation is inspired by the fact that pedestrians' future trajectories are often influenced by pedestrians in front. A latent variable predictor is proposed to estimate latent variable distributions from observed and ground-truth trajectories. Moreover, the gap between these two distributions is minimized during training. Therefore, the latent variable predictor can estimate the latent variable from observed trajectories to approximate that estimated from ground-truth trajectories. We compare the performance of TPPO with related methods on several public datasets. Results demonstrate that TPPO outperforms state-of-the-art methods with low average and final displacement errors. The ablation study shows that the prediction performance will not dramatically decrease as sampling times decline during tests.


title: Weight distributions of several families of 3-weight binary linear codes

abstract: The linear codes with a few weights have been applied widely in combinatorial designs, secret sharing, association schemes, authentication codes and strongly regular graphs. In this paper, we first correct an erroneous result about the exponential sum $\sum_{x\in \mathbb{F}_{2^{e}}}\chi_1\left(ax^{2^{\alpha}+1}+bx\right).$ Then, using the above exponential sum, we construct several families of binary linear codes of $3$-weight and determine their weight distributions. Moreover, Most of them can be used in secret sharing schemes.


title: Interpretable & Time-Budget-Constrained Contextualization for Re-Ranking

abstract: Search engines operate under a strict time constraint as a fast response is paramount to user satisfaction. Thus, neural re-ranking models have a limited time-budget to re-rank documents. Given the same amount of time, a faster re-ranking model can incorporate more documents than a less efficient one, leading to a higher effectiveness. To utilize this property, we propose TK (Transformer-Kernel): a neural re-ranking model for ad-hoc search using an efficient contextualization mechanism. TK employs a very small number of Transformer layers (up to three) to contextualize query and document word embeddings. To score individual term interactions, we use a document-length enhanced kernel-pooling, which enables users to gain insight into the model. TK offers an optimal ratio between effectiveness and efficiency: under realistic time constraints (max. 200 ms per query) TK achieves the highest effectiveness in comparison to BERT and other re-ranking models. We demonstrate this on three large-scale ranking collections: MSMARCO-Passage, MSMARCO-Document, and TREC CAR. In addition, to gain insight into TK, we perform a clustered query analysis of TK's results, highlighting its strengths and weaknesses on queries with different types of information need and we show how to interpret the cause of ranking differences of two documents by comparing their internal scores.


title: Black hole production at lepton colliders

abstract: Production of black holes has been discussed in a variety of extensions of the Standard Model, and related bounds have been established from data taken at the Large Hadron Collider. We show that, if the Higgs particle has a fully gravitational content via the equivalence principle, enhanced cross-sections of black holes at colliders should be expected within the Standard Model itself. The case of black hole production by precision measurements at electron colliders is discussed. The Coulomb repulsion strongly suppresses the related cross-section with respect to the one based on the hoop conjecture, making the possible production of black holes still unfeasible with current beam technology. At the same time, this suggests the reanalysis of the bounds, based on the hoop conjecture, already determined in hadronic collisions for extra-dimensional models.


title: Bloch's formula for 0-cycles with modulus and higher dimensional Class Field Theory

abstract: We prove Bloch's formula for the Chow group of 0-cycles with modulus on a smooth quasi-projective surface over a field. We use this formula to give a simple proof of the rank one case of a conjecture of Deligne and Drinfeld on lisse $\overline{\mathbb{Q}}_{\ell}$-sheaves. This was originally solved by Kerz and Saito in characteristic $\neq 2$.


title: Does the Eigenstate Thermalization Hypothesis Imply Thermalization?

abstract: Eigenstate thermalization hypothesis (ETH) is discussed. We show that one common formulation of ETH does not necessarily imply thermalization of an observable of isolated many body quantum system. To get thermalization one has to postulate the canonical or microcanonical distribution in the ETH-ansatz. More generally, any other average can be postulated in the generalized ETH-ansatz which leads to a corresponding equilibration condition.


title: Observational and theoretical constraints on the formation and early evolution of the first dust grains in galaxies at 5 < z < 10

abstract: The first generation of stars were born a few hundred million years after the big bang. These stars synthesized elements heavier than H and He, that are later expelled into the interstellar medium, initiating the rise of metals. Within this enriched medium, the first dust grains formed. This event is cosmological crucial for molecule formation as dust plays a major role by cooling low-metallicity star-forming clouds which can fragment to create lower mass stars. Collecting information on these first dust grains is difficult because of the negative alliance of large distances and low dust masses. We combine the observational information from galaxies at redshifts 5 < z < 10 to constrain their dust emission and theoretically understand the first evolutionary phases of the dust cycle. Spectral energy distributions (SEDs) are fitted with CIGALE and the physical parameters and their evolution are modelled. From this SED fitting, we build a dust emission template for this population of galaxies in the epoch of reionization. Our new models explain why some early galaxies are observed and others are not. We follow in time the formation of the first grains by supernovae later destroyed by other supernova blasts and expelled in the circumgalactic and intergalactic media. We have found evidence for the first dust grains formed in the universe. But, above all, this letter underlines the need to collect more data and to develop new facilities to further constrain the dust cycle in galaxies in the epoch of reionization.


title: RGISTools: Downloading, Customizing, and Processing Time Series of Remote Sensing Data in R

abstract: There is a large number of data archives and web services offering free access to multispectral satellite imagery. Images from multiple sources are increasingly combined to improve the spatio-temporal coverage of measurements while achieving more accurate results. Archives and web services differ in their protocols, formats, and data standards, which are barriers to combine datasets. Here, we present RGISTools, an R package to create time-series of multispectral satellite images from multiple platforms in a harmonized and standardized way. We first provide an overview of the package functionalities, namely downloading, customizing, and processing multispectral satellite imagery for a region and time period of interest as well as a recent statistical method for gap-filling and smoothing series of images, called interpolation of the mean anomalies. We further show the capabilities of the package through a case study that combines Landsat-8 and Sentinel-2 satellite optical imagery to estimate the level of a water reservoir in Northern Spain. We expect RGISTools to foster research on data fusion and spatio-temporal modelling using satellite images from multiple programs.


title: Impact of Convective Boundary Mixing on the TP-AGB

abstract: The treatment of convective boundaries remains an important source of uncertainty within stellar evolution, with drastic implications for the thermally-pulsing stars on the Asymptotic Giant Branch (AGB). Various sources are taken as motivation for the incorporation of convective boundary mixing during this phase, from s-process nucleosynthesis to hydrodynamical models. In spite of the considerable evidence in favour of the existence of convective boundary mixing on the pre-AGB evolution, this mixing is not universally included in models of TP-AGB stars. The aim of this investigation is to ascertain the extent of convective boundary mixing, which is compatible with observations when considering full evolutionary models. Additionally, we investigate a theoretical argument that has been made that momentum-driven overshooting at the base of the pulse-driven convection zone should be negligible. We show that, while the argument holds,estimations based on the picture of turbulent entrainment suggest that mixing is possible at both convective boundaries. We demonstrate that additional mixing at convective boundaries during core-burning phases prior to the thermally-pulsing Asymptotic Giant Branch has an impact on the later evolution, changing the mass range at which the third dredge-up and hot-bottom burning occur, and thus also the final surface composition. In addition, an effort has been made to constrain the efficiency of convective boundary mixing at the different convective boundaries, using observational constraints. Our study suggests a strong tension between different constraints that makes it impossible to reproduce all observables simultaneously within the framework of an exponentially decaying overshooting. This result calls for a reassessment of both the models of convective boundary mixing and the observational constraints.


title: Rapid Adaptation of BERT for Information Extraction on Domain-Specific Business Documents

abstract: Techniques for automatically extracting important content elements from business documents such as contracts, statements, and filings have the potential to make business operations more efficient. This problem can be formulated as a sequence labeling task, and we demonstrate the adaption of BERT to two types of business documents: regulatory filings and property lease agreements. There are aspects of this problem that make it easier than "standard" information extraction tasks and other aspects that make it more difficult, but on balance we find that modest amounts of annotated data (less than 100 documents) are sufficient to achieve reasonable accuracy. We integrate our models into an end-to-end cloud platform that provides both an easy-to-use annotation interface as well as an inference interface that allows users to upload documents and inspect model outputs.


title: If I Hear You Correctly: Building and Evaluating Interview Chatbots with Active Listening Skills

abstract: Interview chatbots engage users in a text-based conversation to draw out their views and opinions. It is, however, challenging to build effective interview chatbots that can handle user free-text responses to open-ended questions and deliver engaging user experience. As the first step, we are investigating the feasibility and effectiveness of using publicly available, practical AI technologies to build effective interview chatbots. To demonstrate feasibility, we built a prototype scoped to enable interview chatbots with a subset of active listening skills - the abilities to comprehend a user's input and respond properly. To evaluate the effectiveness of our prototype, we compared the performance of interview chatbots with or without active listening skills on four common interview topics in a live evaluation with 206 users. Our work presents practical design implications for building effective interview chatbots, hybrid chatbot platforms, and empathetic chatbots beyond interview tasks.


title: Z production in pPb collisions and charmonium production in PbPb ultra-peripheral collisions at LHCb

abstract: This article presents recent results of the Z boson production in the proton-lead collisions and the charmonium production in the ultra-peripheral lead-lead collisions at the LHC collected during 2013 to 2018 by the LHCb detector. The potential heavy ion programs in these topics are also discussed.


title: Five Texture Zeros for Dirac Neutrino Mass Matrices

abstract: In this work we propose new five textures zeros for the mass matrices in the lepton sector in order to predict values for the neutrino masses. In our approach we go beyond the Standard Model by assuming Dirac masses for the neutrinos, feature which allows us to make a theoretical prediction for the lightest neutrino mass in the normal ordering. The textures analyzed have enough free parameters to adjust the $V_{\text{\tiny PMNS}}$ mixing matrix including the CP-violating phase, the neutrino mass squared differences $\delta m_{21}^2,\,\delta m_{31}^2$, and the three charged lepton masses. In order to have reliable results two different approaches are used: the first one is based on a least-squares analysis to fit the lepton masses and the mixing parameters to their corresponding experimental values, for this case the best fit for the lightest neutrino mass is $(3.9\pm^{0.6}_{0.8})\times 10^{-3}$ eV.; the second approach is just algebraic, based on the weak basis transformation method, in this case the lightest neutrino mass consistent with the experimental values and the restrictions coming from the five texture zeros of the mass matrices is equal to $(3.5\pm0.9)\times 10^{-3}$ eV.


title: Completing Simple Valuations in K-categories

abstract: We prove that Keimel and Lawson's K-completion Kc of the simple valuation monad Vs defines a monad Kc o Vs on each K-category A. We also characterize the Eilenberg-Moore algebras of Kc o Vs as the weakly locally convex K-cones, and its algebra morphisms as the continuous linear maps. In addition, we explicitly describe the distributive law of Vs over Kc, which allows us to show that the K-completion of any locally convex (respectively, weakly locally convex, locally linear) topological cone is a locally convex (respectively, weakly locally convex, locally linear) K-cone. We also give an example - the Cantor tree with a top - that shows the dcpo-completion of the simple valuations is not the D-completion of the simple valuations in general, where D is the category of monotone convergence spaces and continuous maps.


title: Osmotic stress and pore nucleation in charged biological nanoshells and capsids

abstract: A model system is proposed to investigate the chemical equilibrium and mechanical stability of biological spherical-like nanoshells in contact with an aqueous solution with added dissociated electrolyte of a given concentration. The ionic chemical equilibrium across the permeable shell is investigated in the framework of an accurate Density Functional Theory (DFT) that incorporates electrostatic and hardcore correlations beyond the traditional mean-field ({\it e. g.}, Poisson-Boltzmann) limit. The accuracy of the theory is tested by a direct comparison with Monte Carlo (MC) simulations. A simple analytical expression is then deduced which clearly highlights the entropic, electrostatic, and self-energy contributions to the osmotic stress over the shell in terms of the calculated ionic profiles. By invoking a continuum mean-field elastic approach to account for the shell surface stress upon osmotic stretching, the mechanical equilibrium properties of the shell under a wide variety of ionic strengths and surface charges are investigated. The model is further coupled to a continuum mechanical approach similar in structure to a Classical Nucleation Theory (CNT) to address the question of mechanical stability of the shells against a pore nucleation. This allows us to construct a phase diagram which delimits the mechanical stability of capsids for different ionic strengths and shell surface charges.


title: On existence of some special pair of primitive elements over finite fields

abstract: In this paper we generalize the results of Sharma, Awasthi and Gupta (see \cite{SAG}). We work over a field of any characteristic with $q = p^k$ elements and we give a sufficient condition for the existence of a primitive element $\alpha \in \mathbb{F}_{p^k}$ such that $f(\alpha)$ is also primitive in $\mathbb{F}_{p^k}$, where $f(x) \in \mathbb{F}_{p^k}(x)$ is a quotient of polynomials with some restrictions. We explicitly determine the values of $k$ for which such a pair exists for $p=2,3,5$ and $7$.


title: Diffusive interaction of multiple surface nanobubbles and nanodroplets: shrinkage, growth, and coarsening

abstract: Surface nanobubbles are nanoscopic spherical-cap shaped gaseous domains on immersed substrates which are stable, even for days. After the stability of a {\it single} surface nanobubble has been theoretically explained, i.e. contact line pinning and gas oversaturation are required to stabilize them against diffusive dissolution [Lohse and Zhang, Phys.\ Rev.\ E 91, 031003 (R) (2015)], here we focus on the {\it collective} diffusive interaction of {\it multiple} nanobubbles. For that purpose we develop a finite difference scheme for the diffusion equation with the appropriate boundary conditions and with the immersed boundary method used to represent the growing or shrinking bubbles. After validation of the scheme against the exact results of Epstein and Plesset for a bulk bubble [J. Chem. Phys. 18, 1505 (1950)] and of Lohse and Zhang for a surface bubble, the framework of these simulations is used to describe the coarsening process of competitively growing nanobubbles. The coarsening process for such diffusively interacting nanobubbles slows down with advancing time and thus increasing bubble distance. The present results for surface nanobubbles are also applicable for immersed surface nanodroplets, for which better controlled experimental results of the coarsening process exist.


title: Prediction of head motion from speech waveforms with a canonical-correlation-constrained autoencoder

abstract: This study investigates the direct use of speech waveforms to predict head motion for speech-driven head-motion synthesis, whereas the use of spectral features such as MFCC as basic input features together with additional features such as energy and F0 is common in the literature. We show that, rather than combining different features that originate from waveforms, it is more effective to use waveforms directly predicting corresponding head motion. The challenge with the waveform-based approach is that waveforms contain a large amount of information irrelevant to predict head motion, which hinders the training of neural networks. To overcome the problem, we propose a canonical-correlation-constrained autoencoder (CCCAE), where hidden layers are trained to not only minimise the error but also maximise the canonical correlation with head motion. Compared with an MFCC-based system, the proposed system shows comparable performance in objective evaluation, and better performance in subject evaluation.


title: Non-locality and time-dependent boundary conditions: a Klein-Gordon perspective

abstract: The dynamics of a particle in an expanding cavity is investigated in the Klein-Gordon framework in a regime in which the single particle picture remains valid. The cavity expansion represents a time-dependent boundary condition for the relativistic wavefunction. We show that this expansion induces a non-local effect on the current density throughout the cavity. Our results indicate that a relativistic treatment still contains apparently spurious effects traditionally associated with the unbounded velocities inherent to non-relativistic solutions obtained from the Schroedinger equation. Possible reasons for this behaviour are discussed.


title: Iterative algorithm with structured diagonal Hessian approximation for solving nonlinear least squares problems

abstract: Nonlinear least-squares problems are a special class of unconstrained optimization problems in which their gradient and Hessian have special structures. In this paper, we exploit these structures and proposed a matrix-free algorithm with a diagonal Hessian approximation for solving nonlinear least-squares problems. We devise appropriate safeguarding strategies to ensure the Hessian matrix is positive definite throughout the iteration process. The proposed algorithm generates descent direction and is globally convergent. Preliminary numerical experiments show that the proposed method is competitive with a recently developed similar method.


title: CBR: Controlled Burst Recording

abstract: Collecting traces from software running in the field is both useful and challenging. Traces may indeed help revealing unexpected usage scenarios, detecting and reproducing failures, and building behavioral models that reflect how the software is actually used. On the other hand, recording traces is an intrusive activity that may annoy users, negatively affecting the usability of the applications, if not properly designed. In this paper we address field monitoring by introducing Controlled Burst Recording, a monitoring solution that can collect comprehensive runtime data without compromising the quality of the user experience. The technique encodes the knowledge extracted from the monitored application as a finite state model that both represents the sequences of operations that can be executed by the users and the corresponding internal computations that might be activated by each operation. Our initial assessment with information extracted from ArgoUML shows that Controlled Burst Recording can reconstruct behavioral information more effectively than competing sampling techniques, with a low impact on the system response time.


title: $\epsilon$-shotgun: $\epsilon$-greedy Batch Bayesian Optimisation

abstract: Bayesian optimisation is a popular, surrogate model-based approach for optimising expensive black-box functions. Given a surrogate model, the next location to expensively evaluate is chosen via maximisation of a cheap-to-query acquisition function. We present an $\epsilon$-greedy procedure for Bayesian optimisation in batch settings in which the black-box function can be evaluated multiple times in parallel. Our $\epsilon$-shotgun algorithm leverages the model's prediction, uncertainty, and the approximated rate of change of the landscape to determine the spread of batch solutions to be distributed around a putative location. The initial target location is selected either in an exploitative fashion on the mean prediction, or -- with probability $\epsilon$ -- from elsewhere in the design space. This results in locations that are more densely sampled in regions where the function is changing rapidly and in locations predicted to be good (i.e close to predicted optima), with more scattered samples in regions where the function is flatter and/or of poorer quality. We empirically evaluate the $\epsilon$-shotgun methods on a range of synthetic functions and two real-world problems, finding that they perform at least as well as state-of-the-art batch methods and in many cases exceed their performance.


title: Adhesive and non-adhesive contact of a rigid indenter and a thin elastic layer with surface tension

abstract: We consider an adhesive contact between a thin soft layer on a rigid substrate and a rigid cylindrical indenter ("line contact") with account of the surface tension of the layer. First, it is shown that the boundary condition for the surface outside the contact area is given by the constant contact angle - as in the case of fluids in contact with solid surfaces. In the approximation of thin layer and under usual assumptions of small indentation and small inclination angles of the surface, the problem is solved analytically.


title: Pseudo-differential extension for graded nilpotent Lie groups

abstract: Classical pseudo-differential operators of order zero on a graded nilpotent Lie group $G$ form a $^*$-subalgebra of the bounded operators on $L^2(G)$. We show that its $C^*$-closure is an extension of a noncommutative algebra of principal symbols by compact operators. As a new approach, we use the generalized fixed point algebra of an $\mathbb{R}_{>0}$-action on a certain ideal in the $C^*$-algebra of the tangent groupoid of $G$. The action takes the graded structure of $G$ into account. Our construction allows to compute the $K$-theory of the algebra of symbols.


title: Escape from the Quantum Pigeon Conundrum

abstract: It has recently been argued in Aharonov et. al. (2016) that quantum mechanics violates the Pigeon Counting Principle (PCP) which states that if one distributes three pigeons among two boxes there must be at least two pigeons in one of the boxes. However, this conclusion cannot justified by rigorous theoretical arguments. The issue is further complicated by experimental confirmation of the transition amplitudes predicted in this paper that nevertheless do not support the conclusion of PCP violation. Here we prove via a set of operator identities that the PCP is not violated within quantum mechanics, regardless of interpretation.


title: $W^+W^-$ production at the LHC: NLO QCD corrections to the loop-induced gluon fusion channel

abstract: We compute the NLO QCD corrections to the loop-induced gluon fusion contribution in $W^+W^-$ production at the LHC. We consider the full leptonic process $pp\to \ell^+\ell^{\prime\, -}\nu_{\ell}{\bar\nu}_{\ell^\prime}+X$, by including resonant and non-resonant diagrams, spin correlations and off-shell effects. Quark-gluon partonic channels are included for the first time in the calculation, and our results are combined with NNLO predictions to the quark annihilation channel at the fully differential level. The computed corrections, which are formally of ${\cal O}(\alpha_{\mathrm{S}}^3)$, increase the NNLO cross section by only about 2%, but have an impact on the shapes of kinematical distributions, in part due to the jet veto, which is usually applied to reduce the top-quark background. Our results, supplemented with NLO EW effects, provide the most advanced fixed-order predictions available to date for this process, and are compared with differential ATLAS data at $\sqrt{s}=$ 13 TeV.


title: Wasserstein Exponential Kernels

abstract: In the context of kernel methods, the similarity between data points is encoded by the kernel function which is often defined thanks to the Euclidean distance, a common example being the squared exponential kernel. Recently, other distances relying on optimal transport theory - such as the Wasserstein distance between probability distributions - have shown their practical relevance for different machine learning techniques. In this paper, we study the use of exponential kernels defined thanks to the regularized Wasserstein distance and discuss their positive definiteness. More specifically, we define Wasserstein feature maps and illustrate their interest for supervised learning problems involving shapes and images. Empirically, Wasserstein squared exponential kernels are shown to yield smaller classification errors on small training sets of shapes, compared to analogous classifiers using Euclidean distances.


title: Multivariate normal approximation for traces of random unitary matrices

abstract: In this article, we obtain a super-exponential rate of convergence in total variation between the traces of the first $m$ powers of an $n\times n$ random unitary matrices and a $2m$-dimensional Gaussian random variable. This generalizes previous results in the scalar case to the multivariate setting, and we also give the precise dependence on the dimensions $m$ and $n$ in the estimate with explicit constants. We are especially interested in the regime where $m$ grows with $n$ and our main result basically states that if $m\ll \sqrt{n}$, then the rate of convergence in the Gaussian approximation is $\Gamma(\frac nm+1)^{-1}$ times a correction. We also show that the Gaussian approximation remains valid for all $m\ll n^{2/3}$ without a fast rate of convergence.


title: Shape induced segregation and anomalous particle transport under spherical confinement

abstract: Colloid or nanoparticle mobility under confinement is of central importance to a wide range of physical and biological processes. Here, we introduce a minimal model of particles in a hydrodynamic continuum to examine how particle shape and concentration affect the transport of particles in spherical confinement. Specifically, an immersed boundary-General geometry Ewald-like approach is adopted to simulate the dynamics of spheres and cylinders under the influence of short-and long-range fluctuating hydrodynamic interactions with appropriate non-slip conditions at the confining walls. An efficient $\it{O(N)}$ parallel finite element algorithm is used, thereby allowing simulations at high concentrations, while a Chebyshev polynomial approximation is implemented in order to satisfy the fluctuation-dissipation theorem. A concentration-dependent anomalous diffusion is observed for suspended particles. It is found that introducing cylinders in a background of spheres, i.e. particles with a simple degree of anisotropy, has a pronounced influence on the structure and dynamics of the particles. First, increasing the fraction of cylinders induces a particle segregation effect, where spheres are pushed towards the wall and cylinders remain near the center of the cavity. This segregation leads to lower mobility for the spheres relative to that encountered in a system of pure spheres at the same volume fraction. Second, the diffusive-to-anomalous transition and the degree of anomaly--quantified by the power-law exponent in the mean square displacement vs. time relation-both increase as the fraction of cylinders becomes larger. These findings are of relevance for studies of diffusion in the cytoplasm, where proteins exhibit a distribution of size and shapes that could lead to some of the effects identified in the simulations reported here.


title: Fingerprints of the cosmological constant: Folds in the profiles of the axionic dark matter distribution in a dyon exterior

abstract: We consider the magnetic monopole in the axionic dark matter environment (axionic dyon) in the framework of the Reissner - Nordstr\"om - de Sitter model. Our aim is to study the distribution of the pseudoscalar (axion) and electric fields near the so-called folds, which are characterized by the profiles with the central minimum, the barrier on the left, and the maximum on the right of this minimum. The electric field in the fold-like zones is shown to change the sign twice, i.e., the electric structure of the near zone of the axionic dyon contains the domain similar to a double electric layer. We have shown that the described fold-like structures in the profile of the gravitational potential, and in the profiles of the electric and axion fields can exist, when the value of the dyon mass belongs to the interval enclosed between two critical masses, which depend on the cosmological constant.


title: Locally-Adaptive Nonparametric Online Learning

abstract: One of the main strengths of online algorithms is their ability to adapt to arbitrary data sequences. This is especially important in nonparametric settings, where performance is measured against rich classes of comparator functions that are able to fit complex environments. Although such hard comparators and complex environments may exhibit local regularities, efficient algorithms, which can provably take advantage of these local patterns, are hardly known. We fill this gap by introducing efficient online algorithms (based on a single versatile master algorithm) each adapting to one of the following regularities: (i) local Lipschitzness of the competitor function, (ii) local metric dimension of the instance sequence, (iii) local performance of the predictor across different regions of the instance space. Extending previous approaches, we design algorithms that dynamically grow hierarchical $\epsilon$-nets on the instance space whose prunings correspond to different "locality profiles" for the problem at hand. Using a technique based on tree experts, we simultaneously and efficiently compete against all such prunings, and prove regret bounds each scaling with a quantity associated with a different type of local regularity. When competing against "simple" locality profiles, our technique delivers regret bounds that are significantly better than those proven using the previous approach. On the other hand, the time dependence of our bounds is not worse than that obtained by ignoring any local regularities.


title: Deep Radial-Basis Value Functions for Continuous Control

abstract: A core operation in reinforcement learning (RL) is finding an action that is optimal with respect to a learned value function. This operation is often challenging when the learned value function takes continuous actions as input. We introduce deep radial-basis value functions (RBVFs): value functions learned using a deep network with a radial-basis function (RBF) output layer. We show that the maximum action-value with respect to a deep RBVF can be approximated easily and accurately. Moreover, deep RBVFs can represent any true value function owing to their support for universal function approximation. We extend the standard DQN algorithm to continuous control by endowing the agent with a deep RBVF. We show that the resultant agent, called RBF-DQN, significantly outperforms value-function-only baselines, and is competitive with state-of-the-art actor-critic algorithms.


title: On bounds for the mode and median of the generalized hyperbolic and related distributions

abstract: Except for certain parameter values, a closed form formula for the mode of the generalized hyperbolic (GH) distribution is not available. In this paper, we exploit results from the literature on modified Bessel functions and their ratios to obtain simple but tight two-sided inequalities for the mode of the GH distribution for general parameter values. As a special case, we deduce tight two-sided inequalities for the mode of the variance-gamma (VG) distribution, and through a similar approach we also obtain tight two-sided inequalities for the mode of the McKay Type I distribution. The analogous problem for the median is more challenging, but we conjecture some monotonicity results for the median of the VG and McKay Type I distributions, from we which we conjecture some tight two-sided inequalities for their medians. Numerical experiments support these conjectures and also lead us to a conjectured tight lower bound for the median of the GH distribution.


title: The effect of points fattening on del Pezzo surfaces

abstract: In this paper, we study the fattening effect of points over the complex numbers for del Pezzo surfaces $\mathbb{S}_r$ arising by blowing-up of $\mathbb{P}^2$ at $r$ general points, with $ r \in \{1, \dots, 8 \}$. Basic questions when studying the problem of points fattening on an arbitrary variety are what is the minimal growth of the initial sequence and how are the sets on which this minimal growth happens characterized geometrically. We provide complete answer for del Pezzo surfaces.


title: Polylidar -- Polygons from Triangular Meshes

abstract: This paper presents Polylidar, an efficient algorithm to extract non-convex polygons from 2D point sets, including interior holes. Plane segmented point clouds can be input into Polylidar to extract their polygonal counterpart, thereby reducing map size and improving visualization. The algorithm begins by triangulating the point set and filtering triangles by user configurable parameters such as triangle edge length. Next, connected triangles are extracted into triangular mesh regions representing the shape of the point set. Finally each region is converted to a polygon through a novel boundary following method which accounts for holes. Real-world and synthetic benchmarks are presented to comparatively evaluate Polylidar speed and accuracy. Results show comparable accuracy and more than four times speedup compared to other concave polygon extraction methods.


title: Modelling the solar transition region using an adaptive conduction method

abstract: Modelling the solar Transition Region with the use of an Adaptive Conduction (TRAC) method permits fast and accurate numerical solutions of the field-aligned hydrodynamic equations, capturing the enthalpy exchange between the corona and transition region, when the corona undergoes impulsive heating. The TRAC method eliminates the need for highly resolved numerical grids in the transition region and the commensurate very short time steps that are required for numerical stability. When employed with coarse spatial resolutions, typically achieved in multi-dimensional magnetohydrodynamic codes, the errors at peak density are less than 5% and the computation time is three orders of magnitude faster than fully resolved field-aligned models. This paper presents further examples that demonstrate the versatility and robustness of the method over a range of heating events, including impulsive and quasi-steady footpoint heating. A detailed analytical assessment of the TRAC method is also presented, showing that the approach works through all phases of an impulsive heating event because (i) the total radiative losses and (ii) the total heating when integrated over the transition region are both preserved at all temperatures under the broadening modifications of the method. The results from the numerical simulations complement this conclusion.


title: Weighted Inequalities for One-sided Fractional Minimal Function

abstract: In this article, we present weighted norm inequality for a fractional one-sided minimal function. We prove weighted weak and strong type norm inequalities for the one-sided minimal function on $\mathbb{R}.$ We construct two weight classes to ensure the sufficient and necessary part of the weak and strong type estimates. Later, we establish an equivalence relation between these two weight classes.


title: Automated Deep Abstractions for Stochastic Chemical Reaction Networks

abstract: Predicting stochastic cellular dynamics as emerging from the mechanistic models of molecular interactions is a long-standing challenge in systems biology: low-level chemical reaction network (CRN) models give raise to a highly-dimensional continuous-time Markov chain (CTMC) which is computationally demanding and often prohibitive to analyse in practice. A recently proposed abstraction method uses deep learning to replace this CTMC with a discrete-time continuous-space process, by training a mixture density deep neural network with traces sampled at regular time intervals (which can obtained either by simulating a given CRN or as time-series data from experiment). The major advantage of such abstraction is that it produces a computational model that is dramatically cheaper to execute, while preserving the statistical features of the training data. In general, the abstraction accuracy improves with the amount of training data. However, depending on a CRN, the overall quality of the method -- the efficiency gain and abstraction accuracy -- will also depend on the choice of neural network architecture given by hyper-parameters such as the layer types and connections between them. As a consequence, in practice, the modeller would have to take care of finding the suitable architecture manually, for each given CRN, through a tedious and time-consuming trial-and-error cycle. In this paper, we propose to further automatise deep abstractions for stochastic CRNs, through learning the optimal neural network architecture along with learning the transition kernel of the abstract process. Automated search of the architecture makes the method applicable directly to any given CRN, which is time-saving for deep learning experts and crucial for non-specialists. We implement the method and demonstrate its performance on a number of representative CRNs with multi-modal emergent phenotypes.


title: Dynamic clustering of time series data

abstract: We propose a new method for clustering multivariate time-series data based on Dynamic Linear Models. Whereas usual time-series clustering methods obtain static membership parameters, our proposal allows each time-series to dynamically change their cluster memberships over time. In this context, a mixture model is assumed for the time series and a flexible Dirichlet evolution for mixture weights allows for smooth membership changes over time. Posterior estimates and predictions can be obtained through Gibbs sampling, but a more efficient method for obtaining point estimates is presented, based on Stochastic Expectation-Maximization and Gradient Descent. Finally, two applications illustrate the usefulness of our proposed model to model both univariate and multivariate time-series: World Bank indicators for the renewable energy consumption of EU nations and the famous Gapminder dataset containing life-expectancy and GDP per capita for various countries.


title: A neural network model that learns differences in diagnosis strategies among radiologists has an improved area under the curve for aneurysm status classification in magnetic resonance angiography image series

abstract: Purpose: To construct a neural network model that can learn the different diagnosing strategies of radiologists to better classify aneurysm status in magnetic resonance angiography images. Materials and methods: This retrospective study included 3423 time-of-flight brain magnetic resonance angiography image series (subjects: male 1843 [mean age, 50.2 +/- 11.7 years], female 1580 [50.8 +/- 11.3 years]) recorded from November 2017 through January 2019. The image series were read independently for aneurysm status by one of four board-certified radiologists, who were assisted by an established deep learning-based computer-assisted diagnosis (CAD) system. The constructed neural networks were trained to classify the aneurysm status of zero to five aneurysm-suspicious areas suggested by the CAD system for each image series, and any additional aneurysm areas added by the radiologists, and this classification was compared with the judgment of the annotating radiologist. Image series were randomly allocated to training and testing data in an 8:2 ratio. The accuracy of the classification was compared by receiver operating characteristic analysis between the control model that accepted only image data as input and the proposed model that additionally accepted the information of who the annotating radiologist was. The DeLong test was used to compare areas under the curves (P < 0.05 was considered significant). Results: The area under the curve was larger in the proposed model (0.845) than in the control model (0.793), and the difference was significant (P < 0.0001). Conclusion: The proposed model improved classification accuracy by learning the diagnosis strategies of individual annotating radiologists.


title: Analytical investigations on formations of hairy neutral reflecting shells in the scalar-Gauss-Bonnet gravity

abstract: We study scalarization of spherically symmetric neutral reflecting shells in the scalar-tensor gravity. We consider neutral static massless scalar fields non-minimally coupled to the Gauss-Bonnet invariant. We obtain a relation representing the existence regime of hairy neutral reflecting shells. For parameters unsatisfying this relation, the massless scalar field cannot exist outside the neutral reflecting shell. In the parameter region where this relation holds, we get analytical solutions of scalar field hairs outside neutral reflecting shells.


title: FEA-Net: A Physics-guided Data-driven Model for Efficient Mechanical Response Prediction

abstract: An innovative physics-guided learning algorithm for predicting the mechanical response of materials and structures is proposed in this paper. The key concept of the proposed study is based on the fact that physics models are governed by Partial Differential Equation (PDE), and its loading/ response mapping can be solved using Finite Element Analysis (FEA). Based on this, a special type of deep convolutional neural network (DCNN) is proposed that takes advantage of our prior knowledge in physics to build data-driven models whose architectures are of physics meaning. This type of network is named as FEA-Net and is used to solve the mechanical response under external loading. Thus, the identification of a mechanical system parameters and the computation of its responses are treated as the learning and inference of FEA-Net, respectively. Case studies on multi-physics (e.g., coupled mechanical-thermal analysis) and multi-phase problems (e.g., composite materials with random micro-structures) are used to demonstrate and verify the theoretical and computational advantages of the proposed method.


title: Generation of smoothly-varying infill configurations from a continuous menu of cell patterns and the asymptotic analysis of its mechanical behaviour

abstract: We here introduce a novel scheme for generating smoothly-varying infill graded microstructural (IGM) configurations from a given menu of generating cells. The scheme was originally proposed for essentially improving the variety of describable configurations in a modified asymptotic homogenisation-based topology optimisation framework [1] for fast IGM design. But the proposed scheme, after modification, also demonstrates its unique values in two aspects of applications. First, it provides a fairly simple way of generating an IGM configuration continuously patching any given cell configurations. Second, it tenders a straightforward mean for decorating microstructures on a given manifold. We will further show that the form of topology description function given here effectively offers a platform for unifying most existing approaches for IGM generation. Fuelled by asymptotic analysis of the mechanical behaviour of the resulting IGM configurations, a topology optimisation scheme for compliance minimisation is introduced. We will finally show that, the use of the present scheme helps reduce the compliance value of an optimised structure by nearly a half, if compared with that from the original framework [1].


title: A toolbox of Equation-Free functions in Matlab\Octave for efficient system level simulation

abstract: The `equation-free toolbox' empowers the computer-assisted analysis of complex, multiscale systems. Its aim is to enable you to immediately use microscopic simulators to perform macro-scale system level tasks and analysis, because micro-scale simulations are often the best available description of a system. The methodology bypasses the derivation of macroscopic evolution equations by computing the micro-scale simulator only over short bursts in time on small patches in space, with bursts and patches well-separated in time and space respectively. We introduce the suite of coded equation-free functions in an accessible way, link to more detailed descriptions, discuss their mathematical support, and introduce a novel and efficient algorithm for Projective Integration. Some facets of toolbox development of equation-free functions are then detailed. Download the toolbox functions (https://github.com/uoa1184615/EquationFreeGit) and use to empower efficient and accurate simulation in a wide range of your science and engineering problems.


title: Topology optimization of 2D structures with nonlinearities using deep learning

abstract: The field of optimal design of linear elastic structures has seen many exciting successes that resulted in new architected materials and structural designs. With the availability of cloud computing, including high-performance computing, machine learning, and simulation, searching for optimal nonlinear structures is now within reach. In this study, we develop convolutional neural network models to predict optimized designs for a given set of boundary conditions, loads, and optimization constraints. We have considered the case of materials with a linear elastic response with and without stress constraint. Also, we have considered the case of materials with a hyperelastic response, where material and geometric nonlinearities are involved. For the nonlinear elastic case, the neo-Hookean model is utilized. For this purpose, we generate datasets composed of the optimized designs paired with the corresponding boundary conditions, loads, and constraints, using a topology optimization framework to train and validate the neural network models. The developed models are capable of accurately predicting the optimized designs without requiring an iterative scheme and with negligible inference computational time. The suggested pipeline can be generalized to other nonlinear mechanics scenarios and design domains.


title: Anomalous decay rate of quasinormal modes

abstract: The decay timescales of the quasinormal modes of a massive scalar field have an intriguing behavior: they either grow or decay with increasing angular harmonic numbers $\ell$, depending on whether the mass of the scalar field is small or large. We identify the properties of the effective potential of the scalar field that leads to this behavior and characterize it in detail. If the scalar field is non-minimally coupled, considered here, the scalar quasinormal modes will leak into the gravitational wave signal and will have decaying times that are comparable or smaller than those typical in General Relativity. Hence, these modes could be detectable in the future. Finally, we find that the anomalous behavior in the decay timescales of quasinormal modes is present in a much larger class of models beyond a simple massive scalar field.


title: Bistability in the Unstable Flow of Polymer Solutions Through Porous Media

abstract: Polymer solutions are often injected in porous media for applications such as oil recovery and groundwater remediation. As the fluid navigates the tortuous pore space, elastic stresses build up, causing the flow to become unstable at sufficiently large injection rates. However, it is poorly understood how the spatial and temporal characteristics of this unstable flow depend on pore space geometry. We elucidate this dependence by systematically varying the spacing between pore constrictions in a model porous medium. We find that when the pore spacing is large, unstable eddies form upstream of each pore, similar to observations of an isolated pore. By contrast, when the pore spacing is sufficiently small, the flow exhibits a surprising bistability, stochastically switching between two distinct unstable flow states. We hypothesize that this unusual behavior arises from the interplay between flow-induced polymer elongation and relaxation of polymers as they are advected through the porous medium. Consistent with this idea, we find that the flow state in a given pore persists for long times. Moreover, we find that the flow state is correlated between neighboring pores; however, these correlations do not persist long-range. Our results thus help to elucidate the rich array of flow behaviors that can arise in polymer solution flow through porous media.


title: Generalized Laplacian decomposition of vector fields on fractal surfaces

abstract: We consider the behavior of generalized Laplacian vector fields on a Jordan domain of $\mathbb{R}^{3}$ with fractal boundary. Our approach is based on properties of the Teodorescu transform and suitable extension of the vector fields. Specifically, the present article addresses the decomposition problem of a H\"older continuous vector field on the boundary (also called reconstruction problem) into the sum of two generalized Laplacian vector fields in the domain and in the complement of its closure, respectively. In addition, conditions on a H\"older continuous vector field on the boundary to be the trace of a generalized Laplacian vector field in the domain are also established.


title: Semisimplification of the category of tilting modules for GL_n

abstract: We describe the semisimplification of the monoidal category of tilting modules for the algebraic group GL_n in characteristic p > 0. In particular, we compute the dimensions of the indecomposable tilting modules modulo p.


title: $O(D)-$equivariant fuzzy hyperspheres

abstract: Fuzzy hyperspheres $S^d_\Lambda$ of dimension $d>2$ are constructed here generalizing the procedure adopted in [G. Fiore, F. Pisacane, J. Geom. Phys. 132 (2018), 423-451] for $d=1,2$. The starting point is an ordinary quantum particle in $\mathbb{R}^D$, $D:=d+1$, subject to a rotation invariant potential well $V(r)$ with a very sharp minimum on the sphere of radius $r=1$. The subsequent imposition of a sufficiently low energy cutoff `freezes' the radial excitations, this makes only a finite-dimensional Hilbert subspace $\mathcal{H}_{\Lambda,D}$ accessible and on it the coordinates noncommutative {\it \`a la Snyder}. In addition, the coordinate operators generate the whole algebra of observables $\mathcal{A}_{\Lambda,D}$ which turns out to be realizable through a suitable irreducible vector representation of $Uso(D+1)$. This construction is equivariant not only under $SO(D)$, but under the full orthogonal group $O(D)$, and making the cutoff and the depth of the well grow with a natural number $\Lambda$, the result is a sequence $S^d_{\Lambda}$ of fuzzy spheres converging to $S^d$ as $\Lambda\to\infty$ (where one recovers ordinary quantum mechanics on $S^d$).


title: Rotational and translational dynamics in dense fluids of patchy particles

abstract: We explore the effect of directionality on rotational and translational relaxation in glassy systems of patchy particles. Using molecular dynamics simulations we analyze the impact of two distinct patch geometries, one that enhance local icosahedral structure and one which does not strongly affect local order. We find that in nearly all investigated cases, rotational relaxation takes place on a much faster time scale than translational relaxation. By comparing to a simplified dynamical Monte Carlo model, we illustrate that rotational diffusion can be qualitatively explained as purely local motion within a fixed environment, which is not coupled strongly to the cage-breaking dynamics required for translational relaxation. Nonetheless, icosahedral patch placement has a profound effect on the local structure of the system, resulting in a dramatic slowdown at low temperatures which is strongest at an intermediate "optimal" patch size.


title: Model nuclear energy density functionals derived from ab initio calculations

abstract: We present the first application of a new approach, proposed in [Journal of Physics G: Nuclear and Particle Physics, 43, 04LT01 (2016)] to derive coupling constants of the Skyrme energy density functional (EDF) from ab initio Hamiltonian. By perturbing the ab initio Hamiltonian with several functional generators defining the Skyrme EDF, we create a set of metadata that is then used to constrain the coupling constants of the functional. We use statistical analysis to obtain such an ab initio-equivalent Skyrme EDF. We find that the resulting functional describes properties of atomic nuclei and infinite nuclear matter quite poorly. This may point out to the necessity of building up the ab initio-equivalent functionals from more sophisticated generators. However, we also indicate that the current precision of the ab initio calculations may be insufficient for deriving meaningful nuclear EDFs.


title: An upper bound conjecture for the Yokota invariant

abstract: We conjecture an upper bound on the growth of the Yokota invariant of polyhedral graphs, extending a previous result on the growth of the $6j$-symbol. Using Barrett's Fourier transform we are able to prove this conjecture in a large family of examples. As a consequence of this result, we prove the Turaev-Viro Volume Conjecture for a new infinite family of hyperbolic manifolds.


title: Experimental Implementation of Quantum Walks on IBM Quantum Computers

abstract: The development of universal quantum computers has achieved remarkable success in recent years, culminating with the quantum supremacy reported by Google. Now is possible to implement short-depth quantum circuits with dozens of qubits and to obtain results with significant fidelity. Quantum walks are good candidates to be implemented on the available quantum computers. In this work, we implement discrete-time quantum walks with one and two interacting walkers on cycles, two-dimensional lattices, and complete graphs on IBM quantum computers. We are able to obtain meaningful results using the cycle, the two-dimensional lattice, and the complete graph with 16 nodes each, which require 4-qubit quantum circuits up to depth 100.


title: Bounding tangencies of sections on elliptic surfaces

abstract: Given an elliptic surface $\mathcal{E}\to\mathcal{C}$ over a field $k$ of characteristic zero equipped with zero section $O$ and another section $P$ of infinite order, we give a simple and explicit upper bound on the number of points where $O$ is tangent to a multiple of $P$.