Showing 1–33 of 33 results for author: Landry, W

Complex contagions can outperform simple contagions for network reconstruction with dense networks or saturated dynamics

Authors: Nicholas W. Landry, William Thompson, Laurent Hébert-Dufresne, Jean-Gabriel Young

Abstract: Network scientists often use complex dynamic processes to describe network contagions, but tools for fitting contagion models typically assume simple dynamics. Here, we address this gap by developing a nonparametric method to reconstruct a network and dynamics from a series of node states, using a model that breaks the dichotomy between simple pairwise and complex neighborhood-based contagions. We… ▽ More Network scientists often use complex dynamic processes to describe network contagions, but tools for fitting contagion models typically assume simple dynamics. Here, we address this gap by developing a nonparametric method to reconstruct a network and dynamics from a series of node states, using a model that breaks the dichotomy between simple pairwise and complex neighborhood-based contagions. We then show that a network is more easily reconstructed when observed through the lens of complex contagions if it is dense or the dynamic saturates, and that simple contagions are better otherwise. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 8 pages, 5 figures

The simpliciality of higher-order networks

Authors: Nicholas W. Landry, Jean-Gabriel Young, Nicole Eikmeier

Abstract: Higher-order networks are widely used to describe complex systems in which interactions can involve more than two entities at once. In this paper, we focus on inclusion within higher-order networks, referring to situations where specific entities participate in an interaction, and subsets of those entities also interact with each other. Traditional modeling approaches to higher-order networks tend… ▽ More Higher-order networks are widely used to describe complex systems in which interactions can involve more than two entities at once. In this paper, we focus on inclusion within higher-order networks, referring to situations where specific entities participate in an interaction, and subsets of those entities also interact with each other. Traditional modeling approaches to higher-order networks tend to either not consider inclusion at all (e.g., hypergraph models) or explicitly assume perfect and complete inclusion (e.g., simplicial complex models). To allow for a more nuanced assessment of inclusion in higher-order networks, we introduce the concept of "simpliciality" and several corresponding measures. Contrary to current modeling practice, we show that empirically observed systems rarely lie at either end of the simpliciality spectrum. In addition, we show that generative models fitted to these datasets struggle to capture their inclusion structure. These findings suggest new modeling directions for the field of higher-order network science. △ Less

Submitted 7 March, 2024; v1 submitted 26 August, 2023; originally announced August 2023.

Comments: Main text: 13 pages, 3 figures. Appendix: 8 pages, 1 figure

Journal ref: EPJ Data Sci. 13, 17 (2024)

Filtering higher-order datasets

Authors: Nicholas W. Landry, Ilya Amburg, Mirah Shi, Sinan G. Aksoy

Abstract: Many complex systems often contain interactions between more than two nodes, known as higher-order interactions, which can change the structure of these systems in significant ways. Researchers often assume that all interactions paint a consistent picture of a higher-order dataset's structure. In contrast, the connection patterns of individuals or entities in empirical systems are often stratified… ▽ More Many complex systems often contain interactions between more than two nodes, known as higher-order interactions, which can change the structure of these systems in significant ways. Researchers often assume that all interactions paint a consistent picture of a higher-order dataset's structure. In contrast, the connection patterns of individuals or entities in empirical systems are often stratified by interaction size. Ignoring this fact can aggregate connection patterns that exist only at certain scales of interaction. To isolate these scale-dependent patterns, we present an approach for analyzing higher-order datasets by filtering interactions by their size. We apply this framework to several empirical datasets from three domains to demonstrate that data practitioners can gain valuable information from this approach. △ Less

Submitted 1 November, 2023; v1 submitted 11 May, 2023; originally announced May 2023.

Comments: 22 pages, 17 figures

Opinion disparity in hypergraphs with community structure

Authors: Nicholas W. Landry, Juan G. Restrepo

Abstract: The division of a social group into subgroups with opposing opinions, which we refer to as opinion disparity, is a prevalent phenomenon in society. This phenomenon has been modeled by including mechanisms such as opinion homophily, bounded confidence interactions, and social reinforcement mechanisms. In this paper we study a complementary mechanism for the formation of opinion disparity based on h… ▽ More The division of a social group into subgroups with opposing opinions, which we refer to as opinion disparity, is a prevalent phenomenon in society. This phenomenon has been modeled by including mechanisms such as opinion homophily, bounded confidence interactions, and social reinforcement mechanisms. In this paper we study a complementary mechanism for the formation of opinion disparity based on higher-order interactions, i.e., simultaneous interactions between multiple agents. We present an extension of the planted partition model for uniform hypergraphs as a simple model of community structure and consider the hypergraph SIS model on a hypergraph with two communities where the binary ideology can spread via links (pairwise interactions) and triangles (three-way interactions). We approximate this contagion process with a mean-field model and find that for strong enough community structure, the two communities can hold very different average opinions. We determine the regimes of structural and infectious parameters for which this opinion disparity can exist and find that the existence of these disparities is much more sensitive to the triangle community structure than to the link community structure. We show that the existence and type of opinion disparities are extremely sensitive to differences in the sizes of the two communities. △ Less

Submitted 24 August, 2023; v1 submitted 27 February, 2023; originally announced February 2023.

Comments: 14 pages, 8 figures

On limitations of uniplex networks for modeling multiplex contagion

Authors: Nicholas W. Landry, jimi adams

Abstract: Many network contagion processes are inherently multiplex in nature, yet are often reduced to processes on uniplex networks in analytic practice. We therefore examine how data modeling choices can affect the predictions of contagion processes. We demonstrate that multiplex contagion processes are not simply the union of contagion processes over their constituent uniplex networks. We use multiplex… ▽ More Many network contagion processes are inherently multiplex in nature, yet are often reduced to processes on uniplex networks in analytic practice. We therefore examine how data modeling choices can affect the predictions of contagion processes. We demonstrate that multiplex contagion processes are not simply the union of contagion processes over their constituent uniplex networks. We use multiplex network data from two different contexts -- (1) a behavioral network to represent their potential for infectious disease transmission using a "simple" epidemiological model, and (2) users from online social network sites to represent their potential for information spread using a threshold-based "complex" contagion process. Our results show that contagion on multiplex data is not captured accurately in models developed from the uniplex networks even when they are combined, and that the nature of the differences between the (combined) uniplex and multiplex results depends on the specific spreading process over these networks. △ Less

Submitted 20 January, 2023; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: Main text: 12 pages, 3 figures; Supporting information: 5 pages, 5 figures

Hypergraph assortativity: a dynamical systems perspective

Authors: Nicholas W. Landry, Juan G. Restrepo

Abstract: The largest eigenvalue of the matrix describing a network's contact structure is often important in predicting the behavior of dynamical processes. We extend this notion to hypergraphs and motivate the importance of an analogous eigenvalue, the expansion eigenvalue, for hypergraph dynamical processes. Using a mean-field approach, we derive an approximation to the expansion eigenvalue in terms of t… ▽ More The largest eigenvalue of the matrix describing a network's contact structure is often important in predicting the behavior of dynamical processes. We extend this notion to hypergraphs and motivate the importance of an analogous eigenvalue, the expansion eigenvalue, for hypergraph dynamical processes. Using a mean-field approach, we derive an approximation to the expansion eigenvalue in terms of the degree sequence for uncorrelated hypergraphs. We introduce a generative model for hypergraphs that includes degree assortativity, and use a perturbation approach to derive an approximation to the expansion eigenvalue for assortative hypergraphs. We define the dynamical assortativity, a dynamically sensible definition of assortativity for uniform hypergraphs, and describe how reducing the dynamical assortativity of hypergraphs through preferential rewiring can extinguish epidemics. We validate our results with both synthetic and empirical datasets. △ Less

Submitted 10 May, 2022; v1 submitted 2 September, 2021; originally announced September 2021.

Comments: 11 pages, 4 figures

Journal ref: Chaos 32, 053113 (2022)

Bootstrapping Heisenberg Magnets and their Cubic Instability

Authors: Shai M. Chester, Walter Landry, Junyu Liu, David Poland, David Simmons-Duffin, Ning Su, Alessandro Vichi

Abstract: We study the critical $O(3)$ model using the numerical conformal bootstrap. In particular, we use a recently developed cutting-surface algorithm to efficiently map out the allowed space of CFT data from correlators involving the leading $O(3)$ singlet $s$, vector $φふぁい$, and rank-2 symmetric tensor $t$. We determine their scaling dimensions to be… ▽ More We study the critical $O(3)$ model using the numerical conformal bootstrap. In particular, we use a recently developed cutting-surface algorithm to efficiently map out the allowed space of CFT data from correlators involving the leading $O(3)$ singlet $s$, vector $φふぁい$, and rank-2 symmetric tensor $t$. We determine their scaling dimensions to be $(Δでるた_{s}, Δでるた_φふぁい, Δでるた_{t}) = (0.518942(51), 1.59489(59), 1.20954(23))$, and also bound various OPE coefficients. We additionally introduce a new "tip-finding" algorithm to compute an upper bound on the leading rank-4 symmetric tensor $t_4$, which we find to be relevant with $Δでるた_{t_4} < 2.99056$. The conformal bootstrap thus provides a numerical proof that systems described by the critical $O(3)$ model, such as classical Heisenberg ferromagnets at the Curie transition, are unstable to cubic anisotropy. △ Less

Submitted 19 November, 2021; v1 submitted 30 November, 2020; originally announced November 2020.

Comments: 35 pages, 9 figures, 8 tables. v2: clearing typos

Report number: CALT-TH-2020-053

Journal ref: Phys. Rev. D 104, 105013, Editors' Suggestion, (2021)

blocks_3d: Software for general 3d conformal blocks

Authors: Rajeev S. Erramilli, Luca V. Iliesiu, Petr Kravchuk, Walter Landry, David Poland, David Simmons-Duffin

Abstract: We introduce the software blocks_3d for computing four-point conformal blocks of operators with arbitrary Lorentz representations in 3d CFTs. It uses Zamolodchikov-like recursion relations to numerically compute derivatives of blocks around a crossing-symmetric configuration. It is implemented as a heavily optimized, multithreaded, C++ application. We give performance benchmarks for correlators co… ▽ More We introduce the software blocks_3d for computing four-point conformal blocks of operators with arbitrary Lorentz representations in 3d CFTs. It uses Zamolodchikov-like recursion relations to numerically compute derivatives of blocks around a crossing-symmetric configuration. It is implemented as a heavily optimized, multithreaded, C++ application. We give performance benchmarks for correlators containing scalars, fermions, and stress tensors. As an example application, we recompute bootstrap bounds on four-point functions of fermions and study whether a previously observed sharp jump can be explained using the "fake primary" effect. We conclude that the fake primary effect cannot fully explain the jump and the possible existence of a "dead-end" CFT near the jump merits further study. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: 33 pages + appendices

Report number: CALT-TH 2020-048

Carving out OPE space and precise $O(2)$ model critical exponents

Authors: Shai M. Chester, Walter Landry, Junyu Liu, David Poland, David Simmons-Duffin, Ning Su, Alessandro Vichi

Abstract: We develop new tools for isolating CFTs using the numerical bootstrap. A "cutting surface" algorithm for scanning OPE coefficients makes it possible to find islands in high-dimensional spaces. Together with recent progress in large-scale semidefinite programming, this enables bootstrap studies of much larger systems of correlation functions than was previously practical. We apply these methods to… ▽ More We develop new tools for isolating CFTs using the numerical bootstrap. A "cutting surface" algorithm for scanning OPE coefficients makes it possible to find islands in high-dimensional spaces. Together with recent progress in large-scale semidefinite programming, this enables bootstrap studies of much larger systems of correlation functions than was previously practical. We apply these methods to correlation functions of charge-0, 1, and 2 scalars in the 3d $O(2)$ model, computing new precise values for scaling dimensions and OPE coefficients in this theory. Our new determinations of scaling dimensions are consistent with and improve upon existing Monte Carlo simulations, sharpening the existing decades-old $8σしぐま$ discrepancy between theory and experiment. △ Less

Submitted 3 July, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Comments: 56 pages, 12 figures, 10 tables, 1 algorithm; V2: small corrections to match JHEP version

Report number: CALT-TH-2019-051

Journal ref: JHEP 06 (2020) 142

Scaling the semidefinite program solver SDPB

Authors: Walter Landry, David Simmons-Duffin

Abstract: We present enhancements to SDPB, an open source, parallelized, arbitrary precision semidefinite program solver designed for the conformal bootstrap. The main enhancement is significantly improved performance and scalability using the Elemental library and MPI. The result is a new version of SDPB that runs on multiple nodes with hundreds of cores with excellent scaling, making it practical to solve… ▽ More We present enhancements to SDPB, an open source, parallelized, arbitrary precision semidefinite program solver designed for the conformal bootstrap. The main enhancement is significantly improved performance and scalability using the Elemental library and MPI. The result is a new version of SDPB that runs on multiple nodes with hundreds of cores with excellent scaling, making it practical to solve larger problems. We demonstrate performance on a moderate-size problem in the 3d Ising CFT and a much larger problem in the $O(2)$ Model. △ Less

Submitted 20 September, 2019; originally announced September 2019.

Comments: 13 pages plus references, 2 figures

The Zwicky Transient Facility: Science Objectives

Authors: Matthew J. Graham, S. R. Kulkarni, Eric C. Bellm, Scott M. Adams, Cristina Barbarino, Nadejda Blagorodnova, Dennis Bodewits, Bryce Bolin, Patrick R. Brady, S. Bradley Cenko, Chan-Kao Chang, Michael W. Coughlin, Kishalay De, Gwendolyn Eadie, Tony L. Farnham, Ulrich Feindt, Anna Franckowiak, Christoffer Fremling, Avishay Gal-yam, Suvi Gezari, Shaon Ghosh, Daniel A. Goldstein, V. Zach Golkhou, Ariel Goobar, Anna Y. Q. Ho , et al. (92 additional authors not shown)

Abstract: The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single… ▽ More The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: 26 pages, 7 figures, Published in PASP Focus Issue on the Zwicky Transient Facility

The Zwicky Transient Facility: System Overview, Performance, and First Results

Authors: Eric C. Bellm, Shrinivas R. Kulkarni, Matthew J. Graham, Richard Dekany, Roger M. Smith, Reed Riddle, Frank J. Masci, George Helou, Thomas A. Prince, Scott M. Adams, C. Barbarino, Tom Barlow, James Bauer, Ron Beck, Justin Belicki, Rahul Biswas, Nadejda Blagorodnova, Dennis Bodewits, Bryce Bolin, Valery Brinnel, Tim Brooke, Brian Bue, Mattia Bulla, Rick Burruss, S. Bradley Cenko , et al. (91 additional authors not shown)

Abstract: The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation… ▽ More The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: Published in PASP Focus Issue on the Zwicky Transient Facility (https://dx.doi.org/10.1088/1538-3873/aaecbe). 21 Pages, 12 Figures

Journal ref: Publications of the Astronomical Society of the Pacific, Volume 131, Issue 995, pp. 018002 (2019)

The Zwicky Transient Facility: Data Processing, Products, and Archive

Authors: Frank J. Masci, Russ R. Laher, Ben Rusholme, David L. Shupe, Steven Groom, Jason Surace, Edward Jackson, Serge Monkewitz, Ron Beck, David Flynn, Scott Terek, Walter Landry, Eugean Hacopians, Vandana Desai, Justin Howell, Tim Brooke, David Imel, Stefanie Wachter, Quan-Zhi Ye, Hsing-Wen Lin, S. Bradley Cenko, Virginia Cunningham, Umaa Rebbapragada, Brian Bue, Adam A. Miller , et al. (24 additional authors not shown)

Abstract: The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of ~3760 square degrees/hour to median depths of g ~ 20.8 and r ~ 20.6 mag (AB, 5sigma in 30 sec). We describe the Science Data System that is housed… ▽ More The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of ~3760 square degrees/hour to median depths of g ~ 20.8 and r ~ 20.6 mag (AB, 5sigma in 30 sec). We describe the Science Data System that is housed at IPAC, Caltech. This comprises the data-processing pipelines, alert production system, data archive, and user interfaces for accessing and analyzing the products. The realtime pipeline employs a novel image-differencing algorithm, optimized for the detection of point source transient events. These events are vetted for reliability using a machine-learned classifier and combined with contextual information to generate data-rich alert packets. The packets become available for distribution typically within 13 minutes (95th percentile) of observation. Detected events are also linked to generate candidate moving-object tracks using a novel algorithm. Objects that move fast enough to streak in the individual exposures are also extracted and vetted. The reconstructed astrometric accuracy per science image with respect to Gaia is typically 45 to 85 milliarcsec. This is the RMS per axis on the sky for sources extracted with photometric S/N >= 10. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. Photometric calibration accuracy with respect to Pan-STARRS1 is generally better than 2%. The products support a broad range of scientific applications: fast and young supernovae, rare flux transients, variable stars, eclipsing binaries, variability from active galactic nuclei, counterparts to gravitational wave sources, a more complete census of Type Ia supernovae, and Solar System objects. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: 30 pages, 16 figures, Published in PASP Focus Issue on the Zwicky Transient Facility (doi: 10.1088/1538-3873/aae8ac)

Fast, accurate solutions for curvilinear earthquake faults and anelastic strain

Authors: Walter Landry, Sylvain Barbot

Abstract: Imaging the anelastic deformation within the crust and lithosphere using surface geophysical data remains a significant challenge in part due to the wide range of physical processes operating at different depths and to various levels of localization that they embody. Models of Earth's elastic properties from seismological imaging combined with geodetic modeling may form the basis of comprehensive… ▽ More Imaging the anelastic deformation within the crust and lithosphere using surface geophysical data remains a significant challenge in part due to the wide range of physical processes operating at different depths and to various levels of localization that they embody. Models of Earth's elastic properties from seismological imaging combined with geodetic modeling may form the basis of comprehensive rheological models of Earth's interior. However, representing the structural complexity of faults and shear zones in numerical models of deformation still constitutes a major difficulty. Here, we present numerical techniques for high-precision models of deformation and stress around both curvilinear faults and volumes undergoing anelastic (irreversible) strain in a heterogenous elastic half-space. To that end, we enhance the software Gamra to model triangular and rectangular fault patches and tetrahedral and cuboidal strain volumes. This affords a means of rapid and accurate calculations of elasto-static Green's functions for localized (e.g., faulting) and distributed (e.g., viscoelastic) deformation in Earth's crust and lithosphere. We demonstrate the correctness of the method with analytic tests, and we illustrate its practical performance by solving for coseismic and postseismic deformation following the 2015 Mw 7.8 Gorkha, Nepal earthquake to extremely high precision. △ Less

Submitted 24 February, 2018; originally announced February 2018.

Processing Images from the Zwicky Transient Facility

Authors: Russ R. Laher, Frank J. Masci, Steve Groom, Benjamin Rusholme, David L. Shupe, Ed Jackson, Jason Surace, Dave Flynn, Walter Landry, Scott Terek, George Helou, Ron Beck, Eugean Hacopians, Umaa Rebbapragada, Brian Bue, Roger M. Smith, Richard G. Dekany, Adam A. Miller, S. B. Cenko, Eric Bellm, Maria Patterson, Thomas Kupfer, Lin Yan, Tom Barlow, Matthew Graham , et al. (3 additional authors not shown)

Abstract: The Zwicky Transient Facility is a new robotic-observing program, in which a newly engineered 600-MP digital camera with a pioneeringly large field of view, 47~square degrees, will be installed into the 48-inch Samuel Oschin Telescope at the Palomar Observatory. The camera will generate $\sim 1$~petabyte of raw image data over three years of operations. In parallel related work, new hardware and s… ▽ More The Zwicky Transient Facility is a new robotic-observing program, in which a newly engineered 600-MP digital camera with a pioneeringly large field of view, 47~square degrees, will be installed into the 48-inch Samuel Oschin Telescope at the Palomar Observatory. The camera will generate $\sim 1$~petabyte of raw image data over three years of operations. In parallel related work, new hardware and software systems are being developed to process these data in real time and build a long-term archive for the processed products. The first public release of archived products is planned for early 2019, which will include processed images and astronomical-source catalogs of the northern sky in the $g$ and $r$ bands. Source catalogs based on two different methods will be generated for the archive: aperture photometry and point-spread-function fitting. △ Less

Submitted 16 October, 2017; v1 submitted 4 August, 2017; originally announced August 2017.

Comments: 6 pages, 4 figures, submitted to RTSRE Proceedings (www.rtsre.org)

Gamra: Simple Meshes for Complex Earthquakes

Authors: Walter Landry, Sylvain Barbot

Abstract: The static offsets caused by earthquakes are well described by elastostatic models with a discontinuity in the displacement along the fault. A traditional approach to model this discontinuity is to align the numerical mesh with the fault and solve the equations using finite elements. However, this distorted mesh can be difficult to generate and update. We present a new numerical method, inspired b… ▽ More The static offsets caused by earthquakes are well described by elastostatic models with a discontinuity in the displacement along the fault. A traditional approach to model this discontinuity is to align the numerical mesh with the fault and solve the equations using finite elements. However, this distorted mesh can be difficult to generate and update. We present a new numerical method, inspired by the Immersed Interface Method, for solving the elastostatic equations with embedded discontinuities. This method has been carefully designed so that it can be used on parallel machines on an adapted finite difference grid. We have implemented this method in Gamra, a new code for earth modelling. We demonstrate the correctness of the method with analytic tests, and we demonstrate its practical performance by solving a realistic earthquake model to extremely high precision. △ Less

Submitted 24 February, 2018; v1 submitted 3 March, 2016; originally announced March 2016.

Comments: 38 pages, 20 figures

Journal ref: Computers and Geosciences Volume 90, Part A, May 2016, Pages 49-63

Learning from FITS: Limitations in use in modern astronomical research

Authors: Brian Thomas, Tim Jenness, Frossie Economou, Perry Greenfield, Paul Hirst, David S. Berry, Erik Bray, Norman Gray, Demitri Muna, James Turner, Miguel de Val-Borro, Juande Santander-Vela, David Shupe, John Good, G. Bruce Berriman, Slava Kitaeff, Jonathan Fay, Omar Laurino, Anastasia Alexov, Walter Landry, Joe Masters, Adam Brazier, Reinhold Schaaf, Kevin Edwards, Russell O. Redman , et al. (13 additional authors not shown)

Abstract: The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard, however, is showing its age. Developed in the late 1970s, the FITS authors made a number of implementation choices that, while common at the time, are now seen to limit its utility with modern data.… ▽ More The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard, however, is showing its age. Developed in the late 1970s, the FITS authors made a number of implementation choices that, while common at the time, are now seen to limit its utility with modern data. The authors of the FITS standard could not anticipate the challenges which we are facing today in astronomical computing. Difficulties we now face include, but are not limited to, addressing the need to handle an expanded range of specialized data product types (data models), being more conducive to the networked exchange and storage of data, handling very large datasets, and capturing significantly more complex metadata and data relationships. There are members of the community today who find some or all of these limitations unworkable, and have decided to move ahead with storing data in other formats. If this fragmentation continues, we risk abandoning the advantages of broad interoperability, and ready archivability, that the FITS format provides for astronomy. In this paper we detail some selected important problems which exist within the FITS standard today. These problems may provide insight into deeper underlying issues which reside in the format and we provide a discussion of some lessons learned. It is not our intention here to prescribe specific remedies to these issues; rather, it is to call attention of the FITS and greater astronomical computing communities to these problems in the hope that it will spur action to address them. △ Less

Submitted 10 February, 2015; v1 submitted 3 February, 2015; originally announced February 2015.

Analysis of Granular Flow in a Pebble-Bed Nuclear Reactor

Authors: Chris H. Rycroft, Gary S. Grest, James W. Landry, Martin Z. Bazant

Abstract: Pebble-bed nuclear reactor technology, which is currently being revived around the world, raises fundamental questions about dense granular flow in silos. A typical reactor core is composed of graphite fuel pebbles, which drain very slowly in a continuous refueling process. Pebble flow is poorly understood and not easily accessible to experiments, and yet it has a major impact on reactor physics… ▽ More Pebble-bed nuclear reactor technology, which is currently being revived around the world, raises fundamental questions about dense granular flow in silos. A typical reactor core is composed of graphite fuel pebbles, which drain very slowly in a continuous refueling process. Pebble flow is poorly understood and not easily accessible to experiments, and yet it has a major impact on reactor physics. To address this problem, we perform full-scale, discrete-element simulations in realistic geometries, with up to 440,000 frictional, viscoelastic 6cm-diameter spheres draining in a cylindrical vessel of diameter 3.5m and height 10m with bottom funnels angled at 30 degrees or 60 degrees. We also simulate a bidisperse core with a dynamic central column of smaller graphite moderator pebbles and show that little mixing occurs down to a 1:2 diameter ratio. We analyze the mean velocity, diffusion and mixing, local ordering and porosity (from Voronoi volumes), the residence-time distribution, and the effects of wall friction and discuss implications for reactor design and the basic physics of granular flow. △ Less

Submitted 16 February, 2006; originally announced February 2006.

Comments: 18 pages, 21 figures

Dynamics of Random Packings in Granular Flow

Authors: Chris H. Rycroft, Martin Z. Bazant, Gary S. Grest, James W. Landry

Abstract: We present a multiscale simulation algorithm for amorphous materials, which we illustrate and validate in a canonical case of dense granular flow. Our algorithm is based on the recently proposed Spot Model, where particles in a dense random packing undergo chain-like collective displacements in response to diffusing "spots" of influence, carrying a slight excess of interstitial free volume. We r… ▽ More We present a multiscale simulation algorithm for amorphous materials, which we illustrate and validate in a canonical case of dense granular flow. Our algorithm is based on the recently proposed Spot Model, where particles in a dense random packing undergo chain-like collective displacements in response to diffusing "spots" of influence, carrying a slight excess of interstitial free volume. We reconstruct the microscopic dynamics of particles from the "coarse grained" dynamics of spots by introducing a localized particle relaxation step after each spot-induced block displacement, simply to enforce packing constraints with a (fairly arbitrary) soft-core repulsion. To test the model, we study to what extent it can describe the dynamics of up to 135,000 frictional, viscoelastic spheres in granular drainage simulated by the discrete-element method (DEM). With only five fitting parameters (the radius, volume, diffusivity, drift velocity, and injection rate of spots), we find that the spot simulations are able to largely reproduce not only the mean flow and diffusion, but also some subtle statistics of the flowing packings, such as spatial velocity correlations and many-body structural correlations. The spot simulations run over 100 times faster than DEM and demonstrate the possibility of multiscale modeling for amorphous materials, whenever a suitable model can be devised for the coarse-grained spot dynamics. △ Less

Submitted 16 February, 2006; originally announced February 2006.

Comments: 8 pages, 7 figures

Plug flow and the breakdown of Bagnold scaling in cohesive granular flows

Authors: Robert Brewster, Gary S. Grest, James W. Landry, Alex J. Levine

Abstract: Cohesive granular media flowing down an inclined plane are studied by discrete element simulations. Previous work on cohesionless granular media demonstrated that within the steady flow regime where gravitational energy is balanced by dissipation arising from intergrain forces, the velocity profile in the flow direction scales with depth in a manner consistent with the predictions of Bagnold. He… ▽ More Cohesive granular media flowing down an inclined plane are studied by discrete element simulations. Previous work on cohesionless granular media demonstrated that within the steady flow regime where gravitational energy is balanced by dissipation arising from intergrain forces, the velocity profile in the flow direction scales with depth in a manner consistent with the predictions of Bagnold. Here we demonstrate that this Bagnold scaling does not hold for the analogous steady-flows in cohesive granular media. We develop a generalization of the Bagnold constitutive relation to account for our observation and speculate as to the underlying physical mechanisms responsible for the different constitutive laws for cohesive and noncohesive granular media. △ Less

Submitted 5 September, 2005; originally announced September 2005.

Comments: 8 pages, 10 figures

Jamming in granular hopper flow

Authors: James W. Landry, Gary S. Grest

Abstract: Large-scale three dimensional molecular dynamics simulations of hopper flow are presented. The flow rate of the system is controlled by the width of the aperture at the bottom. As the steady-state flow rate is reduced, the force distribution $P(f)$ changes only slightly, while there is a large change in the impulse distribution $P(i)$. In both cases, the distributions show an increase in small f… ▽ More Large-scale three dimensional molecular dynamics simulations of hopper flow are presented. The flow rate of the system is controlled by the width of the aperture at the bottom. As the steady-state flow rate is reduced, the force distribution $P(f)$ changes only slightly, while there is a large change in the impulse distribution $P(i)$. In both cases, the distributions show an increase in small forces or impulses as the systems approach jamming, the opposite of that seen in previous Lennard-Jones simulations. This occurs dynamically as well for a hopper that transitions from a flowing to a jammed state over time. The final jammed $P(f)$ is quite distinct from a poured packing $P(f)$ in the same geometry. The change in $P(i)$ is a much stronger indicator of the approach to jamming. The formation of a peak or plateau in $P(f)$ at the average force is not a general feature of the approach to jamming. △ Less

Submitted 4 June, 2005; originally announced June 2005.

Comments: 5 pages, 6 figures, submitted to PRE

Granular packings with moving side walls

Authors: J. W. Landry, G. S. Grest

Abstract: The effects of movement of the side walls of a confined granular packing are studied by discrete element, molecular dynamics simulations. The dynamical evolution of the stress is studied as a function of wall movement both in the direction of gravity as well as opposite to it. For all wall velocities explored, the stress in the final state of the system after wall movement is fundamentally diffe… ▽ More The effects of movement of the side walls of a confined granular packing are studied by discrete element, molecular dynamics simulations. The dynamical evolution of the stress is studied as a function of wall movement both in the direction of gravity as well as opposite to it. For all wall velocities explored, the stress in the final state of the system after wall movement is fundamentally different from the original state obtained by pouring particles into the container and letting them settle under the influence of gravity. The original packing possesses a hydrostatic-like region at the top of the container which crosses over to a depth-independent stress. As the walls are moved in the direction opposite to gravity, the saturation stress first reaches a minimum value independent of the wall velocity, then increases to a steady-state value dependent on the wall-velocity. After wall movement ceases and the packing reaches equilibrium, the stress profile fits the classic Janssen form for high wall velocities, while it has some deviations for low wall velocities. The wall movement greatly increases the number of particle-wall and particle-particle forces at the Coulomb criterion. Varying the wall velocity has only small effects on the particle structure of the final packing so long as the walls travel a similar distance. △ Less

Submitted 19 September, 2003; originally announced September 2003.

Comments: 11 pages, 10 figures, some figures in color

Discrete element simulations of stress distributions in silos: crossover from two to three dimensions

Authors: James W. Landry, Gary S. Grest, Stephen J. Plimpton

Abstract: The transition from two-dimensional (2D) to three-dimensional (3D) granular packings is studied using large-scale discrete element computer simulations. We focus on vertical stress profiles and examine how they change with dimensionality from 2D to 3D. We compare results for packings in 2D, quasi-2D packings between flat plates, and 3D packings. Analysis of these packings suggests that the Janss… ▽ More The transition from two-dimensional (2D) to three-dimensional (3D) granular packings is studied using large-scale discrete element computer simulations. We focus on vertical stress profiles and examine how they change with dimensionality from 2D to 3D. We compare results for packings in 2D, quasi-2D packings between flat plates, and 3D packings. Analysis of these packings suggests that the Janssen theory does not fully describe these packings, especially at the top of the piles, where a hydrostatic-like region of vertical stress is visible in all cases. We find that the interior of the packing is far from incipient failure, while in general, the forces at the walls are close to incipient failure. △ Less

Submitted 4 November, 2003; v1 submitted 6 February, 2003; originally announced February 2003.

Comments: 9 pages, 8 figures v2: some reorganization, sections moved around, discussion made clearer

Dynamics of a Complex Quantum Magnet

Authors: J. W. Landry, S. N. Coppersmith

Abstract: We have computed the low energy quantum states and low frequency dynamical susceptibility of complex quantum spin systems in the limit of strong interactions, obtaining exact results for system sizes enormously larger than accessible previously. The ground state is a complex superposition of a substantial fraction of all the classical ground states, and yet the dynamical susceptibility exhibits… ▽ More We have computed the low energy quantum states and low frequency dynamical susceptibility of complex quantum spin systems in the limit of strong interactions, obtaining exact results for system sizes enormously larger than accessible previously. The ground state is a complex superposition of a substantial fraction of all the classical ground states, and yet the dynamical susceptibility exhibits sharp resonances reminiscent of the behavior of single spins. These results show that strongly interacting quantum systems can organize to generate coherent excitations and shed light on recent experiments demonstrating that coherent excitations are present in a disordered spin liquid. The dependence of the energy spectra on system size differs qualitatively from that of the energy spectra of random undirected bipartite graphs with similar statistics, implying that strong interactions are giving rise to these unusual spectral properties. △ Less

Submitted 14 January, 2003; originally announced January 2003.

Comments: 4 pages, 4 figures, submitted to PRL

Confined granular packings: structure, stress, and forces

Authors: J. W. Landry, G. S. Grest, L. E. Silbert, S. J. Plimpton

Abstract: The structure and stresses of static granular packs in cylindrical containers are studied using large-scale discrete element molecular dynamics simulations in three dimensions. We generate packings by both pouring and sedimentation and examine how the final state depends on the method of construction. The vertical stress becomes depth-independent for deep piles and we compare these stress depth-… ▽ More The structure and stresses of static granular packs in cylindrical containers are studied using large-scale discrete element molecular dynamics simulations in three dimensions. We generate packings by both pouring and sedimentation and examine how the final state depends on the method of construction. The vertical stress becomes depth-independent for deep piles and we compare these stress depth-profiles to the classical Janssen theory. The majority of the tangential forces for particle-wall contacts are found to be close to the Coulomb failure criterion, in agreement with the theory of Janssen, while particle-particle contacts in the bulk are far from the Coulomb criterion. In addition, we show that a linear hydrostatic-like region at the top of the packings unexplained by the Janssen theory arises because most of the particle-wall tangential forces in this region are far from the Coulomb yield criterion. The distributions of particle-particle and particle-wall contact forces $P(f)$ exhibit exponential-like decay at large forces in agreement with previous studies. △ Less

Submitted 31 January, 2003; v1 submitted 11 November, 2002; originally announced November 2002.

Comments: 11 pages, 11 figures, submitted to PRE (v2) added new references, fixed typos

Statistics of the contact network in frictional and frictionless granular packings

Authors: Leonardo E. Silbert, Gary S. Grest, James W. Landry

Abstract: Simulated granular packings with different particle friction coefficient mu are examined. The distribution of the particle-particle and particle-wall normal and tangential contact forces P(f) are computed and compared with existing experimental data. Here f equivalent to F/F-bar is the contact force F normalized by the average value F-bar. P(f) exhibits exponential-like decay at large forces, a… ▽ More Simulated granular packings with different particle friction coefficient mu are examined. The distribution of the particle-particle and particle-wall normal and tangential contact forces P(f) are computed and compared with existing experimental data. Here f equivalent to F/F-bar is the contact force F normalized by the average value F-bar. P(f) exhibits exponential-like decay at large forces, a plateau/peak near f = 1, with additional features at forces smaller than the average that depend on mu. Computations of the force-force spatial distribution function and the contact point radial distribution function indicate that correlations between forces are only weakly dependent on friction and decay rapidly beyond approximately three particle diameters. Distributions of the particle-particle contact angles show that the contact network is not isotropic and only weakly dependent on friction. High force-bearing structures, or force chains, do not play a dominant role in these three dimensional, unloaded packings. △ Less

Submitted 23 August, 2002; originally announced August 2002.

Comments: 11 pages, 13 figures, submitted to PRE

Granular flow down a rough inclined plane: transition between thin and thick piles

Authors: Leonardo E. Silbert, James W. Landry, Gary S. Grest

Abstract: The rheology of granular particles in an inclined plane geometry is studied using molecular dynamics simulations. The flow--no-flow boundary is determined for piles of varying heights over a range of inclination angles $θしーた$. Three angles determine the phase diagram: $θしーた_{r}$, the angle of repose, is the angle at which a flowing system comes to rest; $θしーた_{m}$, the maximum angle of stability, is the… ▽ More The rheology of granular particles in an inclined plane geometry is studied using molecular dynamics simulations. The flow--no-flow boundary is determined for piles of varying heights over a range of inclination angles $θしーた$. Three angles determine the phase diagram: $θしーた_{r}$, the angle of repose, is the angle at which a flowing system comes to rest; $θしーた_{m}$, the maximum angle of stability, is the inclination required to induce flow in a static system; and $θしーた_{max}$ is the maximum angle for which stable, steady state flow is observed. In the stable flow region $θしーた_{r}<θしーた<θしーた_{max}$, three flow regimes can be distinguished that depend on how close $θしーた$ is to $θしーた_{r}$: i) $θしーた>>θしーた_{r}$: Bagnold rheology, characterized by a mean particle velocity $v_{x}$ in the direction of flow that scales as $v_{x}\propto h^{3/2}$, for a pile of height $h$, ii) $θしーた\gtrsimθしーた_{r}$: the slow flow regime, characterized by a linear velocity profile with depth, and iii) $θしーた\approxθしーた_{r}$: avalanche flow characterized by a slow underlying creep motion combined with occasional free surface events and large energy fluctuations. We also probe the physics of the initiation and cessation of flow. The results are compared to several recent experimental studies on chute flows and suggest that differences between measured velocity profiles in these experiments may simply be a consequence of how far the system is from jamming. △ Less

Submitted 11 June, 2002; originally announced June 2002.

Comments: 19 pages, 14 figs, submitted to Physics of Fluids

Radiation-balanced simulations for binary inspiral

Authors: John T Whelan, Christopher Beetle, Walter Landry, Richard H Price

Abstract: The late stage of the inspiral of two black holes may have important non-Newtonian effects that are unrelated to radiation reaction. To understand these effects we approximate a slowly inspiralling binary by a stationary solution to Einstein's equations in which the holes orbit eternally. Radiation reaction is nullified by specifying a boundary condition at infinity containing equal amounts of i… ▽ More The late stage of the inspiral of two black holes may have important non-Newtonian effects that are unrelated to radiation reaction. To understand these effects we approximate a slowly inspiralling binary by a stationary solution to Einstein's equations in which the holes orbit eternally. Radiation reaction is nullified by specifying a boundary condition at infinity containing equal amounts of ingoing and outgoing radiation. The computational problem is then converted from an evolution problem with initial data to a boundary value problem. In addition to providing an approximate inspiral waveform via extraction of the outgoing modes, our approximation can give alternative initial data for numerical relativity evolution. We report results on simplified models and on progress in building 3D numerical solutions. △ Less

Submitted 1 October, 2001; originally announced October 2001.

Comments: 7 pages, uses IOP class files, submitted to the proceedings of the 4th Amaldi meeting (which will be published in Classical and Quantum Gravity)

Journal ref: Class.Quant.Grav. 19 (2002) 1285-1290

Ground states of two-dimensional $\pm$J Edwards-Anderson spin glasses

Authors: J. W. Landry, S. N. Coppersmith

Abstract: We present an exact algorithm for finding all the ground states of the two-dimensional Edwards-Anderson $\pm J$ spin glass and characterize its performance. We investigate how the ground states change with increasing system size and and with increasing antiferromagnetic bond ratio $x$. We find that that some system properties have very large and strongly non-Gaussian variations between realizati… ▽ More We present an exact algorithm for finding all the ground states of the two-dimensional Edwards-Anderson $\pm J$ spin glass and characterize its performance. We investigate how the ground states change with increasing system size and and with increasing antiferromagnetic bond ratio $x$. We find that that some system properties have very large and strongly non-Gaussian variations between realizations. △ Less

Submitted 7 September, 2001; v1 submitted 7 September, 2001; originally announced September 2001.

Comments: 15 pages, 21 figures, 2 tables, uses revtex4 macros

Journal ref: Phys. Rev. B, 65, 134404 (2002)

An Efficient Method for Fully Relativistic Simulations of Coalescing Binary Neutron Stars

Authors: Walter Landry, Saul A. Teukolsky

Abstract: The merger of two neutron stars has been proposed as a source of gamma-ray bursts, r-process elements, and detectable gravitational waves. Extracting information from observations of these phenomena requires fully relativistic simulations. Unfortunately, the only demonstrated method for stably evolving neutron stars requires solving elliptic equations at each time step, adding substantially to t… ▽ More The merger of two neutron stars has been proposed as a source of gamma-ray bursts, r-process elements, and detectable gravitational waves. Extracting information from observations of these phenomena requires fully relativistic simulations. Unfortunately, the only demonstrated method for stably evolving neutron stars requires solving elliptic equations at each time step, adding substantially to the computational resources required. In this paper we present a simpler, more efficient method. The key insight is in how we apply numerical diffusion. We perform a number of tests to validate the method and our implementation. We also carry out a very rough simulation of coalescence and extraction of the gravitational waves to show that the method is viable if realistic initial data are provided. △ Less

Submitted 31 May, 2000; v1 submitted 1 December, 1999; originally announced December 1999.

Comments: 10 pages, 24 figures, fixed some typos, added a paragraph about the diffusion method and time stepping

Stable characteristic evolution of generic 3-dimensional single-black-hole spacetimes

Authors: The Binary Black Hole Grand Challenge Alliance, :, R. Gomez, L. Lehner, R. Marsa, J. Winicour, A. Abrahams, A. Anderson, P. Anninos, T. Baumgarte, N. Bishop, S. Brandt J. Browne, K. Camarda, M. Choptuik, R. Correl, G. Cook, C. Evans, L. Finn, G. Fox, T. Haupt, M. Huq, L. Kidder, S. Klasky, P. Laguna, W. Landry , et al. (20 additional authors not shown)

Abstract: We report new results which establish that the accurate 3-dimensional numerical simulation of generic single-black-hole spacetimes has been achieved by characteristic evolution with unlimited long term stability. Our results cover a selection of distorted, moving and spinning single black holes, with evolution times up to 60,000M. We report new results which establish that the accurate 3-dimensional numerical simulation of generic single-black-hole spacetimes has been achieved by characteristic evolution with unlimited long term stability. Our results cover a selection of distorted, moving and spinning single black holes, with evolution times up to 60,000M. △ Less

Submitted 20 January, 1998; originally announced January 1998.

Comments: 4 pages, 3 figures

Journal ref: Phys.Rev.Lett.80:3915-3918,1998

Boosted three-dimensional black-hole evolutions with singularity excision

Authors: The Binary Black Hole Grand Challenge Alliance, :, G. B. Cook, M. F. Huq, S. A. Klasky, M. A. Scheel, A. M. Abrahams, A. Anderson, P. Anninos, T. W. Baumgarte, N. T. Bishop, S. R. Brandt, J. C. Browne, K. Camarda, M. W. Choptuik, C. R. Evans, L. S. Finn, G. C. Fox, R. Gomez, T. Haupt, L. E. Kidder, P. Laguna, W. Landry, L. Lehner, J. Lenaghan , et al. (21 additional authors not shown)

Abstract: Binary black hole interactions provide potentially the strongest source of gravitational radiation for detectors currently under development. We present some results from the Binary Black Hole Grand Challenge Alliance three- dimensional Cauchy evolution module. These constitute essential steps towards modeling such interactions and predicting gravitational radiation waveforms. We report on singl… ▽ More Binary black hole interactions provide potentially the strongest source of gravitational radiation for detectors currently under development. We present some results from the Binary Black Hole Grand Challenge Alliance three- dimensional Cauchy evolution module. These constitute essential steps towards modeling such interactions and predicting gravitational radiation waveforms. We report on single black hole evolutions and the first successful demonstration of a black hole moving freely through a three-dimensional computational grid via a Cauchy evolution: a hole moving ~6M at 0.1c during a total evolution of duration ~60M. △ Less

Submitted 26 November, 1997; originally announced November 1997.

Journal ref: Phys.Rev.lett.80:2512-2516,1998

Gravitational wave extraction and outer boundary conditions by perturbative matching

Authors: The Binary Black Hole Grand Challenge Alliance, :, A. M. Abrahams, L. Rezzolla, M. E. Rupright, A. Anderson, P. Anninos, T. W. Baumgarte, N. T. Bishop, S. R. Brandt, J. C. Browne, K. Camarda, M. W. Choptuik, G. B. Cook, C. R. Evans, L. S. Finn, G. Fox, R. Gomez, T. Haupt, M. F. Huq, L. E. Kidder, S. Klasky, P. Laguna, W. Landry, L. Lehner , et al. (20 additional authors not shown)

Abstract: We present a method for extracting gravitational radiation from a three-dimensional numerical relativity simulation and, using the extracted data, to provide outer boundary conditions. The method treats dynamical gravitational variables as nonspherical perturbations of Schwarzschild geometry. We discuss a code which implements this method and present results of tests which have been performed wi… ▽ More We present a method for extracting gravitational radiation from a three-dimensional numerical relativity simulation and, using the extracted data, to provide outer boundary conditions. The method treats dynamical gravitational variables as nonspherical perturbations of Schwarzschild geometry. We discuss a code which implements this method and present results of tests which have been performed with a three dimensional numerical relativity code. △ Less

Submitted 30 September, 1997; originally announced September 1997.

Journal ref: Phys.Rev.Lett.80:1812-1815,1998