From stripe to checkerboard order on the square lattice in the presence of quenched disorder

We discuss the effects of quenched disorder on a model of charge density wave (CDW) ordering on the square lattice. Our model may be applicable to the cuprate superconductors, where a random electrostatic potential exists in the CuO2 planes as a result of the presence of charged dopants. We argue that the presence of a random potential can affect the unidirectionality of the CDW order, characterized by an Ising order parameter. Coupling to a unidirectional CDW, the random potential can lead to the formation of domains with 90 degree relative orientation, thus tending to restore the rotational symmetry of the underlying lattice. We find that the correlation length of the Ising order can be significantly larger than the CDW correlation length. For a checkerboard CDW on the other hand, disorder generates spatial anisotropies on short length scales and thus some degree of unidirectionality. We quantify these disorder effects and suggest new techniques for analyzing the local density of states (LDOS) data measured in scanning tunneling microscopy experiments.Comment: 10 pages, 11 figures; added referenc

On the stability of viscous free-surface flow supported by a rotating cylinder

Using an adaptive finite-element (FE) scheme developed recently by the authors, we shed new light on the long-standing fundamental problem of the unsteady free-surface Stokes flow exterior to a circular cylinder rotating about its horizontal axis in a vertical gravitational field. For supportable loads, we observe that the steady-state is more readily attained for near-maximal fluid loads on the cylinder than for significantly sub-maximal loads. For the latter, we investigate large-time dynamics by means of a finite-difference approximation to the thin-film equations, which is also used to validate the adaptive FE simulations (applied to the full Stokes equations) for these significantly sub-maximal loads. Conversely, by comparing results of the two methods, we assess the validity of the thin-film approximation as either the load is increased or the rotation rate of the cylinder is decreased. Results are presented on the independent effects of gravity, surface tension and initial film thickness on the decay to steady-state. Finally, new numerical simulations of load shedding are presented

Analysis of and workarounds for element reversal for a finite element-based algorithm for warping triangular and tetrahedral meshes

We consider an algorithm called FEMWARP for warping triangular and tetrahedral finite element meshes that computes the warping using the finite element method itself. The algorithm takes as input a two- or three-dimensional domain defined by a boundary mesh (segments in one dimension or triangles in two dimensions) that has a volume mesh (triangles in two dimensions or tetrahedra in three dimensions) in its interior. It also takes as input a prescribed movement of the boundary mesh. It computes as output updated positions of the vertices of the volume mesh. The first step of the algorithm is to determine from the initial mesh a set of local weights for each interior vertex that describes each interior vertex in terms of the positions of its neighbors. These weights are computed using a finite element stiffness matrix. After a boundary transformation is applied, a linear system of equations based upon the weights is solved to determine the final positions of the interior vertices. The FEMWARP algorithm has been considered in the previous literature (e.g., in a 2001 paper by Baker). FEMWARP has been succesful in computing deformed meshes for certain applications. However, sometimes FEMWARP reverses elements; this is our main concern in this paper. We analyze the causes for this undesirable behavior and propose several techniques to make the method more robust against reversals. The most successful of the proposed methods includes combining FEMWARP with an optimization-based untangler.Comment: Revision of earlier version of paper. Submitted for publication in BIT Numerical Mathematics on 27 April 2010. Accepted for publication on 7 September 2010. Published online on 9 October 2010. The final publication is available at http://www.springerlink.co

Evaluating Elevated Convection with the Downdraft Convective Inhibition

A method for evaluating the penetration of a stable layer by an elevated convective downdraft is discussed. Some controversy exists on the community’s ability to define truly elevated convection from surface-based convection. By comparing the downdraft convective inhibition (DCIN) to the downdraft convective available potential energy (DCAPE), we determine that downdraft penetration potential is progressively enabled as the DCIN is progressively smaller than the DCAPE; inversely as DCIN increases over DCAPE, so does the likelihood of purely elevated convection. Serial vertical soundings and accompanying analyses are provided to support this finding

Incoherent dynamics in the toric code subject to disorder

We numerically study the effects of two forms of quenched disorder on the anyons of the toric code. Firstly, a new class of codes based on random lattices of stabilizer operators is presented, and shown to be superior to the standard square lattice toric code for certain forms of biased noise. It is further argued that these codes are close to optimal, in that they tightly reach the upper bound of error thresholds beyond which no correctable CSS codes can exist. Additionally, we study the classical motion of anyons in toric codes with randomly distributed onsite potentials. In the presence of repulsive long-range interaction between the anyons, a surprising increase with disorder strength of the lifetime of encoded states is reported and explained by an entirely incoherent mechanism. Finally, the coherent transport of the anyons in the presence of both forms of disorder is investigated, and a significant suppression of the anyon motion is found.Comment: 13 pages, 12 figure

Exploiting asynchrony from exact forward recovery for DUE in iterative solvers

This paper presents a method to protect iterative solvers from Detected and Uncorrected Errors (DUE) relying on error detection techniques already available in commodity hardware. Detection operates at the memory page level, which enables the use of simple algorithmic redundancies to correct errors. Such redundancies would be inapplicable under coarse grain error detection, but become very powerful when the hardware is able to precisely detect errors. Relations straightforwardly extracted from the solver allow to recover lost data exactly. This method is free of the overheads of backwards recoveries like checkpointing, and does not compromise mathematical convergence properties of the solver as restarting would do. We apply this recovery to three widely used Krylov subspace methods, CG, GMRES and BiCGStab, and their preconditioned versions. We implement our resilience techniques on CG considering scenarios from small (8 cores) to large (1024 cores) scales, and demonstrate very low overheads compared to state-of-the-art solutions. We deploy our recovery techniques either by overlapping them with algorithmic computations or by forcing them to be in the critical path of the application. A trade-off exists between both approaches depending on the error rate the solver is suffering. Under realistic error rates, overlapping decreases overheads from 5.37% down to 3.59% for a non-preconditioned CG on 8 cores.This work has been partially supported by the European Research Council under the European Union's 7th FP, ERC Advanced Grant 321253, and by the Spanish Ministry of Science and Innovation under grant TIN2012-34557. L. Jaulmes has been partially supported by the Spanish Ministry of Education, Culture and Sports under grant FPU2013/06982. M. Moreto has been partially supported by the Spanish Ministry of Economy and Competitiveness under Juan de la Cierva postdoctoral fellowship JCI-2012-15047. M. Casas has been partially supported by the Secretary for Universities and Research of the Ministry of Economy and Knowledge of the Government of Catalonia and the Co-fund programme of the Marie Curie Actions of the European Union's 7th FP (contract 2013 BP B 00243).Peer ReviewedPostprint (author's final draft

What Hispanic parents do to encourage or discourage 3-5 year old children to be active: a qualitative study

Conference Theme: Promoting healthy eating and physical activity: The latest international researchabstract no. O8.2PURPOSE: Hispanic preschool children may be less active than Anglo-American children. Our aim was to identify what parents do to encourage or discourage physical activity (PA) among Hispanic preschool children. METHODS: Nominal Group Technique (NGT), a structured multi-step group procedure, was used to elicit and prioritize responses from 10 groups of Hispanic parents regarding what parents do to encourage (5 groups) or discourage (5 groups) preschool children to be active. Five groups consisted of parents with low education (less than high school) and 5 with high education (high school or greater) split between the two NGT questions: What …postprin

Evidence on anti-malarial and diagnostic markets in Cambodia to guide malaria elimination strategies and policies

BACKGROUND: Understanding Cambodia's anti-malarial and diagnostic landscape in 2015 is critical for informing and monitoring strategies and policies as Cambodia moves forward with national efforts to eliminate malaria. The aim of this paper is to present timely and key findings on the public and private sector anti-malarial and diagnostic landscape in Cambodia. This evidence can serve as a baseline benchmark for guiding implementation of national strategies as well as other regional initiatives to address malaria elimination activities. METHODS: From August 17th to October 1st, 2015, a cross sectional, nationally-representative malaria outlet survey was conducted in Cambodia. A census of all public and private outlets with potential to distribute malaria testing and/or treatment was conducted among 180 communes. An audit was completed for all anti-malarials, malaria rapid diagnostic tests (RDT) and microscopy. RESULTS: A total of 26,664 outlets were screened, and 1303 outlets were eligible and interviewed. Among all screened outlets in the public sector, 75.9% of public health facilities and 67.7% of community health workers stocked both malaria diagnostic testing and a first-line artemisinin-based combination therapy (ACT). Among anti-malarial-stocking private sector outlets, 64.7% had malaria blood testing available, and 70.9% were stocking a first-line ACT. Market share data illustrate that most of the anti-malarials were sold or distributed through the private sector (58.4%), including itinerant drug vendors (23.4%). First-line ACT accounted for the majority of the market share across the public and private sectors (90.3%). Among private sector outlets stocking any anti-malarial, the proportion of outlets with a first-line ACT or RDT was higher among outlets that had reportedly received one or more forms of 'support' (e.g. reportedly received training in the previous year on malaria diagnosis [RDT and/or microscopy] and/or the national treatment guidelines for malaria) compared to outlets that did not report receiving any support (ACT: 82.1 and 60.6%, respectively; RDT: 78.2 and 64.0%, respectively). CONCLUSION: The results point to high availability and distribution of first-line ACT and widespread availability of malaria diagnosis, especially in the public sector. This suggests that there is a strong foundation for achieving elimination goals in Cambodia. However, key gaps in terms of availability of malaria commodities for case management must be addressed, particularly in the private sector where most people seek treatment. Continued engagement with the private sector will be important to ensure accelerated progress towards malaria elimination

Whirling Hexagons and Defect Chaos in Hexagonal Non-Boussinesq Convection

We study hexagon patterns in non-Boussinesq convection of a thin rotating layer of water. For realistic parameters and boundary conditions we identify various linear instabilities of the pattern. We focus on the dynamics arising from an oscillatory side-band instability that leads to a spatially disordered chaotic state characterized by oscillating (whirling) hexagons. Using triangulation we obtain the distribution functions for the number of pentagonal and heptagonal convection cells. In contrast to the results found for defect chaos in the complex Ginzburg-Landau equation and in inclined-layer convection, the distribution functions can show deviations from a squared Poisson distribution that suggest non-trivial correlations between the defects.Comment: 4 mpg-movies are available at http://www.esam.northwestern.edu/~riecke/lit/lit.html submitted to New J. Physic

Variational tetrahedral meshing

In this paper, a novel Delaunay-based variational approach to isotropic tetrahedral meshing is presented. To achieve both robustness and efficiency, we minimize a simple mesh-dependent energy through global updates of both vertex positions and connectivity. As this energy is known to be the ∠1 distance between an isotropic quadratic function and its linear interpolation on the mesh, our minimization procedure generates well-shaped tetrahedra. Mesh design is controlled through a gradation smoothness parameter and selection of the desired number of vertices. We provide the foundations of our approach by explaining both the underlying variational principle and its geometric interpretation. We demonstrate the quality of the resulting meshes through a series of examples