title: Two impurities in a Bose-Einstein condensate: from Yukawa to Efimov attracted polarons

abstract: The well-known Yukawa and Efimov potentials are two different mediated interaction potentials. The first one arises in quantum field theory from the exchange of virtual particles. The second one is mediated by a real particle resonantly interacting with two other particles. This Letter shows how two impurities immersed in a Bose-Einstein condensate can exhibit both phenomena. For a weak impurity-boson attractive interactionattraction with the condensate, the two impurities form two polarons that interact through a weak Yukawa attraction mediated by virtual excitations. For a resonant attraction with the condensate, the exchanged excitation becomes a real boson and the mediated interaction changes to a strong Efimov attraction that can bind the two polarons. The resulting bipolarons turn into in-medium Efimov trimers made of the two impurities and one boson. Evidence of this physics could be seen in ultracold mixtures of atoms.


title: Spatio-Orientational Decoherence of Nanoparticles

abstract: Motivated by trapping and cooling experiments with non-spherical nanoparticles, we discuss how their combined rotational and translational quantum state is affected by the interaction with a gaseous environment. Based on the quantum master equation in terms of orientation-dependent scattering amplitudes, we evaluate the localization rate for gas collisions off an anisotropic van der Waals-type potential and for photon scattering off an anisotropic dielectric. We also show how pure angular momentum diffusion arises from these open quantum dynamics in the limit of weak anisotropies.


title: Optimizing Synchronization Stability of the Kuramoto Model in Complex Networks and Power Grids

abstract: Maintaining the stability of synchronization state is crucial for the functioning of many natural and artificial systems. In this study, we develop methods to optimize the synchronization stability of the Kuramoto model by minimizing the dominant Lyapunov exponent. With the help of the recently proposed cut-set space approximation of the steady states, we greatly simplify the objective function, and further derive its gradient and Hessian with respect to natural frequencies, which leads to an efficient algorithm with the quasi-Newton's method. The optimized systems are demonstrated to achieve better synchronization stability for the Kuramoto model with or without inertia in certain regimes. Hence our method is applicable in improving the stability of power grids. It is also viable to adjust the coupling strength of each link to improve the stability of the system. Various operational constraints can also be easily integrated into our scope by employing the interior point method in convex optimization. The properties of the optimized networks are also discussed.


title: Existence of positive solution for a system of elliptic equations via bifurcation theory

abstract: In this paper we study the existence of solution for the following class of system of elliptic equations $$ \left\{ \begin{array}{lcl} -\Delta u=\left(a-\int_{\Omega}K(x,y)f(u,v)dy\right)u+bv,\quad \mbox{in} \quad \Omega -\Delta v=\left(d-\int_{\Omega}\Gamma(x,y)g(u,v)dy\right)v+cu,\quad \mbox{in} \quad \Omega u=v=0,\quad \mbox{on} \quad \partial\Omega \end{array} \right. \eqno{(P)} $$ where $\Omega\subset\R^N$ is a smooth bounded domain, $N\geq1$, and $K,\Gamma:\Omega\times\Omega\rightarrow\R$ is a nonnegative function checking some hypotheses and $a,b,c,d\in\R$. The functions $f$ and $g$ satisfy some conditions which permit to use Bifurcation Theory to prove the existence of solution for $(P)$.


title: Radio emissions from double RHESSI TGFs

abstract: A detailed analysis of Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) terrestrial gamma ray flashes (TGFs) is performed in association with World Wide Lightning Location Network (WWLLN) sources and very low frequency (VLF) sferics recorded at Duke University. RHESSI clock offset is evaluated and found to experience changes on the 5 August 2005 and 21 October 2013, based on the analysis of TGF-WWLLN matches. The clock offsets were found for all three periods of observations with standard deviations less than 100 {\mu}s. This result opens the possibility for the precise comparative analyses of RHESSI TGFs with the other types of data (WWLLN, radio measurements, etc.) In case of multiple-peak TGFs, WWLLN detections are observed to be simultaneous with the last TGF peak for all 16 cases of multipeak RHESSI TGFs simultaneous with WWLLN sources. VLF magnetic field sferics were recorded for two of these 16 events at Duke University. These radio measurements also attribute VLF sferics to the second peak of the double TGFs, exhibiting no detectable radio emission during the first TGF peak. Possible scenarios explaining these observations are proposed. Double (multipeak) TGFs could help to distinguish between the VLF radio emission radiated by the recoil currents in the +IC leader channel and the VLF emission from the TGF producing electrons.


title: THz Discrimination of materials: demonstration of a bioinspired apparatus based on metasurfaces selective filters

abstract: We present an apparatus for terahertz fingerprint discrimination of materials designed to be fast, simple, compact and economical in order to be suitable for preliminary on-field analysis. The system working principles, bioinspired by the human vision of colors, are based on the use of microfabricated metamaterials selective filters and of a very compact optics based on metallic ellipsoidal mirrors in air. We experimentally demonstrate the operation of the apparatus in discriminating simple substances such as salt, staple foods and grease in an accurate and reproducible manner. We present the system and the obtained results and discuss issues and possible developments.


title: Projective spaces over $\mathbb{F}_{1^{\ell}}$

abstract: In this essay we study various notions of projective space (and other schemes) over $\mathbb{F}_{1^\ell}$, with $\mathbb{F}_1$ denoting the field with one element. Our leading motivation is the "Hiden Points Principle," which shows a huge deviation between the set of rational points as closed points defined over $\mathbb{F}_{1^\ell}$, and the set of rational points defined as morphisms $\texttt{Spec}(\mathbb{F}_{1^\ell}) \mapsto \mathcal{X}$. We also introduce, in the same vein as Kurokawa [13], schemes of $\mathbb{F}_{1^\ell}$-type, and consider their zeta functions.


title: Black Box Linear Algebra: Extending Wiedemann's Analysis of a Sparse Matrix Preconditioner for Computations over Small Fields

abstract: Wiedemann's paper, introducing his algorithm for sparse and structured matrix computations over arbitrary fields, also presented a pair of matrix preconditioners for computations over small fields. The analysis of the second of these is extended in order to provide more explicit statements of the expected number of nonzero entries in the matrices obtained as well as bounds on the probability that such matrices have maximal rank. This is part of ongoing work to establish that this matrix preconditioner can also be used to bound the number of nontrivial nilpotent blocks in the Jordan normal form of a preconditioned matrix, in such a way that one can also sample uniformly from the null space of the originally given matrix. If successful this will result in a black box algorithm for the type of matrix computation required when using the number field sieve for integer factorization that is provably reliable and - by a small factor - asymptotically more efficient than alternative techniques that make use of other matrix preconditioners or require computations over field extensions.


title: Multi-body Non-rigid Structure-from-Motion

abstract: Conventional structure-from-motion (SFM) research is primarily concerned with the 3D reconstruction of a single, rigidly moving object seen by a static camera, or a static and rigid scene observed by a moving camera --in both cases there are only one relative rigid motion involved. Recent progress have extended SFM to the areas of {multi-body SFM} (where there are {multiple rigid} relative motions in the scene), as well as {non-rigid SFM} (where there is a single non-rigid, deformable object or scene). Along this line of thinking, there is apparently a missing gap of "multi-body non-rigid SFM", in which the task would be to jointly reconstruct and segment multiple 3D structures of the multiple, non-rigid objects or deformable scenes from images. Such a multi-body non-rigid scenario is common in reality (e.g. two persons shaking hands, multi-person social event), and how to solve it represents a natural {next-step} in SFM research. By leveraging recent results of subspace clustering, this paper proposes, for the first time, an effective framework for multi-body NRSFM, which simultaneously reconstructs and segments each 3D trajectory into their respective low-dimensional subspace. Under our formulation, 3D trajectories for each non-rigid structure can be well approximated with a sparse affine combination of other 3D trajectories from the same structure (self-expressiveness). We solve the resultant optimization with the alternating direction method of multipliers (ADMM). We demonstrate the efficacy of the proposed framework through extensive experiments on both synthetic and real data sequences. Our method clearly outperforms other alternative methods, such as first clustering the 2D feature tracks to groups and then doing non-rigid reconstruction in each group or first conducting 3D reconstruction by using single subspace assumption and then clustering the 3D trajectories into groups.


title: Exploiting vibrational strong coupling to make an optical parametric oscillator out of a Raman laser

abstract: When the collective coupling of the rovibrational states in organic molecules and confined electromagnetic modes is sufficiently strong, the system enters into vibrational strong coupling, leading to the formation of hybrid light-matter quasiparticles. In this work we demonstrate theoretically how this hybridization in combination with stimulated Raman scattering can be utilized to widen the capabilities of Raman laser devices. We explore the conditions under which the lasing threshold can be diminished and the system can be transformed into an optical parametric oscillator. Finally, we show how the dramatic reduction of the many final molecular states into two collective excitations can be used to create an all-optical switch with output in the mid-infrared.


title: Critical-angle x-ray transmission grating spectrometer with extended bandpass and resolving power > 10,000

abstract: Several high priority subjects in astrophysics can be addressed by a state-of-the-art soft x-ray grating spectrometer (XGS). An Explorer-scale, large-area (> 1,000 cm2), high resolving power (R > 3,000) XGS is highly feasible based on Critical-Angle Transmission (CAT) gratings, even for telescopes with angular resolution of 5-10 arcsec. Significantly higher performance can be provided by a CAT XGS on an X-ray-Surveyor-type mission. CAT gratings combine the advantages of blazed reflection gratings (high efficiency, use of higher diffraction orders) with those of transmission gratings (low mass, relaxed alignment and temperature requirements, transparent at high energies) with minimal mission resource demands. They are high-efficiency blazed transmission gratings that consist of freestanding, ultra-high aspect-ratio grating bars made from SOI wafers using anisotropic dry and wet etch techniques. Blazing is achieved through reflection off grating bar sidewalls. Silicon is well matched to the soft x-ray band, and existing silicon CAT gratings exceed 30% absolute diffraction efficiency, with clear paths for improvement. CAT gratings coated with heavier elements allow extension of the CAT grating principle to higher energies and larger angles, enabling higher resolving power at shorter wavelengths. We show x-ray data from CAT gratings coated with platinum using atomic layer deposition, and demonstrate blazing to higher energies and much larger blaze angles than possible with silicon. We measure resolving power of a CAT XGS consisting of a Wolter-I focusing mirror pair from GSFC and CAT gratings, performed at the MSFC SLF. Measurement of the Al Ka doublet in 18th order shows resolving power > 10,000, based on preliminary analysis. This demonstrates that currently fabricated CAT gratings are compatible with the most advanced XGS designs for future soft x-ray spectroscopy missions.


title: Empirical study of the influence of social groups in evacuation scenarios

abstract: The effects of social groups on pedestrian dynamics, especially in evacuation scenarios, have attracted some interest recently. However, due to the lack of reliable empirical data, most of the studies focussed on modelling aspects. It was shown that social groups can have a considerable effect, e.g. on evacuation times. In order to test the model predictions we have performed laboratory experiments of evacuations with different types and sizes of the social groups. The experiments have been performed with pupils of different ages. Parameters that have been considered are (1) group size, (2) strength of intra-group interactions, and (3) composition of the groups (young adults, children, and mixtures). For all the experiments high-quality trajectories for all participants have been obtained using the PeTrack software. This allows for a detailed analysis of the group effects. One surprising observation is a decrease of the evacuation time with increasing group size.


title: Mirror and cavity formations by chains of collectively radiating atoms

abstract: We search for mirror and cavity-like features of a linear chain of atoms in which one of the atoms is specially chosen as a probe atom that is initially prepared in its excited state or is continuously driven by a laser field. Short chains are considered, composed of only three and five atoms. The analysis demonstrate the importance of the inter atomic dipole-dipole interaction which may lead to a collective ordering of the emission along some specific directions. We examine the conditions under which the radiative modes available for the emission are only those contained inside a cone centered about the inter atomic axis. Particular interest is in achieving the one-way emission along the inter atomic axis, either into left (backward) or right (forward) direction, which is referred to as a mirror-like behavior of the atomic chain. A direction dependent quantity called the directivity function, which determines how effective the system is in concentrating the radiation into a given direction, is introduced. We show that the function depends crucially on the distance between the atoms and find that there is a threshold for the inter atomic distances above which a strongly directional emission can be achieved. The one-sided emission as a manifestation of the mirror-like behavior and a highly focused emission along the inter atomic axis as a characteristic of a single-mode cavity are demonstrated to occur in the stationary field. Below the threshold the directivity function is spherically symmetric. However, we find that the population can be trapped in one of the atoms, and sometimes in all atoms indicating that at these small distances the system decays to a state for which there are no radiative modes available for emission.


title: Normalized bound states for the nonlinear Schrodinger equation in bounded domains

abstract: Given $\rho>0$, we study the elliptic problem \[ \text{find } (U,\lambda)\in H^1_0(\Omega)\times \mathbb{R} \text{ such that } \begin{cases} -\Delta U+\lambda U=|U|^{p-1}U \int_{\Omega} U^2\, dx=\rho, \end{cases} \] where $\Omega\subset\mathbb{R}^N$ is a bounded domain and $p>1$ is Sobolev-subcritical, searching for conditions (about $\rho$, $N$ and $p$) for the existence of solutions. By the Gagliardo-Nirenberg inequality it follows that, when $p$ is $L^2$-subcritical, i.e. $1<p\leq1+4/N$, the problem admits solution for every $\rho>0$. In the $L^2$-critical and supercritical case, i.e. when $1+4/N \leq p < 2^*-1$, we show that, for any $k\in\mathbb{N}$, the problem admits solutions having Morse index bounded above by $k$ only if $\rho$ is sufficiently small. Next we provide existence results for certain ranges of $\rho$, which can be estimated in terms of the Dirichlet eigenvalues of $-\Delta$ in $H^1_0(\Omega)$, extending to general domains and to changing sign solutions some results obtained in [Noris, Tavares, Verzini, Analysis & PDE, 2014] for positive solutions in the ball.


title: Yamabe type equations on graphs

abstract: Let $G=(V,E)$ be a locally finite graph, $\Omega\subset V$ be a bounded domain, $\Delta$ be the usual graph Laplacian, and $\lambda_1(\Omega)$ be the first eigenvalue of $-\Delta$ with respect to Dirichlet boundary condition. Using the mountain pass theorem due to Ambrosetti-Rabinowitz, we prove that if $\alpha<\lambda_1(\Omega)$, then for any $p>2$, there exists a positive solution to $-\Delta u-\alpha u=|u|^{p-2}u$ in $\Omega^\circ$, $u=0$ on $\partial\Omega$, where $\Omega^\circ$ and $\partial\Omega$ denote the interior and the boundary of $\Omega$ respectively. Also we consider similar problems involving the $p$-Laplacian and poly-Laplacian by the same method. Such problems can be viewed as discrete versions of the Yamabe type equations on Euclidean space or compact Riemannian manifolds.


title: Modelling high-dimensional time series efficiently by means of constrained spatio--temporal models

abstract: Many econometric analyses involve spatio--temporal data. A considerable amount of literature has addressed spatio--temporal models, with Spatial Dynamic Panel Data (SDPD) being widely investigated and applied. In real data applications, checking the validity of the theoretical assumptions underlying the SDPD models is essential but sometimes difficult. At other times, the assumptions are clearly violated. For example, the spatial matrix is assumed to be known but it may actually be unknown and needs to be estimated. In such cases, the performance of the SDPD model's estimator is generally affected. Motivated by such considerations, we propose a new model (called stationary SDPD) and a new estimation procedure based on simple and clear assumptions that can be easily checked with real data. The new model is highly adaptive, and the estimation procedure has a rate of convergence that is not affected by the dimension of the time series (under general assumptions), notwithstanding the relatively high number of parameters to be estimated. The new model may be used to represent a wide class of multivariate time series, not necessarily spatio-temporal. So, it can be used as a valid alternative to vector autoregressive (VAR) models with two immediate advantages: i) a faster rate of convergence of the estimation procedure and ii) the possibility of estimating the model even when the dimension is higher than the time series length, overcoming the curse of dimensionality typical of the VAR models. The simulation study shows that the new estimation procedure performs well compared with the classic alternative procedure, even when the spatial matrix is unknown and therefore estimated.


title: Highly Itinerant Atomic Vacancies in Phosphorene

abstract: Using detailed first-principles calculations, we investigate the hopping rate of vacancies in phosphorene, an emerging elemental 2D material besides graphene. Our work predicts that a direct observation of these mono-vacancies (MVs), showing a highly mobile and anisotropic motion, is possible only at low temperatures around 70 K or below where the thermal activity is greatly suppressed. At room temperature, the motion of a MV is sixteen orders faster than that in graphene, because of the low diffusion barrier of 0.3 eV. Built-in strain associated with the vacancies extends far along the zigzag direction while attenuating rapidly along the armchair direction. We reveal new features of the motion of di-vacancies (DVs) in phosphorene via multiple dissociation-recombination processes of vacancies owing to a small energy cost of ~ 1.05 eV for the splitting of a DV into two MVs. Furthermore, we find that uniaxial tensile strain along the zigzag direction can promote the motion of MVs, while the tensile strain along the armchair direction has the opposite effect. These itinerant features of vacancies, rooted in the unique puckering structure facilitating bond reorganization, enable phosphorene to be a bright new opportunity to broaden the knowledge of the evolution of vacancies, and a proper control of the exceedingly active and anisotropic movement of the vacancies should be critical for applications based on phosphorene.


title: Tachyon inflation in an AdS braneworld with back-reaction

abstract: We analyze the inflationary scenario based on the tachyon field coupled with the radion of the second Randall-Sundrum model (RSII). The tachyon Lagrangian is derived from the dynamics of a 3-brane moving in the five dimensional bulk. The AdS$_5$ geometry of the bulk is extended to include the radion. Using the Hamiltonian formalism we find four nonlinear field equations supplemented by the modified Friedmann equations of the RSII braneworld cosmology. After a suitable rescaling we reduce the parameters of our model to only one free parameter related to the brane tension and the AdS$_5$ curvature. We solve the equations numerically assuming a reasonably wide range of initial conditions determined by physical considerations. Varying the free parameter and initial conditions we confront our results with the Planck 2015 data.


title: Secure Analog Network Coding in Layered Networks

abstract: We consider a class of Gaussian layered networks where a source communicates with a destination through $L$ intermediate relay layers with $N$ nodes in each layer in the presence of a single eavesdropper which can overhear the transmissions of the nodes in any one layer. The problem of maximum secrecy rate achievable with analog network coding for a unicast communication over such layered wireless relay networks with directed links is considered. A relay node performing analog network coding scales and forwards the signals received at its input. The key contribution of this work is a lemma that provides the globally optimal set of scaling factors for the nodes that maximizes the end-to-end secrecy rate for a class of layered networks. We also show that in the high-SNR regime, ANC achieves secrecy rates within a constant gap of the cutset upper bound on the secrecy capacity. To the best of our knowledge, this work offers the first characterization of the performance of secure ANC in multi-layered networks in the presence of an eavesdropper.


title: Linear Quadratic Stochastic Two-Person Nonzero-Sum Differential Games: Open-Loop and Closed-Loop Nash Equilibria

abstract: In this paper, we consider a linear quadratic stochastic two-person nonzero-sum differential game. Open-loop and closed-loop Nash equilibria are introduced. The existence of the former is characterized by the solvability of a system of forward-backward stochastic differential equations, and that of the latter is characterized by the solvability of a system of coupled symmetric Riccati differential equations. Sometimes, open-loop Nash equilibria admit a closed-loop representation, via the solution to a system of non-symmetric Riccati equations, which is different from the outcome of the closed-loop Nash equilibria in general. However, it is found that for the case of zero-sum differential games, the Riccati equation system for the closed-loop representation of open-loop saddle points coincides with that for the closed-loop saddle points, which leads to the conclusion that the closed-loop representation of open-loop saddle points is the outcome of the corresponding closed-loop saddle point as long as both exist. In particular, for linear quadratic optimal control problem, the closed-loop representation of open-loop optimal controls coincides with the outcome of the corresponding closed-loop optimal strategy, provided both exist.


title: Maximizing a Monotone Submodular Function with a Bounded Curvature under a Knapsack Constraint

abstract: We consider the problem of maximizing a monotone submodular function under a knapsack constraint. We show that, for any fixed $\epsilon > 0$, there exists a polynomial-time algorithm with an approximation ratio $1-c/e-\epsilon$, where $c \in [0,1]$ is the (total) curvature of the input function. This approximation ratio is tight up to $\epsilon$ for any $c \in [0,1]$. To the best of our knowledge, this is the first result for a knapsack constraint that incorporates the curvature to obtain an approximation ratio better than $1-1/e$, which is tight for general submodular functions. As an application of our result, we present a polynomial-time algorithm for the budget allocation problem with an improved approximation ratio.


title: Equiangular tight frames and unistochastic matrices

abstract: In this work, we show that a complex equiangular tight frame (ETF) composed by $N$ vectors in dimension $d$ exists if and only if a certain bistochastic matrix, univocally determined by $N$ and $d$, belongs to a special class of unistochastic matrices. This connection allows us to find new complex ETF in infinitely many dimensions and to derive a method to introduce non-trivial free parameters in ETF. We derive a 6-parametric family of complex ETF(6,16), which defines a family of symmetric POVM. Minimal and maximal possible average entanglement of the vectors within this qubit-qutrit family are presented. Furthermore, we propose an efficient numerical procedure to find the unitary matrix underlying a unistochastic matrix, which we apply to find all existing classes of complex ETF containing up to 19 vectors.


title: Relative Cohen-Macaulay filtered modules with a view toward relative Cohen-Macaulay modules

abstract: Let R be a commutative Noetherian ring, a a proper ideal of R and M a finite R-module. It is shown that, if (R;m) is a complete local ring, then under certain conditions a contains a regular element on DR(Hc a(M)), where c = cd(a;M). A non-zerodivisor characterization of relative Cohen-Macaulay modules w.r.t a is given. We introduce the concept of relative Cohen-Macaulay filtered modules w.r.t a and study some basic properties of such modules. In paticular, we provide a non-zerodivisor characterization of relative Cohen-Macaulay filtered modules w.r.t a. Furthermore, a characterization of cohomological dimension filtration of M by the associated prime ideals of its factors is established. As a consequence, we present a cohomological dimension filtration for those modules whose zero submodule has a primary decomposition. Finally, we bring some new results about relative Cohen- Macaulay modules w.r.t a.


title: Prohorov-type local limit theorems on abstract Wiener spaces

abstract: We prove that the density of $\frac{X_1+\cdot\cdot\cdot+X_n-nE[X_1]}{\sqrt{n}}$, where $\{X_n\}_{n\geq 1}$ is a sequence of independent and identically distributed random variables taking values on an abstract Wiener space, converges in $\mathcal{L}^1$ to the density of a certain Gaussian measure which is absolutely continuous with respect to the reference Wiener measure. The crucial feature in our investigation is that we do not require the covariance structure of $\{X_n\}_{n\geq 1}$ to coincide with the one of the Wiener measure. This produces a non trivial (different from the constant function one) limiting object which reflects the different covariance structures involved. The present paper generalizes the results proved in [18] and deepens the connection between local limit theorems on (infinite dimensional) Gaussian spaces and some key tools from the Analysis on the Wiener space, like the Wiener-It\^o chaos decomposition, Ornstein-Uhlenbeck semigroup and Wick product. We also verify and discuss our main assumptions on some examples arising from the applications: dimension independent Berry-Esseen-type bounds and weak solutions of stochastic differential equations.


title: Nano strain-amplifier: making ultra-sensitive piezoresistance in nanowires possible without the need of quantum and surface charge effects

abstract: This paper presents an innovative nano strain-amplifier employed to significantly enhance the sensitivity of piezoresistive strain sensors. Inspired from the dogbone structure, the nano strain-amplifier consists of a nano thin frame released from the substrate, where nanowires were formed at the centre of the frame. Analytical and numerical results indicated that a nano strain-amplifier significantly increases the strain induced into a free standing nanowire, resulting in a large change in their electrical conductance. The proposed structure was demonstrated in p-type cubic silicon carbide nanowires fabricated using a top down process. The experimental data showed that the nano strain-amplifier can enhance the sensitivity of SiC strain sensors at least 5.4 times larger than that of the conventional structures. This result indicates the potential of the proposed strain-amplifier for ultra-sensitive mechanical sensing applications.


title: Should I stay or should I go? A latent threshold approach to large-scale mixture innovation models

abstract: This paper proposes a straightforward algorithm to carry out inference in large time-varying parameter vector autoregressions (TVP-VARs) with mixture innovation components for each coefficient in the system. We significantly decrease the computational burden by approximating the latent indicators that drive the time-variation in the coefficients with a latent threshold process that depends on the absolute size of the shocks. The merits of our approach are illustrated with two applications. First, we forecast the US term structure of interest rates and demonstrate forecast gains of the proposed mixture innovation model relative to other benchmark models. Second, we apply our approach to US macroeconomic data and find significant evidence for time-varying effects of a monetary policy tightening.


title: Infrared Interferometry and AGNs: Parsec-scale Disks and Dusty Outflows

abstract: The "torus" is the central element of the most popular theory unifying various classes of AGNs, but it is usually described as "putative" because it has not been imaged yet. Since it is too small to be resolved with single-dish telescopes, one can only make indirect assumptions about its structure using models. Using infrared interferometry, however, we were able to resolve the circum-nuclear dust distributions for several nearby AGNs and achieved constraints on some further two dozen sources. We discovered circum-nuclear dust on parsec scales in all sources and, in two nearby sources, were able to dissect this dust into two distinct components. The compact component, a very thin disk, appears to be connected to the maser disk and the extended one, which is responsible for most of the mid-IR flux, is oriented perpendicularly to the circum-nuclear gas disks. What may come as a surprise when having in mind the standard unification cartoon actually connects well to observations on larger scales. Optically thin dust in the polar region, perhaps driven by a disk wind, could solve both the scale height problem of the torus and explain the missing anisotropy in the mid-IR - X-ray relation.


title: Bulk Higgs and the 750 GeV diphoton signal

abstract: We consider scenarios of warped extra-dimensions with all matter fields in the bulk and in which both the hierarchy and the flavor puzzles of the Standard Model are addressed. The simplest extra dimensional extension of the Standard Model Higgs sector, i.e a 5D bulk Higgs doublet, can be a natural and simple explanation to the 750 GeV excess of diphotons hinted at the LHC, with the resonance responsible for the signal being the lightest CP odd excitation coming from the Higgs sector. No new matter content is invoked, the only new ingredient being the presence of (positive) brane localized kinetic terms associated to the 5D bulk Higgs, which allow to reduce the mass of the lightest CP odd Higgs excitation to 750 GeV. Production and decay of this resonance can naturally fit the observed signal when the mass scale of the rest of extradimensional resonances is of order 1 TeV.


title: Theory of spin-coherent transport through a defect spin state in a metal/ferromagnet tunnel junction during ferromagnetic resonance

abstract: We describe the coherent interaction between a defect spin at the interface of a ferromagnet and a non-magnetic material, under bias and when the magnetization of the ferromagnetic contact precesses during ferromagnetic resonance. The magnet filters charge carriers by preferentially allowing in parallel spins, which leads to a dynamic spin accumulation on the defect. Local effective fields acting on the defect spin site modify the defect spin's precession, which also modifies the charge current through the defect. This new form of current-detected spin resonance reveals the local environment of the defect, and thus can yield the defect identity.


title: Hamiltonian effective field theory study of the $\mathbf{N^*(1440)}$ resonance in lattice QCD

abstract: We examine the phase shifts and inelasticities associated with the $N^*(1440)$ Roper resonance and connect these infinite-volume observables to the finite-volume spectrum of lattice QCD using Hamiltonian effective field theory. We explore three hypotheses for the structure of the Roper resonance. All three hypotheses are able to describe the scattering data well. In the third hypothesis the Roper resonance couples the low-lying bare basis-state component associated with the ground state nucleon with the virtual meson-baryon contributions. Here the non-trivial superpositions of the meson-baryon scattering states are complemented by bare basis-state components explaining their observation in contemporary lattice QCD calculations. The merit of this scenario lies in its ability to not only describe the observed nucleon energy levels in large-volume lattice QCD simulations but also explain why other low-lying states have been missed in today's lattice QCD results for the nucleon spectrum.


title: Trajectory Optimization Through Contacts and Automatic Gait Discovery for Quadrupeds

abstract: In this work we present a trajectory Optimization framework for whole-body motion planning through contacts. We demonstrate how the proposed approach can be applied to automatically discover different gaits and dynamic motions on a quadruped robot. In contrast to most previous methods, we do not pre-specify contact switches, timings, points or gait patterns, but they are a direct outcome of the optimization. Furthermore, we optimize over the entire dynamics of the robot, which enables the optimizer to fully leverage the capabilities of the robot. To illustrate the spectrum of achievable motions, here we show eight different tasks, which would require very different control structures when solved with state-of-the-art methods. Using our trajectory Optimization approach, we are solving each task with a simple, high level cost function and without any changes in the control structure. Furthermore, we fully integrated our approach with the robot's control and estimation framework such that optimization can be run online. By demonstrating a rough manipulation task with multiple dynamic contact switches, we exemplarily show how optimized trajectories and control inputs can be directly applied to hardware.


title: An NLO+PS generator for $t\bar{t}$ and $Wt$ production and decay including non-resonant and interference effects

abstract: We present a Monte Carlo generator that implements significant theoretical improvements in the simulation of top-quark pair production and decay at the LHC. Spin correlations and off-shell effects in top-decay chains are described in terms of exact matrix elements for $p p \to \ell^+\nu_{\scriptscriptstyle\ell}\, l^-\bar{\nu}_{\scriptscriptstyle l} b \bar{b}$ at NLO QCD, where the leptons $\ell$ and $l$ belong to different families, and $b$ quarks are massive. Thus, the contributions from $t\bar{t}$ and $Wt$ single-top production as well as their quantum interference are fully included. Matrix elements are matched to the Pythia8 parton shower using a recently proposed method that allows for a consistent treatment of resonances in the POWHEG framework. These theoretical improvements are especially important for the interpretation of precision measurements of the top-quark mass, for single-top analyses in the $Wt$ channel, and for $t\bar{t}$ and $Wt$ backgrounds in the presence of jet vetoes or cuts that enhance off-shell effects. The new generator is based on a process-independent interface of the OpenLoops amplitude generator with the POWHEGBOX framework.


title: Spectroscopic investigations of divertor detachment in TCV

abstract: The aim of this work is to provide an understanding of detachment at TCV with emphasis on analysis of the Balmer line emission. A new Divertor Spectroscopy System has been developed for this purpose. Further development of Balmer line analysis techniques has allowed detailed information to be extracted from the three-body recombination contribution to the n=7 Balmer line intensity. During density ramps, the plasma at the target detaches as inferred from a drop in ion current to the target. At the same time the Balmer $6\rightarrow2$ and $7\rightarrow2$ line emission near the target is dominated by recombination. As the core density increases further, the density and recombination rate are rising all along the outer leg to the x-point while remaining highest at the target. Even at the highest core densities accessed (Greenwald fraction 0.7) the peaks in recombination and density may have moved not more than a few cm poloidally away from the target which is different to other, higher density tokamaks, where both the peak in recombination and density continue to move towards the x-point as the core density is increased. The inferred magnitude of recombination is small compared to the target ion current at the time detachment (particle flux drop) starts at the target. However, recombination may be having more localized effects (to a flux tube) which we cannot discern at this time. Later, at the highest densities achieved, the total recombination does reach levels similar to the particle flux.


title: Kazdan-Warner equation on graph

abstract: Let $G=(V,E)$ be a finite graph and $\Delta$ be the usual graph Laplacian. Using the calculus of variations and a method of upper and lower solutions, we give various conditions such that the Kazdan-Warner equation $\Delta u=c-he^u$ has a solution on $V$, where $c$ is a constant, and $h:V\rightarrow\mathbb{R}$ is a function. We also consider similar equations involving higher order derivatives on graph. Our results can be compared with the original manifold case of Kazdan-Warner (Ann. Math., 1974).


title: Electrical conductivity tensor of dense plasma in magnetic fields

abstract: Electrical conductivity of finite-temperature plasma in neutron star crusts is studied for applications in magneto-hydrodynamical description of compact stars. We solve the Boltzmann kinetic equation in relaxation time approximation taking into account the anisotropy of transport due to the magnetic field, the effects of dynamical screening in the scattering matrix element and corre- lations among the nuclei. We show that conductivity has a minimum at a non-zero temperature, a low-temperature decrease and a power-law increase with increasing temperature. Selected numerical results are shown for matter composed of carbon, iron, and heavier nuclei present in the outer crusts of neutron star.


title: Diffusions under a local strong H\"ormander condition. Part I: density estimates

abstract: We study lower and upper bounds for the density of a diffusion process in ${\mathbb{R}}^n$ in a small (but not asymptotic) time, say $\delta$. We assume that the diffusion coefficients $\sigma_1,\ldots,\sigma_d$ may degenerate at the starting time $0$ and point $x_0$ but they satisfy a strong H\"ormander condition involving the first order Lie brackets. The density estimates are written in terms of a norm which accounts for the non-isotropic structure of the problem: in a small time $\delta$, the diffusion process propagates with speed $\sqrt{\delta}$ in the direction of the diffusion vector fields $\sigma_{j}$ and with speed $\delta=\sqrt{\delta}\times \sqrt{\delta}$ in the direction of $[\sigma_{i},\sigma_{j}]$. In the second part of this paper, such estimates will be used in order to study lower and upper bounds for the probability that the diffusion process remains in a tube around a skeleton path up to a fixed time.


title: The gmwm R package: a comprehensive tool for time series analysis from state-space models to robustness

abstract: The gmwm R package for inference on time series models is mainly based on the quantity called wavelet variance which is derived from a wavelet decomposition of a time series. This quantity provides a means to summarize and graphically represent the features of time series in order to identify possible models. Moreover, it is used as a moment condition for model estimation through the generalized method of wavelet moments. Based on the latter method, this package not only provides an alternative method to estimate classical ARMA models but also delivers a general framework for the robust estimation of many time series models as well as a quick and efficient estimation of many linear state-space models.


title: Diffusions under a local strong H\"ormander condition. Part II: tube estimates

abstract: We study lower and upper bounds for the probability that a diffusion process in $\mathbb{R}^n$ remains in a tube around a skeleton path up to a fixed time. We assume that the diffusion coefficients $\sigma_1,\ldots,\sigma_d$ may degenerate but they satisfy a strong H\"ormander condition involving the first order Lie brackets around the skeleton of interest. The tube is written in terms of a norm which accounts for the non-isotropic structure of the problem: in a small time $\delta$, the diffusion process propagates with speed $\sqrt{\delta}$ in the direction of the diffusion vector fields $\sigma_{j}$ and with speed $\delta=\sqrt{\delta}\times \sqrt{\delta}$ in the direction of $[\sigma_{i},\sigma_{j}]$. The proof consists in a concatenation technique which strongly uses the lower and upper bounds for the density proved in the part I.


title: On Distance-$d$ Independent Set and other problems in graphs with few minimal separators

abstract: Fomin and Villanger (STACS 2010) proved that Maximum Independent Set, Feedback Vertex Set, and more generally the problem of finding a maximum induced subgraph of treewith at most a constant $t$, can be solved in polynomial time on graph classes with polynomially many minimal separators. We extend these results in two directions. Let $\Gpoly$ be the class of graphs with at most $\poly(n)$ minimal separators, for some polynomial $\poly$. We show that the odd powers of a graph $G$ have at most as many minimal separators as $G$. Consequently, \textsc{Distance-$d$ Independent Set}, which consists in finding maximum set of vertices at pairwise distance at least $d$, is polynomial on $\Gpoly$, for any even $d$. The problem is NP-hard on chordal graphs for any odd $d \geq 3$. We also provide polynomial algorithms for Connected Vertex Cover and Connected Feedback Vertex Set on subclasses of $\Gpoly$ including chordal and circular-arc graphs, and we discuss variants of independent domination problems.


title: Maximal Orthoplectic Fusion Frames from Mutually Unbiased Bases and Block Designs

abstract: The construction of optimal line packings in real or complex Euclidean spaces has shown to be a tantalizingly difficult task, because it includes the problem of finding maximal sets of equiangular lines. In the regime where equiangular lines are not possible, some optimal packings are known, for example, those achieving the orthoplex bound related to maximal sets of mutually unbiased bases. In this paper, we investigate the packing of subspaces instead of lines and determine the implications of maximality in this context. We leverage the existence of real or complex maximal mutually unbiased bases with a combinatorial design strategy in order to find optimal subspace packings that achieve the orthoplex bound. We also show that maximal sets of mutually unbiased bases convert between coordinate projections associated with certain balanced incomplete block designs and Grassmannian 2-designs. Examples of maximal orthoplectic fusion frames already appeared in the works by Shor, Sloane and by Zauner. They are realized in dimensions that are a power of four in the real case or a power of two in the complex case.


title: A continuous/discontinuous Galerkin solution of the shallow water equations with dynamic viscosity, high-order wetting and drying, and implicit time integration

abstract: The high-order numerical solution of the non-linear shallow water equations (and of hyperbolic systems in general) is susceptible to unphysical Gibbs oscillations that form in the proximity of strong gradients. The solution to this problem is still an active field of research as no general cure has been found yet. In this paper, we tackle this issue by presenting a dynamically adaptive viscosity based on a residual-based sub-grid scale model that has the following properties: $(i)$ it removes the spurious oscillations in the proximity of strong wave fronts while preserving the overall accuracy and sharpness of the solution. This is possible because of the residual-based definition of the dynamic diffusion coefficient. $(ii)$ For coarse grids, it prevents energy from building up at small wave-numbers. $(iii)$ The model has no tunable parameter. Our interest in the shallow water equations is tied to the simulation of coastal inundation, where a careful handling of the transition between dry and wet surfaces is particularly challenging for high-order Galerkin approximations. In this paper, we extend to a unified continuous/discontinuous Galerkin (CG/DG) framework a very simple, yet effective wetting and drying algorithm originally designed for DG [Xing, Zhang, Shu (2010)]. We show its effectiveness for problems in one and two dimensions on domains of increasing characteristic lengths varying from centimeters to kilometers. Finally, to overcome the time-step restriction incurred by the high-order Galerkin approximation, we advance the equations forward in time via a three stage, second order explicit-first-stage, singly diagonally implicit Runge-Kutta (ESDIRK) time integration scheme. Via ESDIRK, we are able to preserve numerical stability for an advective CFL number 10 times larger than its explicit counterpart.


title: Existence of positive solutions to some nonlinear equations on locally finite graphs

abstract: Let $G=(V,E)$ be a locally finite graph, whose measure $\mu(x)$ have positive lower bound, and $\Delta$ be the usual graph Laplacian. Applying the mountain-pass theorem due to Ambrosetti-Rabinowitz, we establish existence results for some nonlinear equations, namely $\Delta u+hu=f(x,u)$, $x\in V$. In particular, we prove that if $h$ and $f$ satisfy certain assumptions, then the above mentioned equation has strictly positive solutions. Also, we consider existence of positive solutions of the perturbed equation $\Delta u+hu=f(x,u)+\epsilon g$. Similar problems have been extensively studied on the Euclidean space as well as on Riemannian manifolds.


title: DiaSys: Improving SoC Insight Through On-Chip Diagnosis

abstract: To find the cause of a functional or non-functional defect (bug) in software running on a multi-processor System-on-Chip (MPSoC), developers need insight into the chip. Tracing systems provide this insight non-intrusively, at the cost of high off-chip bandwidth requirements. This I/O bottleneck limits the observability, a problem becoming more severe as more functionality is integrated on-chip. In this paper, we present DiaSys, an MPSoC diagnosis system with the potential to replace today's tracing systems. Its main idea is to partially execute the analysis of observation data on the chip; in consequence, more information and less data is sent to the attached host PC. With DiaSys, the data analysis is performed by the diagnosis application. Its input are events, which are generated by observation hardware at interesting points in the program execution (like a function call). Its outputs are events with higher information density. The event transformation is modeled as dataflow application. For execution, it is mapped in part to dedicated and distributed on-chip components, and in part to the host PC; the off-chip boundary is transparent to the developer of the diagnosis application. We implement DiaSys as extension to an existing SoC with four tiles and a mesh network running on an FPGA platform. Two usage examples confirm that DiaSys is flexible enough to replace a tracing system, while significantly lowering the off-chip bandwidth requirements. In our examples, the debugging of a race-condition bug, and the creation of a lock contention profile, we see a reduction of trace bandwidth of more than three orders of magnitude, compared to a full trace created by a common tracing system.


title: Geometry of the Fisher-Rao metric on the space of smooth densities on a compact manifold

abstract: It is known that on a closed manifold of dimension greater than one, every smooth weak Riemannian metric on the space of smooth positive densities that is invariant under the action of the diffeomorphism group, is of the form $$ G_\mu(\alpha,\beta)=C_1(\mu(M)) \int_M \frac{\alpha}{\mu}\frac{\beta}{\mu}\,\mu + C_2(\mu(M)) \int_M\alpha \cdot \int_M\beta $$ for some smooth functions $C_1,C_2$ of the total volume $\mu(M)$. Here we determine the geodesics and the curvature of this metric and study geodesic and metric completeness.


title: Approach to a Proof of the Riemann Hypothesis by the Second Mean-Value Theorem of Calculus

abstract: By the second mean-value theorem of calculus (Gauss-Bonnet theorem) we prove that the class of functions ${\mit \Xi}(z)$ with an integral representation of the form $\int_{0}^{+\infty}du\,{\mit \Omega}(u)\,{\rm ch}(uz)$ with a real-valued function ${\mit \Omega}(u) \ge 0$ which is non-increasing and decreases in infinity more rapidly than any exponential functions $\exp\left(-\lambda u\right),\,\lambda >0$ possesses zeros only on the imaginary axis. The Riemann zeta function $\zeta(s)$ as it is known can be related to an entire function $\xi(s)$ with the same non-trivial zeros as $\zeta(s)$. Then after a trivial argument displacement $s\leftrightarrow z=s-\frac{1}{2}$ we relate it to a function ${\mit \Xi}(z)$ with a representation of the form ${\mit \Xi}(z)=\int_{0}^{+\infty}du\,{\mit \Omega}(u)\,{\rm ch}(uz)$ where ${\mit \Omega}(u)$ is rapidly decreasing in infinity and satisfies all requirements necessary for the given proof of the position of its zeros on the imaginary axis $z={\rm i} y$ by the second mean-value theorem. Besides this theorem we apply the Cauchy-Riemann differential equation in an integrated operator form derived in the Appendix B. All this means that we prove a theorem for zeros of ${\mit \Xi}(z)$ on the imaginary axis $z={\rm i} y$ for a whole class of function ${\mit \Omega}(u)$ which includes in this way the proof of the Riemann hypothesis. This whole class includes, in particular, the modified Bessel functions ${\rm I}_{\nu}(z)$ for which it is known that their zeros lie on the imaginary axis and which affirms our conclusions. A class of almost-periodic functions to piece-wise constant nonincreasing functions ${\rm \Omega}(u)$ belongs also to this case. At the end we give shortly an equivalent way of a more formal description of the obtained results using the Mellin transform of functions with its variable substituted by an operator.


title: Optimal Orderings of k-subsets for Star Identification

abstract: Finding the optimal ordering of k-subsets with respect to an objective function is known to be an extremely challenging problem. In this paper we introduce a new objective for this task, rooted in the problem of star identification on spacecrafts: subsets of detected spikes are to be generated in an ordering that minimizes time to detection of a valid star constellation. We carry out an extensive analysis of the combinatorial optimization problem, and propose multiple algorithmic solutions, offering different quality-complexity trade-offs. Three main approaches are investigated: exhaustive search (branch and prune), goal-driven (greedy scene elimination, minimally intersecting subsets), and stateless algorithms which implicitly seek to satisfy the problem's goals (pattern shifting, base unrank). In practical terms, these last algorithms are found to provide satisfactory approximations to the ideal performance levels, at small computational costs.


title: Generalized Optimal Liquidation Problems Across Multiple Trading Venues

abstract: In this paper, we generalize the Almgren-Chriss's market impact model to a more realistic and flexible framework and employ it to derive and analyze some aspects of optimal liquidation problem in a security market. We illustrate how a trader's liquidation strategy alters when multiple venues and extra information are brought into the security market and detected by the trader. This study gives some new insights into the relationship between liquidation strategy and market liquidity, and provides a multi-scale approach to the optimal liquidation problem with randomly varying volatility.


title: Measurement of open heavy-flavour production as a function of charged-particle multiplicity with ALICE at the LHC

abstract: Heavy quarks are produced in the early stages of ultra-relativistic hadron collisions via hard scatterings and are an important tool for studying different aspects of Quantum Chromodynamics (QCD) in hadronic collisions. Charged-particle multiplicity gives information on the global characteristics of the event and could be used to characterize particle production mechanisms. In hadronic collisions at Large Hadron Collider (LHC) energies, there is a significant contribution of multi-parton interactions. The measurement of heavy-flavour yields as a function of charged-particle multiplicity gives insight into the mechanisms influencing their production in hadronic collisions at these energies and it is a tool to test the possible influence of multi-parton interactions. Furthermore, the charged-particle multiplicity dependence of open heavy flavours is used to test the ability of QCD theoretical models to describe the data. In ALICE, heavy-flavour production is measured via the hadronic and semi-leptonic decay channels (electrons at central rapidity and muons at forward rapidity). Charged-particle multiplicity is measured at central and forward rapidity. We will present the results on open heavy-flavour production as a function of the charged-particle multiplicity in pp and pPb collisions. Results will be compared to quarkonia measurements as well as theoretical model calculations.


title: Simultaneous models for commuting holomorphic self-maps of the ball

abstract: We prove that a finite family of commuting holomorphic self-maps of the unit ball $\mathbb{B}^q\subset \mathbb{C}^q$ admits a simultaneous holomorphic conjugacy to a family of commuting automorphisms of a possibly lower dimensional ball, and that such conjugacy satisfies a universal property. As an application we describe when a hyperbolic and a parabolic holomorphic self-map of $\mathbb{B}^q$ can commute.


title: Effect of thin film on the generation of vorticity by surface waves

abstract: Recently a theoretical scheme explaining the vorticity generation by surface waves in liquids was developed [S. Filatov et al., Phys. Rev. Lett. 116, 054501 (2016)]. Here we study how a thin (monomolecular) film presented at the surface of liquid affects the generated vorticity. We demonstrate that the vorticity becomes parametrically larger than for the case with a clean surface and now it depends on viscosity of the liquid. We also discuss the motion of particles passively advected by the generated surface flow. The results can be used in different applications: from the analysis of pollutants' diffusion on the ocean surface till the reconstruction of vorticity based on the particle image velocimetry (PIV) measurements.


title: Local Search for Max-Sum Diversification

abstract: We provide simple and fast polynomial time approximation schemes (PTASs) for several variants of the max-sum diversification problem which, in its most basic form, is as follows: Given n points p_1,...,p_n in R^d and an integer k, select k points such that the average Euclidean distance between these points is maximized. This problem commonly appears in information retrieval and web-search in order to select a diverse set of points from the input. In this context, it has recently received a lot of attention. We present new techniques to analyze natural local search algorithms. This leads to a (1-O(1/k))-approximation for distances of negative type, even subject to any matroid constraint of rank k, in time O(n k^2 log k), when assuming that distance evaluations and calls to the independence oracle are constant time. Negative type distances include as special cases Euclidean distances and many further natural distances. Our result easily transforms into a PTAS and improves on the only previously known PTAS for this setting, which relies on convex optimization techniques in an n-dimensional space and is impractical for large data sets. In contrast, our procedure has an (optimal) linear dependence on n. Using generalized exchange properties of matroid intersection, we show that a PTAS can be obtained for matroid intersection constraints as well. Moreover, our techniques, being based on local search, are conceptually simple and allow for various extensions. In particular, we get asymptotically optimal O(1)-approximations when combining the classic dispersion function with a monotone submodular objective, which is a very common class of functions to measure diversity and relevance. This result leverages recent advances on local search techniques based on proxy functions to obtain optimal approximations for monotone submodular function maximization subject to a matroid constraint.


title: Equations and Rational Points of the Modular Curves $X^+_0(p)$

abstract: Let $p$ be an odd prime number and let $X_0^+(p)$ be the quotient of the classical modular curve $X_0(p)$ by the action of the Atkin-Lehner operator $w_p$. In this paper we show how to compute explicit equations for the canonical model of $X_0^+(p)$. Then we show how to compute the modular parametrization, when it exists, from $X_0^+(p)$ to an isogeny factor $E$ of dimension 1 of its jacobian $J_0^+(p)$. Finally we show how use this map to determine the rational points on $X_0^+(p)$ up to a large fixed height.


title: Low RF-Complexity Technologies to Enable Millimeter-Wave MIMO with Large Antenna Array for 5G Wireless Communications

abstract: Millimeter-wave (mmWave) MIMO with large antenna array has attracted considerable interests from academic and industry communities, as it can provide larger bandwidth and higher spectrum efficiency. However, with hundreds of antennas, the number of radio frequency (RF) chains required by mmWave MIMO is also huge, leading to unaffordable hardware cost and power consumption in practice. In this paper, we investigate low RF-complexity technologies to solve this bottleneck. We first review the evolution of low RF-complexity technologies from microwave frequencies to mmWave frequencies. Then, we discuss two promising low RF-complexity technologies for mmWave MIMO systems in detail, i.e., phased array based hybrid precoding (PAHP) and lens array based hybrid precoding (LAHP), including their principles, advantages, challenges, and recent results. We compare the performance of these two technologies to draw some insights about how they can be deployed in practice. Finally, we conclude this paper and point out some future research directions in this area.


title: Measurement of the cosmogenic activation of germanium detectors in EDELWEISS-III

abstract: We present a measurement of the cosmogenic activation in the germanium cryogenic detectors of the EDELWEISS III direct dark matter search experiment. The decay rates measured in detectors with different exposures to cosmic rays above ground are converted into production rates of different isotopes. The measured production rates in units of nuclei/kg/day are 82 $\pm$ 21 for $^3$H, 2.8 $\pm$ 0.6 for $^{49}$V, 4.6 $\pm$ 0.7 for $^{55}$Fe, and 106 $\pm$ 13 for $^{65}$Zn. These results are the most accurate for these isotopes. A lower limit on the production rate of $^{68}$Ge of 74 nuclei/kg/day is also presented. They are compared to model predictions present in literature and to estimates calculated with the ACTIVIA code.

categories: astro-ph.CO hep-ex nucl-ex physics.ins-det


title: FOXTAIL: Modeling the nonlinear interaction between Alfv\'en eigenmodes and energetic particles in tokamaks

abstract: FOXTAIL is a new hybrid magnetohydrodynamic-kinetic code used to describe interactions between energetic particles and Alfv\'en eigenmodes in tokamaks with realistic geometries. The code simulates the nonlinear dynamics of the amplitudes of individual eigenmodes and of a set of discrete markers in five-dimensional phase space representing the energetic particle distribution. Action-angle coordinates of the equilibrium system are used for efficient tracing of energetic particles, and the particle acceleration by the wave fields of the eigenmodes is Fourier decomposed in the same angles. The eigenmodes are described using temporally constant eigenfunctions with dynamic complex amplitudes. Possible applications of the code are presented, e.g., making a quantitative validity evaluation of the one-dimensional bump-on-tail approximation of the system. Expected effects of the fulfillment of the Chirikov criterion in two-mode scenarios have also been verified.


title: Learning from the New Higgs-like Scalar before It Vanishes

abstract: Motivated by a di-photon anomaly observed by ATLAS and CMS we develop an SFitter analysis for a combined electroweak-Higgs sector, and a scalar portal at the LHC. The theoretical description is based on the linear effective Lagrangian for the Higgs and gauge fields, combined with an additional singlet scalar. The key feature is the extraction of reliable information on the portal structure of the combined scalar potential. For the specific di-photon anomaly we find that the new state might well form such a Higgs portal. To obtain more conclusive results we define and test the connection of the Wilson coefficients in the Higgs and heavy scalar sectors, as suggested by a portal setup.


title: A Priori Assessment of Prediction Confidence for Data-Driven Turbulence Modeling

abstract: Although Reynolds-Averaged Navier-Stokes (RANS) equations are still the dominant tool for engineering design and analysis applications involving turbulent flows, standard RANS models are known to be unreliable in many flows of engineering relevance, including flows with separation, strong pressure gradients or mean flow curvature. With increasing amounts of 3-dimensional experimental data and high fidelity simulation data from Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS), data-driven turbulence modeling has become a promising approach to increase the predictive capability of RANS simulations. Recently, a data-driven turbulence modeling approach via machine learning has been proposed to predict the Reynolds stress anisotropy of a given flow based on high fidelity data from closely related flows. In this work, the closeness of different flows is investigated to assess the prediction confidence a priori. Specifically, the Mahalanobis distance and the kernel density estimation (KDE) technique are used as metrics to quantify the distance between flow data sets in feature space. The flow over periodic hills at Re=10595 is used as the test set and seven flows with different configurations are individually used as training set. The results show that the prediction error of the Reynolds stress anisotropy is positively correlated with Mahalanobis distance and KDE distance, demonstrating that both extrapolation metrics can be used to estimate the prediction confidence a priori. A quantitative comparison using correlation coefficients shows that the Mahalanobis distance is less accurate in estimating the prediction confidence than KDE distance. The extrapolation metrics introduced in this work and the corresponding analysis provide an approach to aid in the choice of the data source and to assess the prediction confidence for data-driven turbulence modeling.


title: DAVE: A Unified Framework for Fast Vehicle Detection and Annotation

abstract: Vehicle detection and annotation for streaming video data with complex scenes is an interesting but challenging task for urban traffic surveillance. In this paper, we present a fast framework of Detection and Annotation for Vehicles (DAVE), which effectively combines vehicle detection and attributes annotation. DAVE consists of two convolutional neural networks (CNNs): a fast vehicle proposal network (FVPN) for vehicle-like objects extraction and an attributes learning network (ALN) aiming to verify each proposal and infer each vehicle's pose, color and type simultaneously. These two nets are jointly optimized so that abundant latent knowledge learned from the ALN can be exploited to guide FVPN training. Once the system is trained, it can achieve efficient vehicle detection and annotation for real-world traffic surveillance data. We evaluate DAVE on a new self-collected UTS dataset and the public PASCAL VOC2007 car and LISA 2010 datasets, with consistent improvements over existing algorithms.


title: Extension and boundedness of operators on Morrey spaces from extrapolation techniques and embeddings

abstract: We prove that operators satisfying the hypotheses of the extrapolation theorem for Muckenhoupt weights are bounded on weighted Morrey spaces. As a consequence, we obtain at once a number of results that have been proved individually for many operators. On the other hand, our theorems provide a variety of new results even for the unweighted case because we do not use any representation formula or pointwise bound of the operator as was assumed by previous authors. To extend the operators to Morrey spaces we show different (continuous) embeddings of (weighted) Morrey spaces into appropriate Muckenhoupt $A_1$ weighted $L_p$ spaces, which enable us to define the operators on the considered Morrey spaces by restriction. In this way, we can avoid the delicate problem of the definition of the operator, often ignored by the authors. In dealing with the extension problem through the embeddings (instead of using duality) one is neither restricted in the parameter range of the $p$'s (in particular $p=1$ is admissible and applies to weak-type inequalities) nor the operator has to be linear. Another remarkable consequence of our results is that vector-valued inequalities in Morrey spaces are automatically deduced. On the other hand, we also obtain $A_\infty$-weighted inequalities with Morrey quasinorms.


title: Spectral Echolocation via the Wave Embedding

abstract: Spectral embedding uses eigenfunctions of the discrete Laplacian on a weighted graph to obtain coordinates for an embedding of an abstract data set into Euclidean space. We propose a new pre-processing step of first using the eigenfunctions to simulate a low-frequency wave moving over the data and using both position as well as change in time of the wave to obtain a refined metric to which classical methods of dimensionality reduction can then applied. This is motivated by the behavior of waves, symmetries of the wave equation and the hunting technique of bats. It is shown to be effective in practice and also works for other partial differential equations -- the method yields improved results even for the classical heat equation.


title: Optical Characterization of the BICEP3 CMB Polarimeter at the South Pole

abstract: BICEP3 is a small-aperture refracting cosmic microwave background (CMB) telescope designed to make sensitive polarization maps in pursuit of a potential B-mode signal from inflationary gravitational waves. It is the latest in the BICEP/Keck Array series of CMB experiments at the South Pole, which has provided the most stringent constraints on inflation to date. For the 2016 observing season, BICEP3 was outfitted with a full suite of 2400 optically coupled detectors operating at 95 GHz. In these proceedings we report on the far field beam performance using calibration data taken during the 2015-2016 summer deployment season in situ with a thermal chopped source. We generate high-fidelity per-detector beam maps, show the array-averaged beam profile, and characterize the differential beam response between co-located, orthogonally polarized detectors which contributes to the leading instrumental systematic in pair differencing experiments. We find that the levels of differential pointing, beamwidth, and ellipticity are similar to or lower than those measured for BICEP2 and Keck Array. The magnitude and distribution of BICEP3's differential beam mismatch - and the level to which temperature-to-polarization leakage may be marginalized over or subtracted in analysis - will inform the design of next-generation CMB experiments with many thousands of detectors.


title: Polarized near-infrared light of the Dusty S-cluster Object (DSO/G2) at the Galactic Center

abstract: We investigate an infrared-excess source called G2 or Dusty S-cluster Object (DSO) moving on a highly eccentric orbit around the Galaxy's central black hole, Sgr A*. We use, for the first time, near-infrared polarimetric imaging data to determine the nature and the properties of the DSO, and obtain an improved K_s-band identification of this source in median polarimetry images of different observing years. The source starts to deviate from the stellar confusion in 2008, and it does not show any flux density variability over the years we analyzed it. We measure the polarization degree and angle of the DSO between 2008 and 2012 and conclude, based on the significance analysis on polarization parameters, that it is an intrinsically polarized source (>20%) with a varying polarization angle as it approaches Sgr A* position. DSO shows a near-infrared excess of K_s-L' > 3 that remains compact close to the pericenter of its orbit. Its observed parameters and the significant polarization obtained in this work show that the DSO might be a dust-enshrouded young star, forming a bow shock as it approaches the super massive black hole. The significantly high measured polarization degree indicates that it has a non-spherical geometry and it can be modelled as a combination of a bow shock with a bipolar wind of the star. We use a 3D radiative transfer model that can reproduce the observed properties of the source such as the total flux density and the polarization degree. We obtain that the change of the polarization angle can be due to an intrinsic change in the source structure. Accretion disc precession of the young star in the gravitational field of the black hole can lead to the change of the bipolar outflow and therefore the polarization angle variation. It might also be the result of the source interaction with the ambient medium.


title: The redshift evolution of massive galaxy clusters in the MACSIS simulations

abstract: We present the MAssive ClusterS and Intercluster Structures (MACSIS) project, a suite of 390 clusters simulated with baryonic physics that yields realistic massive galaxy clusters capable of matching a wide range of observed properties. MACSIS extends the recent BAHAMAS simulation to higher masses, enabling robust predictions for the redshift evolution of cluster properties and an assessment of the effect of selecting only the hottest systems. We study the observable-mass scaling relations and the X-ray luminosity-temperature relation over the complete observed cluster mass range. As expected, we find the slope of these scaling relations and the evolution of their normalization with redshift departs significantly from the self-similar predictions. However, for a sample of hot clusters with core-excised temperatures $k_{\rm{B}}T\geq5\,\rm{keV}$ the normalization and slope of the observable-mass relations and their evolution are significantly closer to self-similar. The exception is the temperature-mass relation, for which the increased importance of non-thermal pressure support and biased X-ray temperatures leads to a greater departure from self-similarity in the hottest systems. As a consequence, these also affect the slope and evolution of the normalization in the luminosity-temperature relation. The median hot gas profiles show good agreement with observational data at $z=0$ and $z=1$, with their evolution again departing significantly from the self-similar prediction. However, selecting a hot sample of clusters yields profiles that evolve significantly closer to the self-similar prediction. In conclusion, our results show that understanding the selection function is vital for robust calibration of cluster properties with mass and redshift.


title: Deterministic strong-field quantum control

abstract: Strong-field quantum-state control is investigated, taking advantage of the full---amplitude and phase---characterization of the interaction between matter and intense ultrashort pulses via transient-absorption spectroscopy. A sequence of intense delayed pulses is used, whose parameters are tailored to steer the system into a desired quantum state. We show how to experimentally enable this optimization by retrieving all quantum features of the light-matter interaction from observable spectra. This provides a full characterization of the action of strong fields on the atomic system, including the dependence upon possibly unknown pulse properties and atomic structures. Precision and robustness of the scheme are tested, in the presence of surrounding atomic levels influencing the system's dynamics.


title: Translating solitons from semi-Riemannian foliations

abstract: We recall the notion of (vertical) translating solitons in a product of a semi-Riemannian manifold $(M,g)$ and the real line. Mainly, we restrict our attention to those which are the graph of a smooth function. When dealing with submersions, we show a criteria to lift (or project) translating solitons from the base manifold to the total space (or viceversa). In particular, manifolds foliated by codimension 1 orbits of a Lie group action give rise to such solitons, up to solving a first-order ordinary differential equation. This gives us explicit criteria under which the graph of a function is a soliton, and we employ them to construct many examples of solitons, both new and old, in a unified way.


title: Moving Charged Particles in Lattice Boltzmann-Based Electrokinetics

abstract: The motion of ionic solutes and charged particles under the influence of an electric field and the ensuing hydrodynamic flow of the underlying solvent is ubiquitous in aqueous colloidal suspensions. The physics of such systems is described by a coupled set of differential equations, along with boundary conditions, collectively referred to as the electrokinetic equations. Capuani et al. [J. Chem. Phys. 121, 973 (2004)] introduced a lattice-based method for solving this system of equations, which builds upon the lattice Boltzmann algorithm for the simulation of hydrodynamic flow and exploits computational locality. However, thus far, a description of how to incorporate moving boundary conditions into the Capuani scheme has been lacking. Moving boundary conditions are needed to simulate multiple arbitrarily-moving colloids. In this paper, we detail how to introduce such a particle coupling scheme, based on an analogue to the moving boundary method for the pure LB solver. The key ingredients in our method are mass and charge conservation for the solute species and a partial-volume smoothing of the solute fluxes to minimize discretization artifacts. We demonstrate our algorithm's effectiveness by simulating the electrophoresis of charged spheres in an external field; for a single sphere we compare to the equivalent electro-osmotic (co-moving) problem. Our method's efficiency and ease of implementation should prove beneficial to future simulations of the dynamics in a wide range of complex nanoscopic and colloidal systems that was previously inaccessible to lattice-based continuum algorithms.


title: Analyzing features learned for Offline Signature Verification using Deep CNNs

abstract: Research on Offline Handwritten Signature Verification explored a large variety of handcrafted feature extractors, ranging from graphology, texture descriptors to interest points. In spite of advancements in the last decades, performance of such systems is still far from optimal when we test the systems against skilled forgeries - signature forgeries that target a particular individual. In previous research, we proposed a formulation of the problem to learn features from data (signature images) in a Writer-Independent format, using Deep Convolutional Neural Networks (CNNs), seeking to improve performance on the task. In this research, we push further the performance of such method, exploring a range of architectures, and obtaining a large improvement in state-of-the-art performance on the GPDS dataset, the largest publicly available dataset on the task. In the GPDS-160 dataset, we obtained an Equal Error Rate of 2.74%, compared to 6.97% in the best result published in literature (that used a combination of multiple classifiers). We also present a visual analysis of the feature space learned by the model, and an analysis of the errors made by the classifier. Our analysis shows that the model is very effective in separating signatures that have a different global appearance, while being particularly vulnerable to forgeries that very closely resemble genuine signatures, even if their line quality is bad, which is the case of slowly-traced forgeries.


title: Spin-Dipole Oscillation and Polarizability of a Binary Bose-Einstein Condensate near the Miscible-Immiscible Phase Transition

abstract: We report on the measurement of the spin-dipole (SD) polarizability and of the frequency of the SD oscillation of a two-component Bose-Einstein condensate of sodium atoms occupying the $|3^2S_{1/2}, F=1, m_F=\pm1\rangle$ hyperfine states. This binary spin-mixture presents the important properties of being, at the same time, fully miscible and rid of the limit set by buoyancy. It is also characterized by a huge enhancement of the SD polarizability and by the consequent softening of the frequency of the SD oscillation, due to the vicinity to the transition to the immiscible phase. The experimental data are successfully compared with the predictions of theory.


title: Electro-optic dual-comb interferometry over 40-nm bandwidth

abstract: Dual-comb interferometry is a measurement technique that uses two laser frequency combs to retrieve complex spectra in a line-by-line basis. This technique can be implemented with electro-optic frequency combs, offering intrinsic mutual coherence, high acquisition speed and flexible repetition-rate operation. A challenge with the operation of this kind of frequency comb in dual-comb interferometry is its limited optical bandwidth. Here, we use coherent spectral broadening and demonstrate electro-optic dual-comb interferometry over the entire telecommunications C band (200 lines covering ~ 40 nm, measured within 10 microseconds at 100 signal-to-noise ratio per spectral line). These results offer new prospects for electro-optic dual-comb interferometry as a suitable technology for high-speed broadband metrology, for example in optical coherence tomography or coherent Raman microscopy.


title: Neural Discourse Modeling of Conversations

abstract: Deep neural networks have shown recent promise in many language-related tasks such as the modeling of conversations. We extend RNN-based sequence to sequence models to capture the long range discourse across many turns of conversation. We perform a sensitivity analysis on how much additional context affects performance, and provide quantitative and qualitative evidence that these models are able to capture discourse relationships across multiple utterances. Our results quantifies how adding an additional RNN layer for modeling discourse improves the quality of output utterances and providing more of the previous conversation as input also improves performance. By searching the generated outputs for specific discourse markers we show how neural discourse models can exhibit increased coherence and cohesion in conversations.


title: Optical dispersion of composite particles consisting of millicharged constituents

abstract: Composite dark matter (DM) comprised of electrically charged constituents can interact with the electromagnetic field via the particle's dipole moment. This interaction results in a dispersive optical index of refraction for the DM medium. We compute this refractive index for atomic dark matter and more strongly bound systems, modeled via a harmonic oscillator potential. The dispersive nature of the index will result in a time lag between high and low energy photons simultaneously emitted from a distant astrophysical observable. This time lag, due to matter dispersion, could confound potential claims of Lorentz invariance violation (LIV) which can also result in such time lags. We compare the relative size of the two effects and determine that the dispersion due to DM is dwarfed by potential LIV effects for energies below the Planck scale.


title: Bell inequalities tailored to maximally entangled states

abstract: Bell inequalities have traditionally been used to demonstrate that quantum theory is nonlocal, in the sense that there exist correlations generated from composite quantum states that cannot be explained by means of local hidden variables. With the advent of device-independent quantum information protocols, Bell inequalities have gained an additional role as certificates of relevant quantum properties. In this work we consider the problem of designing Bell inequalities that are tailored to detect maximally entangled states. We introduce a class of Bell inequalities valid for an arbitrary number of measurements and results, derive analytically their tight classical, non-signalling and quantum bounds and prove that the latter is attained by maximally entangled states. Our inequalities can therefore find an application in device-independent protocols requiring maximally entangled states.


title: Learning from Conditional Distributions via Dual Embeddings

abstract: Many machine learning tasks, such as learning with invariance and policy evaluation in reinforcement learning, can be characterized as problems of learning from conditional distributions. In such problems, each sample $x$ itself is associated with a conditional distribution $p(z|x)$ represented by samples $\{z_i\}_{i=1}^M$, and the goal is to learn a function $f$ that links these conditional distributions to target values $y$. These learning problems become very challenging when we only have limited samples or in the extreme case only one sample from each conditional distribution. Commonly used approaches either assume that $z$ is independent of $x$, or require an overwhelmingly large samples from each conditional distribution. To address these challenges, we propose a novel approach which employs a new min-max reformulation of the learning from conditional distribution problem. With such new reformulation, we only need to deal with the joint distribution $p(z,x)$. We also design an efficient learning algorithm, Embedding-SGD, and establish theoretical sample complexity for such problems. Finally, our numerical experiments on both synthetic and real-world datasets show that the proposed approach can significantly improve over the existing algorithms.


title: Extra dimensions and violations of Lorentz symmetry

abstract: We use experimental limits on Lorentz violation to obtain new constraints on Kaluza-Klein-type theories in which the extra dimensions may be large but do not necessarily have units of length. The associated variation in fundamental quantities such as rest mass must occur slowly, on cosmological scales.


title: On the local monodromy of A-hypergeometric functions and some monodromy invariant subspaces

abstract: We obtain an explicit formula for the characteristic polynomial of the local monodromy of $A$-hypergeometric functions with respect to small loops around a coordinate hyperplane $x_i =0$. This formula is similar to the one obtained by Ando, Esterov and Takeuchi for the local monodromy at infinity. Our proof is combinatorial and can be adapted to provide an alternative proof for the latter formula as well. On the other hand, we also prove that the solution space at a nonsingular point of certain irregular and irreducible $A$--hypergeometric $D$--modules has a nontrivial global monodromy invariant subspace.


title: Global mild solutions for the nonautonomous 2D Navier-Stokes equations with impulse effects

abstract: The present paper deals with existence and uniqueness of global mild solutions for the 2D Navier-Stokes equations with impulses. Using the framework of nonautonomous dynamical systems, we extend previous results considering the 2D Navier-Stokes equations with impulse effects and allowing that the nonlinear terms are explicitly time-dependent. Additionally, we present sufficient conditions to obtain dissipativity (boundedness) for solutions starting in bounded sets.


title: A Counterexample to the Forward Recursion in Fuzzy Critical Path Analysis Under Discrete Fuzzy Sets

abstract: Fuzzy logic is an alternate approach for quantifying uncertainty relating to activity duration. The fuzzy version of the backward recursion has been shown to produce results that incorrectly amplify the level of uncertainty. However, the fuzzy version of the forward recursion has been widely proposed as an approach for determining the fuzzy set of critical path lengths. In this paper, the direct application of the extension principle leads to a proposition that must be satisfied in fuzzy critical path analysis. Using a counterexample it is demonstrated that the fuzzy forward recursion when discrete fuzzy sets are used to represent activity durations produces results that are not consistent with the theory presented. The problem is shown to be the application of the fuzzy maximum. Several methods presented in the literature are described and shown to provide results that are consistent with the extension principle.


title: Temperature coefficients of crystalline-quartz elastic constants over the cryogenic range [4 K, 15 K]

abstract: This paper brings out results of a measurement campaign aiming to determine the temperature coefficients of synthetic quartz elastic constants at liquid helium temperature. The method is based on the relationship between the resonance frequencies of a quartz acoustic cavity and the elastic constants of the material. The temperature coefficients of the elastic constants are extracted from experimental frequency-temperature data collected from a set of resonators of various cut angles, because of the anisotropy of quartz, measured on the very useful cryogenic range [4 K - 15 K]. The knowledge of these temperature coefficients would allow to further design either quartz temperature sensors or conversely frequency-temperature compensated quartz cuts. With extremely low losses, lower than $10^{-9}$ for the best ones, key applications of such devices are ultra-low loss mechanical systems used in many research areas including frequency control and fundamental measurements. The Eulerian formalism is used in this study to identify the temperature coefficients.


title: Chiral Phonons and Electrical Resistivity of Ferromagnetic Metals at Low Temperatures

abstract: Ferromagnetism is an exciting phase of matter exhibiting strongly correlated electron behavior and a standard example of spontaneously broken rotational symmetry: below the Curie temperature, atomic magnets in an isotropic single-domain ferromagnetic metal align along a spontaneously chosen direction. The scattering of conduction electrons from thermal perturbations to this spin order, together with electron-electron collisions, mark the material electrical behavior at low temperatures, where the resistivity varies mostly quadratically with the temperature. Around liquid-helium temperatures however, an interesting phenomenon occurs, giving rise to an extra \emph{linear} contribution to the variation of the electrical resistivity with temperature, whose theoretical explanation has encountered problems for a long time. Here I introduce a spin-flip scattering mechanism of conduction electrons in ferromagnetic metals arising from their interaction with the internal magnetic induction and mediated by chiral modes of the crystal lattice vibrations carrying spin 1. This mechanism is able to explain the above anomaly and give a good account of the spin-lattice relaxation times of iron, cobalt and nickel at room temperatures.


title: p-adic functionals on torsion-free abelian groups

abstract: Let p be a prime. A p-adic functional on a torsion-free abelian group G is a group homomorphism from G to the p-adic integers. The group of all such p-adic functionals is viewed as a p-adic dual group of G, and is studied from the point of view of functional analysis. An analogue of the Hahn--Banach Theorem is proved; this result shows that there are sufficiently many p-adic functionals to be interesting. There is a natural homomorphism from G to its p-adic double dual, and one main result that is proved is that the image of G in this double dual is dense in an appropriate topology. This is used to prove the second main result, which says that the p-adic double dual of G is the same as the pro-p completion of G. The theory of p-adic functionals can then be used to produce a matrix description of G if G has finite rank. This matrix description is related to another matrix description due to Malcev.


title: Spin glass behavior in a random Coulomb antiferromagnet

abstract: We study spin glass behavior in a random Ising Coulomb antiferromagnet in two and three dimensions using Monte Carlo simulations. In two dimensions, we find a transition at zero temperature with critical exponents consistent with those of the Edwards Anderson model, though with large uncertainties. In three dimensions, evidence for a finite-temperature transition, as occurs in the Edwards-Anderson model, is rather weak. This may indicate that the sizes are too small to probe the asymptotic critical behavior, or possibly that the universality class is different from that of the Edwards-Anderson model and has a lower critical dimension equal to three.


title: Active learning in pre-class assignments: Exploring the use of interactive simulations to enhance reading assignments

abstract: Pre-class reading assignments help prepare students for active classes by providing a first exposure to the terms and concepts to be used during class. We investigate if the use of inquiry-oriented PhET-based activities in conjunction with pre-class reading assignments can improve both the preparation of students for in-class learning and student attitudes towards and engagement with pre-class assignments. Over three course modules covering different topics, students were assigned randomly to complete either a textbook-only pre-class assignment or both a textbook pre-class assignment and a PhET-based activity. The assignments helped prepare students for class, as measured by performance on the pre-class quiz relative to a beginning-of-semester pre-test, but no evidence for increased learning due the PhET activity was observed. Students rated the assignments which included PhET as more enjoyable and, for the topic latest in the semester, reported engaging more with the assignments when PhET was included.


title: Automatic Environmental Sound Recognition: Performance versus Computational Cost

abstract: In the context of the Internet of Things (IoT), sound sensing applications are required to run on embedded platforms where notions of product pricing and form factor impose hard constraints on the available computing power. Whereas Automatic Environmental Sound Recognition (AESR) algorithms are most often developed with limited consideration for computational cost, this article seeks which AESR algorithm can make the most of a limited amount of computing power by comparing the sound classification performance em as a function of its computational cost. Results suggest that Deep Neural Networks yield the best ratio of sound classification accuracy across a range of computational costs, while Gaussian Mixture Models offer a reasonable accuracy at a consistently small cost, and Support Vector Machines stand between both in terms of compromise between accuracy and computational cost.


title: A Comparison of Popular Point Configurations on $\mathbb{S}^2$

abstract: There are many ways to generate a set of nodes on the sphere for use in a variety of problems in numerical analysis. We present a survey of quickly generated point sets on $\mathbb{S}^2$, examine their equidistribution properties, separation, covering, and mesh ratio constants and present a new point set, equal area icosahedral points, with low mesh ratio. We analyze numerically the leading order asymptotics for the Riesz and logarithmic potential energy for these configurations with total points $N<50,000$ and present some new conjectures.


title: Autonomous quantum machines and the finite sized Quasi-Ideal clock

abstract: Processes such as quantum computation, or the evolution of quantum cellular automata are typically described by a unitary operation implemented by an external observer. In particular, an interaction is generally turned on for a precise amount of time, using a classical clock. A fully quantum mechanical description of such a device would include a quantum description of the clock whose state is generally disturbed because of the back-reaction on it. Such a description is needed if we wish to consider finite sized autonomous quantum machines requiring no external control. The extent of the back-reaction has implications on how small the device can be, on the length of time the device can run, and is required if we want to understand what a fully quantum mechanical treatment of an observer would look like. Here, we consider the implementation of a unitary by a finite sized device which we call the "Quasi-Ideal clock", and show that the back-reaction on it can be made exponentially small in the device's dimension with only a linear increase in energy. As a result, an autonomous quantum machine need only be of modest size and or energy. We are also able to solve a long-standing open problem by using a finite sized quantum clock to approximate the continuous evolution of an Idealised clock. The result has implications on the equivalence of different paradigms of quantum thermodynamics, some which allow external control and some which only allow autonomous thermal machines.


title: Interplay of magnetic responses in all-dielectric oligomers to realize magnetic Fano resonances

abstract: We study the interplay between collective and individual optically-induced magnetic responses in quadrumers made of identical dielectric nanoparticles. Unlike their plasmonic counterparts, all-dielectric nanoparticle clusters are shown to exhibit multiple dimensions of resonant magnetic responses that can be employed for the realization of anomalous scattering signatures. We focus our analysis on symmetric quadrumers made from silicon nanoparticles and verify our theoretical results in proof-of-concept radio frequency experiments demonstrating the existence of a novel type of magnetic Fano resonance in nanophotonics.


title: Spatial Context based Angular Information Preserving Projection for Hyperspectral Image Classification

abstract: Dimensionality reduction is a crucial preprocessing for hyperspectral data analysis - finding an appropriate subspace is often required for subsequent image classification. In recent work, we proposed supervised angular information based dimensionality reduction methods to find effective subspaces. Since unlabeled data are often more readily available compared to labeled data, we propose an unsupervised projection that finds a lower dimensional subspace where local angular information is preserved. To exploit spatial information from the hyperspectral images, we further extend our unsupervised projection to incorporate spatial contextual information around each pixel in the image. Additionally, we also propose a sparse representation based classifier which is optimized to exploit spatial information during classification - we hence assert that our proposed projection is particularly suitable for classifiers where local similarity and spatial context are both important. Experimental results with two real-world hyperspectral datasets demonstrate that our proposed methods provide a robust classification performance.


title: Accurate, robust and reliable calculations of Poisson-Boltzmann solvation energies

abstract: Developing accurate solvers for the Poisson Boltzmann (PB) model is the first step to make the PB model suitable for implicit solvent simulation. Reducing the grid size influence on the performance of the solver benefits to increasing the speed of solver and providing accurate electrostatics analysis for solvated molecules. In this work, we explore the accurate coarse grid PB solver based on the Green's function treatment of the singular charges, matched interface and boundary (MIB) method for treating the geometric singularities, and posterior electrostatic potential field extension for calculating the reaction field energy. We made our previous PB software, MIBPB, robust and provides almost grid size independent reaction field energy calculation. Large amount of the numerical tests verify the grid size independence merit of the MIBPB software. The advantage of MIBPB software directly make the acceleration of the PB solver from the numerical algorithm instead of utilization of advanced computer architectures. Furthermore, the presented MIBPB software is provided as a free online sever.