title: Stability of stagnation via an expanding accretion shock wave

abstract: Stagnation of a cold plasma streaming to the center or axis of symmetry via an expanding accretion shock wave is ubiquitous in inertial confinement fusion (ICF) and high-energy-density plasma physics, the examples ranging from plasma flows in x-ray-generating Z pinches [Y. Maron et al., Phys. Rev. Lett. 111, 035001 (2013)] to the experiments in support of the recently suggested concept of impact ignition in ICF [H. Azechi et al., Phys. Rev. Lett. 102, 235002 (2009); M. Murakami et al., Nucl. Fusion 54, 054007 (2014)]. Some experimental evidence indicates that stagnation via an expanding shock wave is stable, but its stability has never been studied theoretically. We present such analysis for the stagnation that does not involve a rarefaction wave behind the expanding shock front and is described by the classic ideal-gas Noh solution in spherical and cylindrical geometry. In either case the stagnated flow has been demonstrated to be stable, initial perturbations exhibiting a power-law, oscillatory or monotonic, decay with time for all the eigenmodes. This conclusion has been supported by our simulations done both on a Cartesian grid and on a curvilinear grid in spherical coordinates. Dispersion equation determining the eigenvalues of the problem and explicit formulas for the eigenfunction profiles corresponding to these eigenvalues are presented, making it possible to use the theory for hydro code verification in two and three dimensions.


title: Crossing symmetry and Higher spin towers

abstract: We consider higher spin operators in weakly coupled gauge conformal field theories. Crossing symmetry of mixed scalar correlators relates different higher spin towers and we study the consequences for the spectrum and structure constants of higher spin operators of different twists. Constraints are obtained to all loops in perturbation theory. The large spin contributions to the structure constants can be resummed into a theory-dependent prefactor times a universal factor, whose structure of poles agrees with the one that would be obtained from a Witten diagram supergravity computation, although only crossing symmetry is assumed. Finally, our results provide an all loop expression for the double null limit of mixed correlators, which is in perfect agreement with the correlator/Wilson loop correspondence.


title: A calculus for the moduli space of curves

abstract: This article accompanies my lecture at the 2015 AMS summer institute in algebraic geometry in Salt Lake City. I survey the recent advances in the study of tautological classes on the moduli spaces of curves. After discussing the Faber-Zagier relations on the moduli spaces of nonsingular curves and the kappa rings of the moduli spaces of curves of compact type, I present Pixton's proposal for a complete calculus of tautological classes on the moduli spaces of stable curves. Several open questions are discussed. An effort has been made to condense a great deal of mathematics into as few pages as possible with the hope that the reader will follow through to the end.


title: Feature Selection as a Multiagent Coordination Problem

abstract: Datasets with hundreds to tens of thousands features is the new norm. Feature selection constitutes a central problem in machine learning, where the aim is to derive a representative set of features from which to construct a classification (or prediction) model for a specific task. Our experimental study involves microarray gene expression datasets, these are high-dimensional and noisy datasets that contain genetic data typically used for distinguishing between benign or malicious tissues or classifying different types of cancer. In this paper, we formulate feature selection as a multiagent coordination problem and propose a novel feature selection method using multiagent reinforcement learning. The central idea of the proposed approach is to "assign" a reinforcement learning agent to each feature where each agent learns to control a single feature, we refer to this approach as MARL. Applying this to microarray datasets creates an enormous multiagent coordination problem between thousands of learning agents. To address the scalability challenge we apply a form of reward shaping called CLEAN rewards. We compare in total nine feature selection methods, including state-of-the-art methods, and show that the proposed method using CLEAN rewards can significantly scale-up, thus outperforming the rest of learning-based methods. We further show that a hybrid variant of MARL achieves the best overall performance.


title: Molecular dynamics simulation of nanoindentation on nanocomposite pearlite

abstract: We carry out molecular dynamics simulations of nanoindentation to investigate the effect of cementite size and temperature on the deformation behavior of nanocomposite pearlite composed of alternating ferrite and cementite layers. We find that, instead of the coherent transmission, dislocation propagates by forming a widespread plastic deformation in cementite layer. We also show that increasing temperature enhances the distribution of plastic strain in the ferrite layer, which reduces the stress acting on the cementite layer. Hence, thickening cementite layer or increasing temperature reduces the likelihood of dislocation propagation through the cementite layer. Our finding sheds a light on the mechanism of dislocation blocking by cementite layer in the pearlite.


title: 2D Discrete Fourier Transform with Simultaneous Edge Artifact Removal for Real-Time Applications

abstract: Two-Dimensional (2D) Discrete Fourier Transform (DFT) is a basic and computationally intensive algorithm, with a vast variety of applications. 2D images are, in general, non-periodic, but are assumed to be periodic while calculating their DFTs. This leads to cross-shaped artifacts in the frequency domain due to spectral leakage. These artifacts can have critical consequences if the DFTs are being used for further processing. In this paper we present a novel FPGA-based design to calculate high-throughput 2D DFTs with simultaneous edge artifact removal. Standard approaches for removing these artifacts using apodization functions or mirroring, either involve removing critical frequencies or a surge in computation by increasing image size. We use a periodic-plus-smooth decomposition based artifact removal algorithm optimized for FPGA implementation, while still achieving real-time ($\ge$23 frames per second) performance for a 512$\times$512 size image stream. Our optimization approach leads to a significant decrease in external memory utilization thereby avoiding memory conflicts and simplifies the design. We have tested our design on a PXIe based Xilinx Kintex 7 FPGA system communicating with a host PC which gives us the advantage to further expand the design for industrial applications.


title: Periods of the motivic fundamental groupoid of $\boldsymbol{\mathbb{P}^{1} \diagdown \lbrace 0, \mu_{N}, \infty \rbrace}$

abstract: In this thesis, following F. Brown's point of view, we look at the Hopf algebra structure of motivic cyclotomic multiple zeta values, which are motivic periods of the fundamental groupoid of $\mathbb{P}^{1} \diagdown \lbrace 0, \mu_{N}, \infty \rbrace $. By application of a surjective \textit{period} map (which, under Grothendieck's period conjecture, is an isomorphism), we deduce results (such as generating families, identities, etc.) on cyclotomic multiple zeta values, which are complex numbers. The coaction of this Hopf algebra (explicitly given by a combinatorial formula from A. Goncharov and F. Brown's works) is the dual of the action of a so-called \textit{motivic} Galois group on these specific motivic periods. This entire study was actually motivated by the hope of a Galois theory for periods, which should extend the usual Galois theory for algebraic numbers.


title: Continuous-wave and Transient Characteristics of Phosphorene Microwave Transistors

abstract: Few-layer phosphorene MOSFETs with 0.3-um-long gate and 15-nm-thick Al2O3 gate insulator was found to exhibit a forward-current cutoff frequency of 2 GHz and a maximum oscillation frequency of 8 GHz after de-embedding for the parasitic capacitance associated mainly with the relatively large probe pads. The gate lag and drain lag of the transistor was found to be on the order of 1 us or less, which is consistent with the lack of hysteresis, carrier freeze-out or persistent photoconductivity in DC characteristics. These results confirm that the phosphorene MOSFET can be a viable microwave transistor for both small-signal and large-signal applications.


title: Comparing Convolutional Neural Networks to Traditional Models for Slot Filling

abstract: We address relation classification in the context of slot filling, the task of finding and evaluating fillers like "Steve Jobs" for the slot X in "X founded Apple". We propose a convolutional neural network which splits the input sentence into three parts according to the relation arguments and compare it to state-of-the-art and traditional approaches of relation classification. Finally, we combine different methods and show that the combination is better than individual approaches. We also analyze the effect of genre differences on performance.


title: A Landau-de Gennes theory for hard colloidal rods: defects and tactoids

abstract: We construct a phenomenological Landau-de Gennes theory for hard colloidal rods by performing an order parameter expansion of the chemical-potential dependent grand potential. By fitting the coefficients to known results of Onsager theory, we are not only able to describe the isotropic-nematic phase transition as function of density, including the well-known density jump, but also the isotropic-nematic planar interface. The resulting theory is applied in calculations of the isotropic core size in a radial hedgehog defect, the density dependence of linear defects of hard rods in square confinement, and the formation of a nematic droplet in an isotropic background.


title: Polarization preservation of partially coherent Hermite-Gaussian beams for multiple-degrees-of-freedom free-space communication

abstract: Multiple-degrees-of-freedom free-space communication combining polarization and high-order spatial modes promises high-capacity communication channel. While high-order spatial modes have been widely exploited for dense coding and high-dimensional quantum information processing, the properties of polarization preservation of high-order spatial beams propagating in turbulent atmosphere have not been comprehensively investigated yet. Here we focus on the properties of polarization preservation of partially coherent Hermite-Gaussian beams propagating along different atmospheric turbulence paths. The analytical expressions for the polarization of partially coherent Hermite-Gaussian beams propagating through atmospheric turbulence along different paths have been derived. It is shown that the larger the coherence length is, and the larger the beam order m, n are, the less the polarization is changed. We find that the evolution properties of the polarization in slant-down paths through turbulent atmosphere are similar to the case in free space if the condition zenith angle {\xi}<{\pi}/4 is satisfied. While at a long propagation distance, evolution properties of polarization in horizontal paths of turbulent atmosphere differs much from that in free space and in slant paths. The results may allow one to choose the optimal propagation path in terms of specific applications, which is helpful for future experimental implementation of multiple-degrees-of-freedom free-space communication.


title: Recent Results from SPLASH: Chemical Abundances and Kinematics of Andromeda's Stellar Halo

abstract: Large scale surveys of Andromeda's resolved stellar populations have revolutionized our view of this galaxy over the past decade. The combination of large-scale, contiguous photometric surveys and pointed spectroscopic surveys has been particularly powerful for discovering substructure and disentangling the structural components of Andromeda. The SPLASH (Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo) survey consists of broad- and narrow-band imaging and spectroscopy of red giant branch stars in lines of sight ranging in distance from 2 kpc to more than 200 kpc from Andromeda's center. The SPLASH data reveal a power-law surface brightness profile extending to at least two-thirds of Andromeda's virial radius (Gilbert et al. 2012), a metallicity gradient extending to at least 100 kpc from Andromeda's center (Gilbert et al. 2014), and evidence of a significant population of heated disk stars in Andromeda's inner halo (Dorman et al. 2013). We are also using the velocity distribution of halo stars to measure the tangential motion of Andromeda (Beaton et al., in prep).


title: Dimension transformation formula for conformal maps into the complement of an SLE curve

abstract: We prove a formula relating the Hausdorff dimension of a deterministic Borel subset of $\mathbb R$ and the Hausdorff dimension of its image under a conformal map from the upper half-plane to a complementary connected component of an SLE$_\kappa$ curve for $\kappa \not =4$. Our proof is based on the relationship between SLE and Liouville quantum gravity together with the one-dimensional KPZ formula of Rhodes-Vargas (2011) and the KPZ formula of Gwynne-Holden-Miller (2015). As an intermediate step we prove a KPZ formula which relates the Euclidean dimension of a subset of an SLE$_\kappa$ curve for $\kappa \in (0,4)\cup(4,8)$ and the dimension of the same set with respect to the $\gamma$-quantum natural parameterization of the curve induced by an independent Gaussian free field, $\gamma = \sqrt \kappa \wedge (4/\sqrt\kappa)$.


title: A (Basis for a) Philosophy of a Theory of Fuzzy Computation

abstract: Vagueness is a linguistic phenomenon as well as a property of physical objects. Fuzzy set theory is a mathematical model of vagueness that has been used to define vague models of computation. The prominent model of vague computation is the fuzzy Turing machine. This conceptual computing device gives an idea of what computing under vagueness means, nevertheless, it is not the most natural model. Based on the properties of this and other models of vague computing, it is aimed to formulate a basis for a philosophy of a theory of fuzzy computation.


title: Accelerating Data Regeneration for Distributed Storage Systems with Heterogeneous Link Capacities

abstract: Distributed storage systems provide large-scale reliable data storage services by spreading redundancy across a large group of storage nodes. In such a large system, node failures take place on a regular basis. When a storage node breaks down, a replacement node is expected to regenerate the redundant data as soon as possible in order to maintain the same level of redundancy. Previous results have been mainly focused on the minimization of network traffic in regeneration. However, in practical networks, where link capacities vary in a wide range, minimizing network traffic does not always yield the minimum regeneration time. In this paper, we investigate two approaches to the problem of minimizing regeneration time in networks with heterogeneous link capacities. The first approach is to download different amounts of repair data from the helping nodes according to the link capacities. The second approach generalizes the conventional star-structured regeneration topology to tree-structured topologies so that we can utilize the links between helping nodes with bypassing low-capacity links. Simulation results show that the flexible tree-structured regeneration scheme that combines the advantages of both approaches can achieve a substantial reduction in the regeneration time.


title: Robust coupling of DPG and BEM for a singularly perturbed transmission problem

abstract: We consider a transmission problem consisting of a singularly perturbed reaction diffusion equation on a bounded domain and the Laplacian in the exterior, connected through standard transmission conditions. We establish a DPG scheme coupled with Galerkin boundary elements for its discretization, and prove its robustness for the field variables in so-called balanced norms. Our coupling scheme is the one from [F\"uhrer, Heuer, Karkulik: On the coupling of DPG and BEM, arXiv:1508.00630], adapted to the singularly perturbed case by using the scheme from [Heuer, Karkulik: A robust DPG method for singularly perturbed reaction diffusion problems, arXiv:1509.07560]. Essential feature of our method is that optimal test functions have to be computed only locally. We report on various numerical experiments in two dimensions.


title: Non-Markovianity in atom-surface dispersion forces

abstract: We discuss the failure of the Markov approximation in the description of atom-surface fluctuation-induced interactions, both at equilibrium (Casimir-Polder forces) and out-of-equilibrium (quantum friction). Using general theoretical arguments, we show that the Markov approximation can lead to erroneous predictions of such phenomena with regard to both strength and functional dependencies on system parameters. Our findings highlight the importance of non-Markovian effects in dispersion interactions. In particular, we show that the long-time power-law tails of temporal correlations, and the corresponding low-frequency behavior, of two-time dipole correlations, neglected in the Markovian limit, dramatically affect the prediction of the force.


title: Fifty Years of Pulsar Candidate Selection: From simple filters to a new principled real-time classification approach

abstract: Improving survey specifications are causing an exponential rise in pulsar candidate numbers and data volumes. We study the candidate filters used to mitigate these problems during the past fifty years. We find that some existing methods such as applying constraints on the total number of candidates collected per observation, may have detrimental effects on the success of pulsar searches. Those methods immune to such effects are found to be ill-equipped to deal with the problems associated with increasing data volumes and candidate numbers, motivating the development of new approaches. We therefore present a new method designed for on-line operation. It selects promising candidates using a purpose-built tree-based machine learning classifier, the Gaussian Hellinger Very Fast Decision Tree (GH-VFDT), and a new set of features for describing candidates. The features have been chosen so as to i) maximise the separation between candidates arising from noise and those of probable astrophysical origin, and ii) be as survey-independent as possible. Using these features our new approach can process millions of candidates in seconds (~1 million every 15 seconds), with high levels of pulsar recall (90%+). This technique is therefore applicable to the large volumes of data expected to be produced by the Square Kilometre Array (SKA). Use of this approach has assisted in the discovery of 20 new pulsars in data obtained during the LOFAR Tied-Array All-Sky Survey (LOTAAS).


title: A sharp counterexample to local existence of low regularity solutions to Einstein's equations in wave coordinates

abstract: We are concerned with how regular initial data have to be to ensure local existence for Einstein's equations in wave coordinates. Klainerman-Rodnianski and Smith-Tataru showed that there in general is local existence for data in Sobolev spaces H^s with regularity s>2. We give an example of data in Sobolev spaces with regularity s=2 for which there is no local solution in this space.


title: Stability and Robustness of RBF Interpolation

abstract: We consider how some methods of uniform and nonuniform interpolation by translates of radial basis functions -- specifically the so-called general multiquadrics -- perform in the presence of certain types of noise. These techniques provide some new avenues for interpolation on bounded domains that is different from the existing literature by using fast Fourier transform methods to approximate cardinal functions associated with the RBF.


title: Stringy origin of 4d black hole microstates

abstract: We derive a precise dictionary between micro-state geometries and open string condensates for a large class of excitations of four dimensional BPS black holes realised in terms of D3-branes intersecting on a six-torus. The complete multipole expansion of the supergravity solutions at weak coupling is extracted from string amplitudes involving one massless closed string and multiple open strings insertions on disks with mixed boundary conditions.


title: Indifference to symmetry in Hrushovski's ab initio construction

abstract: Denote Hrushovski's non-collapsed ab initio construction for an $n$-ary relation by $\mathcal{M}_{\not\sim}$ and the analogous construction for a symmetric $n$-ary relation by $\mathcal{M}_{\sim}$. We show that $\mathcal{M}_{\not\sim}$ is isomorphic to a proper reduct of $\mathcal{M}_{\sim}$ and vice versa, and that the combinatorial pregeometries associated with both structures are isomorphic.


title: Coalescence of Euclidean geodesics on the Poisson-Delaunay triangulation

abstract: Let us consider Euclidean first-passage percolation on the Poisson-Delaunay triangulation. We prove almost sure coalescence of any two semi-infinite geodesics with the same asymptotic direction. The proof is based on an adapted Burton-Keane argument and makes use of the concentration property for shortest-path lengths in the considered graphs. Moreover, by considering the specific example of the relative neighborhood graph, we illustrate that our approach extends to further well-known graphs in computational geometry. As an application, we show that the expected number of semi-infinite geodesics starting at a given vertex and leaving a disk of a certain radius grows at most sublinearly in the radius.


title: Resonance effects in the Raman scattering of mono- and few layers MoSe$_2$

abstract: Using resonant Raman scattering spectroscopy with 25 different laser lines, we describe the Raman scattering spectra of mono- and multi-layers 2H-molybdenum diselenide (MoSe$_2$) as well as the different resonances affecting the most pronounced features. For high-energy phonons, both A- and E- symmetry type phonons present resonances with A and B excitons of MoSe$_2$ together with a marked increase of intensity when exciting at higher energy, close to the C exciton energy. We observe symmetry dependent exciton-phonon coupling affecting mainly the low-energy rigid layer phonon modes. The shear mode for multilayer displays a pronounced resonance with the C exciton while the breathing mode has an intensity that grows with the excitation laser energy, indicating a resonance with electronic excitations at energies higher than that of the C exciton.


title: SUSY partners of the truncated oscillator, Painlev\'e transcendents and B\"acklund transformations

abstract: In this work the supersymmetric technique is applied to the truncated oscillator to generate Hamiltonians ruled by second and third-order polynomial Heisenberg algebras, which are connected to the Painlev\'e IV and Painlev\'e V equations respectively. The aforementioned connection is exploited to produce particular solutions to both non-linear differential equations and the B\"acklund transformations relating them.


title: New parametrization of the deuteron wave function and calculations of the tensor polarization

abstract: The coefficients of new analytical forms for the deuteron wave function in coordinate space for NijmI, NijmII, Nijm93, Reid93 and Argonne v18 potentials have been numerically calculated. The obtained wave functions do not contain any superfluous knots. The designed parameters of the deuteron are in good agreement with the experimental and theoretical data. The new parameterization of the deuteron wave function in momentum space is obtained. The tensor polarization t20 calculated based on the wave functions is proportionate to the earlier published results.


title: Coexistence of topological Dirac fermions in the surface and three-dimensional Dirac cone state in the bulk of ZrTe$_{5}$ single crystal

abstract: Although, the long-standing debate on the resistivity anomaly in ZrTe$_{5}$ somewhat comes to an end, the exact topological nature of the electronic band structure remains elusive till today. Theoretical calculations predicted that bulk ZrTe$_{5}$ to be either a weak or a strong three-dimensional (3D) topological insulator. However, the angle resolved photoemission spectroscopy and transport measurements clearly demonstrate 3D Dirac cone state with a small mass gap between the valence band and conduction band in the bulk. From the magnetization and magneto-transport measurements on ZrTe$_{5}$ single crystal, we have detected both the signature of helical spin texture from topological surface state and chiral anomaly associated with the 3D Dirac cone state in the bulk. This implies that ZrTe$_{5}$ is a novel 3D topological insulator having massless Dirac fermionic excitation in its bulk gap state. Whereas, no 3D topological insulator known in material science holds linear band dispersion in its insulating bulk. Apart from the band topology, it is also apparent from the resistivity and Hall measurements that the anomalous peak in the resistivity can be shifted to a much lower temperature ($T$$<$2 K) by controlling impurity and defects.


title: Laser Pulse Driven THz Generation via Resonant Transition Radiation in Inhomogeneous Plasmas

abstract: An intense, short laser pulse propagating across a plasma boundary ponderomotively drives THz radiation. Full format PIC simulations and theoretical analysis are conducted to investigate the properties of this radiation. Simulation results show the THz emission originates in regions of varying density and covers a broad spectrum with maximum frequency close to the maximum plasma frequency. In the case of a sharp vacuum-plasma boundary, the radiation is generated symmetrically at the plasma entrance and exit, and its properties are independent of plasma density when the density exceeds a characteristic value determined by the product of the plasma frequency and the laser pulse duration. For a diffuse vacuum-plasma boundary, the emission from the plasma entrance and exit is asymmetric: increasing and decreasing density ramps enhance and diminish the radiated energy respectively. Enhancements by factors of 50 are found and simulations show that a 1.66 J, 50 fs driver pulse can generate ~400 \mu J of THz radiation in a 1.2 mm increasing density ramp. We present a model that attributes this effect to a plasma resonance process in the density ramp. The results from the model match those of the simulations for ramp lengths less than 600 \mu m. For longer ramps for which simulations are too time consuming the model shows that the amount of radiation reaches a maximum at a ramp length determined by collisional absorption.


title: Mechanism of the fcc-hcp Phase Transformation in Solid Ar

abstract: We present an atomistic description of the {\it fcc}--to--{\it hcp} transformation mechanism in solid argon (Ar) obtained from transition path sampling molecular dynamics simulation. The phase transition pathways collected during the sampling for an 8000--particle system reveal three transition types according to the lattice deformation and relaxation details. In all three transition types, we see a critical accumulation of defects and uniform growth of a less ordered transition state, followed by a homogeneous growth of an ordered phase. Stacking disorder is discussed to describe the transition process and the cooperative motions of atoms in \{111\} planes. We investigate the nucleation with larger system. In a system of 18000--particles, the collective movements of atoms required for this transition are facilitated by the formation and growth of stacking faults. However the enthalpy barrier is still far beyond the thermal fluctuation. The high barrier explains previous experimental observations of the inaccessibility of the bulk transition at low pressure and its sluggishness even at extremely high pressure. The transition mechanism in bulk Ar is different from Ar nanoclusters as the orthorhombic intermediate structure proposed for the latter is not observed in any of our simulations.


title: Unramified Euler sums and Hoffman $\star$ basis

abstract: When looking at how periods of $\pi_{1}^{\mathfrak{m}}(\mathbb{P}^{1}\diagdown \lbrace 0, 1, \infty \rbrace )$, i.e. multiple zeta values, embeds into periods of $\pi_{1}^{\mathfrak{m}}(\mathbb{P}^{1}\diagdown \lbrace 0, \pm 1, \infty \rbrace )$, i.e. Euler sums, an explicit criteria via the coaction $\Delta$ acting on their motivic versions comes out. In this paper, adopting this Galois descent approach, we present a new basis for the space $\mathcal{H}^{1}$ of motivic multiple zeta values via motivic Euler sums. Up to an analytic conjecture, we also prove that the motivic Hoffman star basis $\zeta^{\star, \mathfrak{m}} (2^{a_{1}},3,\cdots,3, 2^{a_{p}}, 3, 2^{b})$ is a basis of $\mathcal{H}^{1}$. Under a general motivic identity that we conjecture, these bases are identical. Other examples of unramified ES with alternating patterns of even and odds are also provided.


title: Gaps, rings, and non-axisymmetric structures in protoplanetary disks - Emission from large grains

abstract: Dust grains with sizes around (sub)mm are expected to couple only weakly to the gas motion in regions beyond 10 au of circumstellar disks. In this work, we investigate the influence of the spatial distribution of such grains on the (sub)mm appearance of magnetized protoplanetary disks. We perform non-ideal global 3D magneto-hydrodynamic stratified disk simulations including particles of different sizes (50 micron to 1 cm), using a Lagrangian particle solver. We calculate the spatial dust temperature distribution, including the dynamically coupled submicron-sized dust grains, and derive ideal continuum re-emission maps of the disk through radiative transfer simulations. Finally, we investigate the feasibility to observe specific structures in the thermal re-emission maps with the Atacama Large Millimeter/submillimeter Array (ALMA). The pressure bump close to the outer edge of the dead-zone leads to particle trapping in ring structures. More specifically, vortices in the disk concentrate the dust and create an inhomogeneous distribution of the solid material in the azimuthal direction. The large-scale disk perturbations are preserved in the (sub)mm re-emission maps. The observable structures are very similar to those expected to result from planet-disk interaction. The larger dust particles increase the brightness contrast between gap and ring structures. We find that rings, gaps and the dust accumulation in the vortex could be traced with ALMA down to a scale of a few astronomical units in circumstellar disks located in nearby star-forming regions. Finally, we present a brief comparison of these structures with those recently found with ALMA in the young circumstellar disks of HL Tau and Oph IRS 48.


title: Linear embedding of free energy minimization

abstract: Exact free energy minimization is a convex optimization problem that is usually approximated with stochastic sampling methods. Deterministic approximations have been less successful because many desirable properties have been difficult to attain. Such properties include the preservation of convexity, lower bounds on free energy, and applicability to systems without subsystem structure. We satisfy all of these properties by embedding free energy minimization into a linear program over energy-resolved expectation values. Numerical results on small systems are encouraging, but a lack of size consistency necessitates further development for large systems.


title: Strength of weak layers in cascading failures on multiplex networks: case of the international trade network

abstract: Many real-world complex systems across natural, social, and economical domains consist of manifold layers to form multiplex networks. The multiple network layers give rise to nonlinear effect for the emergent dynamics of systems. Especially, weak layers that can potentially play significant role in amplifying the vulnerability of multiplex networks might be shadowed in the aggregated single-layer network framework which indiscriminately accumulates all layers. Here we present a simple model of cascading failure on multiplex networks of weight-heterogeneous layers. By simulating the model on the multiplex network of international trades, we found that the multiplex model produces more catastrophic cascading failures which are the result of emergent collective effect of coupling layers, rather than the simple sum thereof. Therefore risks can be systematically underestimated in single-layer network analyses because the impact of weak layers can be overlooked. We anticipate that our simple theoretical study can contribute to further investigation and design of optimal risk-averse real-world complex systems.


title: Fractal symmetries: Ungauging the cubic code

abstract: Gauging is a ubiquitous tool in many-body physics. It allows one to construct highly entangled topological phases of matter from relatively simple phases and to relate certain characteristics of the two. Here we develop a gauging procedure for general submanifold symmetries of Pauli Hamiltonians, including symmetries of fractal type. We show a relation between the pre- and post-gauging models and use this to construct short-range entangled phases with fractal-like symmetries, one of which is mapped to the cubic code by the gauging.


title: On the effect of randomness on planted 3-coloring models

abstract: We present the hosted coloring framework for studying algorithmic and hardness results for the $k$-coloring problem. There is a class ${\cal H}$ of host graphs. One selects a graph $H\in{\cal H}$ and plants in it a balanced $k$-coloring (by partitioning the vertex set into $k$ roughly equal parts, and removing all edges within each part). The resulting graph $G$ is given as input to a polynomial time algorithm that needs to $k$-color $G$ (any legal $k$-coloring would do -- the algorithm is not required to recover the planted $k$-coloring). Earlier planted models correspond to the case that ${\cal H}$ is the class of all $n$-vertex $d$-regular graphs, a member $H\in{\cal H}$ is chosen at random, and then a balanced $k$-coloring is planted at random. Blum and Spencer [1995] designed algorithms for this model when $d=n^{\delta}$ (for $0<\delta\le1$), and Alon and Kahale [1997] managed to do so even when $d$ is a sufficiently large constant. The new aspect in our framework is that it need not involve randomness. In one model within the framework (with $k=3$) $H$ is a $d$ regular spectral expander (meaning that except for the largest eigenvalue of its adjacency matrix, every other eigenvalue has absolute value much smaller than $d$) chosen by an adversary, and the planted 3-coloring is random. We show that the 3-coloring algorithm of Alon and Kahale [1997] can be modified to apply to this case. In another model $H$ is a random $d$-regular graph but the planted balanced $3$-coloring is chosen by an adversary, after seeing $H$. We show that for a certain range of average degrees somewhat below $\sqrt{n}$, finding a 3-coloring is NP-hard. Together these results (and other results that we have) help clarify which aspects of randomness in the planted coloring model are the key to successful 3-coloring algorithms.


title: Redshift-space distortions around voids

abstract: We have derived estimators for the linear growth rate of density fluctuations using the cross-correlation function of voids and haloes in redshift space, both directly and in Fourier form. In linear theory, this cross-correlation contains only monopole and quadrupole terms. At scales greater than the void radius, linear theory is a good match to voids traced out by haloes in N-body simulations; small-scale random velocities are unimportant at these radii, only tending to cause small and often negligible elongation of the redshift-space cross-correlation function near its origin. By extracting the monopole and quadrupole from the cross-correlation function, we measure the linear growth rate without prior knowledge of the void profile or velocity dispersion. We recover the linear growth parameter $\beta$ to 9% precision from an effective volume of 3(Gpc/h)^3 using voids with radius greater than 25Mpc/h. Smaller voids are predominantly sub-voids, which may be more sensitive to the random velocity dispersion; they introduce noise and do not help to improve the measurement. Adding velocity dispersion as a free parameter allows us to use information at radii as small as half of the void radius. The precision on $\beta$ is reduced to approximately 5%. Contrary to the simple redshift-space distortion pattern in overdensities, voids show diverse shapes in redshift space, and can appear either elongated or flattened along the line of sight. This can be explained by the competing amplitudes of the local density contrast, plus the radial velocity profile and its gradient, with the latter two factors being determined by the cumulative density profile of voids. The distortion pattern is therefore determined solely by the void profile and is different for void-in-cloud and void-in-void. This diversity of redshift-space void morphology complicates measurements of the Alcock-Paczynski effect using voids.


title: Negativity in the Generalized Valence Bond Solid State

abstract: Using a graphical presentation of the spin $S$ one dimensional Valence Bond Solid (VBS) state, based on the representation theory of the $SU(2)$ Lie-algebra of spins, we compute the spectrum of a mixed state reduced density matrix. This mixed state of two blocks of spins $A$ and $B$ is obtained by tracing out the spins outside $A$ and $B$, in the pure VBS state density matrix. We find in particular that the negativity of the mixed state is non-zero only for adjacent subsystems. The method introduced here can be generalized to the computation of entanglement properties in Levin-Wen models, that possess a similar algebraic structure to the VBS state in the groundstate.

categories: quant-ph cond-mat.str-el math-ph math.MP


title: Acoustic Scattering from Corners, Edges and Circular Cones

abstract: Consider the time-harmonic acoustic scattering from a bounded penetrable obstacle imbedded in an isotropic homogeneous medium. The obstacle is supposed to possess a circular conic point or an edge point on the boundary in three dimensions and a planar corner point in two dimensions. The opening angles of cones and edges are allowed to be any number in $(0,2\pi)\backslash\{\pi\}$. We prove that such an obstacle scatters any incoming wave non-trivially (i.e., the far field patterns cannot vanish identically), leading to the absence of real non-scattering wavenumbers. Local and global uniqueness results for the inverse problem of recovering the shape of a penetrable scatterers are also obtained using a single incoming wave. Our approach relies on the singularity analysis of the inhomogeneous Laplace equation in a cone.


title: Clumpy wind accretion in supergiant neutron star high mass X-ray binaries

abstract: The accretion of the stellar wind material by a compact object represents the main mechanism powering the X-ray emission in classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. In this work we present the first attempt to simulate the accretion process of a fast and dense massive star wind onto a neutron star, taking into account the effects of the centrifugal and magnetic inhibition of accretion ("gating") due to the spin and magnetic field of the compact object. We made use of a radiative hydrodynamical code to model the non-stationary radiatively driven wind of an O-B supergiant star and then place a neutron star characterized by a fixed magnetic field and spin period at a certain distance from the massive companion. Our calculations follow, as a function of time (on a total time scale of several hours), the transition of the system through all different accretion regimes that are triggered by the intrinsic variations in the density and velocity of the non-stationary wind. The X-ray luminosity released by the system is computed at each time step by taking into account the relevant physical processes occurring in the different accretion regimes. Synthetic lightcurves are derived and qualitatively compared with those observed from classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. Although a number of simplifications are assumed in these calculations, we show that taking into account the effects of the centrifugal and magnetic inhibition of accretion significantly reduces the average X-ray luminosity expected for any neutron star wind-fed binary. The present model calculations suggest that long spin periods and stronger magnetic fields are favoured in order to reproduce the peculiar behavior of supergiant fast X-ray transients in the X-ray domain.


title: Computational Analysis of Sparsity-Exploiting Moment Relaxations of the OPF Problem

abstract: With the potential to find global solutions, significant research interest has focused on convex relaxations of the non-convex OPF problem. Recently, "moment-based" relaxations from the Lasserre hierarchy for polynomial optimization have been shown capable of globally solving a broad class of OPF problems. Global solution of many large-scale test cases is accomplished by exploiting sparsity and selectively applying the computationally intensive higher-order relaxation constraints. Previous work describes an iterative algorithm that indicates the buses for which the higher-order constraints should be enforced. In order to speed computation of the moment relaxations, this paper provides a study of the key parameter in this algorithm as applied to relaxations from both the original Lasserre hierarchy and a recent complex extension of the Lasserre hierarchy.


title: Applying Artifical Neural Networks To Predict Nominal Vehicle Performance

abstract: This paper investigates the use of artificial neural networks (ANNs) to replace traditional algorithms and manual review for identifying anomalies in vehicle run data. The specific data used for this study is from undersea vehicle qualification tests. Such data is highly non-linear, therefore traditional algorithms are not adequate and manual review is time consuming. By using ANNs to predict nominal vehicle performance based solely on information available pre-run, vehicle deviation from expected performance can be automatically identified in the post-run data. Such capability is only now becoming available due to the rapid increase in understanding of ANN framework and available computing power in the past decade. The ANN trained for the purpose of this investigation is relatively simple, to keep the computing requirements within the parameters of a modern desktop PC. This ANN showed potential in predicting vehicle performance, particularly during transient events within the run data. However, there were also several performance cases, such as steady state operation and cases which did not have sufficient training data, where the ANN showed deficiencies. It is expected that as computational power becomes more readily available, ANN understanding matures, and more training data is acquired from real world tests, the performance predictions of the ANN will surpass traditional algorithms and manual human review.


title: gadfly: A pandas-based Framework for Analyzing GADGET Simulation Data

abstract: We present the first public release (v0.1) of the open-source GADGET Dataframe Library: gadfly. The aim of this package is to leverage the capabilities of the broader python scientific computing ecosystem by providing tools for analyzing simulation data from the astrophysical simulation codes GADGET and GIZMO using pandas, a thoroughly documented, open-source library providing high-performance, easy-to-use data structures that is quickly becoming the standard for data analysis in python. Gadfly is a framework for analyzing particle-based simulation data stored in the HDF5 format using pandas DataFrames. The package enables efficient memory management, includes utilities for unit handling, coordinate transformations, and parallel batch processing, and provides highly optimized routines for visualizing smoothed-particle hydrodynamics (SPH) datasets.


title: Distributed Inexact Damped Newton Method: Data Partitioning and Load-Balancing

abstract: In this paper we study inexact dumped Newton method implemented in a distributed environment. We start with an original DiSCO algorithm [Communication-Efficient Distributed Optimization of Self-Concordant Empirical Loss, Yuchen Zhang and Lin Xiao, 2015]. We will show that this algorithm may not scale well and propose an algorithmic modifications which will lead to less communications, better load-balancing and more efficient computation. We perform numerical experiments with an regularized empirical loss minimization instance described by a 273GB dataset.


title: First-principles study of the magnetic ground state in kagome francisites Cu3Bi(SeO3)2O2X (X=Cl, Br)

abstract: We explore magnetic behavior of kagome francisites Cu3Bi(SeO3)2O2X (X = Cl and Br) using first-principles calculations. To this end, we propose an approach based on the Hubbard model in the Wannier functions basis constructed on the level of local-density approximation (LDA). The ground-state spin configuration is determined by a Hartree-Fock solution of the Hubbard model both in zero magnetic field and in applied magnetic fields. Additionally, parameters of an effective spin Hamiltonian are obtained by taking into account the hybridization effects and spin-orbit coupling. We show that only the former approach, the Hartree-Fock solution of the Hubbard model, allows for a complete description of the anisotropic magnetization process. While our calculations confirm that the canted zero-field ground state arises from a competition between ferromagnetic nearest-neighbor and antiferromagnetic next-nearest-neighbor couplings in the kagome planes, weaker anisotropic terms are crucial for fixing spin directions and for the overall magnetization process. We thus show that the Hartree-Fock solution of an electronic Hamiltonian is a viable alternative to the analysis of effective spin Hamiltonians when a magnetic ground state and effects of external field are considered.


title: The effects of three-dimensional defects on one-way surface plasmon propagation for photonic topological insulators comprised of continuous media

abstract: We have investigated the one-way surface plasmon-polariton (SPP) at the interface of a continuous magneto-plasma material and metal. We demonstrated that TM modes inside a continuous material can be assigned non-trivial Chern numbers analogous to those of topological photonic crystals; moreover these Chern numbers can be calculated analytically. This leads to the appearance of topologically protected surface modes propagating at frequencies inside the bandgap of the magneto-plasma. Previous works considered 2D structures; here we consider the effects of 3D defects, and show that, although backward propagation/reflection cannot occur, side scattering does take place and it has significant effect on the propagation of the surface mode. Several different waveguiding geometries are considered for reducing the effects of side-scattering, and we also consider the effects of metal loss.


title: Shocks in nova outflows. II. Synchrotron radio emission

abstract: The discovery of GeV gamma-rays from classical novae indicates that shocks and relativistic particle acceleration are energetically key in these events. Further evidence for shocks comes from thermal keV X-ray emission and an early peak in the radio light curve on a timescale of months with a brightness temperature which is too high to result from freely expanding photo-ionized gas. Paper I developed a one dimensional model for the thermal emission from nova shocks. This work concluded that the shock-powered radio peak cannot be thermal if line cooling operates in the post-shock gas at the rate determined by collisional ionization equilibrium. Here we extend this calculation to include non-thermal synchrotron emission. Applying our model to three classical novae, we constrain the amplification of the magnetic field $\epsilon_B$ and the efficiency $\epsilon_e$ of accelerating relativistic electrons of characteristic Lorentz factor $\gamma \sim 100$. If the shocks are radiative (low velocity $v_{\rm sh} \lesssim 1000$ km s$^{-1}$) and cover a large solid angle of the nova outflow, as likely characterize those producing gamma-rays, then values of $\epsilon_e \sim 0.01-0.1$ are required to achieve the peak radio brightness for $\epsilon_B = 10^{-2}$. Such high efficiencies exclude secondary pairs from pion decay as the source of the radio-emitting particles, instead favoring the direct acceleration of electrons at the shock. If the radio-emitting shocks are instead adiabatic (high velocity), as likely characterize those responsible for the thermal X-rays, then much higher brightness temperatures are possible, allowing the radio-emitting shocks to cover a smaller outflow solid angle.


title: Setting the Stage for a Non-Supersymmetric UV-Complete String Phenomenology

abstract: In this talk, I discuss our recent work concerning the construction of non-supersymmetric heterotic string models which have exponentially suppressed dilaton tadpoles and cosmological constants, and thus greatly enhanced stability properties. The existence of such models opens the door to non-supersymmetric string model-building, and I discuss how semi-realistic string models resembling the Standard Model or any of its unified variants may be constructed within this framework. These models maintain modular invariance and exhibit a misaligned supersymmetry which ensures UV finiteness, even without spacetime supersymmetry. I also discuss the potential implications for phenomenology.


title: (p,q)-Five Brane and (p,q)-String Solutions, their Bound State and its Near Horizon Limit

abstract: We determine (p,q)-string and (p,q)-five brane solutions of type IIB supergravity using SL(2,Z)-symmetry of the full type IIB superstring theory. We also determine SL(2,Z)-transformed solution corresponding to the bound state of NS5-branes and fundamental strings. Then we analyze its near horizon limit and we show that it leads to the AdS(3)xS(3) with mixed fluxes.


title: Embedding by Normalisation

abstract: This paper presents the insight that practical embedding techniques, commonly used for implementing Domain-Specific Languages, correspond to theoretical Normalisation-By-Evaluation (NBE) techniques, commonly used for deriving canonical form of terms with respect to an equational theory. NBE constitutes of four components: a syntactic domain, a semantic domain, and a pair of translations between the two. Embedding also often constitutes of four components: an object language, a host language, encoding of object terms in the host, and extraction of object code from the host. The correspondence is deep in that all four components in embedding and NBE correspond to each other. Based on this correspondence, this paper introduces Embedding-By-Normalisation (EBN) as a principled approach to study and structure embedding. The correspondence is useful in that solutions from NBE can be borrowed to solve problems in embedding. In particular, based on NBE techniques, such as Type-Directed Partial Evaluation, this paper presents a solution to the problem of extracting object code from embedded programs involving sum types, such as conditional expressions, and primitives, such as literals and operations on them.


title: Universal relations between nongaussian fluctuations in heavy-ion collisions

abstract: We show that universality near a critical end point implies a characteristic relation between third- and fourth-order baryon susceptibilities $\chi_3$ and $\chi_4$, resulting in a banana-shaped loop when $\chi_4$ is plotted as a function of $\chi_3$ along a freeze-out line. Including the individual enhancements of $\chi_3$ and $\chi_4$ near a critical point, these features may be a consistent set of observations supporting the interpretation of baryon fluctuations data as arising from criticality.


title: Triangle-hinge models for unoriented membranes

abstract: Triangle-hinge models [arXiv:1503.08812] are introduced to describe worldvolume dynamics of membranes. The Feynman diagrams consist of triangles glued together along hinges and can be restricted to tetrahedral decompositions in a large N limit. In this paper, after clarifying that all the tetrahedra resulting in the original models are orientable, we define a version of triangle-hinge models that can describe the dynamics of unoriented membranes. By regarding each triangle as representing a propagation of an open membrane of disk topology, we introduce a local worldvolume parity transformation which inverts the orientation of triangle, and define unoriented triangle-hinge models by gauging the transformation. Unlike two-dimensional cases, this local transformation generally relates a manifold to a nonmanifold, but still is a well-defined manipulation among tetrahedral decompositions. We further show that matter fields can be introduced in the same way as in the original oriented models. In particular, the models will describe unoriented membranes in a target spacetime by taking matter fields to be the target space coordinates.


title: Multi-Vehicle Collision Avoidance via Hamilton-Jacobi Reachability and Mixed Integer Programming

abstract: Multi-agent differential games are important and useful tools for analyzing many practical problems. With the recent surge of interest in using UAVs for civil purposes, the importance and urgency of developing tractable multi-agent analysis techniques that provide safety and performance guarantees is at an all-time high. Hamilton-Jacobi (HJ) reachability has successfully provided safety guarantees to small-scale systems and is flexible in terms of system dynamics. However, the exponential complexity scaling of HJ reachability prevents its direct application to large scale problems when the number of vehicles is greater than two. In this paper, we overcome the scalability limitations of HJ reachability by using a mixed integer program that exploits the properties of HJ solutions to provide higher-level control logic. Our proposed method provides safety guarantee for three-vehicle systems -- a previously intractable task for HJ reachability -- without incurring significant additional computation cost. Furthermore, our method is scalable beyond three vehicles and performs significantly better by several metrics than an extension of pairwise collision avoidance to multi-vehicle collision avoidance. We demonstrate our proposed method in simulations.


title: Understanding and Improving Convolutional Neural Networks via Concatenated Rectified Linear Units

abstract: Recently, convolutional neural networks (CNNs) have been used as a powerful tool to solve many problems of machine learning and computer vision. In this paper, we aim to provide insight on the property of convolutional neural networks, as well as a generic method to improve the performance of many CNN architectures. Specifically, we first examine existing CNN models and observe an intriguing property that the filters in the lower layers form pairs (i.e., filters with opposite phase). Inspired by our observation, we propose a novel, simple yet effective activation scheme called concatenated ReLU (CRelu) and theoretically analyze its reconstruction property in CNNs. We integrate CRelu into several state-of-the-art CNN architectures and demonstrate improvement in their recognition performance on CIFAR-10/100 and ImageNet datasets with fewer trainable parameters. Our results suggest that better understanding of the properties of CNNs can lead to significant performance improvement with a simple modification.


title: Packing, coding, and ground states

abstract: These are the lecture notes from my 2014 PCMI graduate summer school lectures. In these lectures, we'll study simple models of materials from several different perspectives: geometry (packing problems), information theory (error-correcting codes), and physics (ground states of interacting particle systems). These perspectives each shed light on some of the same problems and phenomena, while highlighting different techniques and connections. One noteworthy phenomenon is the exceptional symmetry that is found in certain special cases, and we'll examine when and why it occurs. The overall theme of the lectures is thus order vs. disorder. How much symmetry can we expect to see in optimal geometric structures?


title: Convergence of loop-erased random walk in the natural parametrization

abstract: Loop-erased random walk, abbreviated LERW, is one of the most well-studied critical lattice models. It is the self-avoiding random walk one gets after erasing the loops from a simple random walk in order or alternatively by considering the branches in a uniformly chosen spanning tree. This paper proves that planar LERW parametrized by renormalized length converges in the lattice size scaling limit to SLE(2) parametrized by 5/4-dimensional Minkowski content. In doing this we also provide a method for proving similar convergence results for other models converging to SLE. Besides the main theorem, several of our results about LERW are of independent interest: for example, two-point estimates, estimates on maximal content, and a "separation lemma".


title: A Probabilistic Angle on One Loop Scalar Integrals

abstract: Recasting the $N$-point one loop scalar integral as a probabilistic problem, allows the derivation of integral recurrence relations as well as exact analytical expressions in the most common cases. $\epsilon$ expansions are derived by writing a formula that relates an $N$-point function in decimal dimension to an $N$-point function in integer dimension. As an example, we give relations for the massive 5-point function in dimension $n=4-2\epsilon$, $n=6-2\epsilon$. The reduction of tensor integrals of rank 2 with $N=5$ is achieved showing the method's potential. Hypergeometric functions are not needed but only integration of arcsine function whose analytical continuation is well known.


title: Fast Reachable Set Approximations via State Decoupling Disturbances

abstract: With the recent surge of interest in using robotics and automation for civil purposes, providing safety and performance guarantees has become extremely important. In the past, differential games have been successfully used for the analysis of safety-critical systems. In particular, the Hamilton-Jacobi (HJ) formulation of differential games provides a flexible way to compute the reachable set, which can characterize the set of states which lead to either desirable or undesirable configurations, depending on the application. While HJ reachability is applicable to many small practical systems, the curse of dimensionality prevents the direct application of HJ reachability to many larger systems. To address computation complexity issues, various efficient computation methods in the literature have been developed for approximating or exactly computing the solution to HJ partial differential equations, but only when the system dynamics are of specific forms. In this paper, we propose a flexible method to trade off optimality with computation complexity in HJ reachability analysis. To achieve this, we propose to simplify system dynamics by treating state variables as disturbances. We prove that the resulting approximation is conservative in the desired direction, and demonstrate our method using a four-dimensional plane model.


title: Fractal scattering of Gaussian solitons in directional couplers with logarithmic nonlinearities

abstract: In this paper we study the interaction of Gaussian solitons in a dispersive and nonlinear media with log-law nonlinearity. The model is described by the coupled logarithmic nonlinear Schr\"odinger equations, which is a nonintegrable system that allows the observation of a very rich scenario in the collision patterns. By employing a variational approach and direct numerical simulations, we observe a fractal-scattering phenomenon from the exit velocities of each soliton as a function of the input velocities. Furthermore, we introduce a linearization model to identify the position of the reflection/transmission window that emerges within the chaotic region. This enable us the possibility of controlling the scattering of solitons as well as the lifetime of bound states.


title: Spin-valley locking in the normal state of a transition-metal dichalocogenide superconductor

abstract: The metallic transition-metal dichalcogenides (TMDCs) are benchmark systems for studying and controlling intertwined electronic orders in solids, with superconductivity developing upon cooling from a charge density wave state. The interplay between such phases is thought to play a critical role in the unconventional superconductivity of cuprates, Fe-based, and heavy-fermion systems, yet even for the more moderately-correlated TMDCs, their nature and origins have proved highly controversial. Here, we study a prototypical example, $2H$-NbSe$_2$, by spin- and angle-resolved photoemission and first-principles theory. We find that the normal state, from which its hallmark collective phases emerge, is characterised by quasiparticles whose spin is locked to their valley pseudospin. This results from a combination of strong spin-orbit interactions and local inversion symmetry breaking. Non-negligible interlayer coupling further drives a rich three-dimensional momentum-dependence of the underlying Fermi surface spin texture. Together, these findings necessitate a fundamental re-investigation of the nature of charge order and superconducting pairing in NbSe$_2$ and related TMDCs.

categories: cond-mat.supr-con cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el


title: Safe Sequential Path Planning Under Disturbances and Imperfect Information

abstract: Multi-UAV systems are safety-critical, and guarantees must be made to ensure no unsafe configurations occur. Hamilton-Jacobi (HJ) reachability is ideal for analyzing such safety-critical systems; however, its direct application is limited to small-scale systems of no more than two vehicles due to an exponentially-scaling computational complexity. Previously, the sequential path planning (SPP) method, which assigns strict priorities to vehicles, was proposed; SPP allows multi-vehicle path planning to be done with a linearly-scaling computational complexity. However, the previous formulation assumed that there are no disturbances, and that every vehicle has perfect knowledge of higher-priority vehicles' positions. In this paper, we make SPP more practical by providing three different methods to account for disturbances in dynamics and imperfect knowledge of higher-priority vehicles' states. Each method has different assumptions about information sharing. We demonstrate our proposed methods in simulations.


title: Causality, initial conditions, and inflationary magnetogenesis

abstract: The post-inflationary evolution of inflation-produced magnetic fields, conventional or not, can change dramatically when two fundamental issues are accounted for. The first is causality, which demands that local physical processes can never affect superhorizon perturbations. The second is the nature of the transition from inflation to reheating and then to the radiation era, which determine the initial conditions at the start of these epochs. Causality implies that inflationary magnetic fields dot not freeze into the matter until they have re-entered the causal horizon. The nature of the cosmological transitions and the associated initial conditions, on the other hand, determine the large-scale magnetic evolution after inflation. Put together, the two can slow down the adiabatic magnetic decay on superhorizon scales throughout the universe's post-inflationary evolution and thus lead to considerably stronger residual magnetic fields. This is "good news" for both the conventional and the non-conventional scenarios of cosmic magnetogenesis. Mechanisms operating outside standard electromagnetism, in particular, do not need to enhance their fields too much during inflation, in order to produce seeds that can feed the galactic dynamo today. In fact, even conventionally produced inflationary magnetic fields might be able to sustain the dynamo.


title: New longitudinal mode and compression of pair ions in plasma

abstract: Positive and negative ions forming so-called pair plasma differing in sign of their charge but asymmetric in mass and temperature support a new acoustic mode where damping of heavier ion dominates. The condition for the excitation of ion sound wave through electron beam induced Cherenkov instability is also investigated. This beam can generate a perturbation in the pair ion plasmas in the presence of electrons when there is number density, temperature and mass difference in the two species of ions. Basic emphasis is on the focusing of ion sound waves and we show how, in the area of localization of wave energy, the density of pair particles increases while electrons are pushed away from that region. Further, this localization of wave is dependent on the shape of the pulse. Considering the example of pancake and bullet shaped pulses, we find that only the former leads to compression of pair ions in the supersonic regime of the focusing region. Here possible existence of regions where pure pair particles can exist may also be speculated which is not only useful from academic point of view but also to mimic the situation of plasma (electro positron asymmetric and symmetric) observed in astrophysical environment.


title: Comparison of adaptive multiresolution and adaptive mesh refinement applied to simulations of the compressible Euler equations

abstract: We present a detailed comparison between two adaptive numerical approaches to solve partial differential equations (PDEs), adaptive multiresolution (MR) and adaptive mesh refinement (AMR). Both discretizations are based on finite volumes in space with second order shock-capturing, and explicit time integration either with or without local time-stepping. The two methods are benchmarked for the compressible Euler equations in Cartesian geometry. As test cases a 2D Riemann problem, Lax-Liu 6, and a 3D ellipsoidally expanding shock wave have been chosen. We compare and assess their computational efficiency in terms of CPU time and memory requirements. We evaluate the accuracy by comparing the results of the adaptive computations with those obtained with the corresponding FV scheme using a regular fine mesh. We find that both approaches yield similar trends for CPU time compression for increasing number of refinement levels. MR exhibits more efficient memory compression than AMR and shows slightly enhanced convergence; however, a larger absolute overhead is measured for the tested codes.


title: Analysis of the Cuban journal Bibliotecas: Anales de Investigacion

abstract: The objective of this article is to describe the academic impact, the editorial process quality, and the editorial and visibility strategies of Bibliotecas. Anales de Investigacion (BAI), a scientific Cuban journal edited by National Library of Cuba Jose Marti. The academic impact is determined through a citation analysis, which considers Google Scholar database as reference source. The bibliometric indicators applied are: citation per year, citation vs. self-citation, citable journals vs. non-citable documents, Hirsch Index, and impact factor. The editorial process quality and the visibility strategies are determined through a self-evaluation which takes into account the SciELO, Scopus, CLASE, Redalyc, Latindex, Dialnet, and ERIH PLUS methodologies. The results reveal an ascending citation line that highlights citing journals from the field of Library and Information Science, Medicine and Health Sciences, and Education. Aspects related content and format have negatively influenced on editorial process quality. Some strategies are proposed to improve scientific visibility through the inclusion in databases, directories, and social and academic networks. In general, this study contributes to the editorial decision taking, an issue that could augment the impact and scientific visibility of BAI.


title: Jet effects in high-multiplicity pp events

abstract: The study of the high-multiplicity pp events has become important because we need to understand the origin of the fluid-like features which have been found in such small systems. In this work we concentrate on the radial flow signatures. To this end, the role of jets in high-multiplicity pp collisions is investigated using PYTHIA 8.


title: Guard Your Daggers and Traces: Properties of Guarded (Co-)recursion

abstract: Motivated by the recent interest in models of guarded (co-)recursion, we study their equational properties. We formulate axioms for guarded fixpoint operators generalizing the axioms of iteration theories of Bloom and \'Esik. Models of these axioms include both standard (e.g., cpo-based) models of iteration theories and models of guarded recursion such as complete metric spaces or the topos of trees studied by Birkedal et al. We show that the standard result on the satisfaction of all Conway axioms by a unique dagger operation generalizes to the guarded setting. We also introduce the notion of guarded trace operator on a category, and we prove that guarded trace and guarded fixpoint operators are in one-to-one correspondence. Our results are intended as first steps leading, hopefully, towards future description of classifying theories for guarded recursion.


title: Phase Retrieval from 1D Fourier Measurements: Convexity, Uniqueness, and Algorithms

abstract: This paper considers phase retrieval from the magnitude of 1D over-sampled Fourier measurements, a classical problem that has challenged researchers in various fields of science and engineering. We show that an optimal vector in a least-squares sense can be found by solving a convex problem, thus establishing a hidden convexity in Fourier phase retrieval. We also show that the standard semidefinite relaxation approach yields the optimal cost function value (albeit not necessarily an optimal solution) in this case. A method is then derived to retrieve an optimal minimum phase solution in polynomial time. Using these results, a new measuring technique is proposed which guarantees uniqueness of the solution, along with an efficient algorithm that can solve large-scale Fourier phase retrieval problems with uniqueness and optimality guarantees.

categories: math.OC cs.IR cs.IT math.IT math.ST stat.AP stat.TH


title: Cooperative coupling of static magnetism and bulk superconductivity in the stripe phase of La_2-xBa_xCuO_4: Pressure- and doping-dependent studies

abstract: Static spin-stripe order and superconductivity were systematically studied in La_2-xBa_xCuO_4 (0.11 - 0.17) at ambient pressure by means of magnetization and muSR experiments. We find that all the investigated La_2-xBa_xCuO_4 samples exhibit static spin-stripe order and that the quasi two-dimensional superconducting (SC) transition temperature T_c1 and the static spin-stripe order temperature T_so have very similar values throughout the phase diagram. Moreover, the magnetic and the SC properties of the x = 0.155 (LBCO-0.155) and x = 0.17 (LBCO-0.17) samples were studied under hydrostatic pressure. As a remarkable result, in these bulk cuprate superconductors the three-dimensional SC transition temperature T_c and T_so nearly coincide at all pressure investigated (0 - 2.3 GPa). We also observed a pressure induced transition from long-range spin stripe order to a disordered magnetic state at 1.6 GPa in LBCO-0.155, coexisting with a SC state with substantial superfluid density. In LBCO-0.17 a disordered magnetic state is present at all p. The present results indicate that static magnetic order and SC pairing correlations develop in a cooperative fashion in La_2-xBa_xCuO_4, and provide a new route of understanding the complex interplay between static magnetism and superconductivity in the stripe phase of cuprates.


title: First Upper Limits on the Radar Cross Section of Cosmic-Ray Induced Extensive Air Showers

abstract: TARA (Telescope Array Radar) is a cosmic ray radar detection experiment colocated with Telescope Array, the conventional surface scintillation detector (SD) and fluorescence telescope detector (FD) near Delta, Utah, U.S.A. The TARA detector combines a 40 kW, 54.1 MHz VHF transmitter and high-gain transmitting antenna which broadcasts the radar carrier over the SD array and within the FD field of view, towards a 250 MS/s DAQ receiver. TARA has been collecting data since 2013 with the primary goal of observing the radar signatures of extensive air showers (EAS). Simulations indicate that echoes are expected to be short in duration (~10 microseconds) and exhibit rapidly changing frequency, with rates on the order of 1 MHz/microsecond. The EAS radar cross-section (RCS) is currently unknown although it is the subject of over 70 years of speculation. A novel signal search technique is described in which the expected radar echo of a particular air shower is used as a matched filter template and compared to waveforms obtained by triggering the radar DAQ using the Telescope Array fluorescence detector. No evidence for the scattering of radio frequency radiation by EAS is obtained to date. We report the first quantitative RCS upper limits using EAS that triggered the Telescope Array Fluorescence Detector.


title: Characteristics of Four Upward-pointing Cosmic-ray-like Events Observed with ANITA

abstract: We report on four radio-detected cosmic-ray (CR) or CR-like events observed with the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload. Two of the four were previously identified as stratospheric CR air showers during the ANITA-I flight. A third stratospheric CR was detected during the ANITA-II flight. Here we report on characteristics these three unusual CR events, which develop nearly horizontally, 20-30~km above the surface of the Earth. In addition, we report on a fourth steeply upward-pointing ANITA-I CR-like radio event which has characteristics consistent with a primary that emerged from the surface of the ice. This suggests a possible $\tau$-lepton decay as the origin of this event, but such an interpretation would require significant suppression of the Standard Model ${\tau}$-neutrino cross section.


title: Highly anisotropic and two-fold symmetric superconducting gap in nematically ordered FeSe$_{0.93}$S$_{0.07}$

abstract: FeSe exhibits a novel ground state in which superconductivity coexists with a nematic order in the absence of any long-range magnetic order. Here we report an angle-resolved photoemission study on the superconducting gap structure in the nematic state of FeSe$_{0.93}$S$_{0.07}$, without the complication caused by Fermi surface reconstruction induced by magnetic order. We found that the superconducting gap shows a pronounced 2-fold anisotropy around the elliptical hole pocket near the Z point of the Brillouin zone, with gap minima at the endpoints of its major axis, while no detectable gap was observed around the zone center and zone corner. The large anisotropy and nodal gap distribution demonstrate the substantial effects of the nematicity on the superconductivity, and thus put strong constraints on the current theories.


title: Moment dynamics for a class of time-triggered stochastic hybrid systems

abstract: Stochastic Hybrid Systems (SHS) constitute an important class of mathematical models that integrate discrete stochastic events with continuous dynamics. The time evolution of statistical moments is generally not closed for SHS, in the sense that the time derivative of the lower-order moments depends on higher-order moments. Here, we identify an important class of SHS where moment dynamics is automatically closed, and hence moments can be computed exactly by solving a system of coupled differential equations. This class is referred to as linear time-triggered SHS (TTSHS), where the state evolves according to a linear dynamical system. Stochastic events occur at discrete times and the intervals between them are independent random variables that follow a general class of probability distributions. Moreover, whenever the event occurs, the state of the SHS changes randomly based on a probability distribution. Our approach relies on embedding a Markov chain based on phase-type processes to model timing of events, and showing that the resulting system has closed moment dynamics. Interestingly, we identify a subclass of linear TTSHS, where the first and second-order moments depend only on the mean time interval between events and invariant of their higher-order statistics. Finally, we discuss applicability of our results to different application areas such as network control systems and systems biology.


title: Noncompact RCD(0,N) spaces with linear volume growth

abstract: Since non-compact RCD(0, N) spaces have at least linear volume growth, we study noncompact RCD(0, N) spaces with linear volume growth in this paper. One of the main results is that the diameter of level sets of a Busemann function grow at most linearly on a noncompact RCD(0, N) space satisfying the linear volume growth condition. Another main result in this paper is a splitting theorem at the noncompact end for a RCD(0, N) space with strongly minimal volume growth. These results generalize some theorems on noncompact manifolds with nonnegative Ricci curvature to non-smooth settings.


title: Scaling transition for nonlinear random fields with long-range dependence

abstract: We obtain a complete description of anisotropic scaling limits and the existence of scaling transition for nonlinear functions (Appell polynomials) of stationary linear random fields on $\mathbb{Z}^2$ with moving average coefficients decaying at possibly different rate in the horizontal and vertical direction. The paper extends recent results on scaling transition for linear random fields in Puplinskait\.e and Surgailis (2016), Puplinskait\.e and Surgailis (2015).


title: Time and direction of arrival detection and filtering for imaging in strongly scattering random media

abstract: We study detection and imaging of small reflectors in heavy clutter, using an array of transducers that emits and receives sound waves. Heavy clutter means that multiple scattering of the waves in the heterogeneous host medium is strong and overwhelms the arrivals from the small reflectors. Building on the adaptive time-frequency filter of [1], we propose a robust method for detecting the direction of arrival of the direct echoes from the small reflectors, and suppressing the unwanted clutter backscatter. This improves the resolution of imaging. We illustrate the performance of the method with realistic numerical simulations in a non-destructive testing setup.


title: Confidence Intervals for Causal Effects with Invalid Instruments using Two-Stage Hard Thresholding with Voting

abstract: A major challenge in instrumental variables (IV) analysis is to find instruments that are valid, or have no direct effect on the outcome and are ignorable. Typically one is unsure whether all of the putative IVs are in fact valid. We propose a general inference procedure in the presence of invalid IVs, called Two-Stage Hard Thresholding (TSHT) with voting. TSHT uses two hard thresholding steps to select strong instruments and generate candidate sets of valid IVs. Voting takes the candidate sets and uses majority and plurality rules to determine the true set of valid IVs. In low dimensions, if the sufficient and necessary identification condition under invalid instruments is met, which is more general than the so-called 50% rule or the majority rule, our proposal (i) correctly selects valid IVs, (ii) consistently estimates the causal effect, (iii) produces valid confidence intervals for the causal effect, and (iv) has oracle-optimal width. In high dimensions, we establish nearly identical results without oracle-optimality. In simulations, our proposal outperforms traditional and recent methods in the invalid IV literature. We also apply our method to re-analyze the causal effect of education on earnings.


title: LIGO GW150914 and GW151226 Gravitational Wave Detection and Generalized Gravitation Theory (MOG)

abstract: The nature of gravitational waves in a generalized gravitation theory is investigated. The linearized field equations and the metric tensor quadrupole moment power and the decrease in radius of an inspiralling binary system of two compact objects are derived. The generalized Kerr metric describing a spinning black hole is determined by its mass $M$ and the spin parameter $a=cS/GM^2$. The LIGO-Virgo collaboration data is fitted with smaller binary black hole masses in agreement with the current electromagnetic, observed X-ray binary upper bound for a black hole mass, $M\lesssim 10M_\odot$.


title: Explicit Non-Malleable Extractors, Multi-Source Extractors and Almost Optimal Privacy Amplification Protocols

abstract: We make progress in the following three problems: 1. Constructing optimal seeded non-malleable extractors; 2. Constructing optimal privacy amplification protocols with an active adversary, for any security parameter; 3. Constructing extractors for independent weak random sources, when the min-entropy is extremely small (i.e., near logarithmic). For the first two problems, the best known non-malleable extractors by Chattopadhyay, Goyal and Li [CGL16], and by Cohen [Coh16a,Coh16b] all require seed length and min-entropy at least $\log^2 (1/\epsilon)$, where $\epsilon$ is the error of the extractor. As a result, the best known explicit privacy amplification protocols with an active adversary, which achieve 2 rounds of communication and optimal entropy loss in [Li15c,CGL16], can only handle security parameter up to $s=\Omega(\sqrt{k})$, where $k$ is the min-entropy of the shared secret weak random source. For larger $s$ the best known protocol with optimal entropy loss in [Li15c] requires $O(s/\sqrt{k})$ rounds of communication. In this paper we give an explicit non-malleable extractor that only requires seed length and min-entropy $\log^{1+o(1)} (n/\epsilon)$, which also yields a 2-round privacy amplification protocol with optimal entropy loss for security parameter up to $s=k^{1-\alpha}$ for any constant $\alpha>0$. For the third problem, previously the best known extractor which supports the smallest min-entropy due to Li [Li13a], requires min-entropy $\log^{2+\delta} n$ and uses $O(1/\delta)$ sources, for any constant $\delta>0$. A very recent result by Cohen and Schulman [CS16] improves this, and constructed explicit extractors that use $O(1/\delta)$ sources for min-entropy $\log^{1+\delta} n$, any constant $\delta>0$. In this paper we further improve their result, and give an explicit extractor that uses $O(1)$ (an absolute constant) sources for min-entropy $\log^{1+o(1)} n$.


title: Dynamics of a double-stranded DNA segment in a shear flow

abstract: We study the dynamics of a double-stranded DNA (dsDNA) segment, as a semiflexible polymer, in a shear flow, the strength of which is customarily expressed in terms of the dimensionless Weissenberg number Wi. Polymer chains in shear flows are well-known to undergo tumbling motion. When the chain lengths are much smaller than the persistence length, one expects a (semiflexible) chain to tumble as a rigid rod. At low Wi, a polymer segment shorter than the persistence length does indeed tumble as a rigid rod. However, for higher Wi the chain does not tumble as a rigid rod, even if the polymer segment is shorter than the persistence length. In particular, from time to time the polymer segment may assume a buckled form, a phenomenon commonly known as Euler buckling. Using a bead-spring Hamiltonian model for extensible dsDNA fragments, we first analyze Euler buckling in terms of the oriented deterministic state (ODS), which is obtained as the steady-state solution of the dynamical equations by turning off the stochastic (thermal) forces at a fixed orientation of the chain. The ODS exhibits symmetry breaking at a critical Weissenberg number Wi$_{\text c}$, analogous to a pitchfork bifurcation in dynamical systems. We then follow up the analysis with simulations and demonstrate symmetry breaking in computer experiments, characterized by a unimodal to bimodal transformation of the probability distribution of the second Rouse mode with increasing Wi. Our simulations reveal that shear can cause strong deformation for a chain that is shorter than its persistence length, similar to recent experimental observations.


title: Finite temperature fermion condensate, charge and current densities in a (2+1)-dimensional conical space

abstract: We evaluate the fermion condensate and the expectation values of the charge and current densities for a massive fermionic field in (2+1)-dimensional conical spacetime with a magnetic flux located at the cone apex. The consideration is done for both irreducible representations of the Clifford algebra. The expectation values are decomposed into the vacuum expectation values and contributions coming from particles and antiparticles. All these contributions are periodic functions of the magnetic flux with the period equal to the flux quantum. Related to the non-invariance of the model under the parity and time-reversal transformations, the fermion condensate and the charge density have indefinite parity with respect to the change of the signs of the magnetic flux and chemical potential. The expectation value of the radial current density vanishes. The azimuthal current density is the same for both the irreducible representations of the Clifford algebra. It is an odd function of the magnetic flux and an even function of the chemical potential. The behavior of the expectation values in various asymptotic regions of the parameters are discussed in detail. In particular, we show that for points near the cone apex the vacuum parts dominate. For a massless field with zero chemical potential the fermion condensate and charge density vanish. Simple expressions are derived for the part in the total charge induced by the planar angle deficit and magnetic flux. Combining the results for separate irreducible representations, we also consider the fermion condensate, charge and current densities in parity and time-reversal symmetric models. Possible applications to graphitic nanocones are discussed.


title: PAC-Bayesian bounds for the Gram matrix and least squares regression with a random design

abstract: The topics dicussed in this paper take their origin inthe estimation of the Gram matrix of a random vector from a sample made of n independent copies. They comprise the estimation of the covariance matrix and the study of least squares regression with a random design. We propose four types of results, based on non-asymptotic PAC-Bayesian generalization bounds: a new robust estimator of the Gram matrix and of the covariance matrix, new results on the empirical Gram matrix, new robust least squares estimators and new results on the ordinary least squares estimator, including its exact rate of convergence under polynomial moment assumptions.


title: Restricted rational Cherednik algebras

abstract: We give an overview of the representation theory of restricted rational Cherednik algebras. These are certain finite-dimensional quotients of rational Cherednik algebras at t=0. Their representation theory is connected to the geometry of the Calogero-Moser space, and there is a lot of evidence that they contain certain information about Hecke algebras even though the precise connection is so far unclear. We outline the basic theory along with some open problems and conjectures, and give explicit results in the cyclic and dihedral cases.


title: Convergence of an approximation for rotationally symmetric two-phase lipid bilayer membranes

abstract: We consider a diffuse interface approximation for the lipid phases of rotationally symmetric two-phase bilayer membranes and rigorously derive its $\Gamma$-limit. In particular, we prove that limit vesicles are $C^1$ across interfaces, which justifies a regularity assumption that is widely made in formal asymptotic and numerical studies. Moreover, a limit membrane may consist of several topological spheres, which are connected at the axis of revolution and resemble complete buds of the vesicle.


title: The merging galaxy cluster A520 --- a broken-up cool core, a dark subcluster, and an X-ray channel

abstract: We present results from a deep Chandra X-ray observation of a merging galaxy cluster A520. A high-resolution gas temperature map, after the subtraction of the cluster-scale emission, reveals a long trail of dense, cool clumps --- apparently the fragments of a cool core that has been completely stripped from the infalling subcluster by ram pressure. In this scenario, we can assume that the clumps are still connected by the magnetic field lines. The observed temperature variations imply that thermal conductivity is suppressed by a factor >100 across the presumed direction of the magnetic field (as found in other clusters), and is also suppressed -along- the field lines by a factor of several. Two massive clumps in the periphery of A520, visible in the weak lensing mass map and the X-ray image, have apparently been completely stripped of gas during the merger, but then re-accreted the surrounding high-entropy gas upon exit from the cluster. An X-ray hydrostatic mass estimate for one of the clumps (that has simple geometry) agrees with the lensing mass. Its current gas mass to total mass ratio is very low, 1.5-3%, which makes it a "dark subcluster". We also found a curious low X-ray brightness channel (likely a low-density sheet in projection) going across the cluster along the direction of an apparent secondary merger. The channel may be caused by plasma depletion in a region of an amplified magnetic field (with plasma $\beta\sim 10-20$). The shock in A520 will be studied in a separate paper.


title: Divisor Goldbach Conjecture and its Partition Number

abstract: Based on the Goldbach conjecture and arithmetic fundamental theorem, the Goldbach conjecture was extended to more general situations, i.e., any positive integer can be written as summation of some specific prime numbers, which depends on the divisible factor of this integer, that is: For any positive integer $n~(n>2)$, if there exists an integer $m$, such that $m|n~( 1 < m < n )$, then $n=\sum_{i=1}^m p_i $, where $ p_i~(i=1,2 ,3...m)$ is prime number. In addition, for more prime summands, the combinatorial counting is also discussed. For some special cases, some brief proofs are given. By the use of computer, the preliminary numerical verification was given, there is no an anti-example to be found.


title: The clustering of massive Primordial Black Holes as Dark Matter: measuring their mass distribution with Advanced LIGO

abstract: The recent detection by Advanced LIGO of gravitational waves (GW) from the merging of a binary black hole system sets new limits on the merging rates of massive primordial black holes (PBH) that could be a significant fraction or even the totality of the dark matter in the Universe. aLIGO opens the way to the determination of the distribution and clustering of such massive PBH. If PBH clusters have a similar density to the one observed in ultra-faint dwarf galaxies, we find merging rates comparable to aLIGO expectations. Massive PBH dark matter predicts the existence of thousands of those dwarf galaxies where star formation is unlikely because of gas accretion onto PBH, which would possibly provide a solution to the missing satellite and too-big-to-fail problems. Finally, we study the possibility of using aLIGO and future GW antennas to measure the abundance and mass distribution of PBH in the range [5 - 200] Msun to 10\% accuracy.


title: Convergence of Ricci flows with bounded scalar curvature

abstract: In this paper we prove convergence and compactness results for Ricci flows with bounded scalar curvature and entropy. More specifically, we show that Ricci flows with bounded scalar curvature converge smoothly away from a singular set of codimension $\geq 4$. We also establish a general form of the Hamilton-Tian Conjecture, which is even true in the Riemannian case. These results are based on a compactness theorem for Ricci flows with bounded scalar curvature, which states that any sequence of such Ricci flows converges, after passing to a subsequence, to a metric space that is smooth away from a set of codimension $\geq 4$. In the course of the proof, we will also establish $L^{p < 2}$-curvature bounds on time-slices of such flows.


title: Structure theory of singular spaces

abstract: In this paper we generalize the theory of Cheeger, Colding and Naber to certain singular spaces that arise as limits of sequences of Riemannian manifolds. This theory will have applications in the analysis of Ricci flows of bounded curvature, which we will describe in a subsequent paper.


title: The sharp threshold for the Duarte model

abstract: The class of critical bootstrap percolation models in two dimensions was recently introduced by Bollob\'as, Smith and Uzzell, and the critical threshold for percolation was determined up to a constant factor for all such models by the authors of this paper. Here we develop and refine the techniques introduced in that paper in order to determine a sharp threshold for the Duarte model. This resolves a question of Mountford from 1995, and is the first result of its type for a model with drift.


title: A Universal Coding Scheme for Remote Generation of Continuous Random Variables

abstract: We consider a setup in which Alice selects a pdf $f$ from a set of prescribed pdfs $\mathscr{P}$ and sends a prefix-free codeword $W$ to Bob in order to allow him to generate a single instance of the random variable $X\sim f$. We describe a universal coding scheme for this setup and establish an upper bound on the expected codeword length when the pdf $f$ is bounded, orthogonally concave (which includes quasiconcave pdf), and has a finite first absolute moment. A dyadic decomposition scheme is used to express the pdf as a mixture of uniform pdfs over hypercubes. Alice randomly selects a hypercube according to its weight, encodes its position and size into $W$, and sends it to Bob who generates $X$ uniformly over the hypercube. Compared to previous results on channel simulation, our coding scheme applies to any continuous distribution and does not require two-way communication or shared randomness. We apply our coding scheme to classical simulation of quantum entanglement and obtain a better bound on the average codeword length than previously known.


title: Constraints on decaying early modified gravity from cosmological observations

abstract: Most of the information on our cosmos stems from either late-time observations or the imprint of early-time inhomogeneities on the cosmic microwave background. We explore to what extent early modifications of gravity, which become significant after recombination but then decay towards the present, can be constrained by current cosmological observations. For the evolution of the gravitational modification, we adopt the decaying mode of a hybrid-metric Palatini $f(\mathcal{R})$ gravity model which is designed to reproduce the standard cosmological background expansion history and due to the decay of the modification is naturally compatible with Solar-System tests. We embed the model in the effective field theory description of Horndeski scalar-tensor gravity with an early-time decoupling of the gravitational modification. Since the quasistatic approximation for the perturbations in the model breaks down at high redshifts, where modifications remain relevant, we introduce a computationally efficient correction to describe the evolution of the scalar field fluctuation in this regime. We compare the decaying early-time modification against geometric probes and recent Planck measurements and find no evidence for such effects in the observations. Current data constrains the scalar field value at $|f_{\mathcal{R}}(z=z_{\rm on})| \lesssim 10^{-2}$ for modifications introduced at redshifts $z_{\rm on}\sim(500-1000)$ with present-day value $|f_{\mathcal{R}0}|\lesssim10^{-8}$. Finally, we comment on constraints that will be achievable with future 21~cm surveys and gravitational wave experiments.


title: Bitcoin Mining Decentralization via Cost Analysis

abstract: Bitcoin mining presents a significant economic incentive for efficient hashing and broadcast of data, both parameters stemming from the Proofs of Work used to advance the network. This incentive has led to the development of Bitcoin specific application specific integrated circuits and centralized mining pools, undermining the decentralized motivations behind Bitcoin's design. In addition, the imminent block reward halving threatens the profitability of mining at any scale. Some work has been done in formal models for miner profitability, but existing models do not account for conditions such as the pricing of off-peak power and diverse investment strategies regarding sunken costs. There is also a lack of formal study of how the profit model changes as mining scales from the individual to the industrial level. Given the lack of analysis of these conditions, there are alternative models for profitable or net zero mining that operate at smaller, and therefore more desirable, scale.