title: Nonlinear Cauchy-Riemann Equations and Liouville Equation For Conformal Metrics

abstract: We introduce the Nonlinear Cauchy-Riemann equations as B\"{a}cklund transformations for several nonlinear and linear partial differential equations. From these equations we treat in details the Laplace and the Liouville equations by deriving general solution for the nonlinear Liouville equation. By M\"{o}bius transformation we relate solutions for the Poincare model of hyperbolic geometry, the Klein model in half-plane and the pseudo-sphere. Conformal form of the constant curvature metrics in these geometries, stereographic projections and special solutions are discussed. Then we introduce the hyperbolic analog of the Riemann sphere, which we call the Riemann pseudosphere. We identify point at infinity on this pseudosphere and show that it can be used in complex analysis as an alternative to usual Riemann sphere to extend the complex plane. Interpretation of symmetric and antipodal points on both, the Riemann sphere and the Riemann pseudo-sphere, are given. By M\"{o}bius transformation and homogenous coordinates, the most general solution of Liouville equation as discussed by Crowdy is derived.


title: How driving rates determine the statistics of driven non-equilibrium systems with stationary distributions

abstract: Sample space reducing (SSR) processes offer a simple analytical way to understand of the origin and ubiquity of power-laws in many path-dependent complex systems. SRR processes show a wide range of applications that range from fragmentation processes, language formation to cascading pro- cesses. Here we argue that they also offer a natural framework to understand stationary distributions of generic driven non-equilibrium systems that are composed of a driving and a relaxing process. We show that the statistics of driven non-equilibrium systems can be derived from the understanding of the nature of the underlying driving process. For constant driving rates exact power-laws emerge with exponents that are related to the driving rate. If driving rates become state-dependent, or if they vary across the life-span of the process, the functional form of the state-dependence determines the statistics. Constant driving rates lead to exact power-laws, a linear state-dependence function yields exponential or Gamma distributions, a quadratic function gives the normal distribution. Logarithmic and power-law state dependence leads to log-normal and stretched exponential distribution functions, respectively. Also Weibull, Gompertz and Tsallis-Pareto distributions arise naturally from simple state-dependent driving rates. We discuss a simple physical example of consecutive elastic collisions that exactly represents a SSR process.


title: Joint Optimization of User Association, Data Delivery Rate and Precoding for Cache-Enabled F-RANs

abstract: This paper considers the downlink of a cache-enabled fog radio access network (F-RAN) with limited fronthaul capacity, where user association (UA), data delivery rate (DDR) and signal precoding are jointly optimized. We formulate a mixed-integer nonlinear programming problem in which the weighted difference of network throughput and total power consumption is maximized, subject to the predefined DDR requirements and the maximum transmit power at each eRRH. To address this challenging problem, we first apply the l0-norm approximation and l1-norm minimization techniques to deal with the UA. After this key step, we arrive at an approximated problem that only involves the joint optimization of DDR and precoding. By using the alternating descent method, we further decompose this problem into a convex subproblem for DDR allocation and a nonconvex subproblem for precoding design. While the former is globally solved by the interior-point method, the latter is solved by a specifically tailored successive convex quadratic programming method. Finally, we propose an iterative algorithm for the original joint optimization that is guaranteed to converge. Importantly, each iteration of the developed algorithm only involves solving simple convex problems. Numerical examples demonstrate that the proposed design significantly improves both throughput and power performances, especially in practical F-RANs with limited fronthaul capacity. Compared to the sole precoder design for a given cache placement, our joint design is shown to improve the throughput by 50% while saving at least half of the total power consumption in the considered examples.


title: The Single Particle Excitation Spectrum of Degenerate Fermi gases in a ring cavity

abstract: By considering a spin-$\frac{1}{2}$ degenerate Fermi gases in a ring cavity where strong interaction between atoms and light gives rise to a superradiance, we find the cavity dissipation could cause a severe broadening in some special cases, breaking down the quasi-particle picture which was constantly assumed in mean field theory studies. This broadening happens when the band gap resonant with polariton excitation energy. Interestingly enough, this broadening is highly spin selective depending on how the fermions are filled and the spectrum becomes asymmetric due to dissipation. Further, a non-monotonous dependence of the maximal broadening of the spectrum against cavity decay rate $\kappa$ is found and the largest broadening emerges at $\kappa$ comparable to recoil energy.


title: Imperfections, impacts, and the singularity of Euler's disk

abstract: The motion of a rigid, spinning disk on a flat surface ends with a dissipation-induced finite-time singularity. The problem of finding the dominant energy absorption mechanism during the last phase of the motion generated a lively debate during the past two decades. Various candidates including air drag and various types of friction have been considered, nevertheless impacts have not been examined until now. We investigate the effect of impacts caused by geometric imperfections of the disk and of the underlying flat surface, through analysing the dynamics of polygonal disks with unilateral point contacts. Similarly to earlier works, we determine the rate of energy absorption under the assumption of a regular pattern of motion analogous to precession-free motion of a rolling disk. In addition, we demonstrate that the asymptotic stability of this motion depends on parameters of the impact model. In the case of instability, the emerging irregular motion is investigated numerically. We conclude that there exists a range of model parameters (small radii of gyration or small restitution coefficients) in which absorption by impacts dominates all preiously investigated mechanisms during the last phase of motion. Nevertheless the parameter values associated with a homogenous disk on a hard surface are typically not in this range, hence the effect of impacts is in that case not dominant.


title: Sums of Palindromes: an Approach via Automata

abstract: Recently, Cilleruelo, Luca, & Baxter proved, for all bases b >= 5, that every natural number is the sum of at most 3 natural numbers whose base-b representation is a palindrome. However, the cases b = 2, 3, 4 were left unresolved. We prove, using a decision procedure based on automata, that every natural number is the sum of at most 4 natural numbers whose base-2 representation is a palindrome. Here the constant 4 is optimal. We obtain similar results for bases 3 and 4, thus completely resolving the problem. We consider some other variations on this problem, and prove similar results. We argue that heavily case-based proofs are a good signal that a decision procedure may help to automate the proof.


title: Optimization Methods for Supervised Machine Learning: From Linear Models to Deep Learning

abstract: The goal of this tutorial is to introduce key models, algorithms, and open questions related to the use of optimization methods for solving problems arising in machine learning. It is written with an INFORMS audience in mind, specifically those readers who are familiar with the basics of optimization algorithms, but less familiar with machine learning. We begin by deriving a formulation of a supervised learning problem and show how it leads to various optimization problems, depending on the context and underlying assumptions. We then discuss some of the distinctive features of these optimization problems, focusing on the examples of logistic regression and the training of deep neural networks. The latter half of the tutorial focuses on optimization algorithms, first for convex logistic regression, for which we discuss the use of first-order methods, the stochastic gradient method, variance reducing stochastic methods, and second-order methods. Finally, we discuss how these approaches can be employed to the training of deep neural networks, emphasizing the difficulties that arise from the complex, nonconvex structure of these models.


title: On Fairness, Diversity and Randomness in Algorithmic Decision Making

abstract: Consider a binary decision making process where a single machine learning classifier replaces a multitude of humans. We raise questions about the resulting loss of diversity in the decision making process. We study the potential benefits of using random classifier ensembles instead of a single classifier in the context of fairness-aware learning and demonstrate various attractive properties: (i) an ensemble of fair classifiers is guaranteed to be fair, for several different measures of fairness, (ii) an ensemble of unfair classifiers can still achieve fair outcomes, and (iii) an ensemble of classifiers can achieve better accuracy-fairness trade-offs than a single classifier. Finally, we introduce notions of distributional fairness to characterize further potential benefits of random classifier ensembles.


title: Storage, Communication, and Load Balancing Trade-off in Distributed Cache Networks

abstract: We consider load balancing in a network of caching servers delivering contents to end users. Randomized load balancing via the so-called power of two choices is a well-known approach in parallel and distributed systems. In this framework, we investigate the tension between storage resources, communication cost, and load balancing performance. To this end, we propose a randomized load balancing scheme which simultaneously considers cache size limitation and proximity in the server redirection process. In contrast to the classical power of two choices setup, since the memory limitation and the proximity constraint cause correlation in the server selection process, we may not benefit from the power of two choices. However, we prove that in certain regimes of problem parameters, our scheme results in the maximum load of order $\Theta(\log\log n)$ (here $n$ is the network size). This is an exponential improvement compared to the scheme which assigns each request to the nearest available replica. Interestingly, the extra communication cost incurred by our proposed scheme, compared to the nearest replica strategy, is small. Furthermore, our extensive simulations show that the trade-off trend does not depend on the network topology and library popularity profile details.


title: Study of charged-particle multiplicities with ALICE

abstract: The multiplicity measurements include the pseudorapidity density, $\rm{d}N_{\rm ch}/\rm{d}{\eta}$, and the probability distribution as a function of the number of charged particles, $\rm{P}(N_{\rm ch})$. ALICE has measured the multiplicities for three collision systems, for proton-proton, proton-lead and lead-lead collisions at Run 1 and 2 at the Large Hadron Collider. A selection of these results will be presented in these proceedings, concluding with an overview of new measurements planned.


title: Scalar field descriptions of two dark energy models

abstract: We give a scalar field description of two dark energy parameterizations, and we analyze in detail its cosmology both at the level of background evolution and at the level of linear perturbations. In particular, we compute the statefinder parameters and the growth index as functions of the red-shift for both dark energy parameterizations, and the comparison with the $\Lambda CDM$ model as well as with a few well-known geometrical dark energy models is shown. In addition, the combination parameter $A=f \sigma_8$ of both models is compared against current data.


title: The effective action of a BPS Alice string

abstract: Recently a BPS Alice string has been found in a $U(1)\times SU(2)$ gauge theory coupled with a charged complex adjoint scalar field arXiv:1703.08971. It is a half BPS state preserving a half of supercharges when embedded into a supersymmetric gauge theory. In this paper, we study zero modes of a BPS Alice string. After presenting $U(1)$ and translational zero modes, we construct the effective action of these modes. In contrast to previous analysis of the conventional Alice string for which only large distance behaviors are known, we can perform calculation exactly in the full space thanks to BPS properties.


title: Layer- and doping-tunable long range ferromagnetic order-ing in two-dimensional CrS2

abstract: Interlayer coupling is of vital importance for manipulating physical properties, e.g. electronic bandgap, in two-dimensional materials. However, tuning magnetic proper-ties in these materials is yet to be addressed. Here, we found a striped antiferromag-netic (sAFM) to ferromagnetic (FM) transition undergoing from monolayer to bilayer and thicker CrS2. This transition is attributed to charge sharing of interlayer S atoms and its resulting charge transfer from Cr d eg to t2g orbitals. The transferred charge reduces a portion of Cr4+ to Cr3+, which enhances the double-exchange mechanism favoring FM rather than sAFM ordering. Doping of electrons shares the same effect with stacking induced charge sharing. Therefore, charge doping could either manip-ulate the magnetic ordering between sAFM and FM or tune CrS2 between p- and n-doped magnetic semiconductors. We also proposed several prototype devices ena-bling to use external electric field for manipulating magnetic orderings of CrS2 layers. These results manifest the role of interlayer coupling in modifying magnetic proper-ties of layered materials


title: On the Geil-Matsumoto Bound and the Length of AG codes

abstract: The Geil-Matsumoto bound conditions the number of rational places of a function field in terms of the Weierstrass semigroup of any of the places. Lewittes' bound preceded the Geil-Matsumoto bound and it only considers the smallest generator of the numerical semigroup. It can be derived from the Geil-Matsumoto bound and so it is weaker. However, for general semigroups the Geil-Matsumoto bound does not have a closed formula and it may be hard to compute, while Lewittes' bound is very simple. We give a closed formula for the Geil-Matsumoto bound for the case when the Weierstrass semigroup has two generators. We first find a solution to the membership problem for semigroups generated by two integers and then apply it to find the above formula. We also study the semigroups for which Lewittes's bound and the Geil-Matsumoto bound coincide. We finally investigate on some simplifications for the computation of the Geil-Matsumoto bound.


title: Formation of Stable Strategic Networks with Desired Topologies

abstract: Many real-world networks such as social networks consist of strategic agents. The topology of these networks often plays a crucial role in determining the ease and speed with which certain information driven tasks can be accomplished. Consequently, growing a stable network having a certain desired topology is of interest. Motivated by this, we study the following important problem: given a certain desired topology, under what conditions would best response link alteration strategies adopted by strategic agents, uniquely lead to formation of a stable network having the given topology. This problem is the inverse of the classical network formation problem where we are concerned with determining stable topologies, given the conditions on the network parameters. We study this interesting inverse problem by proposing (1) a recursive model of network formation and (2) a utility model that captures key determinants of network formation. Building upon these models, we explore relevant topologies such as star graph, complete graph, bipartite Turan graph, and multiple stars with interconnected centers. We derive a set of sufficient conditions under which these topologies uniquely emerge, study their social welfare properties, and investigate the effects of deviating from the derived conditions.


title: Free-energy landscapes in magnetic systems from metadynamics

abstract: Knowledge of free energy barriers separating different states is critically important for assessment of long-term stability of information stored in magnetic devices. This information, however, is not directly accessible by standard simulations of microscopic models because of the ubiquitous time-scale problem, related to the fact that the transitions among different free energy minima have character of rare events. Here we show that by employing the metadynamics algorithm based on suitably chosen collective variables, namely helicity and circulation, it is possible to reliably recover the free energy landscape. We demonstrate the effectiveness of the new approach on the example of vortex nucleation process in magnetic nanodot with lowered spatial symmetry. With the help of reconstructed free energy surfaces (FES) we show the origin of the symmetry broken vortex nucleation, where one polarity of the nucleated vortex core is preferred, even though only in-plane magnetic field is present.


title: SMC Faster R-CNN: Toward a scene-specialized multi-object detector

abstract: Generally, the performance of a generic detector decreases significantly when it is tested on a specific scene due to the large variation between the source training dataset and the samples from the target scene. To solve this problem, we propose a new formalism of transfer learning based on the theory of a Sequential Monte Carlo (SMC) filter to automatically specialize a scene-specific Faster R-CNN detector. The suggested framework uses different strategies based on the SMC filter steps to approximate iteratively the target distribution as a set of samples in order to specialize the Faster R-CNN detector towards a target scene. Moreover, we put forward a likelihood function that combines spatio-temporal information extracted from the target video sequence and the confidence-score given by the output layer of the Faster R-CNN, to favor the selection of target samples associated with the right label. The effectiveness of the suggested framework is demonstrated through experiments on several public traffic datasets. Compared with the state-of-the-art specialization frameworks, the proposed framework presents encouraging results for both single and multi-traffic object detections.


title: Preparation of high crystalline nanoparticles of rare-earth based complex pervoskites and comparison of their structural and magnetic properties with bulk counterparts

abstract: A simple route to prepare Gd$_{0.7}$Sr$_{0.3}$MnO$_3$ nanoparticles by ultrasonication of their bulk powder materials is presented in this article. For comparison, Gd$_{0.7}$Sr$_{0.3}$MnO$_3$ nanoparticles are also prepared by ball milling. The prepared samples are characterized by X-ray diffraction (XRD),field emission scanning electron microscope (FESEM), energy dispersive X-ray (EDX), X-ray photoelectron spectroscope (XPS), and Superconducting Quantum Interference Device (SQUID) magnetometer. XRD Rietveld analysis is carried out extensively for the determination of crystallographic parameters and the amount of crystalline and amorphous phases. FESEM images demonstrate the formation of nanoparticles with average particle size in the range of 50-100 nm for both ultrasonication and 4 hours (h) of ball milling. The bulk materials and nanoparticles synthesized by both ultrasonication and 4 h ball milling exhibit a paramagnetic to spin-glass transition. However, nanoparticles synthesized by 8 h and 12 h ball milling do not reveal any phase transition, rather show an upturn of magnetization at low temperature. The degradation of the magnetic properties in ball milled nanoparticles may be associated with amorphization of the nanoparticles due to ball milling particularly for milling time exceeding 8 h. This investigation demonstrates the potential of ultrasonication as a simple route to prepare high crystalline rare-earth based manganite nanoparticles with improved control compared to the traditional ball milling technique.


title: Method of precision increase by averaging with application to numerical differentiation

abstract: If several independent algorithms for a computer-calculated quantity exist, then one can expect their results (which differ because of numerical errors) to follow approximately Gaussian distribution. The mean of this distribution, interpreted as the value of the quantity of interest, can be determined with better precision than what is the precision provided by a single algorithm. Often, with lack of enough independent algorithms, one can proceed differently: many practical algorithms introduce a bias using a parameter, e.g. a small but finite number to compute a limit or a large but finite number (cutoff) to approximate infinity. One may vary such parameter of a single algorithm and interpret the resulting numbers as generated by several algorithms. A numerical evidence for the validity of this approach is shown for differentiation.


title: 6thSense: A Context-aware Sensor-based Attack Detector for Smart Devices

abstract: Sensors (e.g., light, gyroscope, accelerotmeter) and sensing enabled applications on a smart device make the applications more user-friendly and efficient. However, the current permission-based sensor management systems of smart devices only focus on certain sensors and any App can get access to other sensors by just accessing the generic sensor API. In this way, attackers can exploit these sensors in numerous ways: they can extract or leak users' sensitive information, transfer malware, or record or steal sensitive information from other nearby devices. In this paper, we propose 6thSense, a context-aware intrusion detection system which enhances the security of smart devices by observing changes in sensor data for different tasks of users and creating a contextual model to distinguish benign and malicious behavior of sensors. 6thSense utilizes three different Machine Learning-based detection mechanisms (i.e., Markov Chain, Naive Bayes, and LMT) to detect malicious behavior associated with sensors. We implemented 6thSense on a sensor-rich Android smart device (i.e., smartphone) and collected data from typical daily activities of 50 real users. Furthermore, we evaluated the performance of 6thSense against three sensor-based threats: (1) a malicious App that can be triggered via a sensor (e.g., light), (2) a malicious App that can leak information via a sensor, and (3) a malicious App that can steal data using sensors. Our extensive evaluations show that the 6thSense framework is an effective and practical approach to defeat growing sensor-based threats with an accuracy above 96% without compromising the normal functionality of the device. Moreover, our framework costs minimal overhead.


title: Weakly intrusive low-rank approximation method for nonlinear parameter-dependent equations

abstract: This paper presents a weakly intrusive strategy for computing a low-rank approximation of the solution of a system of nonlinear parameter-dependent equations. The proposed strategy relies on a Newton-like iterative solver which only requires evaluations of the residual of the parameter-dependent equation and of a preconditioner (such as the differential of the residual) for instances of the parameters independently. The algorithm provides an approximation of the set of solutions associated with a possibly large number of instances of the parameters, with a computational complexity which can be orders of magnitude lower than when using the same Newton-like solver for all instances of the parameters. The reduction of complexity requires efficient strategies for obtaining low-rank approximations of the residual, of the preconditioner, and of the increment at each iteration of the algorithm. For the approximation of the residual and the preconditioner, weakly intrusive variants of the empirical interpolation method are introduced, which require evaluations of entries of the residual and the preconditioner. Then, an approximation of the increment is obtained by using a greedy algorithm for low-rank approximation, and a low-rank approximation of the iterate is finally obtained by using a truncated singular value decomposition. When the preconditioner is the differential of the residual, the proposed algorithm is interpreted as an inexact Newton solver for which a detailed convergence analysis is provided. Numerical examples illustrate the efficiency of the method.


title: A restriction isomorphism for zero-cycles with coefficients in Milnor K-theory

abstract: We prove a restriction isomorphism for Chow groups of zero-cycles with coefficients in Milnor K-theory for smooth projective schemes over excellent henselian discrete valuation rings. Furthermore, we study torsion subgroups of these groups over local and finite fields.


title: Press F1 for help: participatory design for dealing with on-line and real life security of older adults

abstract: In this paper we present the report on the design and development of a platform for the inter-generational exchange of favors. This platform was designed using participatory design approach during a 24-hour hackathon by a team consisting of younger programmers and older adults. The findings of this report show that inter-generational cooperation in which the older adults serve as representatives of the end user, not only improves the design and development of the application, but also provides an effective method for designing and applying solutions aimed at improving the security of older adults while using online and mobile tools.


title: Bayesian inference using intermediate distribution based on coarse multiscale model for time fractional diffusion equation

abstract: In the paper, we present a strategy for accelerating posterior inference for unknown inputs in time fractional diffusion models. In many inference problems, the posterior may be concentrated in a small portion of the entire prior support. It will be much more efficient if we build and simulate a surrogate only over the significant region of the posterior. To this end, we construct a coarse model using Generalized Multiscale Finite Element Method (GMsFEM), and solve a least-squares problem for the coarse model with a regularizing Levenberg-Marquart algorithm. An intermediate distribution is built based on the approximate sampling distribution. For Bayesian inference, we use GMsFEM and least-squares stochastic collocation method to obtain a reduced coarse model based on the intermediate distribution. To increase the sampling speed of Markov chain Monte Carlo, the DREAM$_\text{ZS}$ algorithm is used to explore the surrogate posterior density, which is based on the surrogate likelihood and the intermediate distribution. The proposed method with lower gPC order gives the approximate posterior as accurate as the the surrogate model directly based on the original prior. A few numerical examples for time fractional diffusion equations are carried out to demonstrate the performance of the proposed method with applications of the Bayesian inversion.


title: Note on correlation functions in conformal quantum mechanics

abstract: We suggest a method to compute the correlation functions in conformal quantum mechanics (CFT$_1$) for the fields that transform under a non-local representation of $\mathfrak{sl}(2)$ basing on the invariance properties. Explicit calculations of 2- and 3-point correlation functions are presented.


title: Dy doped BiFeO$_3$ : A Bulk Ceramic with Improved Multiferroic Properties Compared to Nano Counterparts

abstract: The synthesis as well as structural, multiferroic and optical characterization of Dy doped BiFeO$_3$ multiferroic ceramic are presented. Bulk polycrystalline Bi$_{0.9}$Dy$_{0.1}$FeO$_3$ sample is synthesized by solid state reaction, while their nano counterparts are prepared using ultrasonic probe sonication technique. Significant improvement of phase purity in the as synthesized samples is observed after the doping of Dy both in bulk Bi$_{0.9}$Dy$_{0.1}$FeO$_3$ sample and corresponding nanoparticles as evidenced from Rietveld refinement. Magnetization measurements using SQUID magnetometer exhibit enhanced magnetic properties for Dy doped bulk Bi$_{0.9}$Dy$_{0.1}$FeO$_3$ ceramic compared to their nanostructured counterparts as well as undoped BiFeO$_3$. Within the applied field range, saturation polarization is observed for Bi$_{0.9}$Dy$_{0.1}$FeO$_3$ bulk ceramic only. As a result, intrinsic ferroelectric behavior is obtained just for this sample. Optical bandgap measurements reveal lower bandgap for Dy doped bulk Bi$_{0.9}$Dy$_{0.1}$FeO$_3$ ceramic compared to that of corresponding nanoparticles and undoped BiFeO$_3$. The outcome of this investigation demonstrates the potential of Dy as a doping element in BiFeO$_3$ that provides a bulk ceramic material with improved multiferroic and optical properties compared to those of corresponding nanoparticles which involve rigorous synthesis procedure.


title: Levitated optomechanics with a fiber Fabry-Perot interferometer

abstract: In recent years quantum phenomena have been experimentally demonstrated on variety of optomechanical systems ranging from micro-oscillators to photonic crystals. Since single photon couplings are quite small, most experimental approaches rely on the realization of high finesse Fabry-Perot cavities in order to enhance the effective coupling. Here we show that by exploiting a, long path, low finesse fiber Fabry-Perot interferometer ground state cooling can be achieved. We model a 100 m long cavity with a finesse of 10 and analyze the impact of additional noise sources arising from the fiber. As a mechanical oscillator we consider a levitated microdisk but the same approach could be applied to other optomechanical systems.


title: Contracting a Planar Graph Efficiently

abstract: We present a data structure that can maintain a simple planar graph under edge contractions in linear total time. The data structure supports adjacency queries and provides access to neighbor lists in $O(1)$ time. Moreover, it can report all the arising self-loops and parallel edges. By applying the data structure, we can achieve optimal running times for decremental bridge detection, 2-edge connectivity, maximal 3-edge connected components, and the problem of finding a unique perfect matching for a static planar graph. Furthermore, we improve the running times of algorithms for several planar graph problems, including decremental 2-vertex and 3-edge connectivity, and we show that using our data structure in a black-box manner, one obtains conceptually simple optimal algorithms for computing MST and 5-coloring in planar graphs.


title: Achieving DFT accuracy with a machine-learning interatomic potential: thermomechanics and defects in bcc ferromagnetic iron

abstract: We show that the Gaussian Approximation Potential machine learning framework can describe complex magnetic potential energy surfaces, taking ferromagnetic iron as a paradigmatic challenging case. The training database includes total energies, forces, and stresses obtained from density-functional theory in the generalized-gradient approximation, and comprises approximately 150,000 local atomic environments, ranging from pristine and defected bulk configurations to surfaces and generalized stacking faults with different crystallographic orientations. We find the structural, vibrational and thermodynamic properties of the GAP model to be in excellent agreement with those obtained directly from first-principles electronic-structure calculations. There is good transferability to quantities, such as Peierls energy barriers, which are determined to a large extent by atomic configurations that were not part of the training set. We observe the benefit and the need of using highly converged electronic-structure calculations to sample a target potential energy surface. The end result is a systematically improvable potential that can achieve the same accuracy of density-functional theory calculations, but at a fraction of the computational cost.


title: Mapping the hot gas temperature in galaxy clusters using X-ray and Sunyaev-Zel'dovich imaging

abstract: We propose a method to map the temperature distribution of the hot gas in galaxy clusters that uses resolved images of the thermal Sunyaev-Zel'dovich (tSZ) effect in combination with X-ray data. Application to images from the New IRAM KIDs Array (NIKA) and XMM-Newton allows us to measure and determine the spatial distribution of the gas temperature in the merging cluster MACS J0717.5+3745, at $z=0.55$. Despite the complexity of the target object, we find a good morphological agreement between the temperature maps derived from X-ray spectroscopy only -- using XMM-Newton ($T_{\rm XMM}$) and Chandra ($T_{\rm CXO}$) -- and the new gas-mass-weighted tSZ+X-ray imaging method ($T_{\rm SZ}$). We correlate the temperatures from tSZ+X-ray imaging and those from X-ray spectroscopy alone and find that $T_{\rm SZ}$ is higher than $T_{\rm XMM}$ and lower than $T_{\rm CXO}$ by $\sim 10\%$ in both cases. Our results are limited by uncertainties in the geometry of the cluster gas, contamination from kinetic SZ ($\sim 10\%$), and the absolute calibration of the tSZ map ($7\%$). Investigation using a larger sample of clusters would help minimise these effects.


title: Improving Session Recommendation with Recurrent Neural Networks by Exploiting Dwell Time

abstract: Recently, Recurrent Neural Networks (RNNs) have been applied to the task of session-based recommendation. These approaches use RNNs to predict the next item in a user session based on the previ- ously visited items. While some approaches consider additional item properties, we argue that item dwell time can be used as an implicit measure of user interest to improve session-based item recommen- dations. We propose an extension to existing RNN approaches that captures user dwell time in addition to the visited items and show that recommendation performance can be improved. Additionally, we investigate the usefulness of a single validation split for model selection in the case of minor improvements and find that in our case the best model is not selected and a fold-like study with different validation sets is necessary to ensure the selection of the best model.


title: Relativistic Three-Dimensional Lippmann-Schwinger Cross Sections for Space Radiation Applications

abstract: Radiation transport codes require accurate nuclear cross sections to compute particle fluences inside shielding materials. The Tripathi semi-empirical reaction cross section, which includes over 60 parameters tuned to nucleon-nucleus (NA) and nucleus-nucleus (AA) data, has been used in many of the world's best-known transport codes. Although this parameterization fits well to reaction cross section data, the predictive capability of any parameterization is questionable when it is used beyond the range of the data to which it was tuned. Using uncertainty analysis, it is shown that a relativistic Three-Dimensional Lippmann-Schwinger (LS3D) equation model based on Multiple Scattering Theory (MST) that uses 5 parameterizations---3 fundamental parameterizations to nucleon-nucleon (NN) data and 2 nuclear charge density parameterizations---predicts NA and AA reaction cross sections as well as the Tripathi cross section parameterization for reactions in which the kinetic energy of the projectile in the laboratory frame ($T_{\rm Lab}$) is greater than 220 MeV/n. The relativistic LS3D model has the additional advantage of being able to predict highly accurate total and elastic cross sections. Consequently, it is recommended that the relativistic LS3D model be used for space radiation applications in which $T_{\rm Lab}>$ 220 MeV/n.


title: Conservation in two-particle self-consistent extensions of dynamical-mean-field-theory

abstract: Extensions of dynamical-mean-field-theory (DMFT) make use of quantum impurity models as non-perturbative and exactly solvable reference systems which are essential to treat the strong electronic correlations. Through the introduction of retarded interactions on the impurity, these approximations can be made two-particle self-consistent. This is of interest for the Hubbard model, because it allows to suppress the antiferromagnetic phase transition in two-dimensions in accordance with the Mermin-Wagner theorem, and to include the effects of bosonic fluctuations. For a physically sound description of the latter, the approximation should be conserving. In this paper we show that the mutual requirements of two-particle self-consistency and conservation lead to fundamental problems. For an approximation that is two-particle self-consistent in the charge- and longitudinal spin channel, the double occupancy of the lattice and the impurity are no longer consistent when computed from single-particle properties. For the case of self-consistency in the charge- and longitudinal as well as transversal spin channels, these requirements are even mutually exclusive so that no conserving approximation can exist. We illustrate these findings for a two-particle self-consistent and conserving DMFT approximation.


title: Probabilistic Active Learning of Functions in Structural Causal Models

abstract: We consider the problem of learning the functions computing children from parents in a Structural Causal Model once the underlying causal graph has been identified. This is in some sense the second step after causal discovery. Taking a probabilistic approach to estimating these functions, we derive a natural myopic active learning scheme that identifies the intervention which is optimally informative about all of the unknown functions jointly, given previously observed data. We test the derived algorithms on simple examples, to demonstrate that they produce a structured exploration policy that significantly improves on unstructured base-lines.


title: Connecting macroscopic dynamics with microscopic properties in active microtubule network contraction

abstract: The cellular cytoskeleton is an active material, driven out of equilibrium by molecular motor proteins. It is not understood how the collective behaviors of cytoskeletal networks emerge from the properties of the network's constituent motor proteins and filaments. Here we present experimental results on networks of stabilized microtubules in Xenopus oocyte extracts, which undergo spontaneous bulk contraction driven by the motor protein dynein, and investigate the effects of varying the initial microtubule density and length distribution. We find that networks contract to a similar final density, irrespective of the length of microtubules or their initial density, but that the contraction timescale varies with the average microtubule length. To gain insight into why this microscopic property influences the macroscopic network contraction time, we developed simulations where microtubules and motors are explicitly represented. The simulations qualitatively recapitulate the variation of contraction timescale with microtubule length, and allowed stress contributions from different sources to be estimated and decoupled.


title: Cold collisions in a molecular synchrotron

abstract: We study collisions between neutral, deuterated ammonia molecules (ND$_3$) stored in a 50 cm diameter synchrotron and argon atoms in co-propagating supersonic beams. The advantages of using a synchrotron in collision studies are twofold: (i) By storing ammonia molecules many round-trips, the sensitivity to collisions is greatly enhanced; (ii) The collision partners move in the same direction as the stored molecules, resulting in low collision energies. We tune the collision energy in three different ways: by varying the velocity of the stored ammonia packets, by varying the temperature of the pulsed valve that releases the argon atoms, and by varying the timing between the supersonic argon beam and the stored ammonia packets. These give consistent results. We determine the relative, total, integrated cross-section for $\mathrm{ND}_3+\mathrm{Ar}$ collisions in the energy range of 40-140 cm$^{-1}$, with a resolution of 5-10 cm$^{-1}$ and an uncertainty of 7-15%. Our measurements are in good agreement with theoretical scattering calculations.


title: Is it ethical to avoid error analysis?

abstract: Machine learning algorithms tend to create more accurate models with the availability of large datasets. In some cases, highly accurate models can hide the presence of bias in the data. There are several studies published that tackle the development of discriminatory-aware machine learning algorithms. We center on the further evaluation of machine learning models by doing error analysis, to understand under what conditions the model is not working as expected. We focus on the ethical implications of avoiding error analysis, from a falsification of results and discrimination perspective. Finally, we show different ways to approach error analysis in non-interpretable machine learning algorithms such as deep learning.


title: The role of quantum coherence in dimer and trimer excitation energy transfer

abstract: Recent progress in resource theory of quantum coherence has resulted in measures to quantify coherence in quantum systems. Especially, the l1-norm and relative entropy of coherence have been shown to be proper quantifiers of coherence and have been used to investigate coherence properties in different operational tasks. Since long-lasting quantum coherence has been experimentally confirmed in a number of photosynthetic complexes, it has been debated if and how coherence is connected to the known efficiency of population transfer in such systems. In this study, we investigate quantitatively the relationship between coherence, as quantified by l1-norm and relative entropy of coherence, and efficiency, as quantified by fidelity, for population transfer between end-sites in a network of two-level quantum systems. In particular, we use the coherence averaged over the duration of the population transfer in order to carry out a quantitative comparision between coherence and fidelity. Our results show that although coherence is a necessary requirement for population transfer, there is no unique relation between coherence and the efficiency of the transfer process.


title: Towards Understanding Generalization of Deep Learning: Perspective of Loss Landscapes

abstract: It is widely observed that deep learning models with learned parameters generalize well, even with much more model parameters than the number of training samples. We systematically investigate the underlying reasons why deep neural networks often generalize well, and reveal the difference between the minima (with the same training error) that generalize well and those they don't. We show that it is the characteristics the landscape of the loss function that explains the good generalization capability. For the landscape of loss function for deep networks, the volume of basin of attraction of good minima dominates over that of poor minima, which guarantees optimization methods with random initialization to converge to good minima. We theoretically justify our findings through analyzing 2-layer neural networks; and show that the low-complexity solutions have a small norm of Hessian matrix with respect to model parameters. For deeper networks, extensive numerical evidence helps to support our arguments.


title: Bridging the Gap between Probabilistic and Deterministic Models: A Simulation Study on a Variational Bayes Predictive Coding Recurrent Neural Network Model

abstract: The current paper proposes a novel variational Bayes predictive coding RNN model, which can learn to generate fluctuated temporal patterns from exemplars. The model learns to maximize the lower bound of the weighted sum of the regularization and reconstruction error terms. We examined how this weighting can affect development of different types of information processing while learning fluctuated temporal patterns. Simulation results show that strong weighting of the reconstruction term causes the development of deterministic chaos for imitating the randomness observed in target sequences, while strong weighting of the regularization term causes the development of stochastic dynamics imitating probabilistic processes observed in targets. Moreover, results indicate that the most generalized learning emerges between these two extremes. The paper concludes with implications in terms of the underlying neuronal mechanisms for autism spectrum disorder and for free action.


title: A selectional auto-encoder approach for document image binarization

abstract: Binarization plays a key role in the automatic information retrieval from document images. This process is usually performed in the first stages of documents analysis systems, and serves as a basis for subsequent steps. Hence it has to be robust in order to allow the full analysis workflow to be successful. Several methods for document image binarization have been proposed so far, most of which are based on hand-crafted image processing strategies. Recently, Convolutional Neural Networks have shown an amazing performance in many disparate duties related to computer vision. In this paper we discuss the use of convolutional auto-encoders devoted to learning an end-to-end map from an input image to its selectional output, in which activations indicate the likelihood of pixels to be either foreground or background. Once trained, documents can therefore be binarized by parsing them through the model and applying a threshold. This approach has proven to outperform existing binarization strategies in a number of document domains.


title: Comprehensive study of the dynamics of a classical Kitaev spin liquid

abstract: We study the spin-$S$ Kitaev model in the classical ($S \to \infty$) limit using Monte Carlo simulations combined with semi-classical spin dynamics. We discuss differences and similarities in the dynamical structure factors of the spin-$1/2$ and the classical Kitaev liquids. Interestingly, the low-temperature and low-energy spectrum of the classical model exhibits a finite energy peak, which is the precursor of the one produced by the Majorana modes of the $S=1/2$ model. The classical peak is spectrally narrowed compared to the quantum result and can be explained by magnon excitations within fluctuating one-dimensional manifolds (loops). Hence the difference from the classical limit to the quantum limit can be understood by the fractionalization of magnons propagating in one-dimensional manifolds. Moreover, we show that the momentum space distribution of the low-energy spectral weight of the $S=1/2$ model follows the momentum space distribution of zero modes of the classical model.


title: Inverse Design of Conjugated Polymers from Computed Electronic Structure Properties: Model Chemistries of Polythiophenes

abstract: Using a diverse set of 100 oligothiophenes, extrapolated oligomer and polymer HOMO, LUMO and HOMO-LUMO gaps are shown to be accurately estimated from the computed values of the trimer HOMO, trimer LUMO and trimer HOMO-LUMO gaps. These simple approximations can be improved further through easy to calculate properties of monomers or small oligomers. Polymer reorganization energies, related to the hole transport, are also shown to be predicted accurately from small oligomers. Correlations between the HOMO slope with the HOMO polymer energy and HOMO-LUMO gap slope with the polymer HOMO-LUMO gap energy suggest that the degree of delocalization reflected in the slope is higher in monomers and polymers with electron-rich, less negative HOMO eigenvalues. These statistical models can serve as efficient screens for the computational prediction of polythiophene properties.


title: A null test of General Relativity: New limits on Local Position Invariance and the variation of fundamental constants

abstract: We compare the long-term fractional frequency variation of four hydrogen masers that are part of an ensemble of clocks comprising the National Institute of Standards and Technology,(NIST), Boulder, timescale with the fractional frequencies of primary frequency standards operated by leading metrology laboratories in the United States, France, Germany, Italy and the United Kingdom for a period extending more than 14 years. The measure of the assumed variation of non-gravitational interaction,(LPI parameter, $\beta$)---within the atoms of H and Cs---over time as the earth orbits the sun, has been constrained to $\beta=(2.2 \pm 2.5)\times 10^{-7}$, a factor of two improvement over previous estimates. Using our results together with the previous best estimates of $\beta$ based on Rb vs. Cs, and Rb vs. H comparisons, we impose the most stringent limits to date on the dimensionless coupling constants that relate the variation of fundamental constants such as the fine-structure constant and the scaled quark mass with strong(QCD) interaction to the variation in the local gravitational potential. For any metric theory of gravity $\beta=0$.


title: Higher Order Tur\'an Inequalities for the Partition Function

abstract: The Tur\'{a}n inequalities and the higher order Tur\'{a}n inequalities arise in the study of Maclaurin coefficients of an entire function in the Laguerre-P\'{o}lya class. A real sequence $\{a_{n}\}$ is said to satisfy the Tur\'{a}n inequalities if for $n\geq 1$, $a_n^2-a_{n-1}a_{n+1}\geq 0$. It is said to satisfy the higher order Tur\'{a}n inequalities if for $n\geq 1$, $4(a_{n}^2-a_{n-1}a_{n+1})(a_{n+1}^2-a_{n}a_{n+2})-(a_{n}a_{n+1}-a_{n-1}a_{n+2})^2\geq 0$. A sequence satisfying the Tur\'an inequalities is also called log-concave. For the partition function $p(n)$, DeSalvo and Pak showed that for $n>25$, the sequence $\{ p(n)\}_{n> 25}$ is log-concave, that is, $p(n)^2-p(n-1)p(n+1)>0$ for $n> 25$. It was conjectured by Chen that $p(n)$ satisfies the higher order Tur\'{a}n inequalities for $n\geq 95$. In this paper, we prove this conjecture by using the Hardy-Ramanujan-Rademacher formula to derive an upper bound and a lower bound for $p(n+1)p(n-1)/p(n)^2$. Consequently, for $n\geq 95$, the Jensen polynomials $g_{3,n-1}(x)=p(n-1)+3p(n)x+3p(n+1)x^2+p(n+2)x^3$ have only real zeros. We conjecture that for any positive integer $m\geq 4$ there exists an integer $N(m)$ such that for $n\geq N(m) $, the polynomials $\sum_{k=0}^m {m\choose k}p(n+k)x^k$ have only real zeros. This conjecture was independently posed by Ono.


title: Uniquenss of Some Weak Solutions for 2D Viscous Primitive Equations

abstract: First, a new sufficient condition for uniqueness of weak solutions is proved for the system of 2D viscous Primitive Equations. Second, global existence and uniqueness are established for several classes of weak solutions with partial initial regularity, including but not limited to those weak solutions with initial horizontal regularity, rather than vertical regularity. Our results and analyses for the problem with phys- ical boundary conditions can be extended to those with other typical boundary conditions. Most of the results were not available before, even for the periodic case.


title: Photoinduced Hund excitons in the breakdown of a two-orbital Mott insulator

abstract: We study the photoinduced breakdown of a two-orbital Mott insulator and resulting metallic state. Using time-dependent density matrix renormalization group, we scrutinize the real-time dynamics of the half-filled two-orbital Hubbard model interacting with a resonant radiation field pulse. The breakdown, caused by production of doublon-holon pairs, is enhanced by Hund's exchange, which dynamically activates large orbital fluctuations. The melting of the Mott insulator is accompanied by a high to low spin transition with a concomitant reduction of antiferromagnetic spin fluctuations. Most notably, the overall time response is driven by the photogeneration of excitons with orbital character that are stabilized by Hund's coupling. These unconventional "Hund excitons" correspond to bound spin-singlet orbital-triplet doublon-holon pairs. We study exciton properties such as bandwidth, binding potential, and size within a semiclassical approach. The photometallic state results from a coexistence of Hund excitons and doublon-holon plasma.


title: Existence result for impulsive coupled systems on the half-line

abstract: This work considers a second order impulsive coupled system of differential equations with generalized jump conditions in half-line, which can depend on the impulses of the unknown functions and their first derivatives. The arguments apply the fixed point theory, Green's functions technique, L^11-Carath\'eodory functions and sequences and Schauder's fixed point theorem. The method is based on Carath\'eodory concept of functions and sequences, together with the equiconvergence on infinity and on each impulsive moment, and it allows to consider coupled fully nonlinearities and very general impulsive functions.


title: Robustness of Majorana edge modes and topological order -- exact results for the symmetric interacting Kitaev chain with disorder

abstract: We investigate the robustness of Majorana edge modes under disorder and interactions. We exploit a recently found mapping of the interacting Kitaev chain in the symmetric region ($\mu = 0$, $t = \Delta$) to free fermions. Extending the exact solution to the disordered case allows us to calculate analytically the topological phase boundary for all interaction and disorder strengths, which has been thought to be only accessible numerically. We discover a regime in which moderate disorder in the interaction matrix elements enhances topological order well into the strongly interacting regime $U > t$. We also derive the explicit form of the many-body Majorana edge wave function revealing how it is dressed by many-particle fluctuations from interactions. The qualitative features of our analytical results are valid beyond the fine-tuned integrable point as expected from the robustness of topological order and as corroborated here by an exact diagonalization study of small systems.


title: On collapsing Calabi-Yau fibrations

abstract: We develop some techniques to study the adiabatic limiting behaviour of Calabi-Yau metrics on the total space of a fibration, and obtain strong control near the singular fibres by imposing restrictions on the singularity types. We prove a uniform lower bound on the metric up to the singular fibre, under fairly general hypotheses. Assuming a result in pluripotential theory, we prove a uniform fibre diameter bound for a Lefschetz K3 fibred Calabi-Yau 3-fold, which reduces the study of the collapsing metric to a locally non-collapsed situation, and we identify the Gromov-Hausdorff limit of the rescaled neighbourhood of the singular fibre.


title: Quantitative Normal Approximation of Linear Statistics of $\beta$-ensembles

abstract: We present a new approach, inspired by Stein's method, to prove a central limit theorem (CLT) for linear statistics of $\beta$-ensembles in the one-cut regime. Compared with the previous proofs, our result requires less regularity on the potential and provides a rate of convergence in the quadratic Kantorovich or Wasserstein-2 distance. The rate depends both on the regularity of the potential and the test functions, and we prove that it is optimal in the case of the Gaussian Unitary Ensemble (GUE) for certain polynomial test functions. The method relies on a general normal approximation result of independent interest which is valid for a large class of Gibbs-type distributions. In the context of $\beta$-ensembles, this leads to a multi-dimensional CLT for a sequence of linear statistics which are approximate eigenfunctions of the infinitesimal generator of Dyson Brownian motion once the various error terms are controlled using the rigidity results of Bourgade, Erd\H{o}s and Yau.


title: Lattice polarization effects on the screened Coulomb interaction $W$ of the GW approximation

abstract: In polar insulators where longitudinal and transverse optical phonon modes differ substantially, the electron-phonon coupling affects the energy-band structure primarily through the long-range Fr\"ohlich contribution to the Fan term. This diagram has the same structure as the $GW$ self-energy where $W$ originates from the electron part of the screened coulomb interaction. The two can be conveniently combined by combining electron and lattice contributions to the polarizability. Both contributions are nonanalytic at the origin, and diverge as $1/q^2$ so that the predominant contribution comes from a small region around $q{=}0$. Here we adopt a simple estimate for the Fr\"ohlich contribution by assuming that the entire phonon part can be attributed to a small volume of $q$ near $q{=}0$. We estimate the magnitude for $\mathbf{q}{\rightarrow}0$ from a generalized Lyddane-Sachs-Teller relation, and the radius from the inverse of the polaron length scale. The gap correction is shown to agree with Fr\"ohlich's simple estimate $-\alpha_P\omega_L/2$ of the polaron effect.


title: A review and prospects for Nb3Sn superconductor development

abstract: Nb3Sn superconductors have significant applications in constructing high-field (> 10 T) magnets. This article briefly reviews development of Nb3Sn superconductor and proposes prospects for further improvement. It is shown that significant improvement of critical current density (Jc) is needed for future accelerator magnets. After a brief review of the development of Nb3Sn superconductors, the factors controlling Jc are summarized and correlated with their microstructure and chemistry. The non-matrix Jc of Nb3Sn conductors is mainly determined by three factors: the fraction of current-carrying Nb3Sn phase in the non-matrix area, the upper critical field Bc2, and the flux-line pinning capacity. Then prospects to improve the three factors are discussed respectively. An analytic model was developed to show how the ratios of precursors determine the phase fractions after heat treatment, based on which it is predicted that the limit of current-carrying Nb3Sn fraction in subelements is ~65%. Then, since Bc2 is largely determined by the Nb3Sn stoichiometry, a thermodynamic/kinetic theory was presented to show what essentially determines the Sn content of Nb3Sn conductors. This theory explains the influences of Sn sources and Ti addition on stoichiometry and growth rate of Nb3Sn layers. Next, to improve flux pinning, previous efforts in this community to introduce additional pinning centers (APC) to Nb3Sn wires are reviewed, and an internal oxidation technique is described. Finally, prospects for further improvement of non-matrix Jc of Nb3Sn conductors are discussed, and it is seen that the only opportunity for further significantly improving Jc lies in improving the flux pinning.


title: Graphene-Complex-oxide Nanoscale Device Concepts

abstract: The integration of graphene with complex-oxide heterostructures such as LaAlO$_3$/SrTiO$_3$ offers the opportunity to combine the multifunctional properties of an oxide interface with the electronic properties of graphene. The ability to control interface conduction through graphene and understanding how it affects the intrinsic properties of an oxide interface are critical to the technological development of novel multifunctional devices. Here we demonstrate several device archetypes in which electron transport at an oxide interface is modulated using a patterned graphene top gate. Nanoscale devices are fabricated at the oxide interface by conductive atomic force microscope (c-AFM) lithography, and transport measurements are performed as a function of the graphene gate voltage. Experiments are performed with devices written adjacent to or directly underneath the graphene gate. Unique capabilities of this approach include the ability to create highly flexible device configurations, the ability to modulate carrier density at the oxide interface, and the ability to control electron transport up to the single-electron-tunneling regime, while maintaining intrinsic transport properties of the oxide interface. Our results facilitate the design of a variety of nanoscale devices that combine unique transport properties of these two intimately coupled two-dimensional electron systems.


title: Line Hermitian Grassmann Codes and their Parameters

abstract: In this paper we introduce and study line Hermitian Grassmann codes as those subcodes of the Grassmann codes associated to the $2$-Grassmannian of a Hermitian polar space defined over a finite field of square order. In particular, we determine their parameters and characterize the words of minimum weight.


title: A faster dual algorithm for the Euclidean minimum covering ball problem

abstract: Dearing and Zeck presented a dual algorithm for the problem of the minimum covering ball in $\mathbb{R}^n$. Each iteration of their algorithm has a computational complexity of at least $\mathcal O(n^3)$. In this paper we propose a modification to their algorithm that, together with an implementation that uses updates to the QR factorization of a suitable matrix, achieves a $\mathcal O(n^2)$ iteration.


title: From the GKLS equation to the theory of solar and fuel cells

abstract: The mathematically sound theory of quantum open systems, formulated in the 70-ties of XX century and highlighted by the discovery of Gorini-Kossakowski-Lindblad-Sudarshan (GKLS) equation, found a wide range of applications in various branches of physics and chemistry, notably in the field of quantum information and quantum thermodynamics. However, it took 40 years before this formalism has been applied to explain correctly the operation principles of long existing energy transducers like photovoltaic, thermoelectric and fuel cells. This long path is briefly reviewed from the personal author's perspective. Finally, the new, fully quantum model of chemical engine based on GKLS equation and applicable to fuel cells or replicators is outlined. The model illustrates the difficulty with an entirely quantum operational definition of work, comparable to the Problem of Quantum Measurement.


title: The cohomology of free loop spaces of homogeneous spaces

abstract: The free loops space $\Lambda X$ of a space $X$ has become an important object of study particularly in the case when $X$ is a manifold.The study of free loop spaces is motivated in particular by two main examples. The first is their relation to geometrically distinct periodic geodesics on a manifold, originally studied by Gromoll and Meyer in $1969$. More recently the study of string topology and in particular the Chas-Sullivan loop product has been an active area of research. A complete flag manifold is the quotient of a Lie group by its maximal torus and is one of the nicer examples of a homogeneous space. Both the cohomology and Chas-Sullivan product structure are understood for spaces $S^n$, $\mathbb{C}P^n$ and most simple Lie groups. Hence studying the topology of the free loops space on homogeneous space is a natural next step. In the thesis we compute the differentials in the integral Leray-Serre spectral sequence associated to the free loops space fibrations in the cases of $SU(n+1)/T^n$ and $Sp(n)/T^n$. Study in detail the structure of the third page of the spectral sequence in the case of $SU(n)$ and give the module structure of $H^*(\Lambda(SU(3)/T^2);\mathbb{Z})$ and $H^*(\Lambda(Sp(2)/T^2);\mathbb{Z})$.


title: Hilbert isometries and maximal deviation preserving maps on JB-algebras

abstract: In this paper we characterize the surjective linear variation norm isometries on JB-algebras. Variation norm isometries are precisely the maps that preserve the maximal deviation, the quantum analogue of the standard deviation, which plays an important role in quantum statistics. Consequently, we characterize the Hilbert's metric isometries on cones in JB-algebras.


title: How biased is your model? Concentration Inequalities, Information and Model Bias

abstract: We derive tight and computable bounds on the bias of statistical estimators, or more generally of quantities of interest, when evaluated on a baseline model P rather than on the typically unknown true model Q. Our proposed method combines the scalable information inequality derived by P. Dupuis, K.Chowdhary, the authors and their collaborators together with classical concentration inequalities (such as Bennett's and Hoeffding-Azuma inequalities). Our bounds are expressed in terms of the Kullback-Leibler divergence R(Q||P) of model Q with respect to P and the moment generating function for the statistical estimator under P. Furthermore, concentration inequalities, i.e. bounds on moment generating functions, provide tight and computationally inexpensive model bias bounds for quantities of interest. Finally, they allow us to derive rigorous confidence bands for statistical estimators that account for model bias and are valid for an arbitrary amount of data.


title: Clock Monte Carlo methods

abstract: We propose the clock Monte Carlo technique for sampling each successive chain step in constant time. It is built on a recently proposed factorized transition filter and its core features include its O(1) computational complexity and its generality. We elaborate how it leads to the clock factorized Metropolis (clock FMet) method, and discuss its application in other update schemes. By grouping interaction terms into boxes of tunable sizes, we further formulate a variant of the clock FMet algorithm, with the limiting case of a single box reducing to the standard Metropolis method. A theoretical analysis shows that an overall acceleration of ${\rm O}(N^\kappa)$ ($0 \! \leq \! \kappa \! \leq \! 1$) can be achieved compared to the Metropolis method, where $N$ is the system size and the $\kappa$ value depends on the nature of the energy extensivity. As a systematic test, we simulate long-range O$(n)$ spin models in a wide parameter regime: for $n \! = \! 1,2,3$, with disordered algebraically decaying or oscillatory Ruderman-Kittel-Kasuya-Yoshida-type interactions and with and without external fields, and in spatial dimensions from $d \! = \! 1, 2, 3$ to mean-field. The O(1) computational complexity is demonstrated, and the expected acceleration is confirmed. Its flexibility and its independence from the interaction range guarantee that the clock method would find decisive applications in systems with many interaction terms.


title: Anisotropic superfluidity of two-dimensional excitons in a periodic potential

abstract: We study anisotropies of helicity modulus, excitation spectrum, sound velocity and angle-resolved luminescence spectrum in a two-dimensional system of interacting excitons in a periodic potential. Analytical expressions for anisotropic corrections to the quantities characterizing superfluidity are obtained. We consider particularly the case of dipolar excitons in quantum wells. For GaAs/AlGaAs heterostructures as well as MoS$_2$/hBN/MoS$_2$ and MoSe$_2$/hBN/WSe$_2$ transition metal dichalcogenide bilayers estimates of the magnitude of the predicted effects are given. We also present a method to control superfluid motion and to determine the helicity modulus in generic dipolar systems.


title: A smoothed particle hydrodynamics method for evaporating multiphase flows

abstract: Smoothed particle hydrodynamics (SPH) method has been increasingly used for simulating fluid flows, however its ability to simulate evaporating flow requires significant improvements. This paper proposes an SPH method for evaporating multiphase flows. The present SPH method can simulate the heat and mass transfers across the liquid-gas interfaces. The conservation equations of mass, momentum and energy were reformulated based on SPH, then were used to govern the fluid flow and heat transfer in both the liquid and gas phases. The continuity equation of the vapor species was employed to simulate the vapor mass fraction in the gas phase. The vapor mass fraction at the interface was predicted by the Clausius-Clapeyron correlation. A new evaporation rate was derived to predict the mass transfer from the liquid phase to the gas phase at the interface. Because of the mass transfer across the liquid-gas interface, the mass of an SPH particle was allowed to change. New particle splitting and merging techniques were developed to avoid large mass difference between SPH particles of the same phase. The proposed method was tested by simulating three problems, including the Stefan problem, evaporation of a static drop, and evaporation of a drop impacting on a hot surface. For the Stefan problem, the SPH results of the evaporation rate at the interface agreed well with the analytical solution. For drop evaporation, the SPH result was compared with the result predicted by a level-set method from literature. In the case of drop impact on a hot surface, the evolution of the shape of the drop, temperature, and vapor mass fraction were predicted.


title: Non-degeneracy of Gauss curvature equation with negative conic singularity

abstract: We study the Gauss curvature equation with negative singularities. For a local mean field type equation with only one negative index we prove a uniqueness property. For a global equation with one or two negative indexes we prove the non-degeneracy of the linearized equations.


title: $q$-Poincar\'e supersymmetry in $AdS_5/CFT_4$

abstract: We consider the exact S-matrix governing the planar spectral problem for strings on $AdS_5\times S^5$ and $\mathcal N=4$ super Yang-Mills, and we show that it is invariant under a novel "boost" symmetry, which acts as a differentiation with respect to the particle momentum. This generator leads us also to reinterpret the usual centrally extended $\mathfrak{psu}(2|2)$ symmetry, and to conclude that the S-matrix is invariant under a $q$-Poincar\'e supersymmetry algebra, where the deformation parameter is related to the 't Hooft coupling. We determine the two-particle action (coproduct) that turns out to be non-local, and study the property of the new symmetry under crossing transformations. We look at both the strong-coupling (large tension in the string theory) and weak-coupling (spin-chain description of the gauge theory) limits; in the former regime we calculate the cobracket utilising the universal classical r-matrix of Beisert and Spill. In the eventuality that the boost has higher partners, we also construct a quantum affine version of 2D Poincar\'e symmetry, by contraction of the quantum affine algebra $U_q(\widehat{\mathfrak{sl}_2})$ in Drinfeld's second realisation.


title: Color-opponent mechanisms for local hue encoding in a hierarchical framework

abstract: A biologically plausible computational model for color representation is introduced. We present a mechanistic hierarchical model of neurons that not only successfully encodes local hue, but also explicitly reveals how the contributions of each visual cortical layer participating in the process can lead to a hue representation. Our proposed model benefits from studies on the visual cortex and builds a network of single-opponent and hue-selective neurons. Local hue encoding is achieved through gradually increasing nonlinearity in terms of cone inputs to single-opponent cells. We demonstrate that our model's single-opponent neurons have wide tuning curves, while the hue-selective neurons in our model V4 layer exhibit narrower tunings, resembling those in V4 of the primate visual system. Our simulation experiments suggest that neurons in V4 or later layers have the capacity of encoding unique hues. Moreover, with a few examples, we present the possibility of spanning the infinite space of physical hues by combining the hue-selective neurons in our model.


title: Challenges to Keeping the Computer Industry Centered in the US

abstract: It is undeniable that the worldwide computer industry's center is the US, specifically in Silicon Valley. Much of the reason for the success of Silicon Valley had to do with Moore's Law: the observation by Intel co-founder Gordon Moore that the number of transistors on a microchip doubled at a rate of approximately every two years. According to the International Technology Roadmap for Semiconductors, Moore's Law will end in 2021. How can we rethink computing technology to restart the historic explosive performance growth? Since 2012, the IEEE Rebooting Computing Initiative (IEEE RCI) has been working with industry and the US government to find new computing approaches to answer this question. In parallel, the CCC has held a number of workshops addressing similar questions. This whitepaper summarizes some of the IEEE RCI and CCC findings. The challenge for the US is to lead this new era of computing. Our international competitors are not sitting still: China has invested significantly in a variety of approaches such as neuromorphic computing, chip fabrication facilities, computer architecture, and high-performance simulation and data analytics computing, for example. We must act now, otherwise, the center of the computer industry will move from Silicon Valley and likely move off shore entirely.


title: SafetyNets: Verifiable Execution of Deep Neural Networks on an Untrusted Cloud

abstract: Inference using deep neural networks is often outsourced to the cloud since it is a computationally demanding task. However, this raises a fundamental issue of trust. How can a client be sure that the cloud has performed inference correctly? A lazy cloud provider might use a simpler but less accurate model to reduce its own computational load, or worse, maliciously modify the inference results sent to the client. We propose SafetyNets, a framework that enables an untrusted server (the cloud) to provide a client with a short mathematical proof of the correctness of inference tasks that they perform on behalf of the client. Specifically, SafetyNets develops and implements a specialized interactive proof (IP) protocol for verifiable execution of a class of deep neural networks, i.e., those that can be represented as arithmetic circuits. Our empirical results on three- and four-layer deep neural networks demonstrate the run-time costs of SafetyNets for both the client and server are low. SafetyNets detects any incorrect computations of the neural network by the untrusted server with high probability, while achieving state-of-the-art accuracy on the MNIST digit recognition (99.4%) and TIMIT speech recognition tasks (75.22%).


title: A formalization of convex polyhedra based on the simplex method

abstract: We present a formalization of convex polyhedra in the proof assistant Coq. The cornerstone of our work is a complete implementation of the simplex method, together with the proof of its correctness and termination. This allows us to define the basic predicates over polyhedra in an effective way (i.e., as programs), and relate them with the corresponding usual logical counterparts. To this end, we make an extensive use of the Boolean reflection methodology. The benefit of this approach is that we can easily derive the proof of several fundamental results on polyhedra, such as Farkas' Lemma, the duality theorem of linear programming, and Minkowski's Theorem.


title: The entropic approach to causal correlations

abstract: The existence of a global causal order between events places constraints on the correlations that parties may share. Such "causal correlations" have been the focus of recent attention, driven by the realization that some extensions of quantum mechanics may violate so-called causal inequalities. In this paper we study causal correlations from an entropic perspective, and we show how to use this framework to derive entropic causal inequalities. We consider two different ways to derive such inequalities. Firstly, we consider a method based on the causal Bayesian networks describing the causal relations between the parties. In contrast to the Bell-nonlocality scenario, where this method has previously been shown to be ineffective, we show that it leads to several interesting entropic causal inequalities. Secondly, we consider an alternative method based on counterfactual variables that has previously been used to derive entropic Bell inequalities. We compare the inequalities obtained via these two methods and discuss their violation by noncausal correlations. As an application of our approach, we derive bounds on the quantity of information - which is more naturally expressed in the entropic framework - that parties can communicate when operating in a definite causal order.


title: Lifelong Learning in Costly Feature Spaces

abstract: An important long-term goal in machine learning systems is to build learning agents that, like humans, can learn many tasks over their lifetime, and moreover use information from these tasks to improve their ability to do so efficiently. In this work, our goal is to provide new theoretical insights into the potential of this paradigm. In particular, we propose a lifelong learning framework that adheres to a novel notion of resource efficiency that is critical in many real-world domains where feature evaluations are costly. That is, our learner aims to reuse information from previously learned related tasks to learn future tasks in a feature-efficient manner. Furthermore, we consider novel combinatorial ways in which learning tasks can relate. Specifically, we design lifelong learning algorithms for two structurally different and widely used families of target functions: decision trees/lists and monomials/polynomials. We also provide strong feature-efficiency guarantees for these algorithms; in fact, we show that in order to learn future targets, we need only slightly more feature evaluations per training example than what is needed to predict on an arbitrary example using those targets. We also provide algorithms with guarantees in an agnostic model where not all the targets are related to each other. Finally, we also provide lower bounds on the performance of a lifelong learner in these models, which are in fact tight under some conditions.


title: Nuclear penalized multinomial regression with an application to predicting at bat outcomes in baseball

abstract: We propose the nuclear norm penalty as an alternative to the ridge penalty for regularized multinomial regression. This convex relaxation of reduced-rank multinomial regression has the advantage of leveraging underlying structure among the response categories to make better predictions. We apply our method, nuclear penalized multinomial regression (NPMR), to Major League Baseball play-by-play data to predict outcome probabilities based on batter-pitcher matchups. The interpretation of the results meshes well with subject-area expertise and also suggests a novel understanding of what differentiates players.


title: Community Detection by $L_0$-penalized Graph Laplacian

abstract: Community detection in network analysis aims at partitioning nodes in a network into $K$ disjoint communities. Most currently available algorithms assume that $K$ is known, but choosing a correct $K$ is generally very difficult for real networks. In addition, many real networks contain outlier nodes not belonging to any community, but currently very few algorithm can handle networks with outliers. In this paper, we propose a novel model free tightness criterion and an efficient algorithm to maximize this criterion for community detection. This tightness criterion is closely related with the graph Laplacian with $L_0$ penalty. Unlike most community detection methods, our method does not require a known $K$ and can properly detect communities in networks with outliers. Both theoretical and numerical properties of the method are analyzed. The theoretical result guarantees that, under the degree corrected stochastic block model, even for networks with outliers, the maximizer of the tightness criterion can extract communities with small misclassification rates even when the number of communities grows to infinity as the network size grows. Simulation study shows that the proposed method can recover true communities more accurately than other methods. Applications to a college football data and a yeast protein-protein interaction data also reveal that the proposed method performs significantly better.


title: A Model for Attribute Based Role-Role Assignment (ARRA)

abstract: Administrative Role Based Access Control (ARBAC) models specify how to manage user-role assignments (URA), permission-role assignments (PRA), and role-role assignments (RRA). Many approaches have been proposed in the literature for URA, PRA, and RRA. In this paper, we propose a model for attribute-based role-role assignment (ARRA), a novel way to unify prior RRA approaches. We leverage the idea that attributes of various RBAC entities such as admin users and regular roles can be used to administer RRA in a highly flexible manner. We demonstrate that ARRA can express and unify prior RRA models.


title: Signal Reconstruction from Interferometric Measurements under Sensing Constraints

abstract: This paper develops a unifying framework for signal reconstruction from interferometric measurements that is broadly applicable to various applications of interferometry. In this framework, the problem of signal reconstruction in interferometry amounts to one of basis analysis. Its applicability is shown to extend beyond conventional temporal interferometry, which leverages the relative delay between the two arms of an interferometer, to arbitrary degrees of freedom of the input signal. This allows for reconstruction of signals supported in other domains (e.g., spatial) with no modification to the underlying structure except for replacing the standard temporal delay with a generalized delay, that is, a practically realizable unitary transformation for which the basis elements are eigenfunctions. Under the proposed model, the interferometric measurements are shown to be linear in the basis coefficients, thereby enabling efficient and fast recovery of the desired information. While the corresponding linear transformation has only a limited number of degrees of freedom set by the structure of the interferometer giving rise to a highly constrained sensing structure, we show that the problem of signal recovery from such measurements can still be carried out compressively. This signifies significant reduction in sample complexity without introducing any additional randomization as is typically done in prior work leveraging compressive sensing techniques. We provide performance guarantees under constrained sensing by proving that the transformation satisfies sufficient conditions for successful reconstruction of sparse signals using concentration arguments. We showcase the effectiveness of the proposed approach using simulation results, as well as actual experimental results in the context of optical modal analysis of spatial beams.


title: DataLair: Efficient Block Storage with Plausible Deniability against Multi-Snapshot Adversaries

abstract: Sensitive information is present on our phones, disks, watches and computers. Its protection is essential. Plausible deniability of stored data allows individuals to deny that their device contains a piece of sensitive information. This constitutes a key tool in the fight against oppressive governments and censorship. Unfortunately, existing solutions, such as the now defunct TrueCrypt, can defend only against an adversary that can access a users device at most once (single-snapshot adversary). Recent solutions have traded significant performance overheads for the ability to handle more powerful adversaries, that are able to access the device at multiple points in time (multi-snapshot adversary). In this paper we show that this sacrifice is not necessary. We introduce and build DataLair, a practical plausible deniability mechanism. When compared with existing approaches, DataLair is two orders of magnitude faster for public data accesses, and 5 times faster for hidden data accesses. An important component in DataLair is a new write-only ORAM construction which improves on the complexity of the state of the art write-only ORAM by a factor of O(logN ), where N denotes the underlying storage disk size.


title: CO2 packing polymorphism under pressure: mechanism and thermodynamics of the I-III polymorphic transition

abstract: In this work we describe the thermodynamics and mechanism of CO$_2$ polymorphic transitions under pressure from form I to form III combining standard molecular dynamics, well-tempered metadynamics and committor analysis. We find that the phase transformation takes place through a concerted rearrangement of CO$_2$ molecules, which unfolds via an anisotropic expansion of the CO$_2$ supercell. Furthermore, at high pressures we find that defected form I configurations are thermodynamically more stable with respect to form I without structural defects. Our computational approach shows the capability of simultaneously providing an extensive sampling of the configurational space, estimates of the thermodynamic stability and a suitable description of a complex, collective polymorphic transition mechanism.


title: Computing Biomolecular System Steady-states

abstract: A new approach to computing the equilibria and steady-states of biomolecular systems modelled by bond graphs is presented. The approach is illustrated using a model of a biomolecular cycle representing a membrane transporter and a model of the mitochondrial electron transport chain.


title: Extreme magnification of a star at redshift 1.5 by a galaxy-cluster lens

abstract: Galaxy-cluster gravitational lenses can magnify background galaxies by a total factor of up to ~50. Here we report an image of an individual star at redshift z=1.49 (dubbed "MACS J1149 Lensed Star 1 (LS1)") magnified by >2000. A separate image, detected briefly 0.26 arcseconds from LS1, is likely a counterimage of the first star demagnified for multiple years by a >~3 solar-mass object in the cluster. For reasonable assumptions about the lensing system, microlensing fluctuations in the stars' light curves can yield evidence about the mass function of intracluster stars and compact objects, including binary fractions and specific stellar evolution and supernova models. Dark-matter subhalos or massive compact objects may help to account for the two images' long-term brightness ratio.


title: Ultracold heteronuclear three-body systems: How diabaticity limits the universality of recombination into shallow dimers

abstract: The mass-imbalanced three-body recombination process that forms a shallow dimer is shown to possess a rich Efimov-St\"uckelberg landscape, with corresponding spectra that differ fundamentally from the homonuclear case. A semi-analytical treatment of the three-body recombination predicts an unusual spectra with intertwined resonance peaks and minima, and yields in-depth insight into the behavior of the corresponding Efimov spectra. In particular, the patterns of the Efimov-St\"uckelberg landscape are shown to depend inherently on the degree of diabaticity of the three-body collisions, which strongly affects the universality of the heteronuclear Efimov states.


title: Dark matter under the microscope: Constraining compact dark matter with caustic crossing events

abstract: A galaxy cluster acts as a cosmic telescope over background galaxies but also as a cosmic microscope of the lens imperfections. The diverging magnification of lensing caustics enhances the microlensing effect of substructure present within the lensing mass. Fine-scale structure can be accessed as a moving background source brightens and disappears when crossing these caustics. The recent recognition of a distant lensed star near the Einstein radius of the galaxy cluster MACSJ1149.5+2223 (Kelly et al. 2017) allows the rare opportunity to reach subsolar mass microlensing through a super-critical column of cluster matter. Here we compare these observations with high-resolution ray-tracing simulations that include stellar microlensing set by the observed intracluster starlight and also primordial black holes that may be responsible for the recently observed LIGO events. We explore different scenarios with microlenses from the intracluster medium and black holes, including primordial ones, and examine strategies to exploit these unique alignments. We find that the best constraints on the fraction of compact dark matter in the small-mass regime can be obtained in regions of the cluster where the intracluster medium plays a negligible role. This new lensing phenomenon should be widespread and can be detected within modest-redshift lensed galaxies so that the luminosity distance is not prohibitive for detecting individual magnified stars. Continuous {\it Hubble Space Telescope} monitoring of several such optimal arcs will be rewarded by an unprecedented mass spectrum of compact objects that can contribute to uncovering the nature of dark matter.


title: Stellar Chemical Clues As To The Rarity of Exoplanetary Tectonics

abstract: Earth's tectonic processes regulate the formation of continental crust, control its unique deep water and carbon cycles, and are vital to its surface habitability. A major driver of steady-state plate tectonics on Earth is the sinking of the cold subducting plate into the underlying mantle. This sinking is the result of the combined effects of the thermal contraction of the lithosphere and of metamorphic transitions within the basaltic oceanic crust and lithospheric mantle. The latter of these effects is dependent on the bulk composition of the planet, e.g., the major, terrestrial planet-building elements Mg, Si, Fe, Ca, Al, and Na, which vary in abundance across the Galaxy. We present thermodynamic phase-equilibria calculations of planetary differentiation to calculate both melt composition and mantle mineralogy, and show that a planet's refractory and moderately-volatile elemental abundances control a terrestrial planet's likelihood to produce mantle-derived, melt-extracted crusts that sink. Those planets forming with a higher concentration of Si and Na abundances are less likely to undergo sustained tectonics compared to the Earth. We find only 1/3 of the range of stellar compositions observed in the Galaxy is likely to host planets able to sustain density-driven tectonics compared to the Sun/Earth. Systems outside of this compositional range are less likely to produce planets able to tectonically regulate their climate and may be inhospitable to life as we know it.


title: Bolt: Accelerated Data Mining with Fast Vector Compression

abstract: Vectors of data are at the heart of machine learning and data mining. Recently, vector quantization methods have shown great promise in reducing both the time and space costs of operating on vectors. We introduce a vector quantization algorithm that can compress vectors over 12x faster than existing techniques while also accelerating approximate vector operations such as distance and dot product computations by up to 10x. Because it can encode over 2GB of vectors per second, it makes vector quantization cheap enough to employ in many more circumstances. For example, using our technique to compute approximate dot products in a nested loop can multiply matrices faster than a state-of-the-art BLAS implementation, even when our algorithm must first compress the matrices. In addition to showing the above speedups, we demonstrate that our approach can accelerate nearest neighbor search and maximum inner product search by over 100x compared to floating point operations and up to 10x compared to other vector quantization methods. Our approximate Euclidean distance and dot product computations are not only faster than those of related algorithms with slower encodings, but also faster than Hamming distance computations, which have direct hardware support on the tested platforms. We also assess the errors of our algorithm's approximate distances and dot products, and find that it is competitive with existing, slower vector quantization algorithms.


title: Cache Oblivious Algorithms for Computing the Triplet Distance Between Trees

abstract: We study the problem of computing the triplet distance between two rooted unordered trees with $n$ labeled leafs. Introduced by Dobson 1975, the triplet distance is the number of leaf triples that induce different topologies in the two trees. The current theoretically best algorithm is an $\mathrm{O}(n \log n)$ time algorithm by Brodal et al. (SODA 2013). Recently Jansson and Rajaby proposed a new algorithm that, while slower in theory, requiring $\mathrm{O}(n \log^3 n)$ time, in practice it outperforms the theoretically faster $\mathrm{O}(n \log n)$ algorithm. Both algorithms do not scale to external memory. We present two cache oblivious algorithms that combine the best of both worlds. The first algorithm is for the case when the two input trees are binary trees and the second a generalized algorithm for two input trees of arbitrary degree. Analyzed in the RAM model, both algorithms require $\mathrm{O}(n \log n)$ time, and in the cache oblivious model $\mathrm{O}(\frac{n}{B} \log_{2} \frac{n}{M})$ I/Os. Their relative simplicity and the fact that they scale to external memory makes them achieve the best practical performance. We note that these are the first algorithms that scale to external memory, both in theory and practice, for this problem.


title: Probabilistic estimation of the rank 1 cross approximation accuracy

abstract: In the construction of low-rank matrix approximation and maximum element search it is effective to use maxvol algorithm. Nevertheless, even in the case of rank 1 approximation the algorithm does not always converge to the maximum matrix element, and it is unclear how often close to the maximum element can be found. In this article it is shown that with a certain degree of randomness in the matrix and proper selection of the starting column, the algorithm with high probability in a few steps converges to an element, which module differs little from the maximum. It is also shown that with more severe restrictions on the error matrix no restrictions on the starting column need to be introduced.


title: A Comparison of Cosmological Parameters Determined from CMB Temperature Power Spectra from the South Pole Telescope and the Planck Satellite

abstract: The Planck cosmic microwave background (CMB) temperature data are best fit with a LCDM model that is in mild tension with constraints from other cosmological probes. The South Pole Telescope (SPT) 2540 $\text{deg}^2$ SPT-SZ survey offers measurements on sub-degree angular scales (multipoles $650 \leq \ell \leq 2500$) with sufficient precision to use as an independent check of the Planck data. Here we build on the recent joint analysis of the SPT-SZ and Planck data in \citet{hou17} by comparing LCDM parameter estimates using the temperature power spectrum from both data sets in the SPT-SZ survey region. We also restrict the multipole range used in parameter fitting to focus on modes measured well by both SPT and Planck, thereby greatly reducing sample variance as a driver of parameter differences and creating a stringent test for systematic errors. We find no evidence of systematic errors from such tests. When we expand the maximum multipole of SPT data used, we see low-significance shifts in the angular scale of the sound horizon and the physical baryon and cold dark matter densities, with a resulting trend to higher Hubble constant. When we compare SPT and Planck data on the SPT-SZ sky patch to Planck full-sky data but keep the multipole range restricted, we find differences in the parameters $n_s$ and $A_se^{-2\tau}$. We perform further checks, investigating instrumental effects and modeling assumptions, and we find no evidence that the effects investigated are responsible for any of the parameter shifts. Taken together, these tests reveal no evidence for systematic errors in SPT or Planck data in the overlapping sky coverage and multipole range and, at most, weak evidence for a breakdown of LCDM or systematic errors influencing either the Planck data outside the SPT-SZ survey area or the SPT data at $\ell >2000$.


title: Generation of turbulence through frontogenesis in sheared stratified flows

abstract: The large-scale structures in the ocean and the atmosphere are in geostrophic balance, and a conduit must be found to channel the energy to the small scales where it can be dissipated. In turbulence this takes the form of an energy cascade, whereas one possible mechanism in a balanced flow at large scales is through the formation of fronts, a common occurrence in geophysical dynamics. We show in this paper that an iconic configuration in laboratory and numerical experiments for the study of turbulence, that of the Taylor-Green or von K\'arm\'an swirling flow, can be suitably adapted to the case of fluids with large aspect ratios, leading to the creation of an imposed large-scale vertical shear. To this effect we use direct numerical simulations of the Boussinesq equations without net rotation and with no small-scale modeling, and with this idealized Taylor-Green set-up. Various grid spacings are used, up to $2048^2\times 256$ spatial points. The grids are always isotropic, with box aspect ratios of either $1:4$ or $1:8$. We find that when shear and stratification are comparable, the imposed shear layer resulting from the forcing leads to the formation of multiple fronts and filaments which destabilize and further evolve into a turbulent flow in the bulk, with a sizable amount of dissipation and mixing, and with a cycle of front creation, instability, and development of turbulence. The results depend on the vertical length scales for shear and for stratification, with stronger large-scale gradients being generated when the two length scales are comparable.