Understanding the Formation of Pbse Honeycomb Superstructures by Dynamics Simulations

Total Page:16

File Type:pdf, Size:1020Kb

Understanding the Formation of Pbse Honeycomb Superstructures by Dynamics Simulations PHYSICAL REVIEW X 9, 021015 (2019) Understanding the Formation of PbSe Honeycomb Superstructures by Dynamics Simulations Giuseppe Soligno* and Daniel Vanmaekelbergh Condensed Matter and Interfaces, Debye Institute for Nanomaterials Science, Utrecht University, Princetonplein 1, Utrecht 3584 CC, Netherlands (Received 28 October 2018; revised manuscript received 28 January 2019; published 23 April 2019) Using a coarse-grained molecular dynamics model, we simulate the self-assembly of PbSe nanocrystals (NCs) adsorbed at a flat fluid-fluid interface. The model includes all key forces involved: NC-NC short- range facet-specific attractive and repulsive interactions, entropic effects, and forces due to the NC adsorption at fluid-fluid interfaces. Realistic values are used for the input parameters regulating the various forces. The interface-adsorption parameters are estimated using a recently introduced sharp- interface numerical method which includes capillary deformation effects. We find that the final structure in which the NCs self-assemble is drastically affected by the input values of the parameters of our coarse- grained model. In particular, by slightly tuning just a few parameters of the model, we can induce NC self- assembly into either silicene-honeycomb superstructures, where all NCs have a f111g facet parallel to the fluid-fluid interface plane, or square superstructures, where all NCs have a f100g facet parallel to the interface plane. Both of these nanostructures have been observed experimentally. However, it is still not clear their formation mechanism, and, in particular, which are the factors directing the NC self-assembly into one or another structure. In this work, we identify and quantify such factors, showing illustrative assembled-phase diagrams obtained from our simulations. In addition, with our model, we can study the self-assembly dynamics, simulating how the NCs’ structures evolve from few-NCs aggregates to gradually larger domains. For example, we observe linear chains, where all NCs have a f110g facet parallel to the interface plane as typical precursors of the square superstructure, and zigzag aggregates, where all NCs have a f111g facet parallel to the interface plane as typical precursors of the silicene-honeycomb superstructure. Both of these aggregates have also been observed experimentally. Finally, we show indications that our method can be applied to study defects of the obtained superstructures. DOI: 10.1103/PhysRevX.9.021015 Subject Areas: Computational Physics, Materials Science, Nanophysics I. INTRODUCTION superstructure: (i) synthesis of the NCs, i.e., the super- structure building blocks, from their atomic precursors Nanoparticle self-assembly is an emerging route with [48–60], and (ii) self-assembly of the NCs into ordered tremendous potential to build novel nanostructured materi- structures [61–66]. Typically, NCs are, after their synthesis, als [1,2]. In the past two decades, increasing interest has dispersed in solution and fully covered by organic molecules been devoted toward the formation of 3D [3–20] and 2D (ligands), chemically attached (chemisorbed) at the NC [21–32] nanogeometric materials for optoelectronic appli- surface, which largely screen NC-NC attractive interactions. cations [33–47]. These nanomaterials are formed by the The NC self-assembly can be induced, for example, by self-assembly of semiconductor nanocrystals (NCs) with solvent evaporation (i.e., by increasing the NC density). In the size of a few nanometers, and often with a roughly this case, a close-packed NC superstructure is obtained, spherical shape, into various types of superstructures. whose geometry is dictated mainly by entropic factors, i.e., Two main stages are required for the formation of a the shape and size of the NCs [67–71]. NC-NC attractive interactions can be activated (e.g., by ligand-exchange [72] *Corresponding author. or facet-specific ligand-detachment [73] reactions), inducing [email protected] NC self-assembly in different types of superstructures. A special class of NC superstructures, referred to as NC Published by the American Physical Society under the terms of superlattices, is formed when NC self-assembly is followed the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to by crystal-structure alignment and oriented attachment the author(s) and the published article’s title, journal citation, [74–83] of close-by NCs, resulting in atomically coherent and DOI. nanogeometric materials. The formation process of NC 2160-3308=19=9(2)=021015(32) 021015-1 Published by the American Physical Society GIUSEPPE SOLIGNO and DANIEL VANMAEKELBERGH PHYS. REV. X 9, 021015 (2019) superlattices can be confined in 2D when the NCs are superstructures obtained with our MD model present the adsorbed at a fluid-fluid interface [84]. The NCs self- same kind of defects observed in the experimental super- assemble in ordered monolayers at the fluid-fluid interface, lattices, indicating that our MD model can also be used to and subsequently perform oriented attachment [73,85–88], study the formation and stability of defects in these resulting in 2D atomically coherent nanogeometric materi- structures. als. In this last class of processes, in addition to NC entropy and NC-NC interactions, also the different interactions of the II. COARSE-GRAINED MD MODEL NCs with the different molecules of the two fluids forming First, we briefly illustrate our simulation model (more the interface (i.e., interface-adsorption effects) are involved details are provided in Appendix A). in the superlattice formation mechanism. First experimentally observed a few years ago [73,85,86,88], the PbSe silicene-honeycomb 2D super- A. Outline lattice is expected to exhibit a Dirac-type band structure, We model a PbSe NC using a polybead structure. The with the semiconductor band gap preserved [89–91]. beads forming the NC are disposed to reproduce a Hence, such a material would combine the properties of rhombicuboctahedron of size 6 nm, that is a typical semiconductors with those of graphene, making it very experimental shape [76,92]. We simulate the dynamics interesting for optoelectronic applications. This superlattice of the NCs using the position-Verlet algorithm [94] to is formed by PbSe NCs, approximately having a 5-nm size compute the bead’s motion. Bead-bead pair potentials are and a rhombicuboctahedron shape, disposed as in the used to reproduce NC-NC short-range facet-specific attrac- silicene lattice, i.e., in a periodically buckled honeycomb tive and repulsive interactions. The distance between beads lattice, and all oriented with a f111g facet parallel to the belonging to the same NC is kept fixed by constraint forces superlattice plane. The PbSe NCs, before the self-assembly, [95]. The solvent is treated implicitly by modeling the NC i.e., when still dispersed in solvent, are fully covered by Brownian motion. External potentials are applied to the NC oleic acid ligands, which are strongly attached at the f111g, beads to mimic the interface-adsorption forces experienced f110g facets and weakly bonded at the f100g facets [92]. by NCs at a fluid-fluid interface. The formation process of the PbSe silicene-honeycomb superlattice typically occurs at the solvent-air interface. B. NC-NC interactions During the process, ligands desorb from the NC f100g facets, allowing oriented attachment between PbSe NCs The polybead structure representing each NC is formed by opposite f100g facets [73]. In similar experiments by 144 beads (shell beads) disposed to represent the NC [73,85,87], self-assembly of NCs into square superlattices surface, see Fig. 1, plus seven beads (core beads) used, see was observed, with all NCs oriented with a f100g facet later, for the Brownian motion of the NCs. To model parallel to the superlattice plane, and with oriented attach- van der Waals and electrostatic-chemical atomic inter- ment between PbSe NCs by opposite f100g facets also actions between NCs at almost-contact distance, shell performed. Although the relevant forces involved in the beads of different NCs interact with each other by the formation of these superlattices are known [93], the precise Lennard-Jones pair potential, mechanisms are far from understood. In particular, it is not σ 12 σ 6 clear yet why the PbSe NCs form square superlattices in U ¼ 4ϵ − ; ð1Þ some cases and silicene-honeycomb superlattices in others. LJ r r In this work, we shed light on this, introducing a new coarse-grained molecular dynamics (MD) model to study truncated and shifted in r ¼ 2.4σ [see Eq. (A3)], where r is the self-assembly of NCs at fluid-fluid interfaces. By the center-to-center bead-bead distance, σ ¼ 1 nm, and ϵ presenting simulations with our MD model, we illustrate sets the interaction strength for the bead pair. PbSe NCs are, how the intricate interplay between interface-adsorption after their synthesis, dispersed in solution, and ligands, forces orienting and keeping the NCs at the interface and typically oleic acid molecules, are chemisorbed at the NC attractive forces only between f100g facets of close-by facets. At this stage, ligands largely screen NC-NC attrac- NCs drive the self-assembly of PbSe NCs into either square tive interactions and prevent NC assembly. In the experi- or silicene-honeycomb superstructures. As a confirmation ments of interest [73,85–88], the NC self-assembly is that our model correctly captures the key features of the activated by chemically inducing the detachment of ligands NC self-assembly, in our simulations we observe, in the from the PbSe NC f100g facets, i.e., where the bonding formation of both the silicene-honeycomb and square energy between ligand molecules and NC surface atoms is superstructures, the same few-NCs precursors observed lowest [73,92]. Hence, during the self-assembly, NCs experimentally, that is zigzag aggregates, where all NCs attract (and attach to) each other only by f100g facets, have a f111g facet parallel to the interface plane, and while NC-NC bonds by f110g and f111g facets are linear chains, where all NCs have a f110g facet parallel prevented by ligands.
Recommended publications
  • 1 Lifts of Polytopes
    Lecture 5: Lifts of polytopes and non-negative rank CSE 599S: Entropy optimality, Winter 2016 Instructor: James R. Lee Last updated: January 24, 2016 1 Lifts of polytopes 1.1 Polytopes and inequalities Recall that the convex hull of a subset X n is defined by ⊆ conv X λx + 1 λ x0 : x; x0 X; λ 0; 1 : ( ) f ( − ) 2 2 [ ]g A d-dimensional convex polytope P d is the convex hull of a finite set of points in d: ⊆ P conv x1;:::; xk (f g) d for some x1;:::; xk . 2 Every polytope has a dual representation: It is a closed and bounded set defined by a family of linear inequalities P x d : Ax 6 b f 2 g for some matrix A m d. 2 × Let us define a measure of complexity for P: Define γ P to be the smallest number m such that for some C s d ; y s ; A m d ; b m, we have ( ) 2 × 2 2 × 2 P x d : Cx y and Ax 6 b : f 2 g In other words, this is the minimum number of inequalities needed to describe P. If P is full- dimensional, then this is precisely the number of facets of P (a facet is a maximal proper face of P). Thinking of γ P as a measure of complexity makes sense from the point of view of optimization: Interior point( methods) can efficiently optimize linear functions over P (to arbitrary accuracy) in time that is polynomial in γ P . ( ) 1.2 Lifts of polytopes Many simple polytopes require a large number of inequalities to describe.
    [Show full text]
  • The Minkowski Problem for Simplices
    The Minkowski Problem for Simplices Daniel A. Klain Department of Mathematical Sciences University of Massachusetts Lowell Lowell, MA 01854 USA Daniel [email protected] Abstract The Minkowski existence Theorem for polytopes follows from Cramer’s Rule when attention is limited to the special case of simplices. It is easy to see that a convex polygon in R2 is uniquely determined (up to translation) by the directions and lengths of its edges. This suggests the following (less easily answered) question in higher dimensions: given a collection of proposed facet normals and facet areas, is there a convex polytope in Rn whose facets fit the given data, and, if so, is the resulting polytope unique? This question (along with its answer) is known as the Minkowski problem. For a polytope P in Rn denote by V (P ) the volume of P . If Q is a polytope in Rn having dimension strictly less than n, then denote v(Q) the (n ¡ 1)-dimensional volume of Q. For any u non-zero vector u, let P denote the face of P having u as an outward normal, and let Pu denote the orthogonal projection of P onto the hyperplane u?. The Minkowski problem for polytopes concerns the following specific question: Given a collection u1; : : : ; uk of unit vectors and ®1; : : : ; ®k > 0, under what condition does there exist a polytope P ui having the ui as its facet normals and the ®i as its facet areas; that is, such that v(P ) = ®i for each i? A necessary condition on the facet normals and facet areas is given by the following proposition [BF48, Sch93].
    [Show full text]
  • 15 BASIC PROPERTIES of CONVEX POLYTOPES Martin Henk, J¨Urgenrichter-Gebert, and G¨Unterm
    15 BASIC PROPERTIES OF CONVEX POLYTOPES Martin Henk, J¨urgenRichter-Gebert, and G¨unterM. Ziegler INTRODUCTION Convex polytopes are fundamental geometric objects that have been investigated since antiquity. The beauty of their theory is nowadays complemented by their im- portance for many other mathematical subjects, ranging from integration theory, algebraic topology, and algebraic geometry to linear and combinatorial optimiza- tion. In this chapter we try to give a short introduction, provide a sketch of \what polytopes look like" and \how they behave," with many explicit examples, and briefly state some main results (where further details are given in subsequent chap- ters of this Handbook). We concentrate on two main topics: • Combinatorial properties: faces (vertices, edges, . , facets) of polytopes and their relations, with special treatments of the classes of low-dimensional poly- topes and of polytopes \with few vertices;" • Geometric properties: volume and surface area, mixed volumes, and quer- massintegrals, including explicit formulas for the cases of the regular simplices, cubes, and cross-polytopes. We refer to Gr¨unbaum [Gr¨u67]for a comprehensive view of polytope theory, and to Ziegler [Zie95] respectively to Gruber [Gru07] and Schneider [Sch14] for detailed treatments of the combinatorial and of the convex geometric aspects of polytope theory. 15.1 COMBINATORIAL STRUCTURE GLOSSARY d V-polytope: The convex hull of a finite set X = fx1; : : : ; xng of points in R , n n X i X P = conv(X) := λix λ1; : : : ; λn ≥ 0; λi = 1 : i=1 i=1 H-polytope: The solution set of a finite system of linear inequalities, d T P = P (A; b) := x 2 R j ai x ≤ bi for 1 ≤ i ≤ m ; with the extra condition that the set of solutions is bounded, that is, such that m×d there is a constant N such that jjxjj ≤ N holds for all x 2 P .
    [Show full text]
  • Frequently Asked Questions in Polyhedral Computation
    Frequently Asked Questions in Polyhedral Computation http://www.ifor.math.ethz.ch/~fukuda/polyfaq/polyfaq.html Komei Fukuda Swiss Federal Institute of Technology Lausanne and Zurich, Switzerland [email protected] Version June 18, 2004 Contents 1 What is Polyhedral Computation FAQ? 2 2 Convex Polyhedron 3 2.1 What is convex polytope/polyhedron? . 3 2.2 What are the faces of a convex polytope/polyhedron? . 3 2.3 What is the face lattice of a convex polytope . 4 2.4 What is a dual of a convex polytope? . 4 2.5 What is simplex? . 4 2.6 What is cube/hypercube/cross polytope? . 5 2.7 What is simple/simplicial polytope? . 5 2.8 What is 0-1 polytope? . 5 2.9 What is the best upper bound of the numbers of k-dimensional faces of a d- polytope with n vertices? . 5 2.10 What is convex hull? What is the convex hull problem? . 6 2.11 What is the Minkowski-Weyl theorem for convex polyhedra? . 6 2.12 What is the vertex enumeration problem, and what is the facet enumeration problem? . 7 1 2.13 How can one enumerate all faces of a convex polyhedron? . 7 2.14 What computer models are appropriate for the polyhedral computation? . 8 2.15 How do we measure the complexity of a convex hull algorithm? . 8 2.16 How many facets does the average polytope with n vertices in Rd have? . 9 2.17 How many facets can a 0-1 polytope with n vertices in Rd have? . 10 2.18 How hard is it to verify that an H-polyhedron PH and a V-polyhedron PV are equal? .
    [Show full text]
  • A Note on Solid Coloring of Pure Simplicial Complexes Joseph O'rourke Smith College, [email protected]
    Masthead Logo Smith ScholarWorks Computer Science: Faculty Publications Computer Science 12-17-2010 A Note on Solid Coloring of Pure Simplicial Complexes Joseph O'Rourke Smith College, [email protected] Follow this and additional works at: https://scholarworks.smith.edu/csc_facpubs Part of the Computer Sciences Commons, and the Discrete Mathematics and Combinatorics Commons Recommended Citation O'Rourke, Joseph, "A Note on Solid Coloring of Pure Simplicial Complexes" (2010). Computer Science: Faculty Publications, Smith College, Northampton, MA. https://scholarworks.smith.edu/csc_facpubs/37 This Article has been accepted for inclusion in Computer Science: Faculty Publications by an authorized administrator of Smith ScholarWorks. For more information, please contact [email protected] A Note on Solid Coloring of Pure Simplicial Complexes Joseph O'Rourke∗ December 21, 2010 Abstract We establish a simple generalization of a known result in the plane. d The simplices in any pure simplicial complex in R may be colored with d+1 colors so that no two simplices that share a (d−1)-facet have the same 2 color. In R this says that any planar map all of whose faces are triangles 3 may be 3-colored, and in R it says that tetrahedra in a collection may be \solid 4-colored" so that no two glued face-to-face receive the same color. 1 Introduction The famous 4-color theorem says that the regions of any planar map may be colored with four colors such that no two regions that share a positive-length border receive the same color. A lesser-known special case is that if all the regions are triangles, three colors suffice.
    [Show full text]
  • Linear Programming and Polyhedral Combinatorics Summary of What Was Seen in the Introductory Lectures on Linear Programming and Polyhedral Combinatorics
    Massachusetts Institute of Technology Handout 8 18.433: Combinatorial Optimization March 6th, 2007 Michel X. Goemans Linear Programming and Polyhedral Combinatorics Summary of what was seen in the introductory lectures on linear programming and polyhedral combinatorics. Definition 1 A halfspace in Rn is a set of the form fx 2 Rn : aT x ≤ bg for some vector a 2 Rn and b 2 R. Definition 2 A polyhedron is the intersection of finitely many halfspaces: P = fx 2 Rn : Ax ≤ bg. Definition 3 A polytope is a bounded polyhedron. n Definition 4 If P is a polyhedron in R , the projection Pk of P is defined as fy = (x1; x2; · · · ; xk−1; xk+1; · · · ; xn) : x 2 P for some xkg. We claim that Pk is also a polyhedron and this can be proved by giving an explicit description of Pk in terms of linear inequalities. For this purpose, one uses Fourier-Motzkin elimination. Let P = fx : Ax ≤ bg and let • S+ = fi : aik > 0g, • S− = fi : aik < 0g, • S0 = fi : aik = 0g. T Clearly, any element in Pk must satisfy the inequality ai x ≤ bi for all i 2 S0 (these inequal- ities do not involve xk). Similarly, we can take a linear combination of an inequality in S+ and one in S− to eliminate the coefficient of xk. This shows that the inequalities: aik aljxj − alk akjxj ≤ aikbl − alkbi (1) j ! j ! X X for i 2 S+ and l 2 S− are satisfied by all elements of Pk. Conversely, for any vector (x1; x2; · · · ; xk−1; xk+1; · · · ; xn) satisfying (1) for all i 2 S+ and l 2 S− and also T ai x ≤ bi for all i 2 S0 (2) we can find a value of xk such that the resulting x belongs to P (by looking at the bounds on xk that each constraint imposes, and showing that the largest lower bound is smaller than the smallest upper bound).
    [Show full text]
  • Model Reduction in Geometric Tolerancing by Polytopes Vincent Delos, Santiago Arroyave-Tobón, Denis Teissandier
    Model reduction in geometric tolerancing by polytopes Vincent Delos, Santiago Arroyave-Tobón, Denis Teissandier To cite this version: Vincent Delos, Santiago Arroyave-Tobón, Denis Teissandier. Model reduction in geometric tolerancing by polytopes. Computer-Aided Design, Elsevier, 2018. hal-01741395 HAL Id: hal-01741395 https://hal.archives-ouvertes.fr/hal-01741395 Submitted on 23 Mar 2018 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Model reduction in geometric tolerancing by polytopes Vincent Delos, Santiago Arroyave-Tobon,´ Denis Teissandier Univ. Bordeaux, I2M, UMR 5295, F-33400 Talence, France Abstract There are several models used in mechanical design to study the behavior of mechanical systems involving geometric variations. By simulating defects with sets of constraints it is possible to study simultaneously all the configurations of mechanisms, whether over-constrained or not. Using this method, the accumulation of defects is calculated by summing sets of constraints derived from features (toleranced surfaces and joints) in the tolerance chain. These sets are usually unbounded objects (R6-polyhedra, 3 parameters for the small rotation, 3 for the small translation), due to the unbounded displacements associated with the degrees of freedom of features.
    [Show full text]
  • Convex Hull in Higher Dimensions
    Comp 163: Computational Geometry Professor Diane Souvaine Tufts University, Spring 2005 1991 Scribe: Gabor M. Czako, 1990 Scribe: Sesh Venugopal, 2004 Scrib Higher Dimensional Convex Hull Algorithms 1 Convex hulls 1.0.1 Definitions The following definitions will be used throughout. Define • Sd:Ad-Simplex The simplest convex polytope in Rd.Ad-simplex is always the convex hull of some d + 1 affinely independent points. For example, a line segment is a 1 − simplex i.e., the smallest convex subspace which con- tains two points. A triangle is a 2 − simplex and a tetrahedron is a 3 − simplex. • P: A Simplicial Polytope. A polytope where each facet is a d − 1 simplex. By our assumption, a convex hull in 3-D has triangular facets and in 4-D, a convex hull has tetrahedral facets. •P: number of facets of Pi • F : a facet of P • R:aridgeofF • G:afaceofP d+1 • fj(Pi ): The number of j-faces of P .Notethatfj(Sd)=C(j+1). 1.0.2 An Upper Bound on Time Complexity: Cyclic Polytope d 1 2 d Consider the curve Md in R defined as x(t)=(t ,t , ..., t ), t ∈ R.LetH be the hyperplane defined as a0 + a1x1 + a2x2 + ... + adxd = 0. The intersection 2 d of Md and H is the set of points that satisfy a0 + a1t + a2t + ...adt =0. This polynomial has at most d real roots and therefore d real solutions. → H intersects Md in at most d points. This brings us to the following definition: A convex polytope in Rd is a cyclic polytope if it is the convex hull of a 1 set of at least d +1pointsonMd.
    [Show full text]
  • Math.CO] 14 May 2002 in Npriua,Tecne Ulof Hull Convex Descrip the Simple Particular, Similarly in No Tion
    Fat 4-polytopes and fatter 3-spheres 1, 2, † 3, ‡ David Eppstein, ∗ Greg Kuperberg, and G¨unter M. Ziegler 1Department of Information and Computer Science, University of California, Irvine, CA 92697 2Department of Mathematics, University of California, Davis, CA 95616 3Institut f¨ur Mathematik, MA 6-2, Technische Universit¨at Berlin, D-10623 Berlin, Germany φ We introduce the fatness parameter of a 4-dimensional polytope P, defined as (P)=( f1 + f2)/( f0 + f3). It arises in an important open problem in 4-dimensional combinatorial geometry: Is the fatness of convex 4- polytopes bounded? We describe and analyze a hyperbolic geometry construction that produces 4-polytopes with fatness φ(P) > 5.048, as well as the first infinite family of 2-simple, 2-simplicial 4-polytopes. Moreover, using a construction via finite covering spaces of surfaces, we show that fatness is not bounded for the more general class of strongly regular CW decompositions of the 3-sphere. 1. INTRODUCTION Only the two extreme cases of simplicial and of simple 4- polytopes (or 3-spheres) are well-understood. Their f -vectors F correspond to faces of the convex hull of F , defined by the The characterization of the set 3 of f -vectors of convex 4 valid inequalities f 2 f and f 2 f , and the g-Theorem, 3-dimensional polytopes (from 1906, due to Steinitz [28]) is 2 ≥ 3 1 ≥ 0 well-known and explicit, with a simple proof: An integer vec- proved for 4-polytopes by Barnette [1] and for 3-spheres by Walkup [34], provides complete characterizations of their f - tor ( f0, f1, f2) is the f -vector of a 3-polytope if and only if it satisfies vectors.
    [Show full text]
  • The Strongest Facets of the Acyclic Subgraph Polytope Are Unknown
    The Strongest Facets of the Acyclic Subgraph Polytope Are Unknown Michel X. Goemans 1. and Leslie A. Hall 2.* 1 Dept. of Mathematics, Room 2-382, M.I.T., Cambridge, MA 02139.*** 2 Dept. of Math. Sciences, The Johns Hopkins University, Baltimore, MD 21218. ? Abstract. We consider the acyclic subgraph polytope and define the no- tion of strength of a relaxation as the maximum improvement obtained by using this relaxation instead of the most trivial relaxation of the problem. We show that the strength of a relaxation is the maximum of the strengths of the relaxations obtained by simply adding to the trivial relaxation each valid inequality separately. We also derive from the prob- abilistic method that the maximum strength of any inequality is 2. We then consider all (or almost all) the known valid inequalities for the poly- tope and compute their strength. The surprising observation is that their strength is at most slightly more than 3/2, implying that the strongest inequalities are yet unknown. We then consider a pseudo-random con- struction due to Alon and Spencer based on quadratic residues to obtain new facet-defining inequalities for the polytope. These are also facet- defining for the linear ordering polytope. 1 Introduction Given weights wa on the arcs of a complete directed graph (or digraph) D = (V, A), the acyclic subgraph problem is that of determining a set of arcs of maxi- mum total weight that define an acyclic subgraph. The complement of an acyclic subgraph is called a feedback arc set. For general graphs, the acyclic subgraph problem is NP-hard, even for graphs with unit weights and with total indegree and outdegree of every vertex no more than three [GJ79], although the prob- lem is polynomially solvable for planar graphs as was shown by Lucchesi and Younger [LY78].
    [Show full text]
  • Arxiv:1501.06198V1 [Math.MG]
    EMBEDDED FLEXIBLE SPHERICAL CROSS-POLYTOPES WITH NON-CONSTANT VOLUMES ALEXANDER A. GAIFULLIN Abstract. We construct examples of embedded flexible cross-polytopes in the spheres of all dimensions. These examples are interesting from two points of view. First, in dimensions 4 and higher, they are the first examples of embedded flexible polyhedra. Notice that, unlike in the spheres, in the Euclidean spaces and the Lobachevsky spaces of dimensions 4 and higher, still no example of an embedded flexible polyhedron is known. Second, we show that the volumes of the constructed flexible cross-polytopes are non- constant during the flexion. Hence these cross-polytopes give counterexamples to the Bellows Conjecture for spherical polyhedra. Earlier a counterexample to this conjecture was built only in dimension 3 (Alexandrov, 1997), and was not embedded. For flexible polyhedra in spheres we suggest a weakening of the Bellows Conjecture, which we call the Modified Bellows Conjecture. We show that this conjecture holds for all flexible cross-polytopes of the simplest type among which there are our counterexamples to the usual Bellows Conjecture. By the way, we obtain several geometric results on flexible cross-polytopes of the simplest type. In particular, we write relations on the volumes of their faces of codimensions 1 and 2. To Nicolai Petrovich Dolbilin on the occasion of his seventieth birthday 1. Introduction Let Xn be one of the three n-dimensional spaces of constant curvature, that is, the Euclidean space En or the sphere Sn or the Lobachevsky space Λn. For convenience, we shall always normalize metrics on the sphere Sn and on the Lobachevsky space Λn so that their curvatures are equal to 1 and 1 respectively.
    [Show full text]
  • Smooth Fano Polytopes with Many Vertices
    SMOOTH FANO POLYTOPES WITH MANY VERTICES BENJAMIN ASSARF, MICHAEL JOSWIG, AND ANDREAS PAFFENHOLZ Abstract. The d-dimensional simplicial, terminal, and reflexive polytopes with at least 3d−2 vertices are classified. In particular, it turns out that all of them are smooth Fano polytopes. This improves on previous results of (Casagrande, 2006) and (Øbro, 2008). Smooth Fano polytopes play a role in algebraic geometry and mathematical physics. 1. Introduction A lattice polytope P is a convex polytope whose vertices lie in a lattice N contained in the vector space Rd. Fixing a basis of N describes an isomorphism to Zd. Throughout this paper, we restrict our attention to the standard lattice N = Zd. A d-dimensional lattice polytope P ⊂ Rd is called reflexive if it contains the origin 0 as an interior point and its polar polytope is a lattice polytope in the dual lattice M := Hom(N, Z) =∼ Zd. A lattice polytope P is terminal if 0 and the vertices are the only lattice points in P ∩ Zd. It is simplicial if each face is a simplex. We say that P is a smooth Fano polytope if P ⊆ Rd is simplicial with 0 in the interior and the vertices of each facet form a lattice basis of Zd. The fan where every cone is the non negative linear span over a face is called the face fan. It is dual to the normal fan, which is the collection of all normal cones. In algebraic geometry, reflexive polytopes correspond to Gorenstein toric Fano va- rieties. The toric variety XP of a polytope P is determined by the face fan of P , that is, the fan spanned by all faces of P ; see Ewald [10] or Cox, Little, and Schenck [8] for more details.
    [Show full text]