title: Dynamical topological invariant after a quantum quench

abstract: We show how to define a dynamical topological invariant for general one-dimensional topological systems after a quantum quench. Focusing on two-band topological insulators, we demonstrate that the reduced momentum-time manifold can be viewed as a series of submanifold $S^2$, and thus we are able to define a dynamical topological invariant on each of the sphere. We also unveil the intrinsic relation between the dynamical topological invariant and the difference of topological invariant of the initial and final static Hamiltonian. By considering some concrete examples, we illustrate the calculation of the dynamical topological invariant and its geometrical meaning explicitly.


title: Physical realization of topological quantum walks on IBM-Q and beyond

abstract: We discuss an efficient physical realization of topological quantum walks on a finite lattice. The $N$-point lattice is realized with $\log_2 N$ qubits, and the quantum circuit utilizes a number of quantum gates which is polynomial in the number of qubits. In a certain scaling limit, we show that a large number of steps is implemented with a number of quantum gates which is independent of the number of steps. We ran the quantum algorithm on the IBM-Q five-qubit quantum computer, thus experimentally demonstrating topological features, such as boundary bound states, on a lattice with $N=4$ points.


title: Morse Spectra, Homology Measures and Parametric Packing Problems

abstract: We suggest several mathematical counterparts to the idea of "effective degrees of freedom" and formulate specific questions, much of which are inspired by Larry Guth's results and ideas on the Hermann Weyl kind of asymptotics of the Morse (co)homology spectra of the volume energy function on the spaces of cycles in balls.


title: Compactly-Supported Smooth Interpolators for Shape Modeling with Varying Resolution

abstract: In applications that involve interactive curve and surface modeling, the intuitive manipulation of shapes is crucial. For instance, user interaction is facilitated if a geometrical object can be manipulated through control points that interpolate the shape itself. Additionally, models for shape representation often need to provide local shape control and they need to be able to reproduce common shape primitives such as ellipsoids, spheres, cylinders, or tori. We present a general framework to construct families of compactly-supported interpolators that are piecewise-exponential polynomial. They can be designed to satisfy regularity constraints of any order and they enable one to build parametric deformable shape models by suitable linear combinations of interpolators. They allow to change the resolution of shapes based on the refinability of B-splines. We illustrate their use on examples to construct shape models that involve curves and surfaces with applications to interactive modeling and character design.


title: Asymptotic expansion of the mean-field approximation

abstract: We established and estimate the full asymptotic expansion in integer powers of 1 N of the [ $\sqrt$ N ] first marginals of N-body evolutions lying in a general paradigm containing Kac models and non-relativistic quantum evolution. We prove that the coefficients of the expansion are, at any time, explicitly computable given the knowledge of the linearization on the one-body meanfield kinetic limit equation. Instead of working directly with the corresponding BBGKY-type hierarchy, we follows a method developed in [22] for the meanfield limit, dealing with error terms analogue to the v-functions used in previous works. As a by-product we get that the rate of convergence to the meanfield limit in 1 N is optimal.

categories: quant-ph math-ph math.AP math.DS math.MP


title: A Generalized Algebraic Approach to Optimizing SC-LDPC Codes

abstract: Spatially coupled low-density parity-check (SC-LDPC) codes are sparse graph codes that have recently become of interest due to their capacity-approaching performance on memoryless binary input channels. In this paper, we unify all existing SC-LDPC code construction methods under a new generalized description of SC-LDPC codes based on algebraic lifts of graphs. We present an improved low-complexity counting method for the special case of $(3,3)$-absorbing sets for array-based SC-LDPC codes, which we then use to optimize permutation assignments in SC-LDPC code construction. We show that codes constructed in this way are able to outperform previously published constructions, in terms of the number of dominant absorbing sets and with respect to both standard and windowed decoding.


title: Regularization effects of a noise propagating through a chain of differential equations: an almost sharp result

abstract: We investigate the effects of the propagation of a non-degenerate Brownian noise through a chain of deterministic differential equations whose coefficients are rough and satisfy a weak like H{\"o}rmander structure (i.e. a non-degeneracy condition w.r.t. the components which transmit the noise). In particular we characterize, through suitable counterexamples , almost sharp regularity exponents that ensure that weak well posedness holds for the associated SDE. As a by-product of our approach, we also derive some density estimates of Krylov type for the weak solutions of the considered SDEs.


title: Kinetic modeling of multiple scattering of acoustic waves in randomly heterogeneous flows

abstract: We study the propagation of sound waves in a three-dimensional, infinite ambient flow with weak random fluctuations of the mean particle velocity and speed of sound. We more particularly address the regime where the acoustic wavelengths are comparable to the correlation lengths of the weak inhomogeneities--the so-called weak coupling limit. The analysis is carried on starting from the linearized Euler equations and the convected wave equation with variable density and speed of sound, which can be derived from the nonlinear Euler equations. We use a multi-scale expansion of the Wigner distribution of a velocity potential associated to the waves to derive a radiative transfer equation describing the evolution of the angularly resolved wave action in space/time phase space. The latter experiences convection, refraction and scattering when it propagates through the heterogeneous ambient flow, although the overall wave action is conserved. The convection and refraction phenomena are accounted for by the convective part of the transport equation and depend on the smooth variations of the ambient quantities. The scattering phenomenon is accounted for by the collisional part of the transport equation and depends on the cross-power spectral densities of the fluctuations of the ambient quantities at the wavelength scales. The refraction, phase shift, spectral broadening, and multiple scattering effects of the high-frequency regimes described in various previous publications are thus encompassed by the proposed model. The overall derivation is based on the interpretation of spatial-temporal Wigner transforms in terms of semiclassical operators in their standard quantization.


title: Prototype tests for a highly granular scintillator-based hadronic calorimeter

abstract: Within the CALICE collaboration, several concepts for the hadronic calorimeter of a future lepton collider detector are studied. After having demonstrated the capabilities of the measurement methods in "physics prototypes", the focus now lies on improving their implementation in "technological prototypes", that are scalable to the full linear collider detector. The Analogue Hadronic Calorimeter (AHCAL) concept is a sampling calorimeter of tungsten or steel absorber plates and plastic scintillator tiles read out by silicon photomultipliers (SiPMs) as active components. The front-end electronics is fully integrated into the active layers of the calorimeter and is designed for minimal power consumption (i.e. power pulsing). The versatile electronics enables the prototype to be equipped with different types of scintillator tiles and SiPMs. In recent beam tests, a prototype with $\sim$3700 channels, equipped with several types of scintillator tiles and SiPMs, was exposed to electron, muon and hadron beams. The experience of these beam tests resulted in an optimal detector design with surface-mounted SiPMs suitable for the automated mass assembly. The proceeding will cover topics including the testbeam measurements with the AHCAL technological prototype, the improved detector design and the ongoing development of a large prototype for hadronic showers.


title: Near-misses in Wilf's conjecture

abstract: Let S $\subseteq$ N be a numerical semigroup with multiplicity m, conductor c and minimal generating set P. Let L = S $\cap$ [0, c -- 1] and W(S) = |P||L| -- c. In 1978, Herbert Wilf asked whether W(S) $\ge$ 0 always holds, a question known as Wilf's conjecture and open since then. A related number W0(S), satisfying W0(S) $\le$ W(S), has recently been introduced. We say that S is a near-miss in Wilf's conjecture if W0(S) < 0. Near-misses are very rare. Here we construct infinite families of them, with c = 4m and W0(S) arbitrarily small, and we show that the members of these families still satisfy Wilf's conjecture.


title: Explicit calculations with Eisenstein series

abstract: We find explicit change-of-basis formulas between Eisenstein series attached to cusps, and newform Eisenstein series attached to pairs of primitive Dirichlet characters. As a consequence, we prove a Bruggeman-Kuznetsov formula for newforms of square-free level and trivial nebentypus. We also derive, in explicit form, the Fourier expansion, Hecke eigenvalues, Atkin-Lehner pseudo-eigenvalues, and other data associated to these Eisenstein series, with arbitrary integer weight, level, and nebentypus.


title: An extension of the Polyak convexity principle with application to nonconvex optimization

abstract: The main problem considered in the present paper is to single out classes of convex sets, whose convexity property is preserved under nonlinear smooth transformations. Extending an approach due to B.T. Polyak, the present study focusses on the class of uniformly convex subsets of Banach spaces. As a main result, a quantitative condition linking the modulus of convexity of such kind of set, the regularity behaviour around a point of a nonlinear mapping and the Lipschitz continuity of its derivative is established, which ensures the images of uniformly convex sets to remain uniformly convex. Applications of the resulting convexity principle to the existence of solutions, their characterization and to the Lagrangian duality theory in constrained nonconvex optimization are then discussed.


title: Base Station Diversity Propagation Measurements at 73 GHz Millimeter-Wave for 5G Coordinated Multipoint (CoMP) Analysis

abstract: This paper describes wideband (1 GHz) base station diversity and coordinated multipoint (CoMP)-style large-scale measurements at 73 GHz in an urban microcell open square scenario in downtown Brooklyn, New York on the NYU campus. The measurements consisted of ten random receiver locations at pedestrian level (1.4 meters) and ten random transmitter locations at lamppost level (4.0 meters) that provided 36 individual transmitter-receiver (TX-RX) combinations. For each of the 36 radio links, extensive directional measurements were made to give insights into small-cell base station diversity at millimeter-wave (mmWave) bands. High-gain steerable horn antennas with 7-degree and 15-degree half-power beamwidths (HPBW) were used at the transmitter (TX) and receiver (RX), respectively. For each TX-RX combination, the TX antenna was scanned over a 120-degree sector and the RX antenna was scanned over the entire azimuth plane at the strongest RX elevation plane and two other elevation planes on both sides of the strongest elevation angle, separated by the 15-degree HPBW. Directional and omnidirectional path loss models were derived and match well with the literature. Signal reception probabilities derived from the measurements for one to five base stations that served a single RX location show significant coverage improvement over all potential beamformed RX antenna pointing angles. CDFs for nearest neighbor and Best-N omnidirectional path loss and cell outage probabilities for directional antennas provide insights into coverage and interference for future mmWave small-cells that will exploit macro-diversity and CoMP.


title: Multilevel Modeling with Structured Penalties for Classification from Imaging Genetics data

abstract: In this paper, we propose a framework for automatic classification of patients from multimodal genetic and brain imaging data by optimally combining them. Additive models with unadapted penalties (such as the classical group lasso penalty or $L_1$-multiple kernel learning) treat all modalities in the same manner and can result in undesirable elimination of specific modalities when their contributions are unbalanced. To overcome this limitation, we introduce a multilevel model that combines imaging and genetics and that considers joint effects between these two modalities for diagnosis prediction. Furthermore, we propose a framework allowing to combine several penalties taking into account the structure of the different types of data, such as a group lasso penalty over the genetic modality and a $L_2$-penalty on imaging modalities. Finally , we propose a fast optimization algorithm, based on a proximal gradient method. The model has been evaluated on genetic (single nucleotide polymorphisms-SNP) and imaging (anatomical MRI measures) data from the ADNI database, and compared to additive models. It exhibits good performances in AD diagnosis; and at the same time, reveals relationships between genes, brain regions and the disease status.


title: The Bramson delay in the non-local Fisher-KPP equation

abstract: We consider the non-local Fisher-KPP equation modeling a population with individuals competing with each other for resources with a strength related to their distance, and obtain the asymptotics for the position of the invasion front starting from a localized population. Depending on the behavior of the competition kernel at infinity, the location of the front is either $2t - ({3}/2)\log t + O(1)$, as in the local case, or $2t - O(t^\beta)$ for some explicit $\beta \in (0,1)$. Our main tools here are alocal-in-time Harnack inequality and an analysis of the linearized problem with a suitable moving Dirichlet boundary condition. Our analysis also yields, for any $\beta\in(0,1)$, examples of Fisher-KPP type non-linearities $f\_\beta$ such that the front for the localFisher-KPP equation with reaction term $f\_\beta$ is at $2t - O(t^\beta)$.


title: Multipactor in High Power Radio-Frequency Systems for Nuclear Fusion

abstract: Magnetic confinement fusion researches are an approach to demonstrate the feasibility of nuclear fusion power for energy production. In these experiments, mega-watt range power Radio-Frequency waves, from tens of MHz to hundreds of GHz, are injected into magnetically confined plasmas in order to increase its temperature and to extend its duration. These RF systems are subjected to the magnetic field environment of the experiments, ranging from few tenth of Tesla to Tesla, and various orientations depending of their locations. As these RF systems made of copper, silver or ceramics are located in vacuum environments, they are subject to multipactor discharges. These discharges are generally considered detrimental since they can lead to detuned RF systems, limit the RF power transmission in the plasma and eventually damage RF sources or components. In some case, especially in the MHz range of frequency, multipactor discharges can be desired for vacuum conditioning. This paper presents the various RF systems used in magnetic fusion experiments and review the work performed in the fusion research community on multipactor discharges for high power coaxial (MHz) and rectangular waveguides (GHz) transmission lines, with their practical implications on power delivery into the plasma.


title: Ultrasonic waves in classical gases

abstract: The velocity and absorption coefficient for the plane sound waves in a classical gas are obtained by solving the Boltzmann kinetic equation, which describes the reaction of the single-particle distribution function to a periodic external field. Within the linear response theory, the nonperturbative dispersion equation valid for all sound frequencies is derived and solved numerically. The results are in agreement with the approximate analytical solutions found for both the frequent- and rare-collision regimes.


title: Angular Accuracy of Steerable Feature Detectors

abstract: The detection of landmarks or patterns is of interest for extracting features in biological images. Hence, algorithms for finding these keypoints have been extensively investigated in the literature, and their localization and detection properties are well known. In this paper, we study the complementary topic of local orientation estimation, which has not received similar attention. Simply stated, the problem that we address is the following: estimate the angle of rotation of a pattern with steerable filters centered at the same location, where the image is corrupted by colored isotropic Gaussian noise. For this problem, we use a statistical framework based on the Cram\'{e}r-Rao lower bound (CRLB) that sets a fundamental limit on the accuracy of the corresponding class of estimators. We propose a scheme to measure the performance of estimators based on steerable filters (as a lower bound), while considering the connection to maximum likelihood estimation. Beyond the general results, we analyze the asymptotic behaviour of the lower bound in terms of the order of steerablility and propose an optimal subset of components that minimizes the bound. We define a mechanism for selecting optimal subspaces of the span of the detectors. These are characterized by the most relevant angular frequencies. Finally, we project our template to a basis of steerable functions and experimentally show that the prediction accuracy achieves the predicted CRLB. As an extension, we also consider steerable wavelet detectors.


title: Radiofrequency and Mechanical Tests of Silver Coated CuCrZr Contacts for the ITER Ion Cyclotron Antenna

abstract: The ITER Ion Cyclotron Resonance Heating (ICRH) system is designed to couple to the plasma 20 MW of RF power from two antennas in the 40-55 MHz frequency range during long pulses of up to 3600 s and under various plasma conditions with Edge Localized Modes. Radio-Frequency (RF) contacts are integrated within the ITER ICRH launcher in order to ensure the RF current continuity and ease the mechanical assembly by allowing the free thermal expansion of the Removable Vacuum Transmission Line coaxial conductors during RF operations or during 250{\textdegree}C baking phases. A material study has been carried out to determine which materials and associated coatings are relevant for RF contacts application in ITER. In parallel, RF tests have been performed with a new prototype of Multi-Contact (r) LA-CUT/0,25/0 contacts made of silver-coated CuCrZr louvers. During these tests on a RF vacuum resonator, currents between 1.2 kA and 1.3 kA peak have been reached a few tens of times in steady-state conditions without any visible damage on the louvers. A final 62MHz pulse ending in a 300s flat top at 1.9kA resulted in severe damage to the contact. In addition, a test bed which performs sliding test cycles has been built in order to reproduce the wear of the contact prototype after 30 000 sliding cycles on a 3 mm stroke at 175{\textdegree}C under vacuum. The silver coating of the louvers is removed after approximately a hundred cycles whilst, to the contrary, damage to the CuCrZr louvers is relatively low.


title: Differential-activity driven instabilities in biphasic active matter

abstract: Active stresses can cause instabilities in contractile gels and living tissues. Here we describe a generic hydrodynamic theory that treats these systems as a mixture of two phases of varying activity and different mechanical properties. We find that differential activity between the phases provides a mechanism causing a demixing instability. We follow the nonlinear evolution of the instability and characterize a phase diagram of the resulting patterns. Our study complements other instability mechanisms in mixtures such as differential growth, shape, motion or adhesion.


title: LinXGBoost: Extension of XGBoost to Generalized Local Linear Models

abstract: XGBoost is often presented as the algorithm that wins every ML competition. Surprisingly, this is true even though predictions are piecewise constant. This might be justified in high dimensional input spaces, but when the number of features is low, a piecewise linear model is likely to perform better. XGBoost was extended into LinXGBoost that stores at each leaf a linear model. This extension, equivalent to piecewise regularized least-squares, is particularly attractive for regression of functions that exhibits jumps or discontinuities. Those functions are notoriously hard to regress. Our extension is compared to the vanilla XGBoost and Random Forest in experiments on both synthetic and real-world data sets.


title: Local-global Galois theory of arithmetic function fields

abstract: We study the relationship between the local and global Galois theory of function fields over a complete discretely valued field. We give necessary and sufficient conditions for local separable extensions to descend to global extensions, and for the local absolute Galois group to inject into the global absolute Galois group. As an application we obtain a local-global principle for the index of a variety over such a function field. In this context we also study algebraic versions of van Kampen's theorem, describing the global absolute Galois group as a pushout of local absolute Galois groups.


title: Learning time-dependent noise to reduce logical errors: Real time error rate estimation in quantum error correction

abstract: Quantum error correction is important to quantum information processing, which allows us to reliably process information encoded in quantum error correction codes. Efficient quantum error correction benefits from the knowledge of error rates. We propose a protocol for monitoring error rates in real time without interrupting the quantum error correction. Any adaptation of the quantum error correction code or its implementation circuit is not required. The protocol can be directly applied to the most advanced quantum error correction techniques, e.g. surface code. A Gaussian processes algorithm is used to estimate and predict error rates based on error correction data in the past. We find that using these estimated error rates, the probability of error correction failures can be significantly reduced by a factor increasing with the code distance.


title: Evaluation of Harmonic Sums with Integrals

abstract: We consider the sums $S(k)=\sum_{n=0}^{\infty}\frac{(-1)^{nk}}{(2n+1)^k}$ and $\zeta(2k)=\sum_{n=1}^{\infty}\frac{1}{n^{2k}}$ with $k$ being a positive integer. We evaluate these sums with multiple integration, a modern technique. First, we start with three different double integrals that have been previously used in the literature to show $S(2)=\pi^2/8,$ which implies Euler's identity $\zeta(2)=\pi^2/6.$ Then, we generalize each integral in order to find the considered sums. The $k$ dimensional analogue of the first integral is the density function of the quotient of $k$ independent, nonnegative Cauchy random variables. In seeking this function, we encounter a special logarithmic integral that we can directly relate to $S(k).$ The $k$ dimensional analogue of the second integral, upon a change of variables, is the volume of a convex polytope, which can be expressed as a probability involving certain pairwise sums of $k$ independent uniform random variables. We use combinatorial arguments to find the volume, which in turn gives new closed formulas for $S(k)$ and $\zeta(2k).$ The $k$ dimensional analogue of the last integral, upon another change of variables, is an integral of the joint density function of $k$ Cauchy random variables over a hyperbolic polytope. This integral can be expressed as a probability involving certain pairwise products of these random variables, and it is equal to the probability from the second generalization. Thus, we specifically highlight the similarities in the combinatorial arguments between the second and third generalizations.


title: On-Chip Cooling by Heating with Superconducting Tunnel Junctions

abstract: Heat management and refrigeration are key concepts for nanoscale devices operating at cryogenic temperatures. The design of an on-chip mesoscopic refrigerator that works thanks to the input heat is presented, thus realizing a solid state implementation of the concept of cooling by heating. The system consists of a circuit featuring a thermoelectric element based on a ferromagnetic insulator-superconductor tunnel junction (N-FI-S) and a series of two normal metal-superconductor tunnel junctions (SINIS). The N-FI-S element converts the incoming heat in a thermovoltage, which is applied to the SINIS, thereby yielding cooling. The cooler's performance is investigated as a function of the input heat current for different bath temperatures. We show that this system can efficiently employ the performance of SINIS refrigeration, with a substantial cooling of the normal metal island. Its scalability and simplicity in the design makes it a promising building block for low-temperature on-chip energy management applications.


title: A quantum light emitting diode for the standard telecom window around 1550 nm

abstract: For the development of long-distance quantum networks, sources of single photons and entangled photon pairs emitting in the low-loss wavelength region around 1550 nm are a crucial building block. Here we show that quantum dot devices based on indium phosphide are capable of electrically injected single photon emission in this wavelength region with multiphoton events suppressed down to 0.11$\pm$0.02. Using the biexciton cascade mechanism, they further produce entangled photons with a fidelity of 87$\pm$4%, sufficient for the application of one-way error correction protocols. The new material allows for entangled photon generation up to an operating temperature of 93 K, reaching a regime accessible by electric coolers. The quantum photon source can be directly integrated with existing long distance quantum communication and cryptography systems and provides a new material platform for developing future quantum network hardware.


title: Beam Management for Millimeter Wave Beamspace MU-MIMO Systems

abstract: Millimeter wave (mmWave) communication has attracted increasing attention as a promising technology for 5G networks. One of the key architectural features of mmWave is the use of massive antenna arrays at both the transmitter and the receiver sides. Therefore, by employing directional beamforming (BF), both mmWave base stations (MBSs) and mmWave users (MUEs) are capable of supporting multi-beam simultaneous transmissions. However, most researches have only considered a single beam, which means that they do not make full potential of mmWave. In this context, in order to improve the performance of short-range indoor mmWave networks with multiple reflections, we investigate the challenges and potential solutions of downlink multi-user multi-beam transmission, which can be described as a high-dimensional (i.e., beamspace) multi-user multiple-input multiple-output (MU-MIMO) technique, including multi-user BF training, simultaneous users' grouping, and multi-user multibeam power allocation. Furthermore, we present the theoretical and numerical results to demonstrate that beamspace MU-MIMO compared with single beam transmission can largely improve the rate performance of mmWave systems.


title: Continuous Adaptation via Meta-Learning in Nonstationary and Competitive Environments

abstract: Ability to continuously learn and adapt from limited experience in nonstationary environments is an important milestone on the path towards general intelligence. In this paper, we cast the problem of continuous adaptation into the learning-to-learn framework. We develop a simple gradient-based meta-learning algorithm suitable for adaptation in dynamically changing and adversarial scenarios. Additionally, we design a new multi-agent competitive environment, RoboSumo, and define iterated adaptation games for testing various aspects of continuous adaptation strategies. We demonstrate that meta-learning enables significantly more efficient adaptation than reactive baselines in the few-shot regime. Our experiments with a population of agents that learn and compete suggest that meta-learners are the fittest.


title: Ambient and high pressure phases of tin sulfide compounds

abstract: Tin sulfides SnxSy are important class of materials that are actively investigated as novel photovoltaic and water splitting materials. First-principles evolutionary crystal structure search is performed with the goal of constructing the complete phase diagram of SnxSy and discover new phases as well as new compounds of varying stoichiometry at ambient conditions and pressures up to 100 GPa. The ambient phase of SnS2 with P3m1 symmetry remains stable up to 30 GPa. Another ambient phase, SnS experiences a series of phase transformations: including \alpha-SnS to \beta-SnS at 9 GPa, followed by \beta-SnS to \gamma-SnS at 40 GPa. \gamma-SnS is a new high-pressure metallic phase with Pm3m space group symmetry stable up to 100 GPa, which becomes a superconductor with a maximum Tc= 9.6 K at 40 GPa. Another new metallic compound, Sn3S4 with I43d space group symmetry, is predicted to be stable at pressures above 15 GPa, which also becomes a superconductor with a maximum Tc=21.8 K at 30 GPa.


title: Constraining the interaction between dark sectors with future HI intensity mapping observations

abstract: We study a model of interacting dark matter and dark energy, in which the two components are coupled. We calculate the predictions for the 21-cm intensity mapping power spectra, and forecast the detectability with future single-dish intensity mapping surveys (BINGO, FAST and SKA-I). Since dark energy is turned on at $z\sim 1$, which falls into the sensitivity range of these radio surveys, the HI intensity mapping technique is an efficient tool to constrain the interaction. By comparing with current constraints on dark sector interactions, we find that future radio surveys will produce tight and reliable constraints on the coupling parameters.


title: One dimensional phase transition problem modelling striped spin orbit coupled Bose-Einstein condensates

abstract: We study the behaviour of a Modica-Mortola phase transition type problem with a non-homogeneous Neumann boundary condition. According to the parameters of the problem, this leads to the existence of either one component occupying most of the condensate with an outer boundary layer containing the other component, or to many interfaces, on a periodic pattern. This is related to the striped behaviour of a two component Bose-Einstein condensate with spin orbit coupling in one dimension. We prove that minimizers of the full Gross-Pitaevskii energy in 1D converge, in the Thomas-Fermi limit of strong intra-component interaction, to those of the simplified Modica-Mortola problem we have studied in the first part.


title: Optimized Frameless ALOHA for Cooperative Base Stations with Overlapped Coverage Areas

abstract: Herein, we consider the problem of cooperative multi-access in the presence of overlapped coverage areas. Assuming a frameless ALOHA transmission scheme, we derive exact analytical throughput expressions for throughput in the aforementioned scenarios as a function of the frame length of the system and for arbitrary average numbers of users transmitting in each slot (target degree). After obtaining these original expressions, we then formulate a utility function whose maximization (obtained, e.g., through genetic algorithms) yields unequal and optimum target degrees to be employed by users in each group in order to maximize the peak throughput of the whole system, while satisfying a given prescribed outage. A comparison of the resulting cooperative multiple base station (BS) multi-access scheme against optimized single-BS frameless ALOHA systems - which presume the perfect isolation of users at each BS and an equal optimum target degree for all users - indicates a significant gain in overall throughput, thereby revealing that a multi-access diversity gain can be reaped by allowing groups of users from different BSs to overlap.


title: NMRCloudQ: A Quantum Cloud Experience on a Nuclear Magnetic Resonance Quantum Computer

abstract: As of today, no one can tell when a universal quantum computer with thousands of logical quantum bits (qubits) will be built. At present, most quantum computer prototypes involve less than ten individually controllable qubits, and only exist in laboratories for the sake of either the great costs of devices or professional maintenance requirements. Moreover, scientists believe that quantum computers will never replace our daily, every-minute use of classical computers, but would rather serve as a substantial addition to the classical ones when tackling some particular problems. Due to the above two reasons, cloud-based quantum computing is anticipated to be the most useful and reachable form for public users to experience with the power of quantum. As initial attempts, IBM Q has launched influential cloud services on a superconducting quantum processor in 2016, but no other platforms has followed up yet. Here, we report our new cloud quantum computing service -- NMRCloudQ (http://nmrcloudq.com/zh-hans/), where nuclear magnetic resonance, one of the pioneer platforms with mature techniques in experimental quantum computing, plays as the role of implementing computing tasks. Our service provides a comprehensive software environment preconfigured with a list of quantum information processing packages, and aims to be freely accessible to either amateurs that look forward to keeping pace with this quantum era or professionals that are interested in carrying out real quantum computing experiments in person. In our current version, four qubits are already usable with in average 1.26% single-qubit gate error rate and 1.77% two-qubit controlled-NOT gate error rate via randomized benchmaking tests. Improved control precisions as well as a new seven-qubit processor are also in preparation and will be available later.


title: Accelerating Energy Games Solvers on Modern Architectures

abstract: Quantitative games, where quantitative objectives are defined on weighted game arenas, provide natural tools for designing faithful models of embedded controllers. Instances of these games that recently gained interest are the so called Energy Games. The fast-known algorithm solves Energy Games in O(EVW) where W is the maximum weight. Starting from a sequential baseline implementation, we investigate the use of massively data computation capabilities supported by modern Graphics Processing Units to solve the `initial credit problem' for Energy Games. We present four different parallel implementations on multi-core CPU and GPU systems. Our solution outperforms the baseline implementation by up to 36x speedup and obtains a faster convergence time on real-world graphs.


title: Stochastic Navier-Stokes equation and advection of a tracer field: One-loop renormalization near $d=4$

abstract: The renormalization group approach and the operator product expansion technique are applied to the model of a tracer field advected by the Navier-Stokes velocity ensemble for a compressible fluid. The model is considered in the vicinity of the specific space dimension $d=4$. The properties of the equal-time structure functions are investigated. The multifractal behaviour of various correlation functions is established. All calculations are performed in the leading one-loop approximation.


title: Decentralized Resource Discovery and Management for Future Manycore Systems

abstract: The next generation of many-core enabled large-scale computing systems relies on thousands of billions of heterogeneous processing cores connected to form a single computing unit. In such large-scale computing environments, resource management is one of the most challenging, and complex issues for efficient resource sharing and utilization, particularly as we move toward Future ManyCore Systems (FMCS). This work proposes a novel resource management scheme for future peta-scale many-core-enabled computing systems, based on hybrid adaptive resource discovery, called ElCore. The proposed architecture contains a set of modules which will dynamically be instantiated on the nodes in the distributed system on demand. Our approach provides flexibility to allocate the required set of resources for various types of processes/applications. It can also be considered as a generic solution (with respect to the general requirements of large scale computing environments) which brings a set of interesting features (such as auto-scaling, multitenancy, multi-dimensional mapping, etc,.) to facilitate its easy adaptation to any distributed technology (such as SOA, Grid and HPC many-core). The achieved evaluation results assured the significant scalability and the high quality resource mapping of the proposed resource discovery and management over highly heterogeneous, hierarchical and dynamic computing environments with respect to several scalability and efficiency aspects while supporting flexible and complex queries with guaranteed discovery results accuracy. The simulation results prove that, using our approach, the mapping between processes and resources can be done with high level of accuracy which potentially leads to a significant enhancement in the overall system performance.


title: Electroporation dynamics for different pore lifetimes based on the standard model

abstract: Standard model of electropration (EP) has long emphasized a single pore lifetime to explain post-pulse transport across cell membranes. However, pore lifetimes estimated from molecular dynamics (MD) models and and those measured from experimental data differ by several orders of magnitude. We hypothesize that a broad distribution of lifetimes may describe the post-pulse behavior. Here, we show that pore distribution, number and size of pores, show interesting behavior in different ranges of pore lifetimes. Interestingly, for large electric fields (greater than 1 kV/cm) and short pore lifetimes (approximately 100 ns), an significant loss in pore number occurs during the pulse. Given the large number of EP applications that apply such fields, this phenomenon may be crucial to post-pulse response of cell membrane to such fields.


title: The shoreline problem for the one-dimensional shallow water and Green-Naghdi equations

abstract: The Green-Naghdi equations are a nonlinear dispersive perturbation of the nonlinear shallow water equations, more precise by one order of approximation. These equations are commonly used for the simulation of coastal flows, and in particular in regions where the water depth vanishes (the shoreline). The local well-posedness of the Green-Naghdi equations (and their justification as an asymptotic model for the water waves equations) has been extensively studied, but always under the assumption that the water depth is bounded from below by a positive constant. The aim of this article is to remove this assumption. The problem then becomes a free-boundary problem since the position of the shoreline is unknown and driven by the solution itself. For the (hyperbolic) nonlinear shallow water equation, this problem is very related to the vacuum problem for a compressible gas. The Green-Naghdi equation include additional nonlinear, dispersive and topography terms with a complex degenerate structure at the boundary. In particular, the degeneracy of the topography terms makes the problem loose its quasilinear structure and become fully nonlinear. Dispersive smoothing also degenerates and its behavior at the boundary can be described by an ODE with regular singularity. These issues require the development of new tools, some of which of independent interest such as the study of the mixed initial boundary value problem for dispersive perturbations of characteristic hyperbolic systems, elliptic regularization with respect to conormal derivatives, or general Hardy-type inequalities.


title: Ternary inorganic compounds containing carbon, nitrogen, and oxygen at high pressures

abstract: Ternary C_{x}N_{y}O_{z} compounds are actively researched as novel high energy density and ultrahard materials. Although some synthesis work has been performed at ambient conditions, very little is known about the high pressure chemistry of of C_{x}N_{y}O_{z} compounds. In this work, first principles variable-composition evolutionary structure prediction calculations are performed with the goal of discovering novel mixed C_{x}N_{y}O_{z} materials at ambient and high pressure conditions. By systematically searching ternary variable composition crystalline materials, the full ternary phase diagram is constructed in the range of pressures from 0 to 100 GPa. The search finds the C_{2}N_{2}O crystal containing extended covalent network of C, N, and O atoms, having space group symmetry Cmc2_{1}, and stable above just 10 GPa. Several other novel metastable (CO)_{x}-(N)_{y} crystalline compounds discovered during the search, including two polymorphs of C_{2}NO_{2} and two polymorphs of C_{3}N_{2}O_{3} crystals are found to be energetically favorable compared to polymeric carbon monoxide (CO) and nitrogen. Predicted new compounds are characterized by their Raman spectra and equations of state.


title: Self-similar profiles for homoenergetic solutions of the Boltzmann equation: particle velocity distribution and entropy

abstract: In this paper we study a class of solutions of the Boltzmann equation which have the form $f\left( x,v,t\right) =g\left( v-L\left( t\right) x,t\right) $ where $L\left( t\right) =A\left( I+tA\right) ^{-1}$ with the matrix $A$ describing a shear flow or a dilatation or a combination of both. These solutions are known as homoenergetic solutions. We prove existence of homoenergetic solutions for a large class of initial data. For different choices for the matrix $A$ and for different homogeneities of the collision kernel, we characterize the long time asymptotics of the velocity distribution for the corresponding homoenergetic solutions. For a large class of choices of $A$ we then prove rigorously the existence of self-similar solutions of the Boltzmann equation. The latter are non Maxwellian distributions and describe far-from-equilibrium flows. For Maxwell molecules we obtain exact formulas for the $H$-function for some of these flows. These formulas show that in some cases, despite being very far from equilibrium, the relationship between density, temperature and entropy is exactly the same as in the equilibrium case. We make conjectures about the asymptotics of homoenergetic solutions that do not have self-similar profiles.


title: Quantized Spectral Compressed Sensing: Cramer-Rao Bounds and Recovery Algorithms

abstract: Efficient estimation of wideband spectrum is of great importance for applications such as cognitive radio. Recently, sub-Nyquist sampling schemes based on compressed sensing have been proposed to greatly reduce the sampling rate. However, the important issue of quantization has not been fully addressed, particularly for high-resolution spectrum and parameter estimation. In this paper, we aim to recover spectrally-sparse signals and the corresponding parameters, such as frequency and amplitudes, from heavy quantizations of their noisy complex-valued random linear measurements, e.g. only the quadrant information. We first characterize the Cramer-Rao bound under Gaussian noise, which highlights the trade-off between sample complexity and bit depth under different signal-to-noise ratios for a fixed budget of bits. Next, we propose a new algorithm based on atomic norm soft thresholding for signal recovery, which is equivalent to proximal mapping of properly designed surrogate signals with respect to the atomic norm that motivates spectral sparsity. The proposed algorithm can be applied to both the single measurement vector case, as well as the multiple measurement vector case. It is shown that under the Gaussian measurement model, the spectral signals can be reconstructed accurately with high probability, as soon as the number of quantized measurements exceeds the order of K log n, where K is the level of spectral sparsity and $n$ is the signal dimension. Finally, numerical simulations are provided to validate the proposed approaches.


title: Touching proteins with virtual bare hands: how to visualize protein-drug complexes and their dynamics in virtual reality

abstract: The ability to precisely visualize the atomic geometry of the interactions between a drug and its protein target in structural models is critical in predicting the correct modifications in previously identified inhibitors to create more effective next generation drugs. It is currently common practice among medicinal chemists while attempting the above to access the information contained in three-dimensional structures by using two-dimensional projections, which can preclude disclosure of useful features. A more precise visualization of the three-dimensional configuration of the atomic geometry in the models can be achieved through the implementation of immersive virtual reality (VR). In this work, we present a freely available software pipeline for visualising protein structures through VR. New customer hardware, such as the HTC Vive and the Oculus Rift utilized in this study, are available at reasonable prices. Moreover, we have combined VR visualization with fast algorithms for simulating intramolecular motions of protein flexibility, in an effort to further improve structure-lead drug design by exposing molecular interactions that might be hidden in the less informative static models.


title: Side-Channel Inference Attacks on Mobile Keypads using Smartwatches

abstract: Smartwatches enable many novel applications and are fast gaining popularity. However, the presence of a diverse set of on-board sensors provides an additional attack surface to malicious software and services on these devices. In this paper, we investigate the feasibility of key press inference attacks on handheld numeric touchpads by using smartwatch motion sensors as a side-channel. We consider different typing scenarios, and propose multiple attack approaches to exploit the characteristics of the observed wrist movements for inferring individual key presses. Experimental evaluation using commercial off-the-shelf smartwatches and smartphones show that key press inference using smartwatch motion sensors is not only fairly accurate, but also comparable with similar attacks using smartphone motion sensors. Additionally, hand movements captured by a combination of both smartwatch and smartphone motion sensors yields better inference accuracy than either device considered individually.


title: The Effects of Protostellar Disk Turbulence on CO Emission Lines: A Comparison Study of Disks with Constant CO Abundance vs. Chemically Evolving Disks

abstract: Turbulence is the leading candidate for angular momentum transport in protoplanetary disks and therefore influences disk lifetimes and planet formation timescales. However, the turbulent properties of protoplanetary disks are poorly constrained observationally. Recent studies have found turbulent speeds smaller than what fully-developed MRI would produce (Flaherty et al. 2015, 2017). However, existing studies assumed a constant CO/H2 ratio of 0.0001 in locations where CO is not frozen-out or photo-dissociated. Our previous studies of evolving disk chemistry indicate that CO is depleted by incorporation into complex organic molecules well inside the freeze-out radius of CO. We consider the effects of this chemical depletion on measurements of turbulence. Simon et al. (2015) suggested that the ratio of the peak line flux to the flux at line center of the CO J=3-2 transition is a reasonable diagnostic of turbulence, so we focus on that metric, while adding some analysis of the more complex effects on spatial distribution. We simulate the emission lines of CO based on chemical evolution models presented in Yu et al. (2016), and find that the peak-to-trough ratio changes as a function of time as CO is destroyed. Specifically, a CO-depleted disk with high turbulent velocity mimics the peak-to-trough ratios of a non-CO-depleted disk with lower turbulent velocity. We suggest that disk observers and modelers take into account the possibility of CO depletion when using line peak-to-trough ratios to constrain the degree of turbulence in disks. Assuming that CO/H2 = 0.0001 at all disk radii can lead to underestimates of turbulent speeds in the disk by at least 0.2 km/s.


title: Dark matter and LHC: Complementarities and limitations

abstract: It is well known that dark matter density measurements, indirect and direct detection experiments, importantly complement the LHC in setting strong constraints on new physics scenarios. Yet, dark matter searches are subject to limitations which need to be considered for realistic analyses. For illustration, we explore the parameter space of the phenomenological MSSM and discuss the interplay of the constraints from dark matter searches and the LHC, and analyse the impact of the astrophysical uncertainties in some detail.


title: Measurement of the production cross-section of a single top quark in association with a $Z$ boson in proton--proton collisions at 13 TeV with the ATLAS detector

abstract: The production of a top quark in association with a $Z$ boson is investigated. The proton--proton collision data collected by the ATLAS experiment at the LHC in 2015 and 2016 at a centre-of-mass energy of $\sqrt{s} = 13 \mathrm{TeV}$ are used, corresponding to an integrated luminosity of $36.1 \mathrm{fb}^{-1}$. Events containing three identified leptons (electrons and/or muons) and two jets, one of which is identified as a $b$-quark jet are selected. The major backgrounds are diboson, $t\bar{t}$ and $Z\mathrm{+jets}$ production. A neural network is used to improve the background rejection and extract the signal. The resulting significance is $4.2\sigma$ in the data and the expected significance is $5.4\sigma$. The measured cross-section for $tZq$ production is $600 \pm 170 \mathrm{(stat.)} \pm 140 \mathrm{(syst.)} \mathrm{fb}$.


title: Reinvestigation of the excited states in the proton emitter $^{151}$Lu: particle-hole excitations across the $N=Z=64$ subshell

abstract: The excited states of the proton emitter $^{151}$Lu were reinvestigated in a recoil-decay tagging experiment at the Accelerator Laboratory of the University of Jyv\"askyl\"a (JYFL). The level scheme built on the ground state of $^{151}$Lu was updated with five new $\gamma$-ray transitions. Large-scale shell model calculations were carried out to interpret the experimental level scheme. It is found that the excitation energies of states above the $27/2^-$ and $23/2^+$ isomeric levels can be sensitive to excitations from $g_{7/2}$ and $d_{5/2}$ to single-particle orbitals above $N=Z=64$.


title: Solute transport within single pores post-pulse

abstract: We present results of studying the post-pulse transport properties of single pores. Single pores exhibit active transport (drift), and not just passive transport (diffusion) at early post-pulse times. In addition, we suggest experimental design methods obtained from model parameters that would enable experimentalists to find localized regions of a few pores.


title: Divisibility of the class numbers of imaginary quadratic fields

abstract: For a given odd integer $n>1$, we provide some families of imaginary quadratic number fields of the form $\mathbb{Q}(\sqrt{x^2-t^n})$ whose ideal class group has a subgroup isomorphic to $\mathbb{Z}/n\mathbb{Z}$.


title: Colored discrete spaces: higher dimensional combinatorial maps and quantum gravity

abstract: In any dimension $D$, the Euclidean Einstein-Hilbert action, which describes gravity in the absence of matter, can be discretized over random discrete spaces obtained by gluing families of polytopes together in all possible ways. In the physical limit of small Newton constant, only the spaces which maximize the mean curvature survive. In two dimensions, this results in a theory of random discrete spheres, which converge in the continuum limit towards the Brownian sphere, a random fractal space interpreted as a quantum random space-time. In this limit, the continuous Liouville theory of $D=2$ quantum gravity is recovered. Previous results in higher dimension regarded triangulations - gluings of tetrahedra or $D$-dimensional generalizations, leading to the continuum random tree, or gluings of simple colored building blocks of small sizes, for which multi-trace matrix model results are recovered. This work aims at providing combinatorial tools which would allow a systematic study of richer building blocks and of the spaces they generate in the continuum. We develop a bijection with stacked two-dimensional discrete surfaces, and detail how it can be used to classify discrete spaces according to their mean curvature and topology. A number of blocks are analyzed, including the new infinite family of bi-pyramids, as well as toroidal and $D$-dimensional generalizations. The relation to random tensor models is detailed. A central concern is the lowest bound on the number of ($D-2$)-cells for any given blocks, or equivalently the right scaling for the associated tensor model to have a well-behaved $1/N$ expansion. We also apply our bijection to the identification of the graphs contributing at any order to the $2n$-point functions of the colored SYK model, and to the enumeration of generalized unicellular maps - spaces obtained from a single building block - according to their mean curvature.


title: Flipping-shuttle oscillations of bright one- and two-dimensional solitons in spin-orbit-coupled Bose-Einstein condensates with Rabi mixing

abstract: We analyze a possibility of macroscopic quantum effects in the form of coupled structural oscillations and shuttle motion of bright two-component spin-orbit-coupled striped (one-dimensional, 1D) and semi-vortex (two-dimensional, 2D) matter-wave solitons, under the action of linear mixing (Rabi coupling) between the components. In 1D, the intrinsic oscillations manifest themselves as flippings between spatially even and odd components of striped solitons, while in 2D the system features periodic transitions between zero-vorticity and vortical components of semi-vortex solitons. The consideration is performed by means of a combination of analytical and numerical methods.


title: Existence of Thin Shell Wormholes using non-linear distributional geometry

abstract: We study for which polynomials $F$ a thin shell wormhole with a continuous metric (connecting two Schwarzschild spacetimes of the same mass) satisfy the null energy condition (NEC) in $F(R)$-gravity. We avoid junction conditions by using the mathematical framework of the Colombeau algebra which describes a generalized framework of the distributional geometry such that one can define multiplications between distributions generalizing the tensor product of smooth tensors. The aim for physics is to motivate a conjecture about the satisfaction of the NEC for suitable quadratic $F$ while the aim for mathematics is to derive a rigorous framework describing this situation. Here the $F(R)$-gravity should be seen as a toy model, important is that the NEC may be satisfied by some form of "microstructure" which does not arise in the classical setting and may have interesting physical meanings.


title: A categorical foundation for structured reversible flowchart languages: Soundness and adequacy

abstract: Structured reversible flowchart languages is a class of imperative reversible programming languages allowing for a simple diagrammatic representation of control flow built from a limited set of control flow structures. This class includes the reversible programming language Janus (without recursion), as well as more recently developed reversible programming languages such as R-CORE and R-WHILE. In the present paper, we develop a categorical foundation for this class of languages based on inverse categories with joins. We generalize the notion of extensivity of restriction categories to one that may be accommodated by inverse categories, and use the resulting decisions to give a reversible representation of predicates and assertions. This leads to a categorical semantics for structured reversible flowcharts, which we show to be computationally sound and adequate, as well as equationally fully abstract with respect to the operational semantics under certain conditions.


title: High-dimensional dynamics of generalization error in neural networks

abstract: We perform an average case analysis of the generalization dynamics of large neural networks trained using gradient descent. We study the practically-relevant "high-dimensional" regime where the number of free parameters in the network is on the order of or even larger than the number of examples in the dataset. Using random matrix theory and exact solutions in linear models, we derive the generalization error and training error dynamics of learning and analyze how they depend on the dimensionality of data and signal to noise ratio of the learning problem. We find that the dynamics of gradient descent learning naturally protect against overtraining and overfitting in large networks. Overtraining is worst at intermediate network sizes, when the effective number of free parameters equals the number of samples, and thus can be reduced by making a network smaller or larger. Additionally, in the high-dimensional regime, low generalization error requires starting with small initial weights. We then turn to non-linear neural networks, and show that making networks very large does not harm their generalization performance. On the contrary, it can in fact reduce overtraining, even without early stopping or regularization of any sort. We identify two novel phenomena underlying this behavior in overcomplete models: first, there is a frozen subspace of the weights in which no learning occurs under gradient descent; and second, the statistical properties of the high-dimensional regime yield better-conditioned input correlations which protect against overtraining. We demonstrate that naive application of worst-case theories such as Rademacher complexity are inaccurate in predicting the generalization performance of deep neural networks, and derive an alternative bound which incorporates the frozen subspace and conditioning effects and qualitatively matches the behavior observed in simulation.


title: Petrovskii elliptic systems in the extended Sobolev scale

abstract: Petrovskii elliptic systems of linear differential equations given on a closed smooth manifold are investigated on the extended Sobolev scale. This scale consists of all Hilbert spaces that are interpolation spaces with respect to the Hilbert Sobolev scale. Theorems on the solvability of the elliptic systems on the extended Sobolev scale are proved. An a priori estimate for solutions is obtained, and their regularity is studied.


title: Transfer, loss and physical processing of water in hit-and-run collisions of planetary embryos

abstract: Collisions between large, similar-sized bodies are believed to shape the final characteristics and composition of terrestrial planets. Their inventories of volatiles such as water, are either delivered or at least significantly modified by such events. Besides the transition from accretion to erosion with increasing impact velocity, similar-sized collisions can also result in hit-and-run outcomes for sufficiently oblique impact angles and large enough projectile-to-target mass ratios. We study volatile transfer and loss focusing on hit-and-run encounters by means of Smooth Particle Hydrodynamics simulations, including all main parameters: impact velocity, impact angle, mass ratio, and also the total colliding mass. We find a broad range of overall water losses, up to 75% in the most energetic hit-and-run events, and confirm the much more severe consequences for the smaller body also for stripping of volatile layers. Transfer of water between projectile and target inventories is found to be mostly rather inefficient, and final water contents are dominated by pre-collision inventories reduced by impact losses, for similar pre-collision water mass fractions. Comparison with our numerical results shows that current collision outcome models are not accurate enough to reliably predict these composition changes in hit-and-run events. To also account for non-mechanical losses we estimate the amount of collisionally vaporized water over a broad range of masses, and find that these contributions are particularly important in collisions of ~Mars-sized bodies, with sufficiently high impact energies, but still relatively low gravity. Our results clearly indicate that the cumulative effect of several (hit-and-run) collisions can efficiently strip protoplanets of their volatile layers, especially the smaller body, as it might be common e.g. for Earth-mass planets in systems with Super-Earths.


title: Hecke modules based on involutions in extended Weyl groups

abstract: Let X be the group of weights of a maximal torus of a simply connected semisimple group over C and let W be the Weyl group. The semidirect product W(Q\otimes X/X) is called the extended Weyl group. There is a natural C(v)-algebra H called the extended Hecke algebra with basis indexed by the extended Weyl group which contains the usual Hecke algebra as a subalgebra. We construct an H-module with basis indexed by the involutions in the extended Weyl group. This generalizes a construction of the author and D.Vogan.


title: Semistable Higgs bundles on Calabi-Yau manifolds

abstract: We provide a partial classification of semistable Higgs bundles over a simply connected Calabi-Yau manifolds. Applications to a conjecture about a special class of semistable Higgs bundles are given. In particular, the conjecture is proved for K3 and Enriques surfaces, and some related classes of surfaces.


title: Identifying Active Sites of the Water-Gas Shift Reaction over Titania Supported Platinum Catalysts under Uncertainty

abstract: A comprehensive uncertainty quantification framework has been developed for integrating computational and experimental kinetic data and to identify active sites and reaction mechanisms in catalysis. Three hypotheses regarding the active site for the water-gas shift reaction on Pt/TiO2 catalysts are tested - Pt(111), an edge interface site, and a corner interface site. Uncertainties associated with DFT calculations and model errors of microkinetic models of the active sites are informed and verified using Bayesian inference and predictive validation. Significant evidence is found for the role of the oxide support in the mechanism. Positive evidence is found in support of the edge interface active site over the corner interface site. For the edge interface site, the CO-promoted redox mechanism is found to be the dominant pathway and only at temperatures above 573 K does the classical redox mechanism contribute significantly to the overall rate. At all reaction conditions, water and surface O-H bond dissociation steps at the Pt/TiO2 interface are the main rate controlling steps.


title: Horizon Surface Gravity in Corotating Black Hole Binaries

abstract: For binary systems of corotating black holes, the zeroth law of black hole mechanics states that the surface gravity is constant over each component of the horizon. Using the approximation of a conformally flat spatial metric, we compute sequences of quasi-equilibrium initial data for corotating black hole binaries with irreducible mass ratios in the range 10:1-1:1. For orbits outside the innermost stable one, the surface gravity is found to be constant on each component of the apparent horizon at the sub-percent level. We compare those numerical results to the analytical predictions from post-Newtonian theory at the fourth (4PN) order and from black hole perturbation theory to linear order in the mass ratio. We find a remarkably good agreement for all mass ratios considered, even in the strong-field regime. In particular, our findings confirm that the domain of validity of black hole perturbative calculations appears to extend well beyond the extreme mass-ratio limit.


title: An inverse problem and a time-like Carleman estimate for parabolic integro-differential equations

abstract: In this article, We investigate an inverse problem of determining the time-dependent source factor in parabolic integro-differential equations from boundary data. We establish the uniqueness and the conditional stability estimate of H\"older type for the inverse source problem in a cylindrical shaped domain. Our methodology is based on the Bukhgeim-Klibanov method by means of the Carleman estimate. Here we also derive the "time-like" Carleman estimate for parabolic integro-differential equations.


title: An Introduction to Rocker: Docker Containers for R

abstract: We describe the Rocker project, which provides a widely-used suite of Docker images with customized R environments for particular tasks. We discuss how this suite is organized, and how these tools can increase portability, scaling, reproducibility, and convenience of R users and developers.


title: Interface conditions for a metamaterial with strong spatial dispersion

abstract: Local constitutive relations, i.e. a weak spatial dispersion, are usually considered in the effective description of metamaterials. However, they are often insufficient and effects due to a nonlocality, i.e. a strong spatial dispersion, are encountered. Recently (K.~Mnasri et al., arXiv:1705.10969), a generic form for a nonlocal constitutive relation has been introduced that could accurately describe the bulk properties of a metamaterial in terms of a dispersion relation. However, the description of functional devices made from such nonlocal metamaterials also requires the identification of suitable interface conditions. In this contribution, we derive the interface conditions for such nonlocal metamaterials.


title: Acoustic scattering theory without large-distance asymptotics

abstract: In conventional acoustic scattering theory, a large-distance asymptotic approximation is employed. In this approximation, a far-field pattern, an asymptotic approximation of the exact result, is used to describe a scattering process. The information of the distance between the target and the observer, however, is lost in the large-distance asymptotic approximation. In this paper, we provide a rigorous theory of acoustic scattering without the large-distance asymptotic approximation. The acoustic scattering treatment developed in this paper provides an improved description for the acoustic wave outside the target. Moreover, as examples, we consider acoustic scattering on a rigid sphere and on a nonrigid sphere. We also illustrate the influence of the near target effect on the angular distribution of outgoing waves. It is shown that for long wavelength acoustic scattering, the near target effect must be reckoned in.


title: Single shot ultrafast laser processing of high-aspect ratio nanochannels using elliptical Bessel beams

abstract: Ultrafast lasers have revolutionized material processing, opening a wealth of new applications in many areas of science. A recent technology that allows the cleaving of transparent materials via non-ablative processes is based on focusing and translating a high-intensity laser beam within a material to induce a well-defined internal stress plane. This then enables material separation without debris generation. Here, we use a non-diffracting beam engineered to have a transverse elliptical spatial profile to generate high aspect ratio elliptical channels in glass of dimension 350 nm x 710 nm, and subsequent cleaved surface uniformity at the sub-micron level.


title: Dynamical regimes of vortex flow in type-II superconductors with parallel twin boundaries

abstract: We explore the dynamics of driven magnetic flux lines in disordered type-II superconductors in the presence of twin boundaries oriented parallel to the direction of the applied magnetic field, using a three-dimensional elastic line model simulated with Langevin molecular dynamics. The lines are driven perpendicular to the planes to model the effect of an electric current applied parallel to the planes and perpendicular to the magnetic field. A study of the long-time non-equilibrium steady states for varying sample orientation and thickness reveals a rich collection of dynamical regimes spanning the depinning crossover region that separates the pinned and moving-lattice states of vortex matter. We observe the emergence of a preferred direction for the ordering of the Abrikosov lattice in the free-flowing vortex regime due to asymmetric pinning by the planar defects. We have performed novel direct measurements of flux line excitations such as half-loops and double kinks to aid the characterization of the topologically rich flux flow profile.


title: Computing diffusivities from particle models out of equilibrium

abstract: A new method is proposed to numerically extract the diffusivity of a (typically nonlinear) diffusion equation from underlying stochastic particle systems. The proposed strategy requires the system to be in local equilibrium and have Gaussian fluctuations but is otherwise allowed to undergo arbitrary out of equilibrium evolutions. This could be potentially relevant for particle data obtained from experimental applications. The key idea underlying the method is that finite, yet large, particle systems formally obey stochastic partial differential equations of gradient flow type satisfying a fluctuation-dissipation relation. The strategy is here applied to three classic particle models, namely independent random walkers, a zero range process and a symmetric simple exclusion process in one space dimension, to allow the comparison with analytic solutions.


title: Elastic metamaterial with simultaneously negative Mass Density, Bulk Modulus and Shear Modulus

abstract: We present a study of elastic metamaterial that possesses multiple local resonances. We demonstrated that the elastic metamaterial can have simultaneously three negative effective parameters, i.e., negative effective mass, effective bulk modulus and effective shear modulus at a certain frequency range. Through the analysis of the resonant field, it has been elucidated that the three negative parameters are induced by dipolar, monopolar and quadrupolar resonance respectively. The dipolar and monopolar resonances result into the negative band for longitudinal waves, while the dipolar and quadrupolar resonances cause the negative band for transverse waves. The two bands have an overlapping frequency regime. A simultaneously negative refraction for both longitudinal waves and transverse waves has been demonstrated in the system.


title: Microscopic quantum theory of nonlinear absorption of strong electromagnetic radiation in doped graphene

abstract: Microscopic quantum theory of nonlinear stimulated scattering of 2D Dirac particles in doped graphene on Coulomb field of impurity ions at the presence of an external strong coherent electromagnetic radiation is developed. We consider high Fermi energies and low frequencies (actually terahertz radiation) to exclude the valence electrons excitations. The Liouville-von Neumann equation for the density matrix is solved analytically, taking into account the interaction of electrons with the scattering potential in the Born approximation. With the help of this solution, the nonlinear inverse-bremsstrahlung absorption rate for a grand canonical ensemble of 2D Dirac fermions is calculated. It is shown that one can achieve the efficient absorption coefficient by this mechanism.


title: Automatic picking of P and S arrivals using a minimal uncertainty wavelet approach

abstract: The so-called "minimum uncertainty" or mu-wavelet approach has been shown to be very effective. Prior work has shown that using special algorithms based on a mu-wavelet decomposition of a seismic signal allows for arrival detection of various waves on a single-component seismogram with good precision. However, the question of interpretation of the type of wave arriving has not been addressed. The problem lies in the fact that the algorithm detects more than two arrivals even on seismograms with little noise, the peaks corresponding to P and S wave arrivals on the indicator function are not always the strongest, and the peak corresponding to the P wave arrival is not always the first peak detected. We propose to use a signal/noise relationship function in a running window as a weight coefficient of the indicator function obtained using the mu-wavelet transformation, which allows to isolate specifically those arrivals which correspond to significant changes in the signal amplitude. To interpret wave arrivals on linear registration systems (well observation systems for monitoring hydrolic fracturing, ultrasound investigations on models and samples of geological material), an algorithm was also developed to isolate arrivals that best correspond to arrivals on all registration points that have observation systems. As a result, we propose an algorithm that allows to automatically isolate arrivals of S and P waves on "frac" event recordings of a well observation system, and provides agreement with manual detection that doesn't exceed four discretes.


title: Sampling the protonation states: pH-dependent UV absorption spectrum of a polypeptide dyad

abstract: When a chromophore interacts with titrable molecular sites, the modeling of its photophysical properties requires to take into account all their possible protonation states. We have developed a multi-scale protocol, based on constant-pH molecular dynamics simulations coupled to QM/MM excitation energy calculations, aimed at sampling both the phase space and protonation state space of a short polypeptide featuring a tyrosine--tryptophan dyad interacting with two aspartic acid residues. We show that such a protocol is accurate enough to reproduce the tyrosine UV absorption spectrum at both acidic and basic pH. Moreover, it is confirmed that UV-induced radical tryptophan is reduced thanks to an electron transfer from tyrosine, ultimately explaining the complex pH-dependent behavior of the peptide spectrum.


title: High-quality trapped modes in all-dielectric metamaterials

abstract: A planar all-dielectric metamaterial made of a double-periodic lattice whose unit cell consists of a single subwavelength dielectric particle having the form of a disk possessing a penetrating hole is considered. The resonant states in the transmitted spectra of the metamaterial are identified considering modes inherent to the individual cylindrical dielectric resonator. A correlation between the asymmetry in particle's geometry, which arises from the off-centered displacement of the hole, and formation of the Mie-type and trapped modes is established. The advantages of using a coaxial-sector notch instead of a round hole for the trapped mode excitation are explained.


title: Large Negative Flattened Dispersion over the S+C+L Band Using Highly Birefringent Photonic Crystal Fiber

abstract: Novel triangular lattice photonic crystal fiber with a very large negative flattened dispersion over the S+C+L band, and very high birefringence is offered here. To investigate different optical properties of the proposed fiber, finite element method (FEM) is deployed. The fiber presents a flattened negative dispersion of $-698.5\pm 5 ps/(nm-km)$ over the wavelength of $1440 nm$ to $1600 nm$. Besides, the proposed photonic crystal fiber (PCF) exhibits a high birefringence of $1.886\times10^{-2}$ at the wavelength of $1550 nm$. Furthermore, the nonlinearity, single modeness, effective area etc. of the proposed PCF are thoroughly discussed. The fiber would have important applications in broadband residual dispersion compensation as well as polarization maintaining applications.


title: Periodic three-body orbits with vanishing angular momentum in the Jacobi-Poincare "strong" potential

abstract: Moore and Montgomery have argued that planar periodic orbits of three bodies moving in the Jacobi-Poincare, or the "strong" pairwise potential $\sum_{i>j}\frac{-1}{r_{ij}^2}$, can have all possible topologies. Here we search systematically for such orbits with vanishing angular momentum and find 24 topologically distinct orbits, 22 of which are new, in a small section of the allowed phase space, with a tendency to overcrowd, due to overlapping initial conditions. The topologies of these 24 orbits belong to three algebraic sequences defined as functions of integer $n=0,1,2, \ldots$. Each sequence extends to $n \to \infty$, but the separation of initial conditions for orbits with $n \geq 10$ becomes practically impossible with a numerical precision of 16 decimal places. Nevertheless, even with a precision of 16 decimals, it is clear that in each sequence both the orbit's initial angle $\phi_n$ and its period $T_n$ approach finite values in the asymptotic limit ($n \to \infty$). Two of three sequences are overlapping in the sense that their initial angles $\phi$ occupy the same segment on the circle and their asymptotic values $\phi_{\infty}$ are (very) close to each other. The actions of these orbits rise linearly with the index $n$ that describes the orbit's topology, which is in agreement with the Newtonian case. We show that this behaviour is consistent with the assumption of analyticity of the action as a function of period.


title: Des premiers travaux de Le Verrier \`a la d\'ecouverte de Neptune

abstract: (in english): Urbain-Jean-Joseph Le Verrier was born in Saint-L\^o on March 11, 1811. He entered the Ecole Polytechnique in 1831, from which he was to emerge 8th two years later. After first devoting himself to chemistry, in 1836 he obtained a position as an astronomy assistant at the Ecole Polytechnique. This choice will decide his future career which culminates with the discovery of Neptune in 1846. Le Verrier wrote more than 200 contributions in the Comptes Rendus de l'Acad\'emie des Sciences. These contributions are very varied: some original articles but also reports on publications published elsewhere, sometimes even simple notes of a single page. The whole set gives a very vivid vision of the development of the science of the XIXth century. At that time, the Comptes Rendus are really a reflection of the debates of the sessions of the Academy. They are published very quickly, and leave a large freedom of speech to the authors. They are therefore a snapshot of the sometimes lively polemics which animated the sessions of the Academy of Sciences. In this limited essay, we will look to the first years of the career of Le Verrier until the discovery of Neptune. (in french) Le Verrier a \'ecrit plus de 200 contributions dans les Comptes Rendus de l'Acad\'emie des Sciences . Ces contributions sont le reflet des d\'ebats des s\'eances de l'Acad\'emie. Ils sont publi\'es tr\`es rapidement, et laissent une tr\`es grande libert\'e de parole aux auteurs. On y retrouve de mani\`ere instantan\'ee les pol\'emiques parfois vives qui ont anim\'e les s\'eances de l'Acad\'emie des sciences. Dans cet essai nous survolerons les premi\`eres ann\'ees de la carri\`ere de Le Verrier jusqu'\`a la d\'ecouverte de Neptune.


title: On the wave dispersion and non-reciprocal power flow in space-time traveling acoustic metamaterials

abstract: This note analytically investigates non-reciprocal wave dispersion in locally resonant acoustic metamaterials. Dispersion relations associated with space-time varying modulations of inertial and stiffness parameters of the base material and the resonant components are derived. It is shown that the resultant dispersion bias onsets intriguing features culminating in a break-up of both acoustic and optic propagation modes and one-way local resonance band gaps. The derived band structures are validated using the full transient displacement response of a finite metamaterial. A mathematical framework is presented to characterize power flow in the modulated acoustic metamaterials to quantify energy transmission patterns associated with the non-reciprocal response. Since local resonance band gaps are size-independent and frequency tunable, the outcome enables the synthesis of a new class of sub-wavelength low-frequency one-way wave guides.


title: Optical properties of spin coated Cu doped ZnO nanocomposite films

abstract: Spin coating technique was used to synthesize pure ZnO and Cu doped ZnO films on amorphous and conducting glass substrates. The doped amount of Cu in ZnO was varied up to 5% in atomic percentage. Speed of spin coating system, coating time, initial chemical solution and annealing conditions were varied to optimize the properties of samples. Transmittance of samples was measured for ZnO doped with 1, 2, 3, 4 and 5% of Cu. Absorbance, reflectance and refractive index were derived from the measured transmittance. Film thickness of each film was calculated using the graphs of refractive index versus wavelength. Film thickness varies in a random manner depending on the amount of ZnO or Cu doped ZnO solutions spread on the substrate. The energy gap of each film was calculated using the graph of square of absorption coefficient time photon energy versus photon energy. The calculated energy gap values of Cu doped ZnO film samples decrease with the Cu concentration in ZnO. This means that the conductivity of ZnO can be increased by adding a trace amount of conducting material such as Cu.


title: Metal-Dielectric Antennas for Efficient Photon Collection from Diamond Color Centers

abstract: A central challenge in quantum technologies based on atom-like defects is the efficient collection of the emitter's fluorescence. Optical antennas are appealing as they offer directional emission together with spontaneous emission rate enhancement across a broad emitter spectrum. In this work, we introduce and optimize metal-dielectric nanoantenna designs recessed into a diamond substrate and aligned with quantum emitters. We analyze trade-offs between external quantum efficiency, collection efficiency, Purcell factor, and overall collected photon rate. This analysis shows that an optimized metal-dielectric hybrid structure can increase the collected photon rate from a nitrogen vacancy center by over two orders of magnitude compared to a bare emitter. As a result, these metal-dielectric antennas should enable single-shot electron spin measurements of NV centers at room temperature.


title: Current Saturation in Nonmetallic Field Emitters

abstract: It has been known for a long time that traditional semiconductor (e.g. intrinsic and doped Si and Ge or binary SiC and GaN) field emitters significantly deviate from Fowler-Nordheim (FN) law and saturate when a large current, on the order of microamperes or more, is attempted to be drawn from them. Many experiments established that the field emission current from carbonic materials, such as carbon nanotubes, amorphous carbon and polycrystalline diamond films, also deviate from FN law and saturate. These findings suggested that the saturation and departure from FN law is a broad and general phenomenon that applies to the class of nonmetallic field emitters. In this letter, we report a universal formula that describes the current saturation effect in nonmetallic field emitters. The formula accounts for material's bulk properties and field emitter geometry.


title: Polymer/2D material nanocomposite manufacturing beyond laboratory frontiers

abstract: Polymer nanocomposites based on 2D materials as fillers are the target in the industrial sector, but the ability to manufacture them on a large scale is very limited, and there is a lack of tools to scale up the manufacturing process of these nanocomposites. Here, for the first time, a systematic and fundamental study showing how 2D materials are inserted into the polymeric matrix in order to obtain nanocomposites using conventional and industrially scalable polymer processing machines leading to large-scale manufacturing are described. Two new strategies were used to insert pre-exfoliated 2D material into the polymer matrix, liquid-phase feeder, and solid-solid deposition. Characterizations were beyond micro and nanoscale, allowing the evaluation of the morphology for millimeter samples size. The methodologies described here are extendable to all thermoplastic polymers and 2D materials providing nanocomposites with suitable morphology to obtain singular properties and also triggering the start of the manufacturing process on a large scale.


title: Broadband, electrically tuneable, third harmonic generation in graphene

abstract: Optical harmonic generation occurs when high intensity light ($>10^{10}$W/m$^{2}$) interacts with a nonlinear material. Electrical control of the nonlinear optical response enables applications such as gate-tunable switches and frequency converters. Graphene displays exceptionally strong-light matter interaction and electrically and broadband tunable third order nonlinear susceptibility. Here we show that the third harmonic generation efficiency in graphene can be tuned by over two orders of magnitude by controlling the Fermi energy and the incident photon energy. This is due to logarithmic resonances in the imaginary part of the nonlinear conductivity arising from multi-photon transitions. Thanks to the linear dispersion of the massless Dirac fermions, ultrabroadband electrical tunability can be achieved, paving the way to electrically-tuneable broadband frequency converters for applications in optical communications and signal processing.


title: Fast and Safe: Accelerated gradient methods with optimality certificates and underestimate sequences

abstract: In this work we introduce the concept of an Underestimate Sequence (UES), which is motivated by Nesterov's estimate sequence. Our definition of a UES utilizes three sequences, one of which is a lower bound (or under-estimator) of the objective function. The question of how to construct an appropriate sequence of lower bounds is addressed, and we present lower bounds for strongly convex smooth functions and for strongly convex composite functions, which adhere to the UES framework. Further, we propose several first order methods for minimizing strongly convex functions in both the smooth and composite cases. The algorithms, based on efficiently updating lower bounds on the objective functions, have natural stopping conditions that provide the user with a certificate of optimality. Convergence of all algorithms is guaranteed through the UES framework, and we show that all presented algorithms converge linearly, with the accelerated variants enjoying the optimal linear rate of convergence.


title: Clustered vacancies in ZnO: Chemical aspects and consequences on physical properties

abstract: Chemical nature of point defects, their segregation, cluster or complex formation in ZnO is an important area of investigation. In this report, 1.2 MeV Ar ion beam is used to incorporate defects in granular ZnO. Evolution of defective state with irradiation fluence 1 x 10^14 and 1 x 10^16 ions/cm2 has been monitored using XPS, PL and Raman spectroscopic study. XPS study shows presence of oxygen vacancies (VO) in the Ar ion irradiated ZnO. Zn(LMM) Auger spectra clearly identifies transition involving metallic zinc in the irradiated samples. Intense PL emission from IZn related shallow donor bound excitons (DBX) is visible in the 10 K spectra for all samples. Although overall PL is largely reduced with irradiation disorder, DBX intensity is increased for the highest fluence irradiated sample. Raman study indicates damage in both zinc and oxygen sub-lattice by energetic ion beam. Representative Raman modes from defect complexes involving VO, IZn and IO are visible after irradiation with intermediate fluence. Further increase of fluence shows, to some extent, a homogenization of disorder. Huge reduction of resistance is also noted for this sample. Certainly, high irradiation fluence induces a qualitative modification of the conventional (and highly resistive) grain boundary (GB) structure of granular ZnO. Low resistive path, involving IZn related shallow donors, across the GB can be presumed to explain resistance reduction. Open volumes (VZn and VO) agglomerate more and more with increasing irradiation fluence and finally get transformed to voids. Results as a whole have been elucidated with a model which emphasizes possible evolution of new defect microstructure that is distinctively different from the GB related disorder. Based on the model, qualitative explanations of commonly observed radiation hardness, colouration and ferromagnetism in disordered ZnO have been put forward.


title: A Novel Asymptotic Solution to the Sommerfeld Radiation Problem: Analytic field expressions and the emergence of the Surface Waves

abstract: The well-known "Sommerfeld radiation problem" of a small -Hertzian- vertical dipole above flat lossy ground is reconsidered. The problem is examined in the spectral domain, through which it is proved to yield relatively simple integral expressions for the received Electromagnetic (EM) field. Then, using the Saddle Point method, novel analytical expressions for the scattered EM field are obtained, including sliding observation angles. As a result, a closed form solution for the subject matter is provided. Also, the necessary conditions for the emergence of the so-called Surface Wave are discussed as well. A complete mathematical formulation is presented, with detailed derivations where necessary.


title: The $\eta^\prime$ meson at the physical point with $N_f=2$ Wilson twisted mass fermions

abstract: We present results for the eta prime meson and the topological susceptibility in two flavour lattice QCD. The results are obtained using Wilson twisted mass fermions at maximal twist with pion masses ranging from 340 MeV down to the physical point. A comparison to literature values is performed giving a handle on discretisation effects.


title: Study of a New Target Design with an Additional Horn for NuMI Beam

abstract: This paper describes the simulation studies done for the standard NOvA target and the proposed new minimal NOvA target design, in forward and reverse horn current for studying the neutrino and antineutrino event yield and the background contamination at the near and far detector. The standard NOvA target is segmented into 48 graphite segments (fins) with two magnetic horns, Horn 1 placed at MCZERO and Horn 2 placed at 19.18 m from MCZERO. The new minimal target design has 24 graphite fins extended into Horn 1. It shows the study for the position optimization of the Horn 2 relative to Horn 1 for the standard and the new minimal target design. The effect of introducing an additional horn (identical to Horn 2) with the new minimal target design, on the neutrino yield and the energy spectra, has also been explored.


title: Hall effect spintronics for gas detection

abstract: We present the concept of magnetic gas detection by the Extraordinary Hall effect (EHE). The technique is compatible with the existing conductometric gas detection technologies and allows simultaneous measurement of two independent parameters: resistivity and magnetization affected by the target gas. Feasibility of the approach is demonstrated by detecting low concentration hydrogen using thin CoPd films as the sensor material. The Hall effect sensitivity of the optimized samples exceeds 240% per 104 ppm at hydrogen concentrations below 0.5% in the hydrogen/nitrogen atmosphere, which is more than two orders of magnitude higher than the sensitivity of the conductance detection.


title: A Stochastic Geometry Model of Backhaul and User Coverage in Urban UAV Networks

abstract: Wireless access points on unmanned aerial vehicles (UAVs) are being considered for mobile service provisioning in commercial networks. To be able to efficiently use these devices in cellular networks it is necessary to first have a qualitative and quantitative understanding of how their design parameters reflect on the service quality experienced by the end user. In this paper we model a network of UAVs operating at a certain height above ground to provide wireless service within coverage areas shaped by their directional antennas, with the UAVs using the existing terrestrial base station network for wireless backhaul. We provide an analytical expression for the coverage probability experienced by a typical user as a function of the UAV parameters. Using our derivations we demonstrate the existence of an optimum UAV height which maximises the end user coverage probability. We then explore a scenario where the UAVs adjust their individual heights to meet their backhaul requirements while at the same time attempting to maximise the coverage probability of the end user on the ground.


title: Representations and evaluation strategies for feasibly approximable functions

abstract: A famous result due to Ko and Friedman (1982) asserts that the problems of integration and maximisation of a univariate real function are computationally hard in a well-defined sense. Yet, both functionals are routinely computed at great speed in practice. We aim to resolve this apparent paradox by studying classes of functions which can be feasibly integrated and maximised, together with representations for these classes of functions which encode the information which is necessary to uniformly compute integral and maximum in polynomial time. The theoretical framework for this is the second-order complexity theory for operators in analysis which was introduced by Kawamura and Cook (2012). The representations we study are based on rigorous approximation by polynomials, piecewise polynomials, and rational functions. We compare these representations with respect to polytime reducibility as well as with respect to their ability to quickly evaluate symbolic expressions in a given language. We show that the representation based on rigorous approximation by piecewise polynomials is polytime equivalent to the representation based on rigorous approximation by rational functions. With this representation, all terms in a certain language, which is expressive enough to contain the maximum and integral of most functions of practical interest, can be evaluated in polynomial time. By contrast, both the representation based on polynomial approximation and the standard representation based on function evaluation, which implicitly underlies the Ko-Friedman result, require exponential time to evaluate certain terms in this language. We confirm our theoretical results by an implementation in Haskell, which provides some evidence that second-order polynomial time computability is similarly closely tied with practical feasibility as its first-order counterpart.


title: Scaling relations for parallel disc Capacitance with Dielectric media

abstract: Total capacitance of a dielectric filled circular disc configuration (radius r and separation d) has been numerically evaluated over relative permittivities ranging from 1 to 80 and aspect ratios ranging from 0.001 to 10. A new perspective is presented by which, the total capacitance, for a given disc radius r, is found to depend only on a simple scaled disc separation for the whole range of aspect ratios and relative permittivities considered. The scaling notion stems from noting that the path length has to be larger than what it would be in free space to have the same potential difference across the charged discs since the electric field strength is reduced within a homogenous dielectric medium by the relative permittivity. This scaling permits translation of practical measurements across a range of aspect ratios and relative permittivities. Further, the scaling feature allows for any computational approach that determines the total capacitance in free space, to be applicable for dielectric media having relative permittivities in the range 1 to 80. It is believed that the scaling would become better with increasing relative permittivities for practically relevant separations.


title: Motor Insurance Accidental Damage Claims Modeling with Factor Collapsing and Bayesian Model Averaging

abstract: Accidental damage is a typical component of motor insurance claim. Modeling of this nature generally involves analysis of past claim history and different characteristics of the insured objects and the policyholders. Generalized linear models (GLMs) have become the industry's standard approach for pricing and modeling risks of this nature. However, the GLM approach utilizes a single "best" model on which loss predictions are based, which ignores the uncertainty among the competing models and variable selection. An additional characteristic of motor insurance data sets is the presence of many categorical variables, within which the number of levels is high. In particular, not all levels of such variables may be statistically significant and rather some subsets of the levels may be merged to give a smaller overall number of levels for improved model parsimony and interpretability. A method is proposed for assessing the optimal manner in which to collapse a factor with many levels into one with a smaller number of levels, then Bayesian model averaging (BMA) is used to blend model predictions from all reasonable models to account for factor collapsing uncertainty. This method will be computationally intensive due to the number of factors being collapsed as well as the possibly large number of levels within factors. Hence a stochastic optimisation is proposed to quickly find the best collapsing cases across the model space.