title: Application of an Extended Random Phase Approximation on Giant Resonances in Light, Medium and Heavy Mass Nuclei

abstract: We present results of the time blocking approximation (TBA) on giant resonances in light, medium and heavy mass nuclei. The TBA is an extension of the widely used random-phase approximation (RPA) adding complex configurations by coupling to phonon excitations. A new method for handling the single-particle continuum is developed and applied in the present calculations. We investigate in detail the dependence of the numerical results on the size of the single particle space and the number of phonons as well as on nuclear matter properties. Our approach is self-consistent, based on an energy-density functional of Skyrme type where we used seven different parameter sets. The numerical results are compared with experimental data.


title: Bounds on LFV Higgs decays in a vector-like lepton model and searching for doubly charged leptons at the LHC

abstract: The Higgs-portal lepton flavor violation is studied in a vectorlike lepton model. For evading the constraints from rare $Z\to \ell^\pm_i \ell^\mp_j$ decays, we introduce two triplet vectorlike leptons, $(1,3)_{-1}$ and $(1,3)_{0}$. The resultant branching ratio for $h\to \mu \tau$ can be up to $10^{-4}$ when the constraints from the invisible $Z$ decays are applied. As a result, the signal strength for $\tau\tau$ channel has a $12\%$ deviation from the standard model prediction, while the muon $g-2$ is two-order of magnitude smaller than the data and $BR(\tau\to \mu \gamma)$ is of order of $10^{-12}$. A predicted doubly charged lepton in $pp$ collisions at $\sqrt{s}=13$ TeV is analyzed and it is found that the interesting production channels are $pp\to (\Psi^{--}_{1} \Psi^{++}_1, \Psi^{\pm\pm}_1 \Psi^\mp_1)$. Both single and pair production cross sections of $\Psi^{++}_1$ are comparable and can be a few hundred fb. The main decay channels for the doubly charged lepton are $\Psi^{\pm\pm} \to \ell^\pm W^\pm$ and for the heavy singly charged lepton are $\Psi^\pm_1 \to \nu W^\pm, \ell^\pm Z$. The numerical analysis is given at $13$ TeV LHC with $100$ fb$^{-1}$ luminosity.


title: Nonlinear interactions for massive spin-2 fields

abstract: We give a basic introduction to ghost-free nonlinear theories involving massive spin-2 fields, focussing on bimetric theory. After motivating the construction of such models from field theoretical considerations, we review the linear theories for massive and massless spin-2 fluctuations propagating on maximally symmetric backgrounds. The structure of general nonlinear spin-2 interactions is explained before we specialise to the ghost-free case. We review the maximally symmetric solutions of bimetric theory, its mass spectrum and the parameter limit which brings the theory close to general relativity. Finally we discuss applications of bimetric theory to cosmology with particular emphasis on the role of the general relativity limit.


title: On the significance of asperity models predictions of rough contact with respect to recent alternative theories

abstract: Recently, it has been shown that while asperity models show correctly qualitative features of rough contact problems (linearity in area-load, negative exponential dependence of load on separation which means also linearity of stiffness with load), the exact value of the coefficients are not precise for the idealized case of Gaussian distribution of heigths. This is due to the intrinsic simplifications, neglecting asperity coalescence and interaction effects. However, the issue of Gaussianity has not been proved or experimentally verified in many cases, and here we show that, for example, assuming a Weibull distribution of asperity heigths, the area-load linear coefficient is not much affected, while the relationships load-separation and therefore also stiffness-load do change largely, particularly when considering bounded distributions of asperity heigths. It is suggested that Gaussianity of surfaces should be further tested in experiments, before applying the most sophisticated rough contact models based on the Gaussian assumption.


title: Antineutrino emission and gamma background characteristics from a thermal research reactor

abstract: The detailed understanding of the antineutrino emission from research reactors is mandatory for any high sensitivity experiments either for fundamental or applied neutrino physics, as well as a good control of the gamma and neutron backgrounds induced by the reactor operation. In this article, the antineutrino emission associated to a thermal research reactor: the OSIRIS reactor located in Saclay, France, is computed in a first part. The calculation is performed with the summation method, which sums all the contributions of the beta decay branches of the fission products, coupled for the first time with a complete core model of the OSIRIS reactor core. The MCNP Utility for Reactor Evolution code was used, allowing to take into account the contributions of all beta decayers in-core. This calculation is representative of the isotopic contributions to the antineutrino flux which can be found at research reactors with a standard 19.75\% enrichment in $^{235}$U. In addition, the required off-equilibrium corrections to be applied to converted antineutrino energy spectra of uranium and plutonium isotopes are provided. In a second part, the gamma energy spectrum emitted at the core level is provided and could be used as an input in the simulation of any reactor antineutrino detector installed at such research facilities. Furthermore, a simulation of the core surrounded by the pool and the concrete shielding of the reactor has been developed in order to propagate the emitted gamma rays and neutrons from the core. The origin of these gamma rays and neutrons is discussed and the associated energy spectrum of the photons transported after the concrete walls is displayed.


title: Quantum Hall effect in momentum space

abstract: We theoretically discuss a momentum-space analog of the quantum Hall effect, which could be observed in topologically nontrivial lattice models subject to an external harmonic trapping potential. In our proposal, the Niu-Thouless-Wu formulation of the quantum Hall effect on a torus is realized in the toroidally shaped Brillouin zone. In this analogy, the position of the trap center in real space controls the magnetic fluxes that are inserted through the holes of the torus in momentum space. We illustrate the momentum-space quantum Hall effect with the noninteracting trapped Harper-Hofstadter model, for which we numerically demonstrate how this effect manifests itself in experimental observables. Extension to the interacting trapped Harper-Hofstadter model is also briefly considered. We finally discuss possible experimental platforms where our proposal for the momentum-space quantum Hall effect could be realized.


title: Microwave-induced resistance oscillations as a classical memory effect

abstract: By numerical simulations and analytical studies, we show that the phenomenon of microwave-induced resistance oscillations can be understood as a classical memory effect caused by re-collisions of electrons with scattering centers after a cyclotron period. We develop a Drude-like approach to magneto-transport in presence of a microwave field, taking account of memory effects, and find an excellent agreement between numerical and analytical results, as well as a qualitative agreement with experiment.


title: Self-consistent Description of Graphene Quantum Amplifier

abstract: High level of dissipation in normal metals makes challenging development of active and passive plasmonic devices. One possible solution to this problem is to use alternative materials. Graphene is a good candidate for plasmonics in near infrared (IR) region. In this paper we develop quantum theory of a graphene plasmon generator. We account for the first time quantum correlations and dissipation effects that allows describing such regimes of quantum plasmonic amplifier as surface plasmon emitting diode and surface plasmon amplifier by stimulated emission of radiation. Switching between these generation types is possible in situ with variance of graphene Fermi-level or gain transition frequency. We provide explicit expressions for dissipation and interaction constants through material parameters and find the generation spectrum and correlation function of second order which predicts laser statistics.


title: Detecting the cosmic web with radio surveys

abstract: We study the challenges to detect the cosmic web at radio wavelengths with state-of-the-art cosmological simulations of extragalactic magnetic fields. The incoming generation of radio surveys operating at low frequency, like LOFAR, SKA-LOW and MWA will have the best chance to detect the large-scale, low surface brightness emission from the shocked cosmic web. The detected radio emission will enable to constrain the average magnetisation level of the gas in filaments and the acceleration efficiency of electrons by strong shocks. In case of detections, through statistical modelling (e.g. correlation functions) it will be possible to discriminate among competing scenarios for the magnetisation of large-scale structures (i.e. astrophysical versus primordial scenarios), making radio surveys an important probe of cosmic magnetogenesis.


title: Differentiation of the Cholesky decomposition

abstract: We review strategies for differentiating matrix-based computations, and derive symbolic and algorithmic update rules for differentiating expressions containing the Cholesky decomposition. We recommend new `blocked' algorithms, based on differentiating the Cholesky algorithm DPOTRF in the LAPACK library, which uses `Level 3' matrix-matrix operations from BLAS, and so is cache-friendly and easy to parallelize. For large matrices, the resulting algorithms are the fastest way to compute Cholesky derivatives, and are an order of magnitude faster than the algorithms in common usage. In some computing environments, symbolically-derived updates are faster for small matrices than those based on differentiating Cholesky algorithms. The symbolic and algorithmic approaches can be combined to get the best of both worlds.


title: Trilateration-based reconstruction of ortho-positronium decays into three photons with the J-PET detector

abstract: This work reports on a new reconstruction algorithm allowing to reconstruct the decays of ortho-positronium atoms into three photons using the places and times of photons recorded in the detector. The method is based on trilateration and allows for a simultaneous reconstruction of both location and time of the decay. Results of resolution tests of the new reconstruction in the J-PET detector based on Monte Carlo simulations are presented, which yield a spatial resolution at the level of 2 cm (FWHM) for X and Y and at the level of 1 cm (FWHM) for Z available with the present resolution of J-PET after application of a kinematic fit. Prospects of employment of this method for studying angular correlations of photons in decays of polarized ortho-positronia for the needs of tests of CP and CPT discrete symmetries are also discussed. The new reconstruction method allows for discrimination of background from random three-photon coincidences as well as for application of a novel method for determination of the linear polarization of ortho-positronium atoms, which is also introduced in this work.


title: A distribution of energies of a charged black hole

abstract: We proposed a gravitational energy momentum tensor in previous paper. By using it, we evaluate an electric and a gravitational energies in and out of an event horizon of a charged black hole. There is no gravitational energy in the event horizon as well as Schwarzschild black hole.


title: Unwinding of the twist-bend nematic structure under an external field

abstract: The recently discovered twist-bend nematic phase, $\rm N_{tb}$, is a non-uniform equilibrium nematic phase that presents a spontaneous bend with a precession of the nematic director, $\textbf{n}$, on a conical helix with a tilt angle $\theta$ and helical pitch $P$. The stability of the $\rm N_{tb}$ phase has been recently demonstrated from the elastic point of view by extending the Frank elastic energy density of the nematic phase to include the symmetry element of helical axis, $\bf t$. In the present communication, we investigate the influence of an external bulk field (magnetic or electric) on the pitch of the $\rm N_{tb}$ phase. For a material with positive magnetic susceptibility anisotropy, when the magnetic field $\bf{H}$ is parallel to the helical axis, $\bf{t}$, the field modifies the coupling parameter between the director and the helical axis, thus shifting the interval of values for which this coupling results in the $\rm N_{tb}$ phase. If $\bf{H}$ is perpendicular to $\bf{t}$, it is shown that the pitch of the nematic director $\bf n$ precession increases with $\bf{H}$, and at a critical value of $H_c$ it becomes infinite as in a cholesteric liquid crystal, that is, one gets a uniform orientation of n. Our results include the behaviour of a cholesteric under field as a limiting case.


title: Is a complete, reduced set necessarily of constant width?

abstract: Is it true that a convex body $K$ being complete and reduced with respect to some gauge body $C$ is necessarily of constant width, that is, satisfies $K-K=\rho(C-C)$ for some $\rho>0$? We prove this implication for several cases including the following: if $K$ is a simplex and or if $K$ possesses a smooth extreme point, then the implication holds. Moreover, we derive several new results on perfect norms.


title: A Note on Kernels, Images, and Cokernels in the Perverse Category

abstract: We discuss the relationship between kernels, images and cokernels of morphisms between perverse sheaves and induced maps on stalk cohomology.


title: 5G 3GPP-like Channel Models for Outdoor Urban Microcellular and Macrocellular Environments

abstract: For the development of new 5G systems to operate in bands up to 100 GHz, there is a need for accurate radio propagation models at these bands that currently are not addressed by existing channel models developed for bands below 6 GHz. This document presents a preliminary overview of 5G channel models for bands up to 100 GHz. These have been derived based on extensive measurement and ray tracing results across a multitude of frequencies from 6 GHz to 100 GHz, and this document describes an initial 3D channel model which includes: 1) typical deployment scenarios for urban microcells (UMi) and urban macrocells (UMa), and 2) a baseline model for incorporating path loss, shadow fading, line of sight probability, penetration and blockage models for the typical scenarios. Various processing methodologies such as clustering and antenna decoupling algorithms are also presented.


title: Zitterbewegung, internal momentum and spin of the circular travelling wave electromagnetic electron

abstract: The study of this paper demonstrates that electron has Dirac delta like internal momentum (u,p_{{\theta}}), going round in a circle of radius equal to half the reduced Compton wavelength of electron with tangential velocity c. The circular momentum p_{{\theta}} and energy u emanate from circular Dirac delta type rotating monochromatic electromagnetic (EM) wave that itself travels in another circle having radius equal to the reduced Compton wavelength of electron. The phenomenon of Zitterbewegung and the spin of electron are the natural consequences of the model. The spin is associated with the internal circulating momentum of electron in terms of four component spinor, which leads to the Dirac equation linking the EM electron model with quantum mechanical theory. Our model accurately explains the experimental results of electron channelling experiment, [P. Catillon et al., Found.Phys. 38, 659 (2008)], in which the momentum resonance is observed at 161.784MeV/c corresponding to Zitterbewegung frequency of 80.874MeV/c electron beam.


title: SHAPE: Linear-Time Camera Pose Estimation With Quadratic Error-Decay

abstract: We propose a novel camera pose estimation or perspective-n-point (PnP) algorithm, based on the idea of consistency regions and half-space intersections. Our algorithm has linear time-complexity and a squared reconstruction error that decreases at least quadratically, as the number of feature point correspondences increase. Inspired by ideas from triangulation and frame quantisation theory, we define consistent reconstruction and then present SHAPE, our proposed consistent pose estimation algorithm. We compare this algorithm with state-of-the-art pose estimation techniques in terms of accuracy and error decay rate. The experimental results verify our hypothesis on the optimal worst-case quadratic decay and demonstrate its promising performance compared to other approaches.


title: Towards Security of Additive Layer Manufacturing

abstract: Additive Layer Manufacturing (ALM), also broadly known as 3D printing, is a new technology to produce 3D objects. As an opposite approach to the conventional subtractive manufacturing process, 3D objects are created by adding thin material layers over layers. Until recently, they have been used, mainly, for plastic models. However, the technology has evolved making it possible to use high-quality printing with metal alloys. Agencies and companies like NASA, ESA, Boeing, Airbus, etc. are investigating various ALM technology application areas. Recently, SpaceX used additive manufacturing to produce engine chambers for the newest Dragon spacecraft. BAE System plans to print on-demand a complete Unmanned Aerial Vehicle (UAV), depending on the operational requirements. Companies expect the implementation of ALM technology will bring a broad variety of technological and economic benefits. This includes, but not limited to, the reduction of the time needed to produce complex parts, reduction of wasted material and thus control of production costs along with minimization of part storage space as companies implement just-in-time and on-demand production solutions. The broad variety of application areas and a high grade of computerization of the manufacturing process will inevitably make ALM an attractive target for various attacks.


title: The local metric dimension of the lexicographic product of graphs

abstract: The metric dimension is quite a well-studied graph parameter. Recently, the adjacency dimension and the local metric dimension have been introduced and studied. In this paper, we give a general formula for the local metric dimension of the lexicographic product $G \circ \mathcal{H}$ of a connected graph $G$ of order $n$ and a family $\mathcal{H}$ composed by $n$ graphs. We show that the local metric dimension of $G \circ \mathcal{H}$ can be expressed in terms of the true twin equivalence classes of $G$ and the local adjacency dimension of the graphs in $\mathcal{H}$.


title: Bipolar Neutrosophic Soft Expert Sets

abstract: In this paper, we define concept of bipolar neutrosophic soft expert set. Also, its basic operations, namely complement, union and intersection. We give examples for these concepts.


title: A micromechanical proof-of-principle experiment for measuring the gravitational force of milligram masses

abstract: This paper addresses a simple question: how small can one make a gravitational source mass and still detect its gravitational coupling to a nearby test mass? We describe an experimental scheme based on micromechanical sensing to observe gravity between milligram-scale source masses, thereby improving the current smallest source mass values by three orders of magnitude and possibly even more. We also discuss the implications of such measurements both for improved precision measurements of Newton's constant and for a new generation of experiments at the interface between quantum physics and gravity.


title: Testing the tetraquark structure for the $X$ resonances in low-lying region

abstract: Assuming four-quark structure for the $X$ resonances in low-lying region, we calculate their masses using the color-spin interaction. In specific, the hyperfine masses of the color-spin interaction are calculated for the possible states in spin-0, spin-1, spin-2 channels. The two states in spin-0 channel as well as the two states in spin-1 channel are diagonalized in order to generate the physical hyperfine masses. By matching the difference in hyperfine masses with the splitting in corresponding hadron masses and using the $X(3872)$ mass as an input, we estimate the masses corresponding to the states $J^{PC}=0^{++}, 1^{+-},2^{++}$. We find the masses of two states in $1^{+-}$ are close to those of $X(3823)$, $X(3900)$, and the mass of the $2^{++}$ state is close to that of $X(3940)$. For them, the discrepancies are about $\sim 10$ MeV. This may suggest that the quantum numbers of the controversial states are $X(3823)=1^{+-}, X(3900)=1^{+-}, X(3940)=2^{++}$. In this work, we use the same inputs parameters, the constituent quark masses and the strength of the color-spin interaction, that have been adopted in the previous work on the $D$ or $B$-meson excited states. There, it was shown that the four-quark structure can be manifested in their excited states. Thus, our results in this work provide a consistent treatment on open- and hidden-charm mesons as far as the four-quark model is concerned.


title: Partial Category Actions on Sets and Topological Spaces

abstract: We introduce (continuous) partial category actions on sets (topological spaces) and show that each such action admits a universal globalization. Thereby, we obtain a simultaneous generalization of corresponding results for groups, by Kellendonk and Lawson, and for monoids, by Megrelishvili and Schroder. We apply this result to the special case of partial groupoid actions where we obtain a sharpening of a result by Gilbert, concerning ordered groupoids, in the sense that mediating functions between universal globalizations always are injective.


title: On the Accuracy of Point Localisation in a Circular Camera-Array

abstract: Although many advances have been made in light-field and camera-array image processing, there is still a lack of thorough analysis of the localisation accuracy of different multi-camera systems. By considering the problem from a frame-quantisation perspective, we are able to quantify the point localisation error of circular camera configurations. Specifically, we obtain closed form expressions bounding the localisation error in terms of the parameters describing the acquisition setup. These theoretical results are independent of the localisation algorithm and thus provide fundamental limits on performance. Furthermore, the new frame-quantisation perspective is general enough to be extended to more complex camera configurations.


title: Measurement of the $B_{s}^{0} \rightarrow D_{s}^{(*)+}D_{s}^{(*)-}$ branching fractions

abstract: The branching fraction of the decay $B_{s}^{0} \rightarrow D_{s}^{(*)+}D_{s}^{(*)-}$ is measured using $pp$ collision data corresponding to an integrated luminosity of $1.0fb^{-1}$, collected using the LHCb detector at a centre-of-mass energy of $7$TeV. It is found to be \begin{align*} {\mathcal{B}}(B_{s}^{0}\rightarrow~D_{s}^{(*)+}D_{s}^{(*)-}) = (3.05 \pm 0.10 \pm 0.20 \pm 0.34)\%, \end{align*} where the uncertainties are statistical, systematic, and due to the normalisation channel, respectively. The branching fractions of the individual decays corresponding to the presence of one or two $D^{*\pm}_{s}$ are also measured. The individual branching fractions are found to be \begin{align*} {\mathcal{B}}(B_{s}^{0}\rightarrow~D_{s}^{*\pm}D_{s}^{\mp}) = (1.35 \pm 0.06 \pm 0.09 \pm 0.15)\%, \newline{\mathcal{B}}(B_{s}^{0}\rightarrow~D_{s}^{*+}D_{s}^{*-}) = (1.27 \pm 0.08 \pm 0.10 \pm 0.14)\%. \end{align*} All three results are the most precise determinations to date.


title: Limits on Fast Radio Bursts and other transient sources at 182 MHz using the Murchison Widefield Array

abstract: We present a survey for transient and variable sources, on timescales from 28 seconds to $\sim$1 year, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of $<4.1 \times 10^{-7}$ deg$^{-2}$ for the shortest timescale considered. At these frequencies, emission from Fast Radio Bursts (FRBs) is expected to be detectable in the shortest timescale images without any corrections for interstellar or intergalactic dispersion. At an FRB limiting flux density of 7980 Jy, we find a rate of $<$82 FRBs per sky per day for dispersion measures $<$700 pc cm$^{-3}$. Assuming a cosmological population of standard candles, our rate limits are consistent with the FRB rates obtained by Thornton et al. (2013) if they have a flat spectral slope. Finally, we conduct an initial variability survey of sources in the field with flux densities $\gtrsim$0.5 Jy and identify no sources with significant variability in their lightcurves. However, we note that substantial further work is required to fully characterise both the short term and low level variability within this field.


title: Contact resistances in trigate and FinFET devices in a Non-Equilibrium Green's Functions approach

abstract: We compute the contact resistances $R_{\rm c}$ in trigate and FinFET devices with widths and heights in the 4 to 24 nm range using a Non-Equilibrium Green's Functions approach. Electron-phonon, surface roughness and Coulomb scattering are taken into account. We show that $R_{\rm c}$ represents a significant part of the total resistance of devices with sub-30 nm gate lengths. The analysis of the quasi-Fermi level profile reveals that the spacers between the heavily doped source/drain and the gate are major contributors to the contact resistance. The conductance is indeed limited by the poor electrostatic control over the carrier density under the spacers. We then disentangle the ballistic and diffusive components of $R_{\rm c}$, and analyze the impact of different design parameters (cross section and doping profile in the contacts) on the electrical performances of the devices. The contact resistance and variability rapidly increase when the cross sectional area of the channel goes below $\simeq 50$ nm$^2$. We also highlight the role of the charges trapped at the interface between silicon and the spacer material.


title: Variational Principle for Gravity with Null and Non-null boundaries: A Unified Boundary Counter-term

abstract: It is common knowledge that the Einstein-Hilbert action does not furnish a well-posed variational principle. The usual solution to this problem is to add an extra boundary term to the action, called a counter-term, so that the variational principle becomes well-posed. When the boundary is spacelike or timelike, the Gibbons-Hawking-York counter-term is the most widely used. For null boundaries, we had proposed a counter-term in a previous paper. In this paper, we extend the previous analysis and propose a counter-term that can be used to eliminate variations of the "off-the-surface" derivatives of the metric on any boundary, regardless of its spacelike, timelike or null nature.


title: Critical behavior of sputter-deposited magnetoelectric antiferromagnetic Cr$_2$O$_3$ films near N\'eel temperature

abstract: Chromium(III) oxide is a classical collinear antiferromagnet with a linear magnetoelectric effect. We are presenting the measurements of the magnetoelectric susceptibility $\alpha$ of a sputter-deposited 500-nm film and a bulk single-crystal of Cr$_\mathrm{2}$O$_\mathrm{3}$. We investigated the magnetic phase-transition and the critical exponent $\beta$ of the sublattice magnetization near N\'eel temperature. For the films, an exponent of 0.49(1) was found below 293 K, and changed to 1.06(4) near the N\'eel temperature of 298 K. For the single-crystal, the exponent was constant at 0.324(4). We investigated the reversal probability of antiferromagnetic domains during magnetoelectric field cooling. For the sputtered films, reversal probability was zero above 298 K and stabilized only below 293 K. We attribute this behavior to formation of grains during film growth, which gives different intergrain and intragrain exchange-coupling energies. For the single-crystal, reversal probability was stabilized immediately at the N\'eel temperature of 307.6 K.


title: Complex Evaluation of Hierarchically-Network Systems

abstract: Methods of complex evaluation based on local, forecasting, aggregated, and interactive evaluation of the state, function quality, and interaction of complex system's objects on the all hierarchical levels is proposed. Examples of analysis of the structural elements of railway transport system are used for illustration of efficiency of proposed approach.


title: The Landau-Lifshitz equation of the ferromagnetic spin chain and Oseen-Frank flow

abstract: In this paper, we consider the Landau-Lifshitz equation of the ferromagnetic spin chain from $\R^2$ to the unit sphere $S^2$ under the general Oseen-Frank energy. We obtain global existence and uniqueness of weak solutions for large energy data; moreover, the number of singular points is finite.


title: Numerical Diagnostics for Systems of Differential Algebraic Equations

abstract: In many commercial and academic settings, numerical solvers fail to achieve their theoretical performance levels due to issues in the system definition, parameterization, and even implementation. We propose a pair of methods for detecting and localizing these convergence rate issues in applications that can be treated as homotopy problems including numerical continuation and the evolution of differential algebraic equations. Both approaches are rooted in dynamical systems theory, in particular, the numerical techniques used to perform bifurcation studies on "black-box" systems, and can be applied across a range of numerical solvers and systems without significant modification. In a general problem, these methods will detect certain classes of convergence anomalies, and indicate the states that are affected by their existence. However in applications like circuit simulation, certain classes of anomalies can be localized to an individual component through the use of a small number of additional computational experiments. To demonstrate the efficacy of the approaches for detection and localization, we apply them to a pair of nonlinear circuits: a diode bridge and a simplified model of a power channel


title: What Mathematical Logic Says about the Foundations of Mathematics

abstract: My purpose is to examine some concepts of mathematical logic, which have been studied by Carlo Cellucci. Today the aim of classical mathematical logic is not to guarantee the certainty of mathematics, but I will argue that logic can help us to explain mathematical activity; the point is to discuss what and in which sense logic can "explain". For example, let's consider the basic concept of an axiomatic system: an axiomatic system can be very useful to organize, to present, and to clarify mathematical knowledge. And, more importantly, logic is a science with its own results: so, axiomatic systems are interesting also because we know several revealing theorems about them. Similarly, I will discuss other topics such as mathematical definitions, and some relationships between mathematical logic and computer science. I will also consider these subjects from an educational point of view: can logical concepts be useful in teaching and learning elementary mathematics?


title: Un problema da discutere: un inganno con le percentuali

abstract: We discuss a non-intuitive situation concerning percentages.


title: Gli angoli alla base di un triangolo isoscele

abstract: This paper deals with the celebrated Euclidean theorem about isosceles triangles, comparing different proofs.


title: Un problema da discutere: una rappresentazione geometrica del teorema del coseno

abstract: We present, discuss and generalize an elegant geometrical proof of the law of cosines, due to Al Cuoco.


title: Graphs of real functions with pathological behaviors

abstract: This paper deals with constructions and properties of unusual function from R to R, as discontinuous additive functions and everywhere surjections.


title: Primitive sets of nonnegative matrices and synchronizing automata

abstract: A set of nonnegative matrices $\mathcal{M}=\{M_1, M_2, \ldots, M_k\}$ is called primitive if there exist indices $i_1, i_2, \ldots, i_m$ such that $M_{i_1} M_{i_2} \ldots M_{i_m}$ is positive (i.e. has all its entries $>0$). The length of the shortest such product is called the exponent of $\mathcal{M}$. The concept of primitive sets of matrices comes up in a number of problems within control theory, non-homogeneous Markov chains, automata theory etc. Recently, connections between synchronizing automata and primitive sets of matrices were established. In the present paper, we significantly strengthen these links by providing equivalence results, both in terms of combinatorial characterization, and computational aspects. We study the maximal exponent among all primitive sets of $n \times n$ matrices, which we denote by $\exp(n)$. We prove that $\lim_{n\rightarrow\infty} \tfrac{\log \exp(n)}{n} = \tfrac{\log 3}{3}$, and moreover, we establish that this bound leads to a resolution of the \v{C}ern\'{y} problem for carefully synchronizing automata. We also study the set of matrices with no zero rows and columns, denoted by $\mathcal{NZ}$, due to its intriguing connections to the \v{C}ern\'{y} conjecture and the recent generalization of Perron-Frobenius theory for this class. We characterize computational complexity of different problems related to the exponent of $\mathcal{NZ}$ matrix sets, and present a quadratic bound on the exponents of sets belonging to a special subclass. Namely, we show that the exponent of a set of matrices having total support is bounded by $2n^2 -5n +5$.


title: The Kelmans-Seymour conjecture II: 2-vertices in $K_4^-$

abstract: We use $K_4^-$ to denote the graph obtained from $K_4$ by removing an edge, and use $TK_5$ to denote a subdivision of $K_5$. Let $G$ be a 5-connected nonplanar graph and $\{x_1,x_2,y_1,y_2\}\subseteq V(G)$ such that $G[\{x_1,x_2,$ $y_1,y_2\}]\cong K_4^-$ with $y_1y_2\notin E(G)$. Let $w_1,w_2,w_3\in N(y_2)-\{x_1,x_2\}$ be distinct. We show that $G$ contains a $TK_5$ in which $y_2$ is not a branch vertex, or $G-y_2$ contains $K_4^-$, or $G$ has a special 5-separation, or $G-\{y_2v:v\notin \{w_1,w_2,w_3,x_1,x_2\}\}$ contains $TK_5$.


title: The swept rule for breaking the latency barrier in time advancing two-dimensional PDEs

abstract: This article describes a method to accelerate parallel, explicit time integration of two-dimensional unsteady PDEs. The method is motivated by our observation that latency, not bandwidth, often limits how fast PDEs can be solved in parallel. The method is called the swept rule of space-time domain decomposition. Compared to conventional, space-only domain decomposition, it communicates similar amount of data, but in fewer messages. The swept rule achieves this by decomposing space and time among computing nodes in ways that exploit the domains of influence and the domain of dependency, making it possible to communicate once per many time steps with no redundant computation. By communicating less often, the swept rule effectively breaks the latency barrier, advancing on average more than one time step per ping-pong latency of the network. The article presents simple theoretical analysis to the performance of the swept rule in two spatial dimensions, and supports the analysis with numerical experiments.


title: Regression of ranked responses when raw responses are censored

abstract: We discuss semiparametric regression when only the ranks of responses are observed. The model is $Y_i = F (\mathbf{x}_i'{\boldsymbol\beta}_0 + \varepsilon_i)$, where $Y_i$ is the unobserved response, $F$ is a monotone increasing function, $\mathbf{x}_i$ is a known $p-$vector of covariates, ${\boldsymbol\beta}_0$ is an unknown $p$-vector of interest, and $\varepsilon_i$ is an error term independent of $\mathbf{x}_i$. We observe $\{(\mathbf{x}_i,R_n(Y_i)) : i = 1,\ldots ,n\}$, where $R_n$ is the ordinal rank function. We explore a novel estimator under Gaussian assumptions. We discuss the literature, apply the method to an Alzheimer's disease biomarker, conduct simulation studies, and prove consistency and asymptotic normality.


title: Lattice-Boltzmann Hydrodynamics of Anisotropic Active Matter

abstract: A plethora of active matter models exist that describe the behavior of self-propelled particles (or swimmers), both with and without hydrodynamics. However, there are few studies that consider shape-anisotropic swimmers and include hydrodynamic interactions. Here, we introduce a simple method to simulate self-propelled colloids interacting hydrodynamically in a viscous medium using the lattice-Boltzmann technique. Our model is based on raspberry-type viscous coupling and a force/counter-force formalism which ensures that the system is force free. We consider several anisotropic shapes and characterize their hydrodynamic multipolar flow field. We demonstrate that shape-anisotropy can lead to the presence of a strong quadrupole and octupole moments, in addition to the principle dipole moment. The ability to simulate and characterize these higher-order moments will prove crucial for understanding the behavior of model swimmers in confining geometries.


title: The Quantum-Classical Boundary for Precision Interferometric Measurements

abstract: Understanding the fundamental limits on the precision to which an optical phase can be estimated is of key interest for many investigative techniques utilized across science and technology. We study the estimation of a fixed optical phase shift due to a sample which has an associated optical loss, and compare phase estimation strategies using classical and non-classical probe states. These comparisons are based on the attainable (quantum) Fisher information calculated per number of photons absorbed or scattered by the sample throughout the sensing process. We find that, for a given number of incident photons upon the unknown phase, non-classical techniques in principle provide less than a 20% reduction in root-mean-square-error (RMSE) in comparison with ideal classical techniques in multi-pass optical setups. Using classical techniques in a new optical setup we analyze, which incorporates additional stages of interference during the sensing process, the achievable reduction in RMSE afforded by non-classical techniques falls to only ~4%. We explain how these conclusions change when non-classical techniques are compared to classical probe states in non-ideal multi-pass optical setups, with additional photon losses due to the measurement apparatus.


title: On a theory of light scattering from a Bose-Einstein condensate

abstract: We consider a quantum theory of elastic light scattering from macroscopic atomic sample existing in the Bose-Einstein condensate (BEC) phase. Following to the second quantized formalism we introduce a set of coupled and closed diagram equations for the polariton propagator contributing to the $T$ -matrix and scattering amplitude. Our approach allows to follow important density correction to the quasi-energy structure caused by static interaction and radiation losses associated with incoherent scattering in the case of near resonance excitation.


title: Multilingual Twitter Sentiment Classification: The Role of Human Annotators

abstract: What are the limits of automated Twitter sentiment classification? We analyze a large set of manually labeled tweets in different languages, use them as training data, and construct automated classification models. It turns out that the quality of classification models depends much more on the quality and size of training data than on the type of the model trained. Experimental results indicate that there is no statistically significant difference between the performance of the top classification models. We quantify the quality of training data by applying various annotator agreement measures, and identify the weakest points of different datasets. We show that the model performance approaches the inter-annotator agreement when the size of the training set is sufficiently large. However, it is crucial to regularly monitor the self- and inter-annotator agreements since this improves the training datasets and consequently the model performance. Finally, we show that there is strong evidence that humans perceive the sentiment classes (negative, neutral, and positive) as ordered.


title: Non-linear modeling of the plasma response to RMPs in ASDEX Upgrade

abstract: The plasma response to Resonant Magnetic Perturbations (RMPs) in ASDEX Upgrade is modeled with the non-linear resistive MHD code JOREK, using input profiles that match those of the experiments as closely as possible. The RMP configuration for which Edge Localized Modes are best mitigated in experiments is related to the largest edge kink response observed near the X-point in modeling. On the edge resonant surfaces $q = m/n$, the coupling between the m + 2 kink component and the m resonant component is found to induce the amplification of the resonant magnetic perturbation. The ergodicity and the 3D-displacement near the X-point induced by the resonant amplification can only partly explain the density pumpout observed in experiments.


title: Stochastic Shortest Path with Energy Constraints in POMDPs

abstract: We consider partially observable Markov decision processes (POMDPs) with a set of target states and positive integer costs associated with every transition. The traditional optimization objective (stochastic shortest path) asks to minimize the expected total cost until the target set is reached. We extend the traditional framework of POMDPs to model energy consumption, which represents a hard constraint. The energy levels may increase and decrease with transitions, and the hard constraint requires that the energy level must remain positive in all steps till the target is reached. First, we present a novel algorithm for solving POMDPs with energy levels, developing on existing POMDP solvers and using RTDP as its main method. Our second contribution is related to policy representation. For larger POMDP instances the policies computed by existing solvers are too large to be understandable. We present an automated procedure based on machine learning techniques that automatically extracts important decisions of the policy allowing us to compute succinct human readable policies. Finally, we show experimentally that our algorithm performs well and computes succinct policies on a number of POMDP instances from the literature that were naturally enhanced with energy levels.


title: Time and Activity Sequence Prediction of Business Process Instances

abstract: The ability to know in advance the trend of running process instances, with respect to different features, such as the expected completion time, would allow business managers to timely counteract to undesired situations, in order to prevent losses. Therefore, the ability to accurately predict future features of running business process instances would be a very helpful aid when managing processes, especially under service level agreement constraints. However, making such accurate forecasts is not easy: many factors may influence the predicted features. Many approaches have been proposed to cope with this problem but all of them assume that the underling process is stationary. However, in real cases this assumption is not always true. In this work we present new methods for predicting the remaining time of running cases. In particular we propose a method, assuming process stationarity, which outperforms the state-of-the-art and two other methods which are able to make predictions even with non-stationary processes. We also describe an approach able to predict the full sequence of activities that a running case is going to take. All these methods are extensively evaluated on two real case studies.


title: New Stability and Exact Observability Conditions for Semilinear Wave Equations

abstract: The problem of estimating the initial state of 1-D wave equations with globally Lipschitz nonlinearities from boundary measurements on a finite interval was solved recently by using the sequence of forward and backward observers, and deriving the upper bound for exact observability time in terms of Linear Matrix Inequalities (LMIs) [5]. In the present paper, we generalize this result to n-D wave equations on a hypercube. This extension includes new LMI-based exponential stability conditions for n-D wave equations, as well as an upper bound on the minimum exact observability time in terms of LMIs. For 1-D wave equations with locally Lipschitz nonlinearities, we find an estimate on the region of initial conditions that are guaranteed to be uniquely recovered from the measurements. The efficiency of the results is illustrated by numerical examples.


title: A new S-type eigenvalue localization set for tensors and its applications

abstract: A new \emph{S}-type eigenvalue localization set for tensors is derived by breaking $N=\{1,2,\cdots,n\}$ into disjoint subsets $S$ and its complement. It is proved that this new set is tighter than those presented by Qi (Journal of Symbolic Computation 40 (2005) 1302-1324), Li et al. (Numer. Linear Algebra Appl. 21 (2014) 39-50) and Li et al. (Linear Algebra Appl. 493 (2016) 469-483). As applications, checkable sufficient conditions for the positive definiteness and the positive semi-definiteness of tensors are proposed. Moreover, based on this new set, we establish a new upper bound for the spectral radius of nonnegative tensors and a lower bound for the minimum \emph{H}-eigenvalue of weakly irreducible strong \emph{M}-tensors in this paper. We demonstrate that these bounds are sharper than those obtained by Li et al. (Numer. Linear Algebra Appl. 21 (2014) 39-50) and He and Huang (J. Inequal. Appl. 114 (2014) 2014). Numerical examples are also given to illustrate this fact.


title: The catastrophic fragmentation of Comet 332P (Ikeya-Murakami)

abstract: We describe 2016 January to April observations of the fragments of 332P/Ikeya-Murakami, a comet earlier observed in a 2010 October outburst (Ishiguro et al 2014). We present photometry of the fragments, and perform simulations to infer the time of breakup. We argue that the eastern-most rapidly brightening fragment ($F4$) best corresponds to the original nucleus, rather than the initial bright fragment $F1$. We compute radial and tangential non-gravitational parameters, $A_1 = (1.5 \pm 0.4) \times 10^{-8}$ AU day$^{-2}$ and $(7.2 \pm 1.9) \times 10^{-9}$ AU day$^{-2}$; both are consistent with zero at the $4\sigma$ level. Monte Carlo simulations indicate that the fragments were emitted on the outbound journey well after the 2010 outburst, with bright fragment $F1$ splitting in mid--2013 and the fainter fragments within months of the 2016 January recovery. Western fragment $F7$ is the oldest, dating from 2011. We suggest that the delayed onset of the splitting is consistent with a self-propagating crystallization of water ice.


title: Bayesian Exploration: Incentivizing Exploration in Bayesian Games

abstract: We consider a ubiquitous scenario in the Internet economy when individual decision-makers (henceforth, agents) both produce and consume information as they make strategic choices in an uncertain environment. This creates a three-way tradeoff between exploration (trying out insufficiently explored alternatives to help others in the future), exploitation (making optimal decisions given the information discovered by other agents), and incentives of the agents (who are myopically interested in exploitation, while preferring the others to explore). We posit a principal who controls the flow of information from agents that came before, and strives to coordinate the agents towards a socially optimal balance between exploration and exploitation, not using any monetary transfers. The goal is to design a recommendation policy for the principal which respects agents' incentives and minimizes a suitable notion of regret. We extend prior work in this direction to allow the agents to interact with one another in a shared environment: at each time step, multiple agents arrive to play a Bayesian game, receive recommendations, choose their actions, receive their payoffs, and then leave the game forever. The agents now face two sources of uncertainty: the actions of the other agents and the parameters of the uncertain game environment. Our main contribution is to show that the principal can achieve constant regret when the utilities are deterministic (where the constant depends on the prior distribution, but not on the time horizon), and logarithmic regret when the utilities are stochastic. As a key technical tool, we introduce the concept of explorable actions, the actions which some incentive-compatible policy can recommend with non-zero probability. We show how the principal can identify (and explore) all explorable actions, and use the revealed information to perform optimally.


title: The role of grain dynamics in determining the onset of sediment transport

abstract: Sediment transport occurs when the nondimensional fluid shear stress $\Theta$ at the bed surface exceeds a minimum value $\Theta_c$. A large collection of data, known as the Shields curve, shows that $\Theta_c$ is primarily a function of the shear Reynolds number ${\rm{Re}}_*$. It is commonly assumed that $\Theta>\Theta_c({\rm{Re}}_*)$ occurs when the ${\rm Re}_*$-dependent fluid forces are too large to maintain static equilibrium for a typical surface grain. A complimentary approach, which remains relatively unexplored, is to identify $\Theta_c({\rm{Re}}_*)$ as the applied shear stress at which grains cannot stop moving. With respect to grain dynamics, ${\rm{Re}}_*$ can be viewed as the viscous time scale for a grain to equilibrate to the fluid flow divided by the typical time for the fluid force to accelerate a grain over the characteristic bed roughness. We performed simulations of granular beds sheared by a model fluid, varying only these two time scales. We find that the critical Shields number $\Theta_c({\rm Re}_*)$ obtained from the model mimics the Shields curve and is insensitive to the grain properties, the model fluid flow, and the form of the drag law. Quantitative discrepancies between the model results and the Shields curve are consistent with previous calculations of lift forces at varying ${\rm Re}_*$. Grains at low ${\rm Re}_*$ find more stable configurations than those at high ${\rm{Re}}_*$ due to differences in the grain reorganization dynamics. Thus, instead of focusing on mechanical equilibrium of a typical grain at the bed surface, $\Theta_c({\rm{Re}}_*)$ may be better described by the stress at which mobile grains cannot find a stable configuration and stop moving.


title: Ultradense Word Embeddings by Orthogonal Transformation

abstract: Embeddings are generic representations that are useful for many NLP tasks. In this paper, we introduce DENSIFIER, a method that learns an orthogonal transformation of the embedding space that focuses the information relevant for a task in an ultradense subspace of a dimensionality that is smaller by a factor of 100 than the original space. We show that ultradense embeddings generated by DENSIFIER reach state of the art on a lexicon creation task in which words are annotated with three types of lexical information - sentiment, concreteness and frequency. On the SemEval2015 10B sentiment analysis task we show that no information is lost when the ultradense subspace is used, but training is an order of magnitude more efficient due to the compactness of the ultradense space.


title: A straightforward method to assess motion blur for different types of displays

abstract: A simulation method based on the liquid crystal response and the human visual system is suitable to characterize motion blur for LCDs but not other display types. We propose a more straightforward and widely applicable method to quantify motion blur based on the width of the moving object. We thus compare various types of displays objectively. A perceptual experiment was conducted to validate the proposed method. We test varying motion velocities for nine commercial displays. We compare the three motion blur evaluation methods (simulation, human perception, and our method) using z-scores. Our comparisons indicate that our method accurately characterizes motion blur for various display types.


title: Collider Tests of (Composite) Diphoton Resonances

abstract: We analyze the Large Hadron Collider sensitivity to new pseudoscalar resonances decaying into diphoton with masses up to scales of few TeVs. We focus on minimal scenarios where the production mechanisms involve either photon or top-mediated gluon fusion, partially motivated by the tantalizing excess around 750 GeV reported by ATLAS and CMS. The two scenarios lead respectively to a narrow and a wide resonance. We first provide a model-independent analysis via effective operators and then introduce minimal models of composite dynamics where the diphoton channel is characterized by their topological sector. The relevant state here is the pseudoscalar associated with the axial anomaly of the new composite dynamics. If the Standard Model top mass is generated via four-fermion operators the coupling of this state to the top remarkably explains the wide-width resonance reported by ATLAS. Beyond the excess, our analysis paves the way to test dynamical electroweak symmetry breaking via topological sectors.


title: Non-dipolar Wilson links for quasi-parton distribution functions

abstract: We propose a modified definition for a quasi-parton distribution function (QPDF) with an equal-time correlator in the large momentum limit, whose two pieces of space-like Wilson links are oriented in orthogonal directions. It is explicitly shown at one-loop level that the linear divergence in the original QPDF with dipolar Wilson links, which complicates its matching to the standard light-cone parton distribution function (LPDF), is removed. The LPDF can then be extracted reliably from Euclidean lattice data for the QPDF with the non-dipolar Wilson links.


title: Group Equivariant Convolutional Networks

abstract: We introduce Group equivariant Convolutional Neural Networks (G-CNNs), a natural generalization of convolutional neural networks that reduces sample complexity by exploiting symmetries. G-CNNs use G-convolutions, a new type of layer that enjoys a substantially higher degree of weight sharing than regular convolution layers. G-convolutions increase the expressive capacity of the network without increasing the number of parameters. Group convolution layers are easy to use and can be implemented with negligible computational overhead for discrete groups generated by translations, reflections and rotations. G-CNNs achieve state of the art results on CIFAR10 and rotated MNIST.


title: Interfacial and morphological features of a twist-bend nematic drop

abstract: In this experimental and theoretical study, we examine the equilibrium shapes of quasitwo-dimensional twist-bend nematic (Ntb) drops formed within a planarly aligned nematic layer of the liquid crystal CB7CB. Initially, at the setting point of the Ntb phase, the drops assume a nonequilibrium cusped elliptical geometry with the major axis orthogonal to the director of the surrounding nematic fluid; this growth is governed principally by anisotropic heat diffusion. The drops attain equilibrium through thermally driven dynamical evolutions close to their melting temperature. They are associated with a characteristic twin-striped morphology that transforms into the familiar focal conic texture as the temperature is lowered. At equilibrium, large millimetric drops are tactoidlike, elongated along the director of the surrounding nematic fluid. This geometry is explained by a mathematical model that features two dimensionless parameters, of which one is the structural cone angle of the Ntb phase and the other, the relative strength of mismatch elastic energy at the drop's interface. Both parameters are extracted from the observations by measuring the aspect ratio of the equilibrium shapes and the inner corner angle of the cusps.


title: An atomically thin matter-wave beamsplitter

abstract: Matter-wave interferometry has become an essential tool in studies on the foundations of quantum physics and for precision measurements. Mechanical gratings have played an important role as coherent beamsplitters for atoms, molecules and clusters since the basic diffraction mechanism is the same for all particles. However, polarizable objects may experience van der Waals shifts when they pass the grating walls and the undesired dephasing may prevent interferometry with massive objects. Here we explore how to minimize this perturbation by reducing the thickness of the diffraction mask to its ultimate physical limit, i.e. the thickness of a single atom. We have fabricated diffraction masks in single-layer and bilayer graphene as well as in 1 nm thin carbonaceous biphenyl membrane. We identify conditions to transform an array of single layer graphene nanoribbons into a grating of carbon nanoscrolls. We show that all these ultra-thin nanomasks can be used for high-contrast quantum diffraction of massive molecules. They can be seen as a nanomechanical answer to the question debated by Bohr and Einstein whether a softly suspended double slit would destroy quantum interference. In agreement with Bohr's reasoning we show that quantum coherence prevails even in the limit of atomically thin gratings.


title: Listen-and-Talk: Protocol Design and Analysis for Full-duplex Cognitive Radio Networks

abstract: In traditional cognitive radio networks, secondary users (SUs) typically access the spectrum of primary users (PUs) by a two-stage "listen-before-talk" (LBT) protocol, i.e., SUs sense the spectrum holes in the first stage before transmitting in the second. However, there exist two major problems: 1) transmission time reduction due to sensing, and 2) sensing accuracy impairment due to data transmission. In this paper, we propose a "listen-and-talk" (LAT) protocol with the help of full-duplex (FD) technique that allows SUs to simultaneously sense and access the vacant spectrum. Spectrum utilization performance is carefully analyzed, with the closed-form spectrum waste ratio and collision ratio with the PU provided. Also, regarding the secondary throughput, we report the existence of a tradeoff between the secondary transmit power and throughput. Based on the power-throughput tradeoff, we derive the analytical local optimal transmit power for SUs to achieve both high throughput and satisfying sensing accuracy. Numerical results are given to verify the proposed protocol and the theoretical results.


title: Nucleation pathway and kinetics of phase-separating active Brownian particles

abstract: Suspensions of purely repulsive but self-propelled Brownian particles might undergo phase separation, a phenomenon that strongly resembles the phase separation of passive particles with attractions. Here we employ computer simulations to study the nucleation kinetics and the microscopic pathway active Brownian disks take in two dimensions when quenched from the homogeneous suspension to propulsion speeds beyond the binodal. We find the same qualitative behavior for the nucleation rate as a function of density as for a passive suspension undergoing liquid-vapor separation, suggesting that the scenario of an effective free energy also extends to the kinetics of phase separation. We study the transition in more detail through a committor analysis and find that transition states are best described by a combination of cluster size and the radial polarization of particles in the cluster.


title: Balloon Flight Test of a Compton Telescope Based on Scintillators with Silicon Photomultiplier Readouts

abstract: We present the results of the first high-altitude balloon flight test of a concept for an advanced Compton telescope making use of modern scintillator materials with silicon photomultiplier (SiPM) readouts. There is a need in the fields of high-energy astronomy and solar physics for new medium-energy gamma-ray (~0.4 - 10 MeV) detectors capable of making sensitive observations. A fast scintillator- based Compton telescope with SiPM readouts is a promising solution to this instrumentation challenge, since the fast response of the scintillators permits the rejection of background via time-of-flight (ToF) discrimination. The Solar Compton Telescope (SolCompT) prototype was designed to demonstrate stable performance of this technology under balloon-flight conditions. The SolCompT instrument was a simple two-element Compton telescope, consisting of an approximately one-inch cylindrical stilbene crystal for a scattering detector and a one-inch cubic LaBr3:Ce crystal for a calorimeter detector. Both scintillator detectors were read out by 2 x 2 arrays of Hamamatsu S11828-3344 MPPC devices. Custom front-end electronics provided optimum signal rise time and linearity, and custom power supplies automatically adjusted the SiPM bias voltage to compensate for temperature-induced gain variations. A tagged calibration source, consisting of ~240 nCi of Co-60 embedded in plastic scintillator, was placed in the field of view and provided a known source of gamma rays to measure in flight. The SolCompT balloon payload was launched on 24 August 2014 from Fort Sumner, NM, and spent ~3.75 hours at a float altitude of ~123,000 feet. The instrument performed well throughout the flight. After correcting for small (~10%) residual gain variations, we measured an in-flight ToF resolution of ~760 ps (FWHM). Advanced scintillators with SiPM readouts continue to show great promise for future gamma-ray instruments.


title: Searching for spectroscopic binaries within transition disk objects

abstract: Transition disks (TDs) are intermediate stage circumstellar disks characterized by an inner gap within the disk structure. To test whether these gaps may have been formed by closely orbiting, previously undetected stellar companions, we collected high-resolution optical spectra of 31 TD objects to search for spectroscopic binaries (SBs). Twenty-four of these objects are in Ophiuchus and seven are within the Coronet, Corona Australis, and Chameleon I star-forming regions. We measured radial velocities for multiple epochs, obtaining a median precision of 400 ms$^{-1}$. We identified double-lined SB SSTc2d J163154.7-250324 in Ophiuchus, which we determined to be composed of a K7($\pm$0.5) and a K9($\pm$0.5) star, with orbital limits of $a<$0.6 AU and $P<$150 days. This results in an SB fraction of 0.04$^{+0.12}_{-0.03}$ in Ophiuchus, which is consistent with other spectroscopic surveys of non-TD objects in the region. This similarity suggests that TDs are not preferentially sculpted by the presence of close binaries and that planet formation around close binaries may take place over similar timescales to that around single stars.


title: Free algebras of discriminator varieties generated by finite algebras are atomic

abstract: We prove that all definable pre-orders are atomic, in a finitely generated free algebra of a discriminator variety of finite similarity type which is generated by its finite members.


title: The emergence of giant multicaloric phenomena near room temperature

abstract: Caloric responses (temperature changes) can be induced in solid-state materials by applying external stimuli such as stress, pressure, and electric and magnetic fields. The magnetic-field-stimulated response is called the magnetocaloric effect, and materials that exhibit this property have long been sought for applications in room temperature magnetic cooling due to their potentially superior efficiency and low impact on the environment. Other solid-state caloric phenomena are less developed, but are likewise under intense investigation. Here we introduce a new material that not only displays giant barocaloric (hydrostatic-pressure-induced) properties, but also a large magnetocaloric response near room temperature. It is unprecedented that two caloric effects of such extreme magnitude occur in the same material and at the same temperature. These effects originate from a magnetostructural transition and a magneto-volume (magnetostriction) effect where the volume change is large enough to force the system from a localized ordered state into an itinerant paramagnetic state.


title: Mapping toric varieties into low dimensional spaces

abstract: A smooth $d$-dimensional projective variety $X$ can always be embedded into $2d+1$-dimensional space. In contrast, a singular variety may require an arbitrary large ambient space. If we relax our requirement and ask only that the map is injective, then any $d$-dimensional projective variety can be mapped injectively to $2d+1$-dimensional projective space. A natural question then arises: what is the minimal $m$ such that a projective variety can be mapped injectively to $m$-dimensional projective space? In this paper we investigate this question for normal toric varieties, with our most complete results being for Segre-Veronese varieties.


title: A parsimonious theory of evidence-based choice

abstract: That an agent's possible evidential states form a Boolean algebra (on which it is natural to define a probability measure) is an assertion that ideally should be proved, rather than assumed, in justifying rational choice as a representation of expected utility. A more parsimonious, axiomatic characterization of evidence is provided here. Two primitive entities are evidential states and a relation, more specific than, between evidential states. The axioms specify that more-specific-than is a partial order, there is a minimally specific e-state, and more-specific-than is a separative order. Choice alternatives are another primitive entity. A plan is an assignment of a choice alternative to each evidential state. In general, plans satisfying a version of the sure-thing principle cannot be rationalized by expected utility. But there is such a rationalization if the evidential structure is a tree.


title: Blockmodels: A R-package for estimating in Latent Block Model and Stochastic Block Model, with various probability functions, with or without covariates

abstract: Analysis of the topology of a graph, regular or bipartite one, can be done by clustering for regular ones or co-clustering for bipartite ones. The Stochastic Block Model and the Latent Block Model are two models, which are very similar for respectively regular and bipartite graphs, based on probabilistic models. Initially developed for binary graphs, these models have been extended to valued networks with optional covariates on the edges. This paper present a implementation of a Variational EM algorithm for Stochastic Block Model and Latent Block Model for some common probability functions, Bernoulli, Gaussian and Poisson, without or with covariates, with some standard flavors, like multivariate extensions. This implementation allow automatic group number exploration and selection via the ICL criterion, and allow analyze networks with thousands of nodes in a reasonable amount of time.


title: Evidence for a $B_s^0 \pi^\pm$ State

abstract: We report evidence of a narrow structure, $X(5568)$, in the decay sequence $X(5568) \rightarrow B_s^0 \pi^{\pm}$, $B_s^0 \rightarrow J/\psi \phi$, $J/\psi\rightarrow \mu^+ \mu^-$, $\phi \rightarrow K^+K^-$. This is evidence for the first instance of a hadronic state with valence quarks of four different flavors. The mass and natural width of this state are measured to be $m = 5567.8 \pm 2.9 {\rm \thinspace (stat)} ^{+0.9}_{-1.9} {\rm \thinspace (syst)}$ MeV/$c^2$ and $\Gamma = 21.9 \pm 6.4 {\rm \thinspace (stat)} ^{+5.0}_{-2.5} {\rm \thinspace (syst)} $ MeV/$c^2$. If the decay is $X(5568) \rightarrow B_s^* \pi^{\pm} \rightarrow B_s^0 \gamma \pi^{\pm}$ with an unseen $\gamma$, $m(X(5568))$ will be shifted up by $m(B_s^*) - m(B_s^0) \sim 49$ MeV/$c^2$. This measurement is based on $10.4~\rm{fb^{-1}}$ of $p \bar p $ collision data at $\sqrt{s}$ = 1.96 TeV collected by the D0 experiment at the Fermilab Tevatron collider.


title: Structure of inactive states of a binary Lennard-Jones mixture

abstract: We study the structure of inactive states in a prototypical model glass, the Kob-Andersen binary Lennard-Jones mixture. These inactive states are obtained by transition path sampling and are at dynamical phase coexistence with an active equilibrium state. Configurations in the inactive states are kinetically stable and are located in deeper basins of the energy landscape than their active counterparts. By analyzing trajectory-to-trajectory fluctuations within the inactive state, we assess correlations between kinetic stability, energy and other structural properties. We show that measures of local order associated to stable local packings and bond-orientational order are weakly correlated with energy and kinetic stability. We discuss what kinds of structural measurement might capture the relevant dynamical features of the inactive state.


title: Geometry-Independent Determination of Radial Density Distributions in Molecular Cloud Cores and Other Astronomical Objects

abstract: We present a geometry-independent method for determining the shapes of radial volume density profiles of astronomical objects whose geometries are unknown, based on a single column density map. Such profiles are often critical to understand the physics and chemistry of molecular cloud cores, in which star formation takes place. The method presented here does not assume any geometry for the object being studied, thus removing a significant source of bias. Instead it exploits contour self-similarity in column density maps which appears to be common in data for astronomical objects. Our method may be applied to many types of astronomical objects and observable quantities so long as they satisfy a limited set of conditions which we describe in detail. We derive the method analytically, test it numerically, and illustrate its utility using 2MASS-derived dust extinction in molecular cloud cores. While not having made an extensive comparison of different density profiles, we find that the overall radial density distribution within molecular cloud cores is adequately described by an attenuated power law.


title: Linear Modeling of a Flexible Substructure Actuated through Piezoelectric Components for Use in Integrated Control/Structure Design

abstract: This study presents a generic TITOP (Two-Input Two-Output Port) model of a substructure actuated with embedded piezoelectric materials as actuators (PEAs), previously modeled with the FE technique. This allows intuitive assembly of actuated flexible substructures in large flexible multi-body systems. The modeling technique is applied to an illustrative example of a flexible beam with bonded piezoelectric strip and vibration attenuation of a chain of flexible beams.


title: Mean-variance risk-averse optimal control of systems governed by PDEs with random parameter fields using quadratic approximations

abstract: We present a method for optimal control of systems governed by partial differential equations (PDEs) with uncertain parameter fields. We consider an objective function that involves the mean and variance of the control objective, leading to a risk-averse optimal control problem. To make the problem tractable, we invoke a quadratic Taylor series approximation of the control objective with respect to the uncertain parameter. This enables deriving explicit expressions for the mean and variance of the control objective in terms of its gradients and Hessians with respect to the uncertain parameter. The risk-averse optimal control problem is then formulated as a PDE-constrained optimization problem with constraints given by the forward and adjoint PDEs defining these gradients and Hessians. The expressions for the mean and variance of the control objective under the quadratic approximation involve the trace of the (preconditioned) Hessian and are thus prohibitive to evaluate. To address this, we employ trace estimators that only require a modest number of Hessian-vector products. We illustrate our approach with two problems: the control of a semilinear elliptic PDE with an uncertain boundary source term, and the control of a linear elliptic PDE with an uncertain coefficient field. For the latter problem, we derive adjoint-based expressions for efficient computation of the gradient of the risk-averse objective with respect to the controls. Our method ensures that the cost of computing the risk-averse objective and its gradient with respect to the control, measured in the number of PDE solves, is independent of the (discretized) parameter and control dimensions, and depends only on the number of random vectors employed in the trace estimation. Finally, we present a comprehensive numerical study of an optimal control problem for fluid flow in a porous medium with uncertain permeability field.


title: Non-Truthful Position Auctions Are More Robust to Misspecification

abstract: In the standard single-dimensional model of position auctions, bidders agree on the relative values of the positions and each of them submits a single bid that is interpreted in terms of these values. Motivated by current practice in sponsored search we consider a situation where the auctioneer uses estimates of the relative values, which may be imprecise, and show that under both complete and incomplete information a non-truthful mechanism is able to support an efficient outcome in equilibrium for a wider range of these estimates than the VCG mechanism. We thus exhibit a property of the VCG mechanism that may help explain the surprising rarity with which it is used even in settings with unit demand, a relative lack of robustness to misspecification of the bidding language. The result for complete information concerns the generalized second-price mechanism and lends additional theoretical support to the use of this mechanism in practice. Particularly interesting from a technical perspective is the result for incomplete information, which is driven by a surprising connection between equilibrium bids in the VCG mechanism and the generalized first-price mechanism.


title: A Griffith-Euler-Bernoulli theory for thin brittle beams derived from nonlinear models in variational fracture mechanics

abstract: We study a planar thin brittle beam subject to elastic deformations and cracks described in terms of a nonlinear Griffith energy functional acting on $SBV$ deformations of the beam. In particular we consider the case in which elastic bulk contributions due to finite bending of the beam are comparable to the surface energy which is necessary to completely break the beam into several large pieces. In the limit of vanishing aspect ratio we rigorously derive an effective Griffith-Euler-Bernoulli functional which acts on piecewise $W^{2,2}$ regular curves representing the midline of the beam. The elastic part of this functional is the classical Euler-Bernoulli functional for thin beams in the bending dominated regime in terms of the curve's curvature. In addition there also emerges a fracture term proportional to the number of discontinuities of the curve and its first derivative.


title: Mean curvature flow of area decreasing maps between Riemann surfaces

abstract: In this article we give a complete description of the evolution of an area decreasing map $f:M\to N$ induced by its mean curvature in the situation where $M$ and $N$ are complete Riemann surfaces with bounded geometry, $M$ being compact, for which their sectional curvatures $\sigma_M$, $\sigma_N$ satisfy $\min\sigma_M\ge\sup\sigma_N$.


title: Optical switching and bistability in a four-level atomic systems

abstract: We explore the coherent control of nonlinear absorption of intense laser fields in four-level atomic systems. For instance, in a four-level ladder system, a coupling field creates electromagnetically induced transparency (EIT) with Aulter-Townes doublet for the probe field while the control field is absent. A large absorption peak appears at resonance as the control field is switched on. We show how such a large absorption leads to optical switching. Further, this large absorption gets diminished and a transparency window appears due to the saturation effects as the strength of the probe field is increased. We further demonstrate that the threshold of the optical bistability can be modified by suitable choices of the coupling and the control fields. In a four-level Y-type configuration, the effect of the control field on saturable absorption (SA) and reverse saturable absorption (RSA) is highlighted in the context of nonlinear absorption of the probe field. We achieve RSA and SA in a simple atomic system just by applying a control field.


title: Conformal Invariance of Graphene Sheets

abstract: Suspended graphene sheets exhibit correlated random deformations that can be studied under the framework of rough surfaces with a Hurst (roughness) exponent $0.72 \pm 0.01$. Here, we show that, independent of the temperature, the iso-height lines at the percolation threshold have a well-defined fractal dimension and are conformally invariant, sharing the same statistical properties as Schramm-Loewner evolution (SLE$_{\kappa}$) curves with $\kappa=2.24\pm0.07$. Interestingly, iso-height lines of other rough surfaces are not necessarily conformally invariant even if they have the same Hurst exponent, e.g. random Gaussian surfaces. We have found that the distribution of the modulus of the Fourier coefficients plays an important role on this property. Our results not only introduce a new universality class and place the study of suspended graphene membranes within the theory of critical phenomena, but also provide hints on the long-standing question about the origin of conformal invariance in iso-height lines of rough surfaces.


title: Forward Modelling of Standing Kink Modes in Coronal Loops II. Applications

abstract: Magnetohydrodynamic waves are believed to play a significant role in coronal heating, and could be used for remote diagnostics of solar plasma. Both the heating and diagnostic applications rely on a correct inversion (or backward modelling) of the observables into the thermal and magnetic structures of the plasma. However, owing to the limited availability of observables, this is an ill-posed issue. Forward Modelling is to establish a plausible mapping of plasma structuring into observables. In this study, we set up forward models of standing kink modes in coronal loops and simulate optically thin emissions in the extreme ultraviolet bandpasses, and then adjust plasma parameters and viewing angles to match three events of transverse loop oscillations observed by the Solar Dynamics Observatory/Atmospheric Imaging Assembly. We demonstrate that forward models could be effectively used to identify the oscillation overtone and polarization, to reproduce the general profile of oscillation amplitude and phase, and to predict multiple harmonic periodicities in the associated emission intensity and loop width variation.


title: Backtesting Lambda Value at Risk

abstract: A new risk measure, the lambda value at risk (Lambda VaR), has been recently proposed from a theoretical point of view as a generalization of the value at risk (VaR). The Lambda VaR appears attractive for its potential ability to solve several problems of the VaR. In this paper we propose three nonparametric backtesting methodologies for the Lambda VaR which exploit different features. Two of these tests directly assess the correctness of the level of coverage predicted by the model. One of these tests is bilateral and provides an asymptotic result. A third test assess the accuracy of the Lambda VaR that depends on the choice of the P&L distribution. However, this test requires the storage of more information. Finally, we perform a backtesting exercise and we compare our results with the ones from Hitaj and Peri (2015)


title: The LHAASO experiment: from Gamma-Ray Astronomy to Cosmic Rays

abstract: LHAASO is expected to be the most sensitive project to face the open problems in Galactic cosmic ray physics through a combined study of photon- and charged particle-induced extensive air showers in the energy range 10$^{11}$ - 10$^{17}$ eV. This new generation multi-component experiment will be able of continuously surveying the gamma-ray sky for steady and transient sources from about 100 GeV to PeV energies, thus opening for the first time the 10$^2$--10$^3$ TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov components) that will be measured in LHAASO will allow the study of the origin, acceleration and propagation of the radiation through a measurement of energy spectrum, elemental composition and anisotropy with unprecedented resolution. The installation of the experiment started at very high altitude in China (Daocheng site, Sichuan province, 4410 m a.s.l.). The commissioning of one fourth of the detector will be implemented in 2018. The completion of the installation is expected by the end of 2021.


title: Hard X-ray view of microlensing events in RX J1131-1231

abstract: RX J1131-1231 is a gravitationally lensed system which includes four images of a quasar lensed by an elliptical galaxy. The flux in the individual images is known to be affected by microlensing effect in the visible and X-ray bands. We study the multi-wavelength properties of RX J1131-1231 over a broad energy range, from optical to hard X-ray, during the periods of the microlensing caustic crossings. We aim to constrain the spatial extent of the X-ray emission region at different energies. We combine the data of the source monitoring in the visible band with the X-ray data of the Burst Alert Telescope (BAT) on board of SWIFT satellite and Chandra X-ray observatory. Inspecting the broad band spectrum and lightcurves of the source we identify several microlensing caustic crossing events, and study the details of variability of the source during these events. The caustic crossings of image A on MJD 55150 and 55500 produce strong variations of the overall X-ray flux from the source. In the soft X-ray band, the caustic crossing events are characterised by the exponential growth / decay of the flux ratio between the source images with the exponentiating time scale T_d~60 d, which corresponds to the caustic transit time of the source of the size ~10^{14} cm. Variability induced by the microlensing dominates the overall variability of the sum of the fluxes of four images of the source which should also be detectable in the hard X-ray lightcurve, if the hard X-ray source is comparably compact. Non-detection of the caustic crossing events in the SWIFT/BAT lightcurve of the source indicates that geometry of the region of emission of X-rays with energies above 15 keV is different from that of the lower energy X-rays, with the hard X-ray emission an order of magnitude larger in size.


title: Security of Quantum Key Distribution

abstract: The security issues facing quantum key distribution (QKD) are explained, herein focusing on those issues that are cryptographic and information theoretic in nature and not those based on physics. The problem of security criteria is addressed. It is demonstrated that an attacker's success probabilities are the fundamental criteria of security that any theoretic security criterion must relate to in order to have operational significance. The errors committed in the prevalent interpretation of the trace distance criterion are analyzed. The security proofs of QKD protocols are discussed and assessed in regard to three main features: their validity, completeness, and adequacy of the achieved numerical security level. Problems are identified in all these features. It appears that the QKD security situation is quite different from the common perception that a QKD-generated key is nearly perfectly secure. Built into our discussion is a simple but complete quantitative description of the information theoretic security of classical key distribution that is also applicable to the quantum situation. In the appendices, we provide a brief outline of the history of some major QKD security proofs, a rather unfavorable comparison of current QKD proven security with that of conventional symmetric key ciphers, and a list of objections and answers concerning some major points of the paper.


title: Minimal dilatation in Penner's construction

abstract: For all orientable closed surfaces, we determine the minimal dilatation among mapping classes arising from Penner's construction. We also discuss generalisations to surfaces with punctures.


title: Cold Milky Way Hi gas in filaments

abstract: We investigate data from the Galactic Effelsberg--Bonn HI Survey (EBHIS), supplemented with data from the third release of the Galactic All Sky Survey (GASS III) observed at Parkes. We explore the all sky distribution of the local Galactic HI gas with $|v_{\rm LSR}| < 25 $ kms$^{-1}$ on angular scales of 11' to 16'. Unsharp masking (USM) is applied to extract small scale features. We find cold filaments that are aligned with polarized dust emission and conclude that the cold neutral medium (CNM) is mostly organized in sheets that are, because of projection effects, observed as filaments. These filaments are associated with dust ridges, aligned with the magnetic field measured on the structures by Planck at 353 GHz. The CNM above latitudes $|b|>20^\circ$ is described by a log-normal distribution, with a median Doppler temperature $T_{\rm D} = 223$ K, derived from observed line widths that include turbulent contributions. The median neutral hydrogen (HI) column density is $N_{\rm HI} \simeq 10^{19.1}\,{\rm cm^{-2}}$. These CNM structures are embedded within a warm neutral medium (WNM) with $N_{\rm HI} \simeq 10^{20} {\rm cm^{-2}}$. Assuming an average distance of 100 pc, we derive for the CNM sheets a thickness of $< 0.3$ pc. Adopting a magnetic field strength of $B_{\rm tot} = (6.0 \pm 1.8)\mu$G, proposed by Heiles & Troland 2005, and assuming that the CNM filaments are confined by magnetic pressure, we estimate a thickness of 0.09 pc. Correspondingly the median volume density is in the range $ 14 < n < 47 {\rm cm^{-3}}$.


title: On occupation times of the first and third quadrants for planar Brownian motion

abstract: An open problem of interest, first infused into the applied probability community in the work of Bingham and Doney in 1988, (see \cite{Bingham}) is stated as follows: find the distribution of the quadrant occupation time of planar Brownian motion. In this short communication, we study an alternate formulation of this longstanding open problem: let $X(t), Y(t), t \geq 0$ be standard Brownian motions starting at $x,y$ respectively. Find the distribution of the total time $T=Leb\{t \in [0,1]: X(t) \times Y(t) >0\}$, when $x=y=0$, i.e., the occupation time of the union of the first and third quadrants. If two adjacent quadrants are used, the problem becomes much easier and the distribution of $T$ follows the arcsine law.


title: Local motifs in GeS$_2$-Ga$_2$S$_3$ glasses

abstract: The structure of (GeS$_2$)$_{0.75}$(Ga$_2$S$_3$)$_{0.25}$ and (GeS$_2$)$_{0.83}$(Ga$_2$S$_3$)$_{0.17}$ glasses was investigated by Raman scattering, high energy X-ray diffraction and extended X-ray absorption fine structure (EXAFS) measurements at the Ga and Ge K-edges. The reverse Monte Carlo simulation technique (RMC) was used to obtain structural models compatible with diffraction and EXAFS datasets. It was found that the coordination number of Ga is close to four. While Ge atoms have only S neighbors, Ga binds to S as well as to Ga atoms showing a violation of chemical ordering in GeS$_2$-Ga$_2$S$_3$ glasses. Analysis of the corner- and edge-sharing between [GeS$_{4/2}$] units revealed that about 30% of germanium atoms participate in the edge-shared tetrahedra.


title: Multicolour Ramsey Numbers of Odd Cycles

abstract: We show that for any positive integer $r$ there exists an integer $k$ and a $k$-colouring of the edges of $K_{2^{k}+1}$ with no monochromatic odd cycle of length less than $r$. This makes progress on a problem of Erd\H{o}s and Graham and answers a question of Chung. We use these colourings to give new lower bounds on the $k$-colour Ramsey number of the odd cycle and prove that, for all odd $r$ and all $k$ sufficiently large, there exists a constant $\epsilon = \epsilon(r) > 0$ such that $R_{k}(C_{r}) > (r-1)(2+\epsilon)^{k-1}$.


title: Classical logic and intuitionistic logic: equivalent formulations in natural deduction, G\"odel-Kolmogorov-Glivenko translation

abstract: This report first shows the equivalence bewteen several formulations of classical logic in intuitionistic logic (tertium non datur, reductio ad absurdum, Pierce's law). Then it establishes the correctness of the G\"odel-Kolmogorov translation, whose restriction to the propositional case is due to Glivenko. This translation maps a formula $F$ of first order logic to a formula $F^{\lnot\lnot}$ in such a way that $F$ is provable in classical logic if and only if $F^{\lnot\lnot}$ is provable in intuitionistic logic. All formal proofs are presented in natural deduction. These questions are well-known proof theoretical facts, but in textbooks, they are often ignored or left to the reader. Because of the combinatorial difficulty of some of the needed formal proofs, we hope that this report may be useful, in particular to students and colleagues from other areas.


title: QED Radiative Corrections to Low-Energy Moller and Bhabha Scattering

abstract: We present a treatment of the next-to-leading-order radiative corrections to unpolarized Moller and Bhabha scattering without resorting to ultra-relativistic approximations. We extend existing soft-photon radiative corrections with new hard-photon bremsstrahlung calculations so that the effect of photon emission is taken into account for any photon energy. This formulation is intended for application in the OLYMPUS experiment and the upcoming DarkLight experiment, but is applicable to a broad range of experiments at energies where QED is a sufficient description.


title: Close packing of rods on spherical surfaces

abstract: We study the optimal packing of short, hard spherocylinders confined to lie tangential to a spherical surface, using simulated annealing and molecular dynamics simulations. For clusters of up to twelve particles, we map out the changes in the geometry of the closest-packed configuration as a function of the aspect ratio $L/D$, where $L$ is the cylinder length and $D$ the diameter of the rods. We find a rich variety of cluster structures. For larger clusters, we find that the best-packed configurations up to around 100 particles are highly dependent on the exact number of particles and aspect ratio. For even larger clusters, we find largely disordered clusters for very short rods ($L/D = 0.25$), while slightly longer rods ($L/D = 0.5$ or $1$) prefer a global baseball-like geometry of smectic-like domains, similar to the behavior of large-scale nematic shells. Our results provide predictions for experimentally realizable systems of colloidal rods trapped at the interface of emulsion droplets.


title: Improving performance in quantum mechanics with explicit incentives to correct mistakes

abstract: An earlier investigation found that the performance of advanced students in a quantum mechanics course did not automatically improve from midterm to final exam on identical problems even when they were provided the correct solutions and their own graded exams. Here, we describe a study, which extended over four years, in which upper-level undergraduate students in a quantum physics course were given four identical problems in both the midterm exam and final exam. Approximately half of the students were given explicit incentives to correct their mistakes in the midterm exam. In particular, they could get back up to 50\% of the points lost on each midterm exam problem. The solutions to the midterm exam problems were provided to all students in both groups but those who corrected their mistakes were provided the solution after they submitted their corrections to the instructor. The performance on the same problems on the final exam suggests that students who were given incentives to correct their mistakes significantly outperformed those who were not given an incentive. The incentive to correct the mistakes had greater impact on the final exam performance of students who had not performed well on the midterm exam.


title: A quantum random walk of a Bose-Einstein condensate in momentum space

abstract: Each step in a quantum random walk is typically understood to have two basic components; a `coin-toss' which produces a random superposition of two states, and a displacement which moves each component of the superposition by different amounts. Here we suggest the realization of a walk in momentum space with a spinor Bose-Einstein condensate subject to a quantum ratchet realized with a pulsed, off-resonant optical lattice. By an appropriate choice of the lattice detuning, we show how the atomic momentum can be entangled with the internal spin states of the atoms. For the coin-toss, we propose to use a microwave pulse to mix these internal states. We present experimental results showing an optimized quantum ratchet, and through a series of simulations, demonstrate how our proposal gives extraordinary control of the quantum walk. This should allow for the investigation of possible biases, and classical-to-quantum dynamics in the presence of natural and engineered noise.


title: Learning Shapes by Convex Composition

abstract: We present a mathematical and algorithmic scheme for learning the principal geometric elements in an image or 3D object. We build on recent work that convexifies the basic problem of finding a combination of a small number shapes that overlap and occlude one another in such a way that they "match" a given scene as closely as possible. This paper derives general sufficient conditions under which this convex shape composition identifies a target composition. From a computational standpoint, we present two different methods for solving the associated optimization programs. The first method simply recasts the problem as a linear program, while the second uses the alternating direction method of multipliers with a series of easily computed proximal operators.


title: A Model of Selective Advantage for the Efficient Inference of Cancer Clonal Evolution

abstract: Recently, there has been a resurgence of interest in rigorous algorithms for the inference of cancer progression from genomic data. The motivations are manifold: (i) growing NGS and single cell data from cancer patients, (ii) need for novel Data Science and Machine Learning algorithms to infer models of cancer progression, and (iii) a desire to understand the temporal and heterogeneous structure of tumor to tame its progression by efficacious therapeutic intervention. This thesis presents a multi-disciplinary effort to model tumor progression involving successive accumulation of genetic alterations, each resulting populations manifesting themselves in a cancer phenotype. The framework presented in this work along with algorithms derived from it, represents a novel approach for inferring cancer progression, whose accuracy and convergence rates surpass the existing techniques. The approach derives its power from several fields including algorithms in machine learning, theory of causality and cancer biology. Furthermore, a modular pipeline to extract ensemble-level progression models from sequenced cancer genomes is proposed. The pipeline combines state-of-the-art techniques for sample stratification, driver selection, identification of fitness-equivalent exclusive alterations and progression model inference. Furthermore, the results are validated by synthetic data with realistic generative models, and empirically interpreted in the context of real cancer datasets; in the later case, biologically significant conclusions are also highlighted. Specifically, it demonstrates the pipeline's ability to reproduce much of the knowledge on colorectal cancer, as well as to suggest novel hypotheses. Lastly, it also proves that the proposed framework can be applied to reconstruct the evolutionary history of cancer clones in single patients, as illustrated by an example from clear cell renal carcinomas.


title: Content Providers Volunteering to Pay Network Providers: Better than Neutrality?

abstract: This paper studies the effects on user welfare of imposing network neutrality, using a game-theoretic model of provider interactions based on a two-sided market framework: we assume that the platform--the last-mile access providers (ISPs)--are monopolists, and consider content providers (CPs) entry decisions. All decisions affect the choices made by users, who are sensitive both to CP and ISP investments (in content creation and quality-of-service, respectively). In a non-neutral regime, CPs and ISPs can charge each other, while such charges are prohibited in the neutral regime. We assume those charges (if any) are chosen by CPs, a direction rarely considered in the literature, where they are assumed fixed by ISPs. Our analysis suggests that, unexpectedly, more CPs enter the market in a non-neutral regime where they pay ISPs, than without such payments. Additionally, in this case ISPs tend to invest more than in the neutral regime. From our results, the best regime in terms of user welfare is parameter dependent, calling for caution in designing neutrality regulations.