title: Benchmarking in cluster analysis: A white paper

abstract: Note: A revised version of this is now published. Please cite and read (it's open access): Van Mechelen, I., Boulesteix, A.-L., Dangl, R., Dean, N., Hennig, C., Leisch, F., Steinley, D., Warrens, M. J. (2023). A white paper on good research practices in benchmarking: The case of cluster analysis. WIREs Data Mining and Knowledge Discovery, e1511. https://doi.org/10.1002/widm.1511 To achieve scientific progress in terms of building a cumulative body of knowledge, careful attention to benchmarking is of the utmost importance. This means that proposals of new methods of data pre-processing, new data-analytic techniques, and new methods of output post-processing, should be extensively and carefully compared with existing alternatives, and that existing methods should be subjected to neutral comparison studies. To date, benchmarking and recommendations for benchmarking have been frequently seen in the context of supervised learning. Unfortunately, there has been a dearth of guidelines for benchmarking in an unsupervised setting, with the area of clustering as an important subdomain. To address this problem, discussion is given to the theoretical conceptual underpinnings of benchmarking in the field of cluster analysis by means of simulated as well as empirical data. Subsequently, the practicalities of how to address benchmarking questions in clustering are dealt with, and foundational recommendations are made.


title: RPC performance vs FE electoronics and detector parameters

abstract: The first Resistive Plate Chambers detectors were developed for cosmic ray experiments, where low rate capability, good time resolution and low cost per unit of area were needed. These same features, except for the low rate capability, were required in the muon spectrometers of the collider experiments like LHC. For this purpose newRPCdetectors with increased rate capability were developed. The rate capability improvement has been achieved thanks to the transition from streamer to saturated avalanche regime, in which the average charge produced in the gas dicharge is smaller. The price to pay for this working mode switch is the transfer of the amplification from the detector to the FE electorinics. The High luminosity LHC and the future colliders will require even greater rate capability in the muon spectrometer compared to the current one. For this reason, further improving of the rate capability is required. The transition from the saturated avalanche to a low saturated avalanche regime moves in this direction and needs a new front end electronics with better signal to noise ratio, because the average induced charge on the electrode is even smaller. The front end electronics design is crucial for the RPC performances. In this paper we discuss the performances of the RPC detector changing the front end design and the detector parameters.


title: Effective dynamics for non-reversible stochastic differential equations: a quantitative study

abstract: Coarse-graining is central to reducing dimensionality in molecular dynamics, and is typically characterized by a mapping which projects the full state of the system to a smaller class of variables. While extensive literature has been devoted to coarse-graining starting from reversible systems, not much is known in the non-reversible setting. In this article, starting with a non-reversible dynamics, we introduce and study an effective dynamics which approximates the (non-closed) projected dynamics. Under fairly weak conditions on the system, we prove error bounds on the trajectorial error between the projected and the effective dynamics. In addition to extending existing results to the non-reversible setting, our error estimates also indicate that the notion of mean force motivated by this effective dynamics is a good one.


title: Collective behavior recognition using compact descriptors

abstract: This paper presents a novel hierarchical approach for collective behavior recognition based solely on ground-plane trajectories. In the first layer of our classifier, we introduce a novel feature called Personal Interaction Descriptor (PID), which combines the spatial distribution of a pair of pedestrians within a temporal window with a pyramidal representation of the relative speed to detect pairwise interactions. These interactions are then combined with higher level features related to the mean speed and shape formed by the pedestrians in the scene, generating a Collective Behavior Descriptor (CBD) that is used to identify collective behaviors in a second stage. In both layers, Random Forests were used as classifiers, since they allow features of different natures to be combined seamlessly. Our experimental results indicate that the proposed method achieves results on par with state of the art techniques with a better balance of class errors. Moreover, we show that our method can generalize well across different camera setups through cross-dataset experiments.


title: $\mathrm{Fe}_{1-x}\mathrm{Ni}_x$ Alloy Nanoparticles Encapsulated inside Carbon Nanotubes: Controlled Synthesis, Structure and Magnetic Properties

abstract: In the present work, different synthesis procedures have been demonstrated to fill carbon nanotubes (CNTs) with $\mathrm{Fe}_{1-x}\mathrm{Ni}_x$ alloy nanoparticles (x = 0.33, 0.5). CNTs act as templates for the encapsulation of magnetic nanoparticles, and provide a protective shield against oxidation as well as prevent nanoparticles agglomeration. By variation of the reaction parameters, the purity of the samples, degree of filling, the composition and size of filling nanoparticles have been tailored and therefore the magnetic properties. The samples were analyzed by scanning electron microscopy (SEM), transmission electron microscopy (TEM), Bright-field (BF) TEM tomography, X-ray powder diffraction, superconducting quantum interference device (SQUID) and thermogravimetric analysis (TGA). The Fe1-xNix-filled CNTs show a huge enhancement in the coercive fields compared to the corresponding bulk materials, which make them excellent candidates for several applications such as magnetic storage devices.


title: Blind Ghost Imaging

abstract: Ghost imaging is an unconventional optical imaging technique that reconstructs the shape of an object combining the measurement of two signals: one that interacted with the object, but without any spatial information, the other containing spatial information, but that never interacted with the object. Ghost imaging is a very flexible technique, that has been generalized to the single-photon regime, to the time domain, to infrared and terahertz frequencies, and many more conditions. Here we demonstrate that ghost imaging can be performed without ever knowing the patterns illuminating the object, but using patterns correlated with them, doesn't matter how weakly. As an experimental proof we exploit the recently discovered correlation between the reflected and transmitted light from a scattering layer, and reconstruct the image of an object hidden behind a scattering layer using only the reflected light, which never interacts with the object. This method opens new perspectives for non-invasive imaging behind or within turbid media.


title: Black Hole Formation and Explosion from Rapidly Rotating Very Massive Stars

abstract: We explore the formation process of a black hole (BH) through the pair-instability collapse of a rotating Population III very massive star in axisymmetric numerical relativity. As the initial condition, we employ a progenitor star which is obtained by evolving a rapidly rotating zero-age main sequence (ZAMS) star with mass $320M_\odot$ until it reaches a pair instability region. We find that for such rapidly rotating model, a fraction of the mass, $\sim 10M_\odot$, forms a torus surrounding the remnant BH of mass $\sim 130M_\odot$ and an outflow is driven by a hydrodynamical effect. We also perform simulations, artificially reducing the initial angular velocity of the progenitor star, and find that only a small or no torus is formed and no outflow is driven. We discuss the possible evolution scenario of the remnant torus for the rapidly rotating model by considering the viscous and recombination effects and show that if the energy of $\sim 10^{52}$ erg is injected from the torus to the envelope, the luminosity and timescale of the explosion could be of the orders of $10^{43}$ erg/s and yrs, respectively. We also point out the possibility for observing gravitational waves associated with the BH formation for the rapidly rotating model by ground-based gravitational-wave detectors.


title: Equilibria in Quantitative Concurrent Games

abstract: Synthesis of finite-state controllers from high-level specifications in multi-agent systems can be reduced to solving multi-player concurrent games over finite graphs. The complexity of solving such games with qualitative objectives for agents, such as reaching a target set, is well understood resulting in tools with applications in robotics. In this paper, we introduce quantitative concurrent graph games, where transitions have separate costs for different agents, and each agent attempts to reach its target set while minimizing its own cost along the path. In this model, a solution to the game corresponds to a set of strategies, one per agent, that forms a Nash equilibrium. We study the problem of computing the set of all Pareto-optimal Nash equilibria, and give a comprehensive analysis of its complexity and related problems such as the price of stability and the price of anarchy. In particular, while checking the existence of a Nash equilibrium is NP-complete in general, with multiple parameters contributing to the computational hardness separately, two-player games with bounded costs on individual transitions admit a polynomial-time solution.


title: A rotation-equivariant convolutional neural network model of primary visual cortex

abstract: Classical models describe primary visual cortex (V1) as a filter bank of orientation-selective linear-nonlinear (LN) or energy models, but these models fail to predict neural responses to natural stimuli accurately. Recent work shows that models based on convolutional neural networks (CNNs) lead to much more accurate predictions, but it remains unclear which features are extracted by V1 neurons beyond orientation selectivity and phase invariance. Here we work towards systematically studying V1 computations by categorizing neurons into groups that perform similar computations. We present a framework to identify common features independent of individual neurons' orientation selectivity by using a rotation-equivariant convolutional neural network, which automatically extracts every feature at multiple different orientations. We fit this model to responses of a population of 6000 neurons to natural images recorded in mouse primary visual cortex using two-photon imaging. We show that our rotation-equivariant network not only outperforms a regular CNN with the same number of feature maps, but also reveals a number of common features shared by many V1 neurons, which deviate from the typical textbook idea of V1 as a bank of Gabor filters. Our findings are a first step towards a powerful new tool to study the nonlinear computations in V1.


title: The Convergence of Sparsified Gradient Methods

abstract: Distributed training of massive machine learning models, in particular deep neural networks, via Stochastic Gradient Descent (SGD) is becoming commonplace. Several families of communication-reduction methods, such as quantization, large-batch methods, and gradient sparsification, have been proposed. To date, gradient sparsification methods - where each node sorts gradients by magnitude, and only communicates a subset of the components, accumulating the rest locally - are known to yield some of the largest practical gains. Such methods can reduce the amount of communication per step by up to three orders of magnitude, while preserving model accuracy. Yet, this family of methods currently has no theoretical justification. This is the question we address in this paper. We prove that, under analytic assumptions, sparsifying gradients by magnitude with local error correction provides convergence guarantees, for both convex and non-convex smooth objectives, for data-parallel SGD. The main insight is that sparsification methods implicitly maintain bounds on the maximum impact of stale updates, thanks to selection by magnitude. Our analysis and empirical validation also reveal that these methods do require analytical conditions to converge well, justifying existing heuristics.


title: Top-Quark Physics at the LHC

abstract: We report on the precision determination of the top-quark mass to next-to-next-to-leading order in QCD in well-defined renormalization schemes using data from the Large Hadron Collider for single-top and top-quark pair production. We also discuss the calibration of the so-called Monte Carlo top-quark mass parameter which is determined from a comparison to events with top-quark decay products. The implications of the measured value of the top-quark mass for conclusions about the stability of the electroweak vacuum state of our Universe are illustrated. At future lepton colliders, we provide for the first time matched exclusive calculations valid both at the top threshold and in the continuum, also fully differentially. In addition, we calculate fully off-shell top-pair production (also with an associated Higgs boson) at next-to-leading order in QCD, which allows to extract the top-Yukawa coupling with an unprecedented precision.


title: H2 chemistry in interstellar ices: the case of CO ice hydrogenation in UV irradiated CO:H2 ice mixtures

abstract: Context. In dense clouds, hydrogenation reactions on icy dust grains are key in the formation of molecules, like formaldehyde, methanol, and complex organic molecules (COMs). These species form through the sequential hydrogenation of CO ice. Although molecular hydrogen (H2) abundances can be four orders of magnitude higher than those of free H-atoms in dense clouds, H2 surface chemistry has been largely ignored; several laboratory studies show that H2 does not actively participate in non-energetic ice chemistry because of the high activation energies required. Aims. For the example of CO ice hydrogenation, we experimentally investigated the potential role of H2 molecules on the surface chemistry when energetic processing (i.e., UV photolysis) is involved. We test whether additional hydrogenation pathways become available upon UV irradiation of a CO:H2 ice mixture and whether this reaction mechanism also applies to other chemical systems. Methods. Ultra-high vacuum (UHV) experiments were performed at 8~20 K. A pre-deposited solid mixture of CO:H2 was irradiated with UV-photons. Reflection absorption infrared spectroscopy (RAIRS) was used as an in situ diagnostic tool. Single reaction steps and possible isotopic effects were studied by comparing results from CO:H2 and CO:D2 ice mixtures. Results. After UV-irradiation of a CO:H2 ice mixture, two photon-induced products, HCO and H2CO, are unambiguously detected. The proposed reaction mechanism involves electronically excited CO in the following reaction steps: CO + h{\nu}->CO*, CO* + H2->HCO + H where newly formed H-atoms are then available for further hydrogenation reactions. The HCO formation yields have a strong temperature dependence for the investigated regime, which is most likely linked to the H2 sticking coefficient. Finally, the astronomical relevance of this photo-induced reaction channel is discussed.


title: Distance and routing labeling schemes for cube-free median graphs

abstract: Distance labeling schemes are schemes that label the vertices of a graph with short labels in such a way that the distance between any two vertices $u$ and $v$ can be determined efficiently by merely inspecting the labels of $u$ and $v$, without using any other information. Similarly, routing labeling schemes label the vertices of a graph in a such a way that given the labels of a source node and a destination node, it is possible to compute efficiently the port number of the edge from the source that heads in the direction of the destination. One of important problems is finding natural classes of graphs admitting distance and/or routing labeling schemes with labels of polylogarithmic size. In this paper, we show that the class of cube-free median graphs on $n$ nodes enjoys distance and routing labeling schemes with labels of $O(\log^3 n)$ bits.


title: Charge-doping-induced variation of BaFe$_2$As$_2$ electronic structure and the emerging physical effects

abstract: We studied the relationship between the charge doping and the correlation, and its effects on the spectral function of the BaFe$_2$As$_2$ compound in the framework of the density functional theory combined with the dynamical mean field theory (DFT+DMFT). The calculated mass enhancements showed that the electronic correlation varies systematically from weak to strong when moving from the heavily electron-doped regime to the heavily hole-doped one. Since the compound has a multi-orbital nature, the correlation is orbital-dependent and it increases as hole-doping increases. The Fe-3d$_{xy}$ (xy) orbital is much more correlated than the other orbitals, because it reaches its half-filled situation and has a narrower energy scale around the Fermi energy. Our findings can be consistently understood as the tendency of the heavily hole-doped BaFe$_2$As$_2$ compound to an orbital-selective Mott phase (OSMP). Moreover, the fact that the superconducting state of the heavily hole-doped BaFe$_2$As$_2$ is an extreme case of such a selective Mottness constrains the non-trivial role of the electronic correlation in iron-pnictide superconductors. In addition, the calculated spectral function shows a behavior that is compatible with experimental results reported for every charge-doped BaFe$_2$As$_2$ compound and clarifies the importance of the characterization of its physical effects on the material.


title: A Weak Reverse Holder Inequality for Caloric Measure

abstract: Following a result of Bennewitz-Lewis for non-doubling harmonic measure, we prove a criterion for non-doubling caloric measure to satisfy a weak reverse Holder inequality on an open set $\Omega$, assuming as a background hypothesis only that the essential boundary of $\Omega$ satisfies an appropriate parabolic version of Ahlfors-David regularity (which entails some backwards in time thickness). We also show that the weak reverse Holder estimate is equivalent to solvability of the initial Dirichlet problem with "lateral" data in $L^p$, for some $p<\infty$, in this setting.


title: Light-Neutrino Exchange and Long-Distance Contributions to $0\nu2\beta$ Decays: An Exploratory Study on $\pi\pi\to ee$

abstract: We present an exploratory lattice QCD calculation of the neutrinoless double beta decay $\pi\pi\to ee$. Under the mechanism of light-neutrino exchange, the decay amplitude involves significant long-distance contributions. The calculation reported here, with pion masses $m_\pi=420$ and 140 MeV, demonstrates that the decay amplitude can be computed from first principles using lattice methods. At unphysical and physical pion masses, we obtain that amplitudes are $24\%$ and $9\%$ smaller than the predication from leading order chiral perturbation theory. Our findings provide the lattice QCD inputs and constraints for effective field theory. A follow-on calculation with fully controlled systematic errors will be possible with adequate computational resources.


title: Kinematics of the Galaxy from OB Stars with Data from the Gaia DR2 Catalogue

abstract: We have selected and analyzed a sample of OB stars with known line-of-sight velocities determined through ground-based observations and with trigonometric parallaxes and proper motions from the Gaia DR2 catalogue. Some of the stars in our sample have distance estimates made from calcium lines. A direct comparison with the trigonometric distance scale has shown that the calcium distance scale should be reduced by 13\%. The following parameters of the Galactic rotation curve have been determined from 495 OB stars with relative parallax errors less than 30\%: $(U,V,W)_\odot=(8.16,11.19,8.55)\pm(0.48,0.56,0.48)$ km s$^{-1}$, $\Omega_0=28.92\pm0.39$ km s$^{-1}$ kpc$^{-1}$, $\Omega^{'}_0=-4.087\pm0.083$ km s$^{-1}$ kpc$^{-2}$ and $\Omega^{''}_0=0.703\pm0.067$ km s$^{-1}$ kpc$^{-3}$, where the circular velocity of the local standard of rest is $V_0=231\pm5$ km s$^{-1}$ (for the adopted $R_0=8.0\pm0.15$ kpc). The parameters of the Galactic spiral density wave have been found from the series of radial, $V_R,$ residual tangential, $\Delta V_{circ}$, and vertical, $W,$ velocities of OB stars by applying a periodogram analysis. The amplitudes of the radial, tangential, and vertical velocity perturbations are $f_R=7.1\pm0.3$ km s$^{-1}$, $f_\theta=6.5\pm0.4$ km s$^{-1}$, and $f_W=4.8\pm0.8$ km s$^{-1}$, respectively; the perturbation wavelengths are $ \lambda_R=3.3\pm0.1$ kpc, $\lambda_\theta=2.3\pm0.2$ kpc, and $ \lambda_W=2.6\pm0.5$ kpc; and the Sun's radial phase in the spiral density wave is $ (\chi_\odot)_R=-135\pm5^\circ$, $(\chi_\odot)_\theta=-123\pm8^\circ$, and $ (\chi_\odot)_W=-132\pm21^\circ$ for the adopted four-armed spiral pattern.


title: Radiative decays of the heavy tensor mesons in light cone QCD sum rules

abstract: The transition form factors of the radiative decays of the heavy tensor mesons to heavy pseudoscalar and heavy vector mesons are calculated in the framework of the light cone QCD sum rules method at the point $Q^2=0$. Using the obtained values of the transition form factors at the point $Q^2=0$ the corresponding decay widths are estimated. The results show that the radiative decays of the heavy--light tensor mesons can be measurable in the future planned experiments at LHCb.


title: The Cosmic Evolution of Magnesium Isotopes

abstract: The abundance of magnesium in the interstellar medium is a powerful probe of star formation processes over cosmological timescales. Magnesium has three stable isotopes, 24Mg, 25Mg, 26Mg, which can be produced both in massive and intermediate-mass (IM) stars with masses between 2 and 8 M_\odot. In this work, we use constraints on the cosmic star formation rate density (SFRD) and explore the role and mass range of intermediate mass stars using the observed isotopic ratios. We compare several models of stellar nucleosynthesis with metallicity-dependent yields and also consider the effect of rotation on the yields massive stars and its consequences on the evolution of the Mg isotopes. We use a cosmic evolution model updated with new observational SFRD data and new reionization constraints coming from 2018 Planck collaboration determinations. We find that the main contribution of 24Mg comes from massive stars whereas 25Mg and 26Mg come from intermediate mass stars. To fit the observational data on magnesium isotopic ratios, an additional intermediate mass SFRD component is preferred. Moreover, the agreement between model and data is further improved when the range of IM masses is narrowed towards higher masses (5-8 M_\odot). While some rotation also improves the fit to data, we can exclude the case where all stars have high rotational velocities due to an over-production of 26Mg.


title: 1/N corrections to $F_1$ and $F_2$ structure functions of vector mesons from holography

abstract: The structure functions $F_1$ and $F_2$ of the hadronic tensor of vector mesons are obtained at order $1/N$ and strong coupling using the gauge/gravity duality. We find that the large $N$ limit and the high energy one do not commute. Thus, by considering the high energy limit first, our results of the first moments of $F_1$ for the rho meson agree well with those from lattice QCD, with an important improvement of the accuracy with respect to the holographic dual calculation in the planar limit.


title: Thermal radiation and dissipative phase transition in a BEC with local loss

abstract: We study the dynamics of an atomic BEC subject to local dissipation in the form of atom losses. We show there is a critical loss rate at which the system undergoes a continuous dissipative phase transition from a homogenous state into a state which contains a sonic horizon. The latter drastically alters the behavior of the system by screening the drain. Dissipation leads to two types of fluctuations. First, fluctuations are generated by particles emitted in the reservoir. Both above and below the critical loss, these result in thermal emission of phonons with a temperature set by the loss rate and the chemical potential. The second type of fluctuation results from scattering on the drain and gives rise to a particular correlation pattern that can be observed in the density-density correlation. Aside from correlations between in an out scattered modes, outgoing particles are correlated with localized modes through a process that is reminiscent of Hawking radiation. Finally, we briefly discuss the dynamics of the system when there are two drains, in which case it is possible to construct a black hole laser.


title: Identifying resonances with wave-packet dynamics

abstract: A new method for the study of resonant behavior - using wave-packet dynamics - is presented, based on the powerful window operator technique. The method is illustrated and quantified by application to the astrophysically-important example of low-energy $^{12}$C + $^{12}$C collisions. For this selected, potential model test case, the technique is shown to provide both resonance energies and widths in agreement with alternative methods, such as complex-energy scattering-matrix pole searches and scattering phase-shift analyses. The method has a more general capability to study resonance phenomena across disciplines, that involve particles temporarily trapped by potential pockets.


title: The Lattice of Profinite Subgroups of Protori

abstract: Compact connected abelian groups, or protori, have intrinsic structural characteristics that present for the entire category. In the case of finite-dimensional torus-free protori, The Resolution Theorem for Compact Abelian Groups sets the stage for demonstrating that the profinite subgroups inducing tori quotients comprise an isogeny class of finitely generated modules over the profinite integers, which is a lattice under intersection (meet) and + (join). The structural results enable the formulation of a universal resolution in the category of protori under morphisms of compact abelian groups. A single profinite subgroup from the lattice in the Resolution Theorem is replaced by the direct limit of the lattice of such subgroups and effects a covering morphism in which the discrete torsion-free Pontryagin dual of the protorus organically emerges as the kernel of the quotient map resolving the protorus. Among other advantages, this enables the study of a finite rank torsion-free abelian group as a canonical subgroup of its dual protorus, with the concomitant topological, analytical, and number-theoretic insights availed by the compact abelian setting.


title: Tuning nonthermal distributions to thermal ones in time-dependent Paul traps

abstract: We study the probability distribution of an atomic ion being laser-cooled in a periodically-driven Paul trap using a Floquet approach to the semiclassical photon scattering dynamics. We show that despite the microscopic nonequilibrium forces, a stationary thermal-like exponential distribution can be obtained in the Hamiltonian action, or equivalently in the number of quanta (phonons) of the motion linearized about the zero of the potential. At the presence of additional stray electric fields, the ion is pushed from the origin of the potential and set into a large-amplitude driven oscillation, and above a threshold amplitude of such "excess micromotion", the action distribution of excitations about the driven oscillation broadens and becomes distinctly nonthermal. We find that by a proper choice of the laser detuning the distribution can be made exponential again, with a mean phonon number close to that of the Doppler cooling limit. We derive a relation allowing to deduce just from the experimentally observable photon scattering rate both the required detuning for optimal cooling and the final mean phonon number. These results are important for quantum information processing and other applications, and in particular the derived approach can be applied to crystals of trapped ions in planar configurations, where the driven motion of ions is unavoidable.


title: Towards the optimal construction of a loss function without spurious local minima for solving quadratic equations

abstract: The problem of finding a vector $x$ which obeys a set of quadratic equations $|a_k^\top x|^2=y_k$, $k=1,\cdots,m$, plays an important role in many applications. In this paper we consider the case when both $x$ and $a_k$ are real-valued vectors of length $n$. A new loss function is constructed for this problem, which combines the smooth quadratic loss function with an activation function. Under the Gaussian measurement model, we establish that with high probability the target solution $x$ is the unique local minimizer (up to a global phase factor) of the new loss function provided $m\gtrsim n$. Moreover, the loss function always has a negative directional curvature around its saddle points.


title: Evolution of galactic magnetic fields

abstract: We study the cosmic evolution of the magnetic fields of a large sample of spiral galaxies in a cosmologically representative volume by employing a semi-analytic galaxy formation model and numerical dynamo solver in tandem. We start by deriving time- and radius-dependent galaxy properties using the galform galaxy formation model, which are then fed into the nonlinear mean-field dynamo equations. These are solved to give the large-scale (mean) field as a function of time and galactocentric radius for a thin disc, assuming axial symmetry. A simple prescription for the evolution of the small-scale (random) magnetic field component is also adopted. We find that, while most massive galaxies are predicted to have large-scale magnetic fields at redshift z=0, a significant fraction of them are expected to contain negligible large-scale field. Our model indicates that, for most of the galaxies containing large-scale magnetic fields today, the mean-field dynamo becomes active at z<3. We compute the radial profiles of pitch angle, and find broad agreement with observational data for nearby galaxies.


title: Consistency and Variation in Kernel Neural Ranking Model

abstract: This paper studies the consistency of the kernel-based neural ranking model K-NRM, a recent state-of-the-art neural IR model, which is important for reproducible research and deployment in the industry. We find that K-NRM has low variance on relevance-based metrics across experimental trials. In spite of this low variance in overall performance, different trials produce different document rankings for individual queries. The main source of variance in our experiments was found to be different latent matching patterns captured by K-NRM. In the IR-customized word embeddings learned by K-NRM, the query-document word pairs follow two different matching patterns that are equally effective, but align word pairs differently in the embedding space. The different latent matching patterns enable a simple yet effective approach to construct ensemble rankers, which improve K-NRM's effectiveness and generalization abilities.


title: Phenomenology of Relativistic $\mathbf{3} \to \mathbf{3}$ Reaction Amplitudes within the Isobar Approximation

abstract: Further progress in hadron spectroscopy necessitates the phenomenological description of three particle reactions. We consider the isobar approximation, where the connected part of the $\mathbf{3}\to\mathbf{3}$ amplitude is first expressed as a sum over initial and final pairs, and then expanded into a truncated partial wave series. The resulting unitarity equation is automatically fulfilled by the $B$-matrix solution, which is an integral equation for the partial wave amplitudes, analogous to the $K$-matrix parameterization used to describe $\mathbf{2}\to\mathbf{2}$ amplitudes. We study the one particle exchange and how its analytic structure impacts rescattering solutions such as the triangle diagram. The analytic structure is compared to other parameterizations discussed in the literature. We briefly discuss the analogies with recent formalisms for extracting $\mathbf{3}\to\mathbf{3}$ scattering amplitudes in lattice QCD.


title: Emergence of anomalous dynamics from the underlying singular continuous spectrum in interacting many-body systems

abstract: We investigate the dynamical properties of an interacting many-body system with a non-trivial energy potential landscape that may induce a singular continuous single-particle energy spectrum. Focusing on the Aubry-Andr\'e model, whose anomalous transport properties in presence of interaction has recently been demonstrated experimentally in an ultracold gas setup, we discuss the anomalous slowing down of the dynamics it exhibits and show that it emerges from the singular-continuous nature of the single-particle excitation spectrum. Our study demonstrates that singular-continuous spectra can be found in interacting systems, unlike previously conjectured by treating the interactions in the mean-field approximation. This, in turns, also highlights the importance of the many-body correlations in giving rise to anomalous dynamics, which, in many-body systems, can result from a non-trivial interplay between geometry and interactions.


title: Packing of Circles on Square Flat Torus as Global Optimization of Mixed Integer Nonlinear problem

abstract: The article demonstrates rather general approach to problems of discrete geometry: treat them as global optimization problems to be solved by one of general purpose solver implementing branch-and-bound algorithm (B&B). This approach may be used for various types of problems, i.e. Tammes problems, Thomson problems, search of minimal potential energy of micro-clusters, etc. Here we consider a problem of densest packing of equal circles in special geometrical object, so called square flat torus $\mathbb{R}^2/\mathbb{Z}^2$ with the induced metric. It is formulated as Mixed-Integer Nonlinear Problem with linear and non-convex quadratic constraints. The open-source B&B-solver SCIP, http://scip.zib.de, and its parallel implementation ParaSCIP, http://ug.zib.de, had been used in computing experiments to find "very good" approximations of optimal arrangements. The main result is a confirmation of the conjecture on optimal packing for N=9 that was published in 2012 by O. Musin and A. Nikitenko. To do that, ParaSCIP took about 2000 CPU*hours (16 hours x 128 CPUs) of cluster HPC4/HPC5, National Research Centre "Kurchatov Institute", http://ckp.nrcki.ru


title: Fast and Scalable Position-Based Layout Synthesis

abstract: The arrangement of objects into a layout can be challenging for non-experts, as is affirmed by the existence of interior design professionals. Recent research into the automation of this task has yielded methods that can synthesize layouts of objects respecting aesthetic and functional constraints that are non-linear and competing. These methods usually adopt a stochastic optimization scheme, which samples from different layout configurations, a process that is slow and inefficient. We introduce an physics-motivated, continuous layout synthesis technique, which results in a significant gain in speed and is readily scalable. We demonstrate our method on a variety of examples and show that it achieves results similar to conventional layout synthesis based on Markov chain Monte Carlo (McMC) state-search, but is faster by at least an order of magnitude and can handle layouts of unprecedented size as well as tightly-packed layouts that can overwhelm McMC.


title: Giant shot noise from Majorana zero modes in topological trijunctions

abstract: The clear-cut experimental identification of Majorana bound states in transport measurements still poses experimental challenges. We here show that the zero-energy Majorana state formed at a junction of three topological superconductor wires is directly responsible for giant shot noise amplitudes, in particular at low voltages and for small contact transparency. The only intrinsic noise limitation comes from the current-induced dephasing rate due to multiple Andreev reflection processes.


title: Conserving approximations in cavity quantum electrodynamics: Implications for density functional theory of electron-photon systems

abstract: By analyzing the many-body problem for non-relativistic electrons strongly coupled to photon modes of a microcavity I derive the exact momentum/force balance equation for cavity quantum electrodynamics. Implications of this equation for the electron self-energy and the exchange-correlation potential of quantum electrodynamic time-dependent density functional (QED-TDDFT) are discussed. In particular I generalize the concept of $\Phi$-derivability to construct approximations which ensure the correct momentum balance. It is shown that a recently proposed optimized effective potential approximation for QED-TDDFT is conserving and its possible improvements are discussed.


title: Transmission from reverse reaction coordinate mappings

abstract: We point out that the transport properties of non-interacting fermionic chains tunnel-coupled to two reservoirs at their ends can be mapped to those of a single quantum dot that is tunnel-coupled to two transformed reservoirs. The parameters of the chain are mapped to additional structure in the spectral densities of the transformed reservoirs. For example, this enables the calculation of the transmission of quantum dot chains by evaluating the known transmission of a single quantum dot together with structured spectral densities. We exemplify this analytically for short chains, which allows to optimize the transmission. In addition, we also demonstrate that the mapping can be performed numerically by computing the transmission of a Su-Schrieffer-Heeger chain.


title: Probability Measures and projections on Quantum Logics

abstract: The present paper is devoted to modelling of a probability measure of logical connectives on a quantum logic (QL), via a $G$-map, which is a special map on it. We follow the work in which the probability of logical conjunction, disjunction and symmetric difference and their negations for non-compatible propositions are studied. We study such a $ G $-map on quantum logics, which is a probability measure of a projection and show, that unlike classical (Boolean) logic, probability measure of projections on a quantum logic are not necessarilly pure projections. We compare properties of a $G$-map on QLs with properties of a probability measure related to logical connectives on a Boolean algebra.


title: Closest-Pair Queries in Fat Rectangles

abstract: In the range closest pair problem, we want to construct a data structure storing a set $S$ of $n$ points in the plane, such that for any axes-parallel query rectangle $R$, the closest pair in the set $R \cap S$ can be reported. The currently best result for this problem is by Xue et al.~(SoCG 2018). Their data structure has size $O(n \log^2 n)$ and query time $O(\log^2 n)$. We show that a data structure of size $O(n \log n)$ can be constructed in $O(n \log n)$ time, such that queries can be answered in $O(\log n + f \log f)$ time, where $f$ is the aspect ratio of $R$. Thus, for fat query rectangles, the query time is $O(\log n)$. This result is obtained by reducing the range closest pair problem to standard range searching problems on the points of $S$.


title: Investigating laser induced phase engineering in MoS2 transistors

abstract: Phase engineering of MoS2 transistors has recently been demonstrated and has led to record low contact resistances. The phase patterning of MoS2 flakes with laser radiation has also been realized via spectroscopic methods, which invites the potential of controlling the metallic and semiconducting phases of MoS2 transistors by simple light exposure. Nevertheless, the fabrication and demonstration of laser patterned MoS2 devices starting from the metallic polymorph has not been demonstrated yet. Here, we study the effects of laser radiation on 1T/1T'-MoS2 transistors with the prospect of driving an in-situ phase transition to the 2H-polymorph through light exposure. We find that although the Raman peaks of 2H-MoS2 become more prominent and the ones from the 1T/1T' phase fade after the laser exposure, the semiconducting properties of the laser patterned devices are not fully restored and the laser treatment ultimately leads to degradation of the transport channel.


title: Real Harmonic Analysis on the Special Orthogonal Group

abstract: This paper presents theoretical analysis and software implementation for real harmonics analysis on the special orthogonal group. Noncommutative harmonic analysis for complex-valued functions on the special orthogonal group has been studied extensively. However, it is customary to treat real harmonic analysis as a special case of complex harmonic analysis, and there have been limited results developed specifically for real-valued functions. Here, we develop a set of explicit formulas for real-valued irreducible unitary representations on the special orthogonal group, and provide several operational properties, such as derivatives, sampling, and Clebsch-Gordon coefficients. Furthermore, we implement both of complex and real harmonics analysis on the special orthogonal group into an open source software package that utilizes parallel processing through the OpenMP library. The efficacy of the presented results are illustrated by benchmark studies and an application to spherical shape matching.


title: On dimensional reduction of 4d N=1 Lagrangians for Argyres-Douglas theories

abstract: Recently, it was found that certain 4d $\mathcal{N}=1$ Lagrangians experience supersymmetry enhancement at their IR fixed point, thereby giving a Lagrangian description for a plethora of Argyres-Douglas theories. A generic feature of these Lagrangians is that a number of gauge invariant operators decouple (as free fields) along the RG-flow. These decoupled operators can be naturally taken into account from the beginning itself by introducing additional gauge singlets (sometimes called `flipping fields') that couple to the decoupled operators via appropriate superpotential terms. It has also been checked that upon dimensionally reducing to 3d, the $(A_1,A_{2n-1})$ type Lagrangians only produce the expected behavior when flipping fields are included in the Lagrangian. In this paper we further investigate the role of flipping fields and find an example where the expected necessity of including the flipping fields in the dimensionally reduced Lagrangians seems to get violated. In the process we find two new dual Lagrangians for the so called 3d $T[SU(2)]$ theory.


title: Physics Informed Topology Learning in Networks of Linear Dynamical Systems

abstract: Learning influence pathways of a network of dynamically related processes from observations is of considerable importance in many disciplines. In this article, influence networks of agents which interact dynamically via linear dependencies are considered. An algorithm for the reconstruction of the topology of interaction based on multivariate Wiener filtering is analyzed. It is shown that for a vast and important class of interactions, that respect flow conservation, the topology of the interactions can be exactly recovered. The class of problems where reconstruction is guaranteed to be exact includes power distribution networks, dynamic thermal networks and consensus networks. The efficacy of the approach is illustrated through simulation and experiments on consensus networks, IEEE power distribution networks and thermal dynamics of buildings.


title: Deep Learning and Holographic QCD

abstract: We apply the relation between deep learning (DL) and the AdS/CFT correspondence to a holographic model of QCD. Using a lattice QCD data of the chiral condensate at a finite temperature as our training data, the deep learning procedure holographically determines an emergent bulk metric as neural network weights. The emergent bulk metric is found to have both a black hole horizon and a finite-height IR wall, so shares both the confining and deconfining phases, signaling the cross-over thermal phase transition of QCD. In fact, a quark antiquark potential holographically calculated by the emergent bulk metric turns out to possess both the linear confining part and the Debye screening part, as is often observed in lattice QCD. From this we argue the discrepancy between the chiral symmetry breaking and the quark confinement in the holographic QCD. The DL method is shown to provide a novel data-driven holographic modeling of QCD, and sheds light on the mechanism of emergence of the bulk geometries in the AdS/CFT correspondence.


title: Cross-Layer Effects on Training Neural Algorithms for Video Streaming

abstract: Nowadays Dynamic Adaptive Streaming over HTTP (DASH) is the most prevalent solution on the Internet for multimedia streaming and responsible for the majority of global traffic. DASH uses adaptive bit rate (ABR) algorithms, which select the video quality considering performance metrics such as throughput and playout buffer level. Pensieve is a system that allows to train ABR algorithms using reinforcement learning within a simulated network environment and is outperforming existing approaches in terms of achieved performance. In this paper, we demonstrate that the performance of the trained ABR algorithms depends on the implementation of the simulated environment used to train the neural network. We also show that the used congestion control algorithm impacts the algorithms' performance due to cross-layer effects.


title: Model-free Study of Ordinary Least Squares Linear Regression

abstract: Ordinary least squares (OLS) linear regression is one of the most basic statistical techniques for data analysis. In the main stream literature and the statistical education, the study of linear regression is typically restricted to the case where the covariates are fixed, errors are mean zero Gaussians with variance independent of the (fixed) covariates. Even though OLS has been studied under misspecification from as early as the 1960's, the implications have not yet caught up with the main stream literature and applied sciences. The present article is an attempt at a unified viewpoint that makes the various implications of misspecification stand out.


title: Theories of truth for countable languages which conform to classical logic

abstract: Every countable language which conforms to classical logic is shown to have an extension which conforms to classical logic, and has a definitional theory of truth. That extension has a semantical theory of truth, if every sentence of the object language is valuated by its meaning either as true or as false. These theories contain both a truth predicate and a non-truth predicate. Theories are equivalent when the sentences of the object language are valuated by their meanings.


title: Identifying Long Term Voltage Stability Caused by Distribution Systems vs Transmission Systems

abstract: Monitoring the long term voltage stability of the power grid is necessary to ensure its secure operation. This paper presents a new phasor based methodology that distinguishes between long term voltage stability caused by distribution systems versus transmission systems. From a conceptual understanding of a simplified system, a Transmission-Distribution Distinguishing Index (TDDI) is proposed to distinguish between the two scenarios. A methodology to calculate the TDDI for multi-bus systems using quasi-steady state phasor measurements is described and validating results are presented for the IEEE 9 Bus system with a load replaced by various distribution feeders. The results verify that the TDDI can indeed be used to distinguish between transmission limited and distribution limited systems. This information can be utilized by the operator to effectively choose controls in distribution and transmission systems to improve the system margin.


title: Novel Sparse Recovery Algorithms for 3D Debris Localization using Rotating Point Spread Function Imagery

abstract: An optical imager that exploits off-center image rotation to encode both the lateral and depth coordinates of point sources in a single snapshot can perform 3D localization and tracking of space debris. When actively illuminated, unresolved space debris, which can be regarded as a swarm of point sources, can scatter a fraction of laser irradiance back into the imaging sensor. Determining the source locations and fluxes is a large-scale sparse 3D inverse problem, for which we have developed efficient and effective algorithms based on sparse recovery using non-convex optimization. Numerical simulations illustrate the efficiency and stability of the algorithms.


title: Towards a classification of Lindenmayer systems

abstract: In this paper we will attempt to classify Lindenmayer systems based on properties of sets of rules and the kind of strings those rules generate. This classification will be referred to as a parametrization of the L-space: the L-space is the phase space in which all possible L-developments are represented. This space is infinite, because there is no halting algorithm for L-grammars; but it is also subjected to hard conditions, because there are grammars and developments which are not possible states of an L-system: a very well-known example is the space of normal grammars. Just as the space of normal grammars is parametrized into Regular, Context-Free, Context-Sensitive, and Unrestricted (with proper containment relations holding among them; see Chomsky, 1959: Theorem 1), we contend here that the L-space is a very rich landscape of grammars which cluster into kinds that are not mutually translatable.


title: Anisotropic thermal conductivity measurement using a new Asymmetric-Beam Time-Domain Thermoreflectance (AB-TDTR) method

abstract: Anisotropic thermal properties are of both fundamental and practical interests, but remain challenging to characterize using conventional methods. In this work, a new metrology based on asymmetric beam time-domain thermoreflectance (AB-TDTR) is developed to measure three-dimensional anisotropic thermal transport by extending the conventional TDTR technique. Using an elliptical laser beam with controlled elliptical ratio and spot size, the experimental signals can be exploited to be dominantly sensitive to measure thermal conductivity along the cross-plane or any specific in-plane directions. An analytic solution for a multi-layer system is derived for the AB-TDTR signal in response to the periodical pulse, elliptical laser beam, and heating geometry, to extract the anisotropic thermal conductivity from experimental measurement. Examples with experimental data are given for various materials with in-plane thermal conductivity from 5 W/mK to 2000W/mK, including isotropic materials (silicon, boron phosphide, boron nitride), transversely isotropic materials (graphite, quartz, sapphire) and transversely anisotropic materials (black phosphorus). Furthermore, a detailed sensitivity analysis is conducted to guide the optimal setting of experimental configurations for different materials. The developed AB-TDTR metrology provides a new approach to accurately measure anisotropic thermal phenomena for rational materials design and thermal applications.


title: On the asymptotic stability of the time--fractional Lengyel--Epstein system

abstract: This paper concerns a time fractional version of the conventional Lengyel--Epstein CIMA reaction model. We define the invariant regions of the system and establish sufficient conditions for the unique equilibrium's local and global asymptotic stability. Numerical results are presented to illustrate the effect of the fractional order on system dynamics.


title: Convergence rate of Markov chains and hybrid numerical schemes to jump-diffusions with application to the Bates model

abstract: We study the rate of weak convergence of Markov chains to diffusion processes under suitable but quite general assumptions. We give an example in the financial framework, applying the convergence analysis to a multiple jumps tree approximation of the CIR process. Then, we combine the Markov chain approach with other numerical techniques in order to handle the different components in jump-diffusion coupled models. We study the speed of convergence of this hybrid approach and we provide an example in finance, applying our results to a tree-finite difference approximation in the Heston or Bates model.


title: Chaotic mixing using micro-rotors in a confined domain

abstract: In this work we study chaotic mixing induced by point micro-rotors in a bounded two dimensional Stokes flow. The dynamics of the pair of rotors, modeled as rotlets, are non Hamiltonian in the bounded domain and produce chaotic advection of fluid tracers in subsets of the domain. A complete parametric investigation of the fluid mixing as a function of the initial locations of the rotlets is performed based on pseudo phase portraits. The mixing of fluid tracers as a function of relative positions of micro-rotors is studied using finite time entropy and locational entropy. The finite time locational entropy is used to identify regions of the fluid that produce good vs poor mixing and this is visualized by the stretching and folding of blobs of tracer particles. Unlike the case of the classic blinking vortex dynamics, the velocity field of the flow modeled using rotlets inside a circular boundary is smooth in time and satisfies the no slip boundary condition. This makes the considered model a more realistic case for studies of mixing in microfluidic devices using magnetic actuated microspheres.


title: Modification of polarization through de-Gaussification

abstract: We analyze the polarization of a quantum radiation field under a de-Gaussification process. Specifically, we consider the addition of photons to a two-mode thermal state to get mixed non-Gaussian and nonclassical states which are still diagonal in the Fock basis. Stokes-operator-based degrees of polarization and two distance-type measures defined with Hilbert-Schmidt and Bures metrics are investigated. For a better insight we here introduce a polarization degree based on the relative entropy. Polarization of the thermal states is fully investigated and simple closed expressions are found for all the defined degrees. The evaluated degrees for photon-added states are then compared to the corresponding ones for the two-mode thermal states they originate from. We present interesting findings which tell us that some popular degrees of polarization are not fully consistent. However, the most solidly defined degrees, which are based on the Bures metric and relative entropy, clearly indicate an enhancement of polarization through de-Gaussification. This conclusion is supported by the behavior of the degrees of polarization of the Fock states, which are finally discussed as a limit case.


title: Real-time 3D Pose Estimation with a Monocular Camera Using Deep Learning and Object Priors On an Autonomous Racecar

abstract: We propose a complete pipeline that allows object detection and simultaneously estimate the pose of these multiple object instances using just a single image. A novel "keypoint regression" scheme with a cross-ratio term is introduced that exploits prior information about the object's shape and size to regress and find specific feature points. Further, a priori 3D information about the object is used to match 2D-3D correspondences and accurately estimate object positions up to a distance of 15m. A detailed discussion of the results and an in-depth analysis of the pipeline is presented. The pipeline runs efficiently on a low-powered Jetson TX2 and is deployed as part of the perception pipeline on a real-time autonomous vehicle cruising at a top speed of 54 km/hr.


title: Rabi oscillations of a quantum dot exciton coupled to acoustic phonons: coherence and population readout

abstract: While the advanced coherent control of qubits is now routinely carried out in low frequency (GHz) systems like single spins, it is far more challenging to achieve for two-level systems in the optical domain. This is because the latter evolve typically in the THz range, calling for tools of ultrafast, coherent, nonlinear optics. Using four-wave mixing micro-spectroscopy, we here measure the optically driven dynamics of a single exciton quantum state confined in a semiconductor quantum dot. In a combined experimental and theoretical approach, we reveal the intrinsic Rabi oscillation dynamics by monitoring both central exciton quantities, i.e., its occupation and the microscopic coherence, as resolved by the four-wave mixing technique. In the frequency domain this oscillation generates the Autler-Townes splitting of the light-exciton dressed states, directly seen in the four-wave mixing spectra. We further demonstrate that the coupling to acoustic phonons strongly influences the FWM dynamics on the picosecond timescale, because it leads to transitions between the dressed states.


title: On the speciality of Tortkara algebras

abstract: A criterion for elements of free Zinbiel algebras to be Lie or Jordan is established. This criterion is used in studying speciality problems of Tortkara algebras. We construct a base of free special Tortkara algebras. Furthermore, we prove analogues of classical Cohn's and Shirshov's theorems for Tortkara algebras.


title: Superconducting spin properties of Majorana nanowires and the associated superconducting anomalous Hall effect

abstract: It is difficult to unambiguously confirm the existence of Majorana zero modes (MZMs) due to the absence of smoking-gun signatures in charge transport measurements. Recent studies suggest that the spin degree of freedom of MZMs may provide an alternative detection method. We study the spin properties of the superconducting state in Majorana nanowires and the associated unconventional Josephson effect with realistic experimental parameters taken from [Zhang {\it et al.}, Nature 556, 74 (2018)]. For a superconducting thin film with in-plane polarized spin-triplet pairing, an out-of-plane electric field can generate a supercurrent perpendicular to both the superconducting spin polarization and the electric field, so we name this phenomena as superconducting anomalous Hall effect (ScAHE). In a Majorana nanowire, the regime with finite polarized spin-triplet pairing almost coincides with the chiral regime, which includes the topological regime. We further study the effects of polarized spin-triplet pairing in two types of Josephson junctions. One dramatic finding is that SOC can induce an anomalous supercurrent at zero phase difference only in the U-shape junction, a basic ingredient of scalable topological quantum computation. This can be viewed as a consequence of the ScAHE. Our work reveals that the spin degree of freedom is indeed helpful for detecting MZMs.


title: A Mechanical Approach to One-Dimensional Interacting Gas

abstract: Traditional derivations of the van der Waals equation typically use standard recipes involving ensemble averages of statistical mechanics. In this work, we study a box of weakly interacting gas particles in one-dimension from a purely mechanical point of view. This has the merit that it not only reproduces the van der Waals equation but also tells us some extra interesting physics not immediately clear from a pure statistical mechanical approach. For example, we find that the traditional handwaving interpretation of the van der Waals equation adopting mean field approximation is actually incorrect. In this investigation of one-dimensional interacting gas, we demonstrate the possibility taking a mechanical point of view and having deeper understanding for the physics of leading order effect of particle-particle interaction, for weakly interacting N-body systems that are usually studied in the framework of statistical mechanics or kinetic theory.


title: A beam-beam monitoring detector for the MPD experiment at NICA

abstract: The Multi-Purpose Detector (MPD) is to be installed at the Nuclotron Ion Collider fAcility (NICA) of the Joint Institute for Nuclear Research (JINR). Its main goal is to study the phase diagram of the strongly interacting matter produced in heavy-ion collisions. These studies, while providing insight into the physics of heavy-ion collisions, are relevant for improving our understanding of the evolution of the early Universe and the formation of neutron stars. In order to extend the MPD trigger capabilities, we propose to include a high granularity beam-beam monitoring detector (BE-BE) to provide a level-0 trigger signal with an expected time resolution of 30 ps. This new detector will improve the determination of the reaction plane by the MPD experiment, a key measurement for flow studies that provides physics insight into the early stages of the reaction. In this work, we use simulated Au+Au collisions at NICA energies to show the potential of such a detector to determine the event plane resolution, providing further redundancy to the detectors originally considered for this purpose namely, the Fast Forward Detector (FFD) and the Hadron Calorimeter (HCAL). We also show our results for the time resolution studies of two prototype cells carried out at the T10 beam line at the CERN PS complex.


title: An Adaptive Tabu Search Algorithm for Market Clearing Problem in Turkish Day-Ahead Market

abstract: In this study, we focus on the market clearing problem of Turkish day-ahead electricity market. We propose a mathematical model by extending the variety of bid types for different price regions. The commercial solvers may not find any feasible solution for the proposed problem in some instances within the given time limits. Hence, we design an adaptive tabu search (ATS) algorithm to solve the problem. ATS discretizes continuous search space arising from the flow variables. Our method has adaptive radius and it achieves backtracking by a commercial solver. Then, we compare the performance of ATS with a heuristic decomposition method from the literature by using synthetic data sets. We evaluate the performances of the algorithms with respect to their solution times and surplus differences. ATS performs better in most of the sets.


title: Conceptual design study for heat exhaust management in the ARC fusion pilot plant

abstract: The ARC pilot plant conceptual design study has been extended beyond its initial scope [B. N. Sorbom et al., FED 100 (2015) 378] to explore options for managing ~525 MW of fusion power generated in a compact, high field (B_0 = 9.2 T) tokamak that is approximately the size of JET (R_0 = 3.3 m). Taking advantage of ARC's novel design - demountable high temperature superconductor toroidal field (TF) magnets, poloidal magnetic field coils located inside the TF, and vacuum vessel (VV) immersed in molten salt FLiBe blanket - this follow-on study has identified innovative and potentially robust power exhaust management solutions.


title: Absence of Buckling in Nerve Fiber

abstract: In this study we give a geometrical model which employs the smoothness of nerve fibers as differentiable curves. We show that a nerve fiber may encounter large curvature due to the possible helicial bending and hence it could cause the fiber to buckle. However, its membrane structure provides a mechanism, entirely geometrical to avoid it. To overcome the challenge of emerging helix we project it into a plane.


title: Locally Polarized Wave Propagation through Metamaterials' Crystallinity

abstract: Wave propagation control is of fundamental interest in many areas of Physics. It can be achieved with wavelength-scaled photonic crystals, hence avoiding low frequency applications. By contrast, metamaterials are structured on a deep-subwavelength scale, and therefore usually described through homogenization, neglecting the unit-cell structuration. Here, we show with microwaves that, by considering their inherent crystallinity, we can induce wave propagation carrying angular momenta within a subwavelength-scaled collection of wires. Then, inspired by the Quantum Valley-Hall Effect in condensed matter physics, we exploit this bulk circular polarization to create modes propagating along particular interfaces. The latter also carry an edge angular momentum whose conservation during the propagation allows wave routing by design in specific directions. This experimental study not only evidences that crystalline metamaterials are a straightforward tabletop platform to emulate exciting solid-state physics phenomena at the macroscopic scale, but it also opens the door to crystalline polarized subwavelength waveguides.


title: Left-handed Band in an Electromagnetic Metamaterial Inducedby Sub-wavelength Multiple Scattering

abstract: Due to the deep sub-wavelength unit cell in metamaterials, the quasi-static approximation is usually employed to describe the propagation. By making pairs of resonators, we highlight that multiple scattering also occurs at this scale and results in the existence of a dipolar resonance, which leads to a negative index of refraction when we consider several resonators. We experimentally verify the possibility of obtaining a negative index of refraction in periodic metamaterials in two different ways, and eventually demonstrates a superlensing effect in both cases.


title: Obladi: Oblivious Serializable Transactions in the Cloud

abstract: This paper presents the design and implementation of Obladi, the first system to provide ACID transactions while also hiding access patterns. Obladi uses as its building block oblivious RAM, but turns the demands of supporting transactions into a performance opportunity. By executing transactions within epochs and delaying commit decisions until an epoch ends, Obladi reduces the amortized bandwidth costs of oblivious storage and increases overall system throughput. These performance gains, combined with new oblivious mechanisms for concurrency control and recovery, allow Obladi to execute OLTP workloads with reasonable throughput: it comes within 5x to 12x of a non-oblivious baseline on the TPC-C, SmallBank, and FreeHealth applications. Latency overheads, however, are higher (70x on TPC-C).


title: Modeling and Loop Shaping of Single-Joint Amplification Exoskeleton with Contact Sensing and Series Elastic Actuation

abstract: In this paper we consider a class of exoskeletons designed to amplify the strength of humans through feedback of sensed human-robot interactions and actuator forces. We define an amplification error signal based on a reference amplification rate, and design a linear feedback compensator to attenuate this error. Since the human operator is an integral part of the system, we design the compensator to be robust to both a realistic variation in human impedance and a large variation in load impedance. We demonstrate our strategy on a one-degree of freedom amplification exoskeleton connected to a human arm, following a three dimensional matrix of experimentation: slow or fast human motion; light or extreme exoskeleton load; and soft or clenched human arm impedances. We demonstrate that a slightly aggressive controller results in a borderline stable system---but only for soft human musculoeskeletal behavior and a heavy load. This class of exoskeleton systems is interesting because it can both amplify a human's interaction forces --- so long as the human contacts the environment through the exoskeleton --- and attenuate the operator's perception of the exoskeleton's reflected dynamics at frequencies within the bandwidth of the control.


title: The first observation of effect of oscillation in Neutrino-4 experiment on search for sterile neutrino (continuation)

abstract: We present the results of the Neutrino-4 experiment on search for a sterile neutrino. The experiment has been carried out on the SM-3 reactor having a compact active zone of $42\times42\times35\textrm{cm}^3$ and operating on the highly enriched uranium-235 at 90 MW thermal power. We report the results of the Neutrino-4 experiment of measurements of reactor antineutrino flux and spectrum dependence on the distance in the range 6-12 meters from the center of the reactor core. Using the measured spectrum and the distance dependence of antineutrino flux, we performed the model independent analysis of restrictions on the oscillation parameters $\Delta m^2_{14}$ and $\sin^2 2\theta_{14}$. The method of coherent addition of results of measurements is proposed. It allows us to directly observe the effect of oscillations. We observed the oscillation effect at CL $3.5\sigma$ in the vicinity of $\Delta m^2_{14} \approx 7.26\textrm{eV}^2$ and $\sin^2 2\theta_{14} \approx 0.38$. Combining the result of the Neutrino-4 experiment and the result of the gallium anomaly effect we obtained value $\sin^2 2\theta_{14} \approx 0.35 \pm 0.07 (5\sigma)$. The analysis of systematics effects is presented. Comparison with results of other experiments is presented. Future prospect of the experiment is discussed. It is necessary to notice that obtained values $\sin^2 2\theta_{14} \approx 0.35 \pm 0.07 (5\sigma)$ and $\Delta m^2_{14} \approx (7.3 \pm 0.7)\textrm{eV}^2$ allow make assessment on the mass of a neutrino: $m_{\beta} \approx 0.8\textrm{eV}$.


title: Dropout Distillation for Efficiently Estimating Model Confidence

abstract: We propose an efficient way to output better calibrated uncertainty scores from neural networks. The Distilled Dropout Network (DDN) makes standard (non-Bayesian) neural networks more introspective by adding a new training loss which prevents them from being overconfident. Our method is more efficient than Bayesian neural networks or model ensembles which, despite providing more reliable uncertainty scores, are more cumbersome to train and slower to test. We evaluate DDN on the the task of image classification on the CIFAR-10 dataset and show that our calibration results are competitive even when compared to 100 Monte Carlo samples from a dropout network while they also increase the classification accuracy. We also propose better calibration within the state of the art Faster R-CNN object detection framework and show, using the COCO dataset, that DDN helps train better calibrated object detectors.


title: Two-field Cosmological $\alpha$-attractors with Noether Symmetry

abstract: We study Noether symmetries in two-field cosmological $\alpha$-attractors, investigating the case when the scalar manifold is an elementary hyperbolic surface. This encompasses and generalizes the case of the Poincare disk. We solve the conditions for the existence of a `separated' Noether symmetry and find the form of the scalar potential compatible with such, for any elementary hyperbolic surface. For this class of symmetries, we find that the $\alpha$-parameter must have a fixed value. Using those Noether symmetries, we also obtain many exact solutions of the equations of motion of these models, which were studied previously with numerical methods.


title: Visualization of correlations in hybrid quantum systems

abstract: In this work we construct Wigner functions for hybrid continuous and discrete variable quantum systems. We demonstrate new capabilities in the visualization of the interactions and correlations within hybrid quantum systems. Specifically, we show how to clearly distinguish signatures that arise due to quantum and classical correlations in an entangled Bell-cat state. We further show how correlations are manifested in different types of interaction, leading to a deeper understanding of how quantum information is shared between two subsystems. Understanding the nature of the correlations between systems is central to harnessing quantum effects for information processing; the methods presented here reveal the nature of these correlations, allowing a clear visualization of the quantum information present in hybrid quantum systems. The methods presented here could be viewed as a form of quantum state spectroscopy.


title: A novel active learning framework for classification: using weighted rank aggregation to achieve multiple query criteria

abstract: Multiple query criteria active learning (MQCAL) methods have a higher potential performance than conventional active learning methods in which only one criterion is deployed for sample selection. A central issue related to MQCAL methods concerns the development of an integration criteria strategy (ICS) that makes full use of all criteria. The conventional ICS adopted in relevant research all facilitate the desired effects, but several limitations still must be addressed. For instance, some of the strategies are not sufficiently scalable during the design process, and the number and type of criteria involved are dictated. Thus, it is challenging for the user to integrate other criteria into the original process unless modifications are made to the algorithm. Other strategies are too dependent on empirical parameters, which can only be acquired by experience or cross-validation and thus lack generality; additionally, these strategies are counter to the intention of active learning, as samples need to be labeled in the validation set before the active learning process can begin. To address these limitations, we propose a novel MQCAL method for classification tasks that employs a third strategy via weighted rank aggregation. The proposed method serves as a heuristic means to select high-value samples of high scalability and generality and is implemented through a three-step process: (1) the transformation of the sample selection to sample ranking and scoring, (2) the computation of the self-adaptive weights of each criterion, and (3) the weighted aggregation of each sample rank list. Ultimately, the sample at the top of the aggregated ranking list is the most comprehensively valuable and must be labeled. Several experiments generating 257 wins, 194 ties and 49 losses against other state-of-the-art MQCALs are conducted to verify that the proposed method can achieve superior results.


title: Some Nontrivial Properties of a Formula for Compound Interest

abstract: We analyze the classical model of compound interest with a constant per-period payment and interest rate. We examine the outstanding balance function as well as the periodic payment function and show that the outstanding balance function is not generally concave in the interest rate, but instead may be initially convex on its domain and then concave.

categories: econ.GN q-fin.EC q-fin.GN q-fin.MF q-fin.PM


title: An Adaptive Algorithm Employing Continuous Linear Functionals

abstract: Automatic algorithms attempt to provide approximate solutions that differ from exact solutions by no more than a user-specified error tolerance. This paper describes an automatic, adaptive algorithm for approximating the solution to a general linear problem on Hilbert spaces. The algorithm employs continuous linear functionals of the input function, specifically Fourier coefficients. We assume that the Fourier coefficients of the solution decay sufficiently fast, but do not require the decay rate to be known a priori. We also assume that the Fourier coefficients decay steadily, although not necessarily monotonically. Under these assumptions, our adaptive algorithm is shown to produce an approximate solution satisfying the desired error tolerance, without prior knowledge of the norm of the function to be approximated. Moreover, the computational cost of our algorithm is shown to be essentially no worse than that of the optimal algorithm. We provide a numerical experiment to illustrate our algorithm.


title: FEM models and Fourier decomposition of three thick septa

abstract: The CEBAF recirculating linear accelerator uses 26 septa: 8 thin, 17 thick and one Lambertson. We report here on detailed FEM models of one thick current-sheet septum, the Lambertson, and a new thick-septum concept which contains features of both. The last has higher-order field components within a factor of two of the current-sheet septum and much lower current density in the conductor were it used as a drop-in replacement in CEBAF. There are accelerator layouts in which the new concept is superior to the current sheet septum in both harmonic content and current density.


title: Irreducible unirational and uniruled components of moduli spaces of polarized Enriques surfaces

abstract: We give an explicit description of the irreducible components of the moduli spaces of polarized Enriques surfaces in terms of decompositions of the polarization as an effective sum of isotropic classes. We prove that infinitely many of these components are unirational (resp. uniruled). In particular, this applies to components of arbitrarily large genus $g$ and $\phi$-invariant of the polarization.


title: Polarized laser-wakefield-accelerated kiloampere electron beams

abstract: High-flux polarized particle beams are of critical importance for the investigation of spin-dependent processes, such as in searches of physics beyond the Standard Model, as well as for scrutinizing the structure of solids and surfaces in material science. Here we demonstrate that kiloampere polarized electron beams can be produced via laser-wakefield acceleration from a gas target. A simple theoretical model for determining the electron beam polarization is presented and supported with self-consistent three-dimensional particle-in-cell simulations that incorporate the spin dynamics. By appropriately choosing the laser and gas parameters, we show that the depolarization of electrons induced by the laser-wakefield-acceleration process can be as low as 10%. Compared to currently available sources of polarized electron beams, the flux is increased by four orders of magnitude.


title: Weak order and descents for monotone triangles

abstract: Monotone triangles are a rich extension of permutations that biject with alternating sign matrices. The notions of weak order and descent sets for permutations are generalized here to monotone triangles, and shown to enjoy many analogous properties. It is shown that any linear extension of the weak order gives rise to a shelling order on a poset, recently introduced by Terwilliger, whose maximal chains biject with monotone triangles; among these shellings are a family of EL-shellings. The weak order turns out to encode an action of the 0-Hecke monoid of type A on the monotone triangles, generalizing the usual bubble-sorting action on permutations. It also leads to a notion of descent set for monotone triangles, having another natural property: the surjective algebra map from the Malvenuto- Reutenauer Hopf algebra of permutations into quasisymmetric functions extends in a natural way to an algebra map out of the recently-defined Cheballah-Giraudo-Maurice algebra of alternating sign matrices.


title: Scalar Arithmetic Multiple Data: Customizable Precision for Deep Neural Networks

abstract: Quantization of weights and activations in Deep Neural Networks (DNNs) is a powerful technique for network compression, and has enjoyed significant attention and success. However, much of the inference-time benefit of quantization is accessible only through the use of customized hardware accelerators or by providing an FPGA implementation of quantized arithmetic. Building on prior work, we show how to construct arbitrary bit-precise signed and unsigned integer operations using a software technique which logically \emph{embeds} a vector architecture with custom bit-width lanes in universally available fixed-width scalar arithmetic. We evaluate our approach on a high-end Intel Haswell processor, and an embedded ARM processor. Our approach yields very fast implementations of bit-precise custom DNN operations, which often match or exceed the performance of operations quantized to the sizes supported in native arithmetic. At the strongest level of quantization, our approach yields a maximum speedup of $\thicksim6\times$ on the Intel platform, and $\thicksim10\times$ on the ARM platform versus quantization to native 8-bit integers.


title: Interpolative Decomposition Butterfly Factorization

abstract: This paper introduces a "kernel-independent" interpolative decomposition butterfly factorization (IDBF) as a data-sparse approximation for matrices that satisfy a complementary low-rank property. The IDBF can be constructed in $O(N\log N)$ operations for an $N\times N$ matrix via hierarchical interpolative decompositions (IDs), if matrix entries can be sampled individually and each sample takes $O(1)$ operations. The resulting factorization is a product of $O(\log N)$ sparse matrices, each with $O(N)$ non-zero entries. Hence, it can be applied to a vector rapidly in $O(N\log N)$ operations. IDBF is a general framework for nearly optimal fast matvec useful in a wide range of applications, e.g., special function transformation, Fourier integral operators, high-frequency wave computation. Numerical results are provided to demonstrate the effectiveness of the butterfly factorization and its construction algorithms.


title: Generalized Whitney Topologies are Baire

abstract: In this paper we show that certain generalizations of the $C^r$-Whitney topology, which include the H\"older-Whitney and Sobolev-Whitney topologies on smooth manifolds, satisfy the Baire property, to wit, the countable intersection of open and dense sets is dense.


title: Topological Insulator in an Atomic Liquid

abstract: We demonstrate theoretically an atomic liquid phase that supports topologically nontrivial electronic structure. A minimum two-orbital model of liquid topological insulator in two dimensions is constructed within the framework of tight-binding molecular dynamics. As temperature approaches zero, our simulations show that the atoms crystallize into a triangular lattice with nontrivial band topology at high densities. Thermal fluctuations at finite temperatures melt the lattice, giving rise to a liquid state which inherits the nontrivial topology from the crystalline phase. The electronic structure of the resultant atomic liquid is characterized by a nonzero Bott index. Our work broadens the notion of topological materials, and points to a new systematic approach for searching topological phases in amorphous and liquid systems.


title: Ultrafast holography enabled by quantum interference of ultrashort electrons

abstract: Holography relies on the interference between a known reference and a signal of interest to reconstruct both the amplitude and phase of that signal. Commonly performed with photons and electrons, it finds numerous applications in imaging, cryptography and arts. With electrons, the extension of holography to the ultrafast time domain remains a challenge, although it would yield the highest possible combined spatio-temporal resolution. Here, we show that holograms of local electromagnetic fields can be obtained with combined attosecond/nanometer resolution in an ultrafast transmission electron microscope (UEM). Unlike conventional holography, where the signal and the reference are spatially separated and then recombined to interfere, in our method we use electromagnetic fields to split an electron wave function in a quantum coherent superposition of different energy states. In the image plane, spatial modulation of the electron-energy distribution reflects the phase relation between reference and signal fields, which we map via energy-filtered UEM. Beyond imaging applications, this approach allows implementing optically-controlled and spatially-resolved quantum measurements in parallel, providing an efficient and versatile tool for the exploration of electron quantum optics.

categories: cond-mat.mes-hall physics.acc-ph physics.optics quant-ph


title: $0\nu\beta\beta$ in left-right theories with Higgs doublets and gauge coupling unification

abstract: We consider a version of Left-Right Symmetric Model in which the scalar sector consists of a Higgs bidoublet ($\Phi$) with $B-L=0$, Higgs doublets ($H_{L,R}$) with $B-L=1$ and a charged scalar ($\delta^+$) with $B-L=2$ leading to radiatively generated Majorana masses for neutrinos and thereby, leads to new physics contributions to neutrinoless double beta decay ($0\nu \beta \beta$). We show that such a novel framework can be embedded in a non-SUSY $SO(10)$ GUT leading to successful gauge coupling unification at around $10^{16}$ GeV with the scale of left-right symmetry breaking around $10^{10}$ GeV. The model can also be extended to have left-right symmetry breaking at TeV scale, enabling detection of $W_R, Z_R$ bosons in LHC and future collider searches. In the context of neutrinoless double beta decay, this model can saturate the present bound from GERDA and KamLAND-Zen experiments. Also, we briefly explain how keV-MeV range RH neutrino arising from our model can saturate various astrophysical and cosmological constraints and can be considered as warm Dark Matter (DM) candidate to address various cosmological issues. We also discuss on left-right theories with Higgs doublets without having scalar bidoublet leading to fermion masses and mixings by inclusion of vector like fermions.


title: Reoptimization of Parameterized Problems

abstract: Parameterized complexity allows us to analyze the time complexity of problems with respect to a natural parameter depending on the problem. Reoptimization looks for solutions or approximations for problem instances when given solutions to neighboring instances. We try to combine both techniques, in order to better classify the complexity of problems in the parameterized setting. Specifically, we see that some problems in the class of compositional problems, which do not have polynomial kernels under standard complexity-theoretic assumptions, do have polynomial kernels under reoptimization for some local modifications. Moreover, we find that the reoptimization version of Vertex Cover has a polynomial kernel of size 2k using crown decomposition. Finally, in a negative result, we prove that the reoptimization version of Connected Vertex Cover does not have a Turing kernelization unless Set Cover has a polynomial kernel


title: On the Impact of Spillover Losses in 28 GHz Rotman Lens Arrays for 5G Applications

abstract: This work demonstrates the sensitivity of lens antenna arrays operating at millimeter-wave (mmWave) frequencies. Considering a Rotman lens array in receive mode, our investigation focuses on its most imperative defect: aberration of electromagnetic (EM) energy. Aberration leads to spillover of electric fields to neighboring ports, reducing the lens' ability to focus the EM energy to a desired port. With full EM simulations, we design a 28 GHz, 13 beam and 13 array port Rotman lens array to characterize its performance with the aforementioned impairment. Our findings show that the impact of aberration is more pronounced when the beam angles are close to the array end-fire. More critically, the corresponding impact of aberration on the desired signal and interference powers is also investigated for an uplink multiuser cellular system operating at 28 GHz. The presented results can be used as a reference to re-calibrate our expectations for Rotman lens arrays at mmWave frequencies.


title: Higher-Spin Quartic Vertices in AdS5 : a Rectangular Limit

abstract: We determine the form of a significantly large class of gauge-invariant quartic vertices for symmetric higher-spin fields (in Vasiliev formalism) in AdS5 space in the rectangular limit, (defined in the paper) using the higher-spin vertex operators in superstring theory, which construction is based on the generators of the higher-spin symmetry algebra, enveloping the AdS isometries. In this limit, a particular simplification is that the quartic vertices do not receive corrections from worldsheet renormalization group flows of the cubic terms, so their structure is fixed directly by the 4-point amplitudes in superstring theory.


title: Acoustic Probing for Estimating the Storage Time and Firmness of Tomatoes and Mandarin Oranges

abstract: This paper introduces an acoustic probing technique to estimate the storage time and firmness of fruits; we emit an acoustic signal to fruit from a small speaker and capture the reflected signal with a tiny microphone. We collect reflected signals for fruits with various storage times and firmness conditions, using them to train regressors for estimation. To evaluate the feasibility of our acoustic probing, we performed experiments; we prepared 162 tomatoes and 153 mandarin oranges, collected their reflected signals using our developed device and measured their firmness with a fruit firmness tester, for a period of 35 days for tomatoes and 60 days for mandarin oranges. We performed cross validation by using this data set. The average estimation errors of storage time and firmness for tomatoes were 0.89 days and 9.47 g/mm2. Those for mandarin oranges were 1.67 days and 15.67 g/mm2. The estimation of storage time was sufficiently accurate for casual users to select fruits in their favorite condition at home. In the experiments, we tested four different acoustic probes and found that sweep signals provide highly accurate estimation results.


title: Kernel based low-rank sparse model for single image super-resolution

abstract: Self-similarity learning has been recognized as a promising method for single image super-resolution (SR) to produce high-resolution (HR) image in recent years. The performance of learning based SR reconstruction, however, highly depends on learned representation coeffcients. Due to the degradation of input image, conventional sparse coding is prone to produce unfaithful representation coeffcients. To this end, we propose a novel kernel based low-rank sparse model with self-similarity learning for single image SR which incorporates nonlocalsimilarity prior to enforce similar patches having similar representation weights. We perform a gradual magnification scheme, using self-examples extracted from the degraded input image and up-scaled versions. To exploit nonlocal-similarity, we concatenate the vectorized input patch and its nonlocal neighbors at different locations into a data matrix which consists of similar components. Then we map the nonlocal data matrix into a high-dimensional feature space by kernel method to capture their nonlinear structures. Under the assumption that the sparse coeffcients for the nonlocal data in the kernel space should be low-rank, we impose low-rank constraint on sparse coding to share similarities among representation coeffcients and remove outliers in order that stable weights for SR reconstruction can be obtained. Experimental results demonstrate the advantage of our proposed method in both visual quality and reconstruction error.


title: Molecular dust precursors in envelopes of oxygen-rich AGB stars and red supergiants

abstract: Condensation of circumstellar dust begins with formation of molecular clusters close to the stellar photosphere. These clusters are predicted to act as condensation cores at lower temperatures and allow efficient dust formation farther away from the star. Recent observations of metal oxides, such as AlO, AlOH, TiO, and TiO$_2$, whose emission can be traced at high angular resolutions with ALMA, have allowed first observational studies of the condensation process in oxygen-rich stars. We are now in the era when depletion of gas-phase species into dust can be observed directly. I review the most recent observations that allow us to identify gas species involved in the formation of inorganic dust of AGB stars and red supergiants. I also discuss challenges we face in interpreting the observations, especially those related to non-equilibrium gas excitation and the high complexity of stellar atmospheres in the dust-formation zone.


title: Industrial and Tramp Ship Routing Problems: Closing the Gap for Real-Scale Instances

abstract: Recent studies in maritime logistics have introduced a general ship routing problem and a benchmark suite based on real shipping segments, considering pickups and deliveries, cargo selection, ship-dependent starting locations, travel times and costs, time windows, and incompatibility constraints, among other features. Together, these characteristics pose considerable challenges for exact and heuristic methods, and some cases with as few as 18 cargoes remain unsolved. To face this challenge, we propose an exact branch-and-price (B&P) algorithm and a hybrid metaheuristic. Our exact method generates elementary routes, but exploits decremental state-space relaxation to speed up column generation, heuristic strong branching, as well as advanced preprocessing and route enumeration techniques. Our metaheuristic is a sophisticated extension of the unified hybrid genetic search. It exploits a set-partitioning phase and uses problem-tailored variation operators to efficiently handle all the problem characteristics. As shown in our experimental analyses, the B&P optimally solves 239/240 existing instances within one hour. Scalability experiments on even larger problems demonstrate that it can optimally solve problems with around 60 ships and 200 cargoes (i.e., 400 pickup and delivery services) and find optimality gaps below 1.04% on the largest cases with up to 260 cargoes. The hybrid metaheuristic outperforms all previous heuristics and produces near-optimal solutions within minutes. These results are noteworthy, since these instances are comparable in size with the largest problems routinely solved by shipping companies.


title: Fast QR decomposition of HODLR matrices

abstract: The efficient and accurate QR decomposition for matrices with hierarchical low-rank structures, such as HODLR and hierarchical matrices, has been challenging. Existing structure-exploiting algorithms are prone to numerical instability as they proceed indirectly, via Cholesky decompositions or a block Gram-Schmidt procedure. For a highly ill-conditioned matrix, such approaches either break down in finite-precision arithmetic or result in significant loss of orthogonality. Although these issues can sometimes be addressed by regularization and iterative refinement, it would be more desirable to have an algorithm that avoids these detours and is numerically robust to ill-conditioning. In this work, we propose such an algorithm for HODLR matrices. It achieves accuracy by utilizing Householder reflectors. It achieves efficiency by utilizing fast operations in the HODLR format in combination with compact WY representations and the recursive QR decomposition by Elmroth and Gustavson. Numerical experiments demonstrate that our newly proposed algorithm is robust to ill-conditioning and capable of achieving numerical orthogonality down to the level of roundoff error.


title: Inflation as a spontaneous symmetry breaking of Weyl symmetry

abstract: In this paper we study a novel realization of inflation, based on Weyl invariant gravity with torsion. We show that requiring the classical action for the scalar field to be Weyl invariant introduces a dilaton which induces a non trivial modification of the field space geometry of the scalar sector, which allows for inflationary phase that begins at the conformal point of the inflaton $\psi$, i.e. $\langle\psi\rangle=0$. Since the model is Weyl invariant, the inflaton condensation models a process of spontaneous Weyl symmetry breaking. For a wide range of parameters the spectral observables of the model are in good agreement with the CMB measurements, such that the scalar spectral index and the tensor-to-scalar ratio approximately agree with those of Starobinsky's inflation, i.e. $n_s\simeq 0.96-0.97$ and $r \approx 3\times 10^{-3}$. The simplest version of our model contains two scalar degrees of freedom, one of them being an exactly flat direction. If that degree is excited early on in inflation and if inflation lasts for about 60 e-foldings, we find that the Unverse undergoes a short period of kination that predates inflation. Such a period strongly suppresses the amplitude of large scale CMB temperature fluctuations providing thus an elegant explanation for the lack of power in the lowest CMB multipoles.


title: Large-deviation theory for a Brownian particle on a ring: a WKB approach

abstract: We study the large deviation function of the displacement of a Brownian particle confined on a ring. In the zero noise limit this large deviation function has a cusp at zero velocity given by the Freidlin-Wentzell theory. We develop a WKB approach to analyse how this cusp is rounded in the weak noise limit.


title: Property testing and expansion in cubical complexes

abstract: We consider expansion and property testing in the language of incidence geometry, covering both simplicial and cubical complexes in any dimension. We develop a general method for passing from an explicit description of the cohomology group, which need not be trivial, to a testability proof with linear ratio between errors. The method is demonstrated by testing functions on 2-cells in cubical complexes to be induced from the edges.