Showing 1–43 of 43 results for author: Peterson, L

Dynamical Control of Excitons in Atomically Thin Semiconductors

Authors: Eric L. Peterson, Trond I. Andersen, Giovanni Scuri, Andrew Y. Joe, Andrés M. Mier Valdivia, Xiaoling Liu, Alexander A. Zibrov, Bumho Kim, Takashi Taniguchi, Kenji Watanabe, James Hone, Valentin Walther, Hongkun Park, Philip Kim, Mikhail D. Lukin

Abstract: Excitons in transition metal dichalcogenides (TMDs) have emerged as a promising platform for novel applications ranging from optoelectronic devices to quantum optics and solid state quantum simulators. While much progress has been made towards characterizing and controlling excitons in TMDs, manipulating their properties during the course of their lifetime - a key requirement for many optoelectron… ▽ More Excitons in transition metal dichalcogenides (TMDs) have emerged as a promising platform for novel applications ranging from optoelectronic devices to quantum optics and solid state quantum simulators. While much progress has been made towards characterizing and controlling excitons in TMDs, manipulating their properties during the course of their lifetime - a key requirement for many optoelectronic device and information processing modalities - remains an outstanding challenge. Here we combine long-lived interlayer excitons in angle-aligned MoSe$_2$/WSe$_2$ heterostructures with fast electrical control to realize dynamical control schemes, in which exciton properties are not predetermined at the time of excitation but can be dynamically manipulated during their lifetime. Leveraging the out-of-plane exciton dipole moment, we use electric fields to demonstrate dynamical control over the exciton emission wavelength. Moreover, employing a patterned gate geometry, we demonstrate rapid local sample doping and toggling of the radiative decay rate through exciton-charge interactions during the exciton lifetime. Spatially mapping the exciton response reveals charge redistribution, offering a novel probe of electronic transport in twisted TMD heterostructures. Our results establish the feasibility of dynamical exciton control schemes, unlocking new directions for exciton-based information processing and optoelectronic devices, and the realization of excitonic phenomena in TMDs. △ Less

Submitted 17 July, 2024; v1 submitted 15 July, 2024; originally announced July 2024.

Comments: 37 pages, 4 figures in main text, 6 figures in supplemental materials; (v2) corrected funding acknowledgements

Constrained Bayesian optimization with merit functions

Authors: J. Wang, C. G. Petra, J. L. Peterson

Abstract: Bayesian optimization is a powerful optimization tool for problems where native first-order derivatives are unavailable. Recently, constrained Bayesian optimization (CBO) has been applied to many engineering applications where constraints are essential. However, several obstacles remain with current CBO algorithms that could prevent a wider adoption. We propose CBO algorithms using merit functions… ▽ More Bayesian optimization is a powerful optimization tool for problems where native first-order derivatives are unavailable. Recently, constrained Bayesian optimization (CBO) has been applied to many engineering applications where constraints are essential. However, several obstacles remain with current CBO algorithms that could prevent a wider adoption. We propose CBO algorithms using merit functions, such as the penalty merit function, in acquisition functions, inspired by nonlinear optimization methods, e.g., sequential quadratic programming. Merit functions measure the potential progress of both the objective and constraint functions, thus increasing algorithmic efficiency and allowing infeasible initial samples. The acquisition functions with merit functions are relaxed to have closed forms, making its implementation readily available wherever Bayesian optimization is. We further propose a unified CBO algorithm that can be seen as extension to the popular expected constrained improvement (ECI) approach. We demonstrate the effectiveness and efficiency of the proposed algorithms through numerical experiments on synthetic problems and a practical data-driven engineering design problem in the field of plasma physics. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Structure of Periodic Orbit Families in the Hill Restricted 4-Body Problem

Authors: Gavin M. Brown, Luke T. Peterson, Damennick B. Henry, Daniel J. Scheeres

Abstract: The Hill Restricted 4-Body Problem (HR4BP) is a coherent time-periodic model that can be used to represent motion in the Sun-Earth-Moon (SEM) system. Periodic orbits were computed in this model to better understand the periodic orbit family structures that exist in these types of systems. First, periodic orbits in the Circular Restricted 3-Body Problem (CR3BP) representation of the Earth-Moon (EM)… ▽ More The Hill Restricted 4-Body Problem (HR4BP) is a coherent time-periodic model that can be used to represent motion in the Sun-Earth-Moon (SEM) system. Periodic orbits were computed in this model to better understand the periodic orbit family structures that exist in these types of systems. First, periodic orbits in the Circular Restricted 3-Body Problem (CR3BP) representation of the Earth-Moon (EM) system were identified. A Melnikov-type function was used to identify a set of candidate points on the EM CR3BP periodic orbits to start a continuation algorithm. A pseudo-arclength continuation scheme was then used to obtain the corresponding periodic orbit families in the HR4BP when including the effect of the Sun. Bifurcation points were identified in the computed families to obtain additional orbit families. △ Less

Submitted 29 February, 2024; originally announced February 2024.

Comments: Submitted to SIAM Journal on Applied Dynamical Systems

MSC Class: 37N05 (Primary) 34C25; 70K60; 37M20 (Secondary)

Dynamics Around the Earth-Moon Triangular Points in the Hill Restricted 4-Body Problem

Authors: Luke T. Peterson, Gavin Brown, Àngel Jorba, Daniel Scheeres

Abstract: This paper investigates the motion of a small particle moving near the triangular points of the Earth-Moon system. The dynamics are modeled in the Hill restricted 4-body problem (HR4BP), which includes the effect of the Earth and Moon as in the circular restricted 3-body problem (CR3BP), as well as the direct and indirect effect of the Sun as a periodic time-dependent perturbation of the CR3BP. Du… ▽ More This paper investigates the motion of a small particle moving near the triangular points of the Earth-Moon system. The dynamics are modeled in the Hill restricted 4-body problem (HR4BP), which includes the effect of the Earth and Moon as in the circular restricted 3-body problem (CR3BP), as well as the direct and indirect effect of the Sun as a periodic time-dependent perturbation of the CR3BP. Due to the periodic perturbation, the triangular points of the CR3BP are no longer equilibrium solutions; rather, the triangular points are replaced by periodic orbits with the same period as the perturbation. Additionally, there is a 2:1 resonant periodic orbit that persists from the CR3BP into the HR4BP. In this work, we investigate the dynamics around these invariant objects by computing families of 2-dimensional invariant tori and their linear normal behavior. We identify bifurcations and relationships between families. Mechanisms for transport between Earth, L4, and Moon are discussed. Comparisons are made between the results presented here and in the bicircular problem (BCP). △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: 37 pages, 26 figures

Learning reduced-order Quadratic-Linear models in Process Engineering using Operator Inference

Authors: Ion Victor Gosea, Luisa Peterson, Pawan Goyal, Jens Bremer, Kai Sundmacher, Peter Benner

Abstract: In this work, we address the challenge of efficiently modeling dynamical systems in process engineering. We use reduced-order model learning, specifically operator inference. This is a non-intrusive, data-driven method for learning dynamical systems from time-domain data. The application in our study is carbon dioxide methanation, an important reaction within the Power-to-X framework, to demonstra… ▽ More In this work, we address the challenge of efficiently modeling dynamical systems in process engineering. We use reduced-order model learning, specifically operator inference. This is a non-intrusive, data-driven method for learning dynamical systems from time-domain data. The application in our study is carbon dioxide methanation, an important reaction within the Power-to-X framework, to demonstrate its potential. The numerical results show the ability of the reduced-order models constructed with operator inference to provide a reduced yet accurate surrogate solution. This represents an important milestone towards the implementation of fast and reliable digital twin architectures. △ Less

Submitted 30 July, 2024; v1 submitted 27 February, 2024; originally announced February 2024.

Comments: 10 pages, 3 figures

A multifidelity Bayesian optimization method for inertial confinement fusion design

Authors: J. Wang, N. Chiang, A. Gillette, J. L. Peterson

Abstract: Due to their cost, experiments for inertial confinement fusion (ICF) heavily rely on numerical simulations to guide design. As simulation technology progresses, so too can the fidelity of models used to plan for new experiments. However, these high-fidelity models are by themselves insufficient for optimal experimental design, because their computational cost remains too high to efficiently and ef… ▽ More Due to their cost, experiments for inertial confinement fusion (ICF) heavily rely on numerical simulations to guide design. As simulation technology progresses, so too can the fidelity of models used to plan for new experiments. However, these high-fidelity models are by themselves insufficient for optimal experimental design, because their computational cost remains too high to efficiently and effectively explore the numerous parameters required to describe a typical experiment. Traditionally, ICF design has relied on low-fidelity modeling to initially identify potentially interesting design regions, which are then subsequently explored via selected high-fidelity modeling. In this paper, we demonstrate that this two-step approach can be insufficient: even for simple design problems, a two-step optimization strategy can lead high-fidelity searching towards incorrect regions and consequently waste computational resources on parameter regimes far away from the true optimal solution. We reveal that a primary cause of this behavior in ICF design problems is the presence of low-fidelity optima in distinct regions of the parameter space from high-fidelity optima. To address this issue, we propose an iterative multifidelity Bayesian optimization method based on Gaussian Process Regression that leverages both low- and high-fidelity modelings. We demonstrate, using both two- and eight-dimensional ICF test problems, that our algorithm can effectively utilize low-fidelity modeling for exploration, while automatically refining promising designs with high-fidelity models. This approach proves to be more efficient than relying solely on high-fidelity modeling for optimization. △ Less

Submitted 15 December, 2023; originally announced December 2023.

MSC Class: 65K10

Journal ref: Physics of Plasmas 31 (3) (2024) 032706

Machine Learning-Based Multi-Objective Design Exploration Of Flexible Disc Elements

Authors: Gehendra Sharma, Sungkwang Mun, Nayeon Lee, Luke Peterson, Daniela Tellkamp, Anand Balu Nellippallil

Abstract: Design exploration is an important step in the engineering design process. This involves the search for design/s that meet the specified design criteria and accomplishes the predefined objective/s. In recent years, machine learning-based approaches have been widely used in engineering design problems. This paper showcases Artificial Neural Network (ANN) architecture applied to an engineering desig… ▽ More Design exploration is an important step in the engineering design process. This involves the search for design/s that meet the specified design criteria and accomplishes the predefined objective/s. In recent years, machine learning-based approaches have been widely used in engineering design problems. This paper showcases Artificial Neural Network (ANN) architecture applied to an engineering design problem to explore and identify improved design solutions. The case problem of this study is the design of flexible disc elements used in disc couplings. We are required to improve the design of the disc elements by lowering the mass and stress without lowering the torque transmission and misalignment capability. To accomplish this objective, we employ ANN coupled with genetic algorithm in the design exploration step to identify designs that meet the specified criteria (torque and misalignment) while having minimum mass and stress. The results are comparable to the optimized results obtained from the traditional response surface method. This can have huge advantage when we are evaluating conceptual designs against multiple conflicting requirements. △ Less

Submitted 14 April, 2023; originally announced April 2023.

Estimating Longitudinal Causal Effects with Unobserved Noncompliance Using a Semi-Parametric G-computation Algorithm

Authors: Ross L Peterson, David M Vock, Joseph S Koopmeiners

Abstract: Participant noncompliance, in which participants do not follow their assigned treatment protocol, often obscures the causal relationship between treatment and treatment effect in randomized trials. In the longitudinal setting, the G-computation algorithm can adjust for confounding to estimate causal effects. Typically, G-computation assumes that both 1) compliance is observed; and 2) the densities… ▽ More Participant noncompliance, in which participants do not follow their assigned treatment protocol, often obscures the causal relationship between treatment and treatment effect in randomized trials. In the longitudinal setting, the G-computation algorithm can adjust for confounding to estimate causal effects. Typically, G-computation assumes that both 1) compliance is observed; and 2) the densities of the confounders can be correctly specified. We aim to develop a G-computation estimator in the setting where both assumptions are violated. For 1), in place of unobserved compliance, we substitute in probability weights derived from modeling a biomarker associated with compliance. For 2), we fit semiparametric models using predictive mean matching. Specifically, we parametrically specify only the conditional mean of the confounders, and then use predictive mean matching to randomly generate confounder data for G-computation. In both the simulation and application, we compare multiple causal estimators already established in the literature with those derived from our method. For the simulation, we generated data across different sample sizes and levels of confounding. For the application, we apply our method to a trial that sought to evaluate the effect of cigarettes with low nicotine on cigarette consumption (Center for the Evaluation of Nicotine in Cigarettes Project 2 - CENIC-P2). △ Less

Submitted 7 February, 2023; originally announced February 2023.

Nonlinear dynamics of asymmetric bistable energy harvesters

Authors: João Pedro Norenberg, Roberto Luo, Vinicius Goncaalves Lopes, João Victor L. L. Peterson, Americo Cunha Jr

Abstract: The paper investigates asymmetries effects over a nonlinear vibration energy harvester dynamics. The asymmetric system performance is compared with symmetric ones. Different asymmetry levels on restoring force and gravity action are investigated from a system-sloping angle variation. Bifurcation diagrams and basins of attraction are used to examine the local and global characteristics underlying d… ▽ More The paper investigates asymmetries effects over a nonlinear vibration energy harvester dynamics. The asymmetric system performance is compared with symmetric ones. Different asymmetry levels on restoring force and gravity action are investigated from a system-sloping angle variation. Bifurcation diagrams and basins of attraction are used to examine the local and global characteristics underlying dynamical systems under different excitation energy. The results show the adverse effects of asymmetries on system dynamics. They also reveal ways to overcome them by canceling asymmetric influence from optimal sloping angle values and improving asymmetric system performance over symmetrical ones. This comprehensive numerical study provides novel valuable insights into asymmetrical energy harvester dynamics, a wide and still less explored topic. △ Less

Submitted 9 June, 2023; v1 submitted 20 August, 2022; originally announced September 2022.

Report number: Volume 257, 1 November 2023, 108542 MSC Class: 37N15 ACM Class: I.6.3

Journal ref: International Journal of Mechanical Sciences 2023

General Relativistic Implicit Monte Carlo Radiation-Hydrodynamics

Authors: Nathaniel Roth, Peter Anninos, Peter B. Robinson, J. Luc Peterson, Brooke Polak, Tymothy K. Mangan, Kyle Beyer

Abstract: We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++: an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of the radiation-hydrodynamics equations, but generalized for both Newtonian and relativistic regimes. A multiple reference frame approach is used to geodesically t… ▽ More We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++: an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of the radiation-hydrodynamics equations, but generalized for both Newtonian and relativistic regimes. A multiple reference frame approach is used to geodesically transport photon packets (and solve the hydrodynamics equations) in the coordinate frame, while radiation-matter interactions are handled either in the fluid or electron frames then communicated via Lorentz boosts and orthonormal tetrad bases attached to the fluid. We describe a method for constructing estimators of radiation moments using path-weighting that generalizes to arbitrary coordinate systems in flat or curved spacetime. Absorption, emission, scattering, and relativistic Comptonization are among the matter interactions considered in this report. We discuss our formulations and numerical methods, and validate our models against a suite of radiation and coupled radiation-hydrodynamics test problems in both flat and curved spacetimes. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: Accepted for publication in ApJS

2022 Review of Data-Driven Plasma Science

Authors: Rushil Anirudh, Rick Archibald, M. Salman Asif, Markus M. Becker, Sadruddin Benkadda, Peer-Timo Bremer, Rick H. S. Budé, C. S. Chang, Lei Chen, R. M. Churchill, Jonathan Citrin, Jim A Gaffney, Ana Gainaru, Walter Gekelman, Tom Gibbs, Satoshi Hamaguchi, Christian Hill, Kelli Humbird, Sören Jalas, Satoru Kawaguchi, Gon-Ho Kim, Manuel Kirchen, Scott Klasky, John L. Kline, Karl Krushelnick , et al. (38 additional authors not shown)

Abstract: Data science and technology offer transformative tools and methods to science. This review article highlights latest development and progress in the interdisciplinary field of data-driven plasma science (DDPS). A large amount of data and machine learning algorithms go hand in hand. Most plasma data, whether experimental, observational or computational, are generated or collected by machines today.… ▽ More Data science and technology offer transformative tools and methods to science. This review article highlights latest development and progress in the interdisciplinary field of data-driven plasma science (DDPS). A large amount of data and machine learning algorithms go hand in hand. Most plasma data, whether experimental, observational or computational, are generated or collected by machines today. It is now becoming impractical for humans to analyze all the data manually. Therefore, it is imperative to train machines to analyze and interpret (eventually) such data as intelligently as humans but far more efficiently in quantity. Despite the recent impressive progress in applications of data science to plasma science and technology, the emerging field of DDPS is still in its infancy. Fueled by some of the most challenging problems such as fusion energy, plasma processing of materials, and fundamental understanding of the universe through observable plasma phenomena, it is expected that DDPS continues to benefit significantly from the interdisciplinary marriage between plasma science and data science into the foreseeable future. △ Less

Submitted 31 May, 2022; originally announced May 2022.

Comments: 112 pages (including 700+ references), 44 figures, submitted to IEEE Transactions on Plasma Science as a part of the IEEE Golden Anniversary Special Issue

Report number: Los Alamos Report number LA-UR-22-24834

Journal ref: IEEE Transactions on Plasma Science 51, 1750 - 1838 (2023)

Transfer learning driven design optimization for inertial confinement fusion

Authors: K. D. Humbird, J. L. Peterson

Abstract: Transfer learning is a promising approach to creating predictive models that incorporate simulation and experimental data into a common framework. In this technique, a neural network is first trained on a large database of simulations, then partially retrained on sparse sets of experimental data to adjust predictions to be more consistent with reality. Previously, this technique has been used to c… ▽ More Transfer learning is a promising approach to creating predictive models that incorporate simulation and experimental data into a common framework. In this technique, a neural network is first trained on a large database of simulations, then partially retrained on sparse sets of experimental data to adjust predictions to be more consistent with reality. Previously, this technique has been used to create predictive models of Omega and NIF inertial confinement fusion (ICF) experiments that are more accurate than simulations alone. In this work, we conduct a transfer learning driven hypothetical ICF campaign in which the goal is to maximize experimental neutron yield via Bayesian optimization. The transfer learning model achieves yields within 5% of the maximum achievable yield in a modest-sized design space in fewer than 20 experiments. Furthermore, we demonstrate that this method is more efficient at optimizing designs than traditional model calibration techniques commonly employed in ICF design. Such an approach to ICF design could enable robust optimization of experimental performance under uncertainty. △ Less

Submitted 26 May, 2022; originally announced May 2022.

The data-driven future of high energy density physics

Authors: Peter W. Hatfield, Jim A. Gaffney, Gemma J. Anderson, Suzanne Ali, Luca Antonelli, Suzan Başeğmez du Pree, Jonathan Citrin, Marta Fajardo, Patrick Knapp, Brendan Kettle, Bogdan Kustowski, Michael J. MacDonald, Derek Mariscal, Madison E. Martin, Taisuke Nagayama, Charlotte A. J. Palmer, J. Luc Peterson, Steven Rose, J J Ruby, Carl Shneider, Matt J. V. Streeter, Will Trickey, Ben Williams

Abstract: The study of plasma physics under conditions of extreme temperatures, densities and electromagnetic field strengths is significant for our understanding of astrophysics, nuclear fusion and fundamental physics. These extreme physical systems are strongly non-linear and very difficult to understand theoretically or optimize experimentally. Here, we argue that machine learning models and data-driven… ▽ More The study of plasma physics under conditions of extreme temperatures, densities and electromagnetic field strengths is significant for our understanding of astrophysics, nuclear fusion and fundamental physics. These extreme physical systems are strongly non-linear and very difficult to understand theoretically or optimize experimentally. Here, we argue that machine learning models and data-driven methods are in the process of reshaping our exploration of these extreme systems that have hitherto proven far too non-linear for human researchers. From a fundamental perspective, our understanding can be helped by the way in which machine learning models can rapidly discover complex interactions in large data sets. From a practical point of view, the newest generation of extreme physics facilities can perform experiments multiple times a second (as opposed to ~daily), moving away from human-based control towards automatic control based on real-time interpretation of diagnostic data and updates of the physics model. To make the most of these emerging opportunities, we advance proposals for the community in terms of research design, training, best practices, and support for synthetic diagnostics and data analysis. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: 14 pages, 4 figures. This work was the result of a meeting at the Lorentz Center, University of Leiden, 13th-17th January 2020. This is a preprint of Hatfield et al., Nature, 593, 7859, 351-361 (2021) https://www.nature.com/articles/s41586-021-03382-w

Journal ref: Nature, 593, 7859, 351-361, 2021

Experiments conducted in the burning plasma regime with inertial fusion implosions

Authors: J. S. Ross, J. E. Ralph, A. B. Zylstra, A. L. Kritcher, H. F. Robey, C. V. Young, O. A. Hurricane, D. A. Callahan, K. L. Baker, D. T. Casey, T. Doeppner, L. Divol, M. Hohenberger, S. Le Pape, A. Pak, P. K. Patel, R. Tommasini, S. J. Ali, P. A. Amendt, L. J. Atherton, B. Bachmann, D. Bailey, L. R. Benedetti, L. Berzak Hopkins, R. Betti , et al. (127 additional authors not shown)

Abstract: An experimental program is currently underway at the National Ignition Facility (NIF) to compress deuterium and tritium (DT) fuel to densities and temperatures sufficient to achieve fusion and energy gain. The primary approach being investigated is indirect drive inertial confinement fusion (ICF), where a high-Z radiation cavity (a hohlraum) is heated by lasers, converting the incident energy into… ▽ More An experimental program is currently underway at the National Ignition Facility (NIF) to compress deuterium and tritium (DT) fuel to densities and temperatures sufficient to achieve fusion and energy gain. The primary approach being investigated is indirect drive inertial confinement fusion (ICF), where a high-Z radiation cavity (a hohlraum) is heated by lasers, converting the incident energy into x-ray radiation which in turn drives the DT fuel filled capsule causing it to implode. Previous experiments reported DT fuel gain exceeding unity [O.A. Hurricane et al., Nature 506, 343 (2014)] and then exceeding the kinetic energy of the imploding fuel [S. Le Pape et al., Phys. Rev. Lett. 120, 245003 (2018)]. We report on recent experiments that have achieved record fusion neutron yields on NIF, greater than 100 kJ with momentary fusion powers exceeding 1PW, and have for the first time entered the burning plasma regime where fusion alpha-heating of the fuel exceeds the energy delivered to the fuel via compression. This was accomplished by increasing the size of the high-density carbon (HDC) capsule, increasing energy coupling, while controlling symmetry and implosion design parameters. Two tactics were successful in controlling the radiation flux symmetry and therefore the implosion symmetry: transferring energy between laser cones via plasma waves, and changing the shape of the hohlraum. In conducting these experiments, we controlled for known sources of degradation. Herein we show how these experiments were performed to produce record performance, and demonstrate the data fidelity leading us to conclude that these shots have entered the burning plasma regime. △ Less

Submitted 8 November, 2021; originally announced November 2021.

A Community Roadmap for Scientific Workflows Research and Development

Authors: Rafael Ferreira da Silva, Henri Casanova, Kyle Chard, Ilkay Altintas, Rosa M Badia, Bartosz Balis, Tainã Coleman, Frederik Coppens, Frank Di Natale, Bjoern Enders, Thomas Fahringer, Rosa Filgueira, Grigori Fursin, Daniel Garijo, Carole Goble, Dorran Howell, Shantenu Jha, Daniel S. Katz, Daniel Laney, Ulf Leser, Maciej Malawski, Kshitij Mehta, Loïc Pottier, Jonathan Ozik, J. Luc Peterson , et al. (4 additional authors not shown)

Abstract: The landscape of workflow systems for scientific applications is notoriously convoluted with hundreds of seemingly equivalent workflow systems, many isolated research claims, and a steep learning curve. To address some of these challenges and lay the groundwork for transforming workflows research and development, the WorkflowsRI and ExaWorks projects partnered to bring the international workflows… ▽ More The landscape of workflow systems for scientific applications is notoriously convoluted with hundreds of seemingly equivalent workflow systems, many isolated research claims, and a steep learning curve. To address some of these challenges and lay the groundwork for transforming workflows research and development, the WorkflowsRI and ExaWorks projects partnered to bring the international workflows community together. This paper reports on discussions and findings from two virtual "Workflows Community Summits" (January and April, 2021). The overarching goals of these workshops were to develop a view of the state of the art, identify crucial research challenges in the workflows community, articulate a vision for potential community efforts, and discuss technical approaches for realizing this vision. To this end, participants identified six broad themes: FAIR computational workflows; AI workflows; exascale challenges; APIs, interoperability, reuse, and standards; training and education; and building a workflows community. We summarize discussions and recommendations for each of these themes. △ Less

Submitted 8 October, 2021; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: arXiv admin note: substantial text overlap with arXiv:2103.09181

The nonlinear dynamics of a bistable energy harvesting system with colored noise disturbances

Authors: Vinicius Gonçalves Lopes, João Victor L. L. Peterson, Americo Cunha Jr

Abstract: This paper deals with the nonlinear stochastic dynamics of a piezoelectric energy harvesting system subjected to a harmonic external excitation disturbed by Gaussian colored noise. A parametric analysis is conducted, where the effects of the standard deviation and the correlation time of colored noise on the system response are investigated. The numerical results suggest a strong influence of nois… ▽ More This paper deals with the nonlinear stochastic dynamics of a piezoelectric energy harvesting system subjected to a harmonic external excitation disturbed by Gaussian colored noise. A parametric analysis is conducted, where the effects of the standard deviation and the correlation time of colored noise on the system response are investigated. The numerical results suggest a strong influence of noise on the system response for higher values of correlation time and standard deviation, and a low (noise level independent) influence for low values of correlation time. △ Less

Submitted 14 July, 2021; originally announced July 2021.

MSC Class: 37N15 ACM Class: I.6.6

Journal ref: Journal of Computational Interdisciplinary Sciences, vol. 10, pp. 125, 2019

A conformally invariant Yang-Mills type energy and equation on 6-manifolds

Authors: A. Rod Gover, Lawrence J. Peterson, Callum Sleigh

Abstract: We define a conformally invariant action S on gauge connections on a closed pseudo-Riemannian manifold M of dimension 6. At leading order this is quadratic in the gauge connection. The Euler-Lagrange equations of S, with respect to variation of the gauge connection, provide a higher-order conformally invariant analogue of the (source-free) Yang-Mills equations. For any gauge connection A on M, w… ▽ More We define a conformally invariant action S on gauge connections on a closed pseudo-Riemannian manifold M of dimension 6. At leading order this is quadratic in the gauge connection. The Euler-Lagrange equations of S, with respect to variation of the gauge connection, provide a higher-order conformally invariant analogue of the (source-free) Yang-Mills equations. For any gauge connection A on M, we define S(A) by first defining a Lagrangian density associated to A. This is not conformally invariant but has a conformal transformation analogous to a Q-curvature. Integrating this density provides the conformally invariant action. In the special case that we apply S to the conformal Cartan-tractor connection, the functional gradient recovers the natural conformal curvature invariant called the Fefferman-Graham obstruction tensor. So in this case the Euler-Lagrange equations are exactly the "obstruction-flat" condition for 6-manifolds. This extends known results for 4-dimensional pseudo-Riemannian manifolds where the Bach tensor is recovered in the Yang-Mills equations of the Cartan-tractor connection. △ Less

Submitted 18 July, 2021; originally announced July 2021.

Comments: 24 pages. Uses tensor.sty

MSC Class: 53C18; 53C07 (Primary) 58E15; 81T13; 53A55 (Secondary)

Probing spin dynamics on diamond surfaces using a single quantum sensor

Authors: Bo L. Dwyer, Lila V. H. Rodgers, Elana K. Urbach, Dolev Bluvstein, Sorawis Sangtawesin, Hengyun Zhou, Yahia Nassab, Mattias Fitzpatrick, Zhiyang Yuan, Kristiaan De Greve, Eric L. Peterson, Jyh-Pin Chou, Adam Gali, V. V. Dobrovitski, Mikhail D. Lukin, Nathalie P. de Leon

Abstract: Understanding the dynamics of a quantum bit's environment is essential for the realization of practical systems for quantum information processing and metrology. We use single nitrogen-vacancy (NV) centers in diamond to study the dynamics of a disordered spin ensemble at the diamond surface. Specifically, we tune the density of "dark" surface spins to interrogate their contribution to the decohere… ▽ More Understanding the dynamics of a quantum bit's environment is essential for the realization of practical systems for quantum information processing and metrology. We use single nitrogen-vacancy (NV) centers in diamond to study the dynamics of a disordered spin ensemble at the diamond surface. Specifically, we tune the density of "dark" surface spins to interrogate their contribution to the decoherence of shallow NV center spin qubits. When the average surface spin spacing exceeds the NV center depth, we find that the surface spin contribution to the NV center free induction decay can be described by a stretched exponential with variable power n. We show that these observations are consistent with a model in which the spatial positions of the surface spins are fixed for each measurement, but some of them reconfigure between measurements. In particular, we observe a depth-dependent critical time associated with a dynamical transition from Gaussian (n=2) decay to n=2/3, and show that this transition arises from the competition between the small decay contributions of many distant spins and strong coupling to a few proximal spins at the surface. These observations demonstrate the potential of a local sensor for understanding complex systems and elucidate pathways for improving and controlling spin qubits at the surface. △ Less

Submitted 23 March, 2021; originally announced March 2021.

Comments: 32 pages, 4 figures main text, 15 figures supplementary information

Cognitive simulation models for inertial confinement fusion: Combining simulation and experimental data

Authors: K. D. Humbird, J. L. Peterson, J. Salmonson, B. K. Spears

Abstract: The design space for inertial confinement fusion (ICF) experiments is vast and experiments are extremely expensive. Researchers rely heavily on computer simulations to explore the design space in search of high-performing implosions. However, ICF multiphysics codes must make simplifying assumptions, and thus deviate from experimental measurements for complex implosions. For more effective design a… ▽ More The design space for inertial confinement fusion (ICF) experiments is vast and experiments are extremely expensive. Researchers rely heavily on computer simulations to explore the design space in search of high-performing implosions. However, ICF multiphysics codes must make simplifying assumptions, and thus deviate from experimental measurements for complex implosions. For more effective design and investigation, simulations require input from past experimental data to better predict future performance. In this work, we describe a cognitive simulation method for combining simulation and experimental data into a common, predictive model. This method leverages a machine learning technique called transfer learning, the process of taking a model trained to solve one task, and partially retraining it on a sparse dataset to solve a different, but related task. In the context of ICF design, neural network models trained on large simulation databases and partially retrained on experimental data, producing models that are far more accurate than simulations alone. We demonstrate improved model performance for a range of ICF experiments at the National Ignition Facility, and predict the outcome of recent experiments with less than ten percent error for several key observables. We discuss how the methods might be used to carry out a data-driven experimental campaign to optimize performance, illustrating the key product -- models that become increasingly accurate as data is acquired. △ Less

Submitted 18 March, 2021; originally announced March 2021.

Workflows Community Summit: Bringing the Scientific Workflows Community Together

Authors: Rafael Ferreira da Silva, Henri Casanova, Kyle Chard, Dan Laney, Dong Ahn, Shantenu Jha, Carole Goble, Lavanya Ramakrishnan, Luc Peterson, Bjoern Enders, Douglas Thain, Ilkay Altintas, Yadu Babuji, Rosa M. Badia, Vivien Bonazzi, Taina Coleman, Michael Crusoe, Ewa Deelman, Frank Di Natale, Paolo Di Tommaso, Thomas Fahringer, Rosa Filgueira, Grigori Fursin, Alex Ganose, Bjorn Gruning , et al. (20 additional authors not shown)

Abstract: Scientific workflows have been used almost universally across scientific domains, and have underpinned some of the most significant discoveries of the past several decades. Many of these workflows have high computational, storage, and/or communication demands, and thus must execute on a wide range of large-scale platforms, from large clouds to upcoming exascale high-performance computing (HPC) pla… ▽ More Scientific workflows have been used almost universally across scientific domains, and have underpinned some of the most significant discoveries of the past several decades. Many of these workflows have high computational, storage, and/or communication demands, and thus must execute on a wide range of large-scale platforms, from large clouds to upcoming exascale high-performance computing (HPC) platforms. These executions must be managed using some software infrastructure. Due to the popularity of workflows, workflow management systems (WMSs) have been developed to provide abstractions for creating and executing workflows conveniently, efficiently, and portably. While these efforts are all worthwhile, there are now hundreds of independent WMSs, many of which are moribund. As a result, the WMS landscape is segmented and presents significant barriers to entry due to the hundreds of seemingly comparable, yet incompatible, systems that exist. As a result, many teams, small and large, still elect to build their own custom workflow solution rather than adopt, or build upon, existing WMSs. This current state of the WMS landscape negatively impacts workflow users, developers, and researchers. The "Workflows Community Summit" was held online on January 13, 2021. The overarching goal of the summit was to develop a view of the state of the art and identify crucial research challenges in the workflow community. Prior to the summit, a survey sent to stakeholders in the workflow community (including both developers of WMSs and users of workflows) helped to identify key challenges in this community that were translated into 6 broad themes for the summit, each of them being the object of a focused discussion led by a volunteer member of the community. This report documents and organizes the wealth of information provided by the participants before, during, and after the summit. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Predicting Opioid Use Disorder from Longitudinal Healthcare Data using Multi-stream Transformer

Authors: Sajjad Fouladvand, Jeffery Talbert, Linda P. Dwoskin, Heather Bush, Amy Lynn Meadows, Lars E. Peterson, Ramakanth Kavuluru, Jin Chen

Abstract: Opioid Use Disorder (OUD) is a public health crisis costing the US billions of dollars annually in healthcare, lost workplace productivity, and crime. Analyzing longitudinal healthcare data is critical in addressing many real-world problems in healthcare. Leveraging the real-world longitudinal healthcare data, we propose a novel multi-stream transformer model called MUPOD for OUD identification. M… ▽ More Opioid Use Disorder (OUD) is a public health crisis costing the US billions of dollars annually in healthcare, lost workplace productivity, and crime. Analyzing longitudinal healthcare data is critical in addressing many real-world problems in healthcare. Leveraging the real-world longitudinal healthcare data, we propose a novel multi-stream transformer model called MUPOD for OUD identification. MUPOD is designed to simultaneously analyze multiple types of healthcare data streams, such as medications and diagnoses, by attending to segments within and across these data streams. Our model tested on the data from 392,492 patients with long-term back pain problems showed significantly better performance than the traditional models and recently developed deep learning models. △ Less

Submitted 7 July, 2021; v1 submitted 15 March, 2021; originally announced March 2021.

Comments: This manuscript has been accepted by AMIA 2021 for oral presentation on November 1, 2021

Computing Research Challenges in Next Generation Wireless Networking

Authors: Elisa Bertino, Daniel Bliss, Daniel Lopresti, Larry Peterson, Henning Schulzrinne

Abstract: By all measures, wireless networking has seen explosive growth over the past decade. Fourth Generation Long Term Evolution (4G LTE) cellular technology has increased the bandwidth available for smartphones, in essence, delivering broadband speeds to mobile devices. The most recent 5G technology is further enhancing the transmission speeds and cell capacity, as well as, reducing latency through the… ▽ More By all measures, wireless networking has seen explosive growth over the past decade. Fourth Generation Long Term Evolution (4G LTE) cellular technology has increased the bandwidth available for smartphones, in essence, delivering broadband speeds to mobile devices. The most recent 5G technology is further enhancing the transmission speeds and cell capacity, as well as, reducing latency through the use of different radio technologies and is expected to provide Internet connections that are an order of magnitude faster than 4G LTE. Technology continues to advance rapidly, however, and the next generation, 6G, is already being envisioned. 6G will make possible a wide range of powerful, new applications including holographic telepresence, telehealth, remote education, ubiquitous robotics and autonomous vehicles, smart cities and communities (IoT), and advanced manufacturing (Industry 4.0, sometimes referred to as the Fourth Industrial Revolution), to name but a few. The advances we will see begin at the hardware level and extend all the way to the top of the software "stack." Artificial Intelligence (AI) will also start playing a greater role in the development and management of wireless networking infrastructure by becoming embedded in applications throughout all levels of the network. The resulting benefits to society will be enormous. At the same time these exciting new wireless capabilities are appearing rapidly on the horizon, a broad range of research challenges loom ahead. These stem from the ever-increasing complexity of the hardware and software systems, along with the need to provide infrastructure that is robust and secure while simultaneously protecting the privacy of users. Here we outline some of those challenges and provide recommendations for the research that needs to be done to address them. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: A Computing Community Consortium (CCC) white paper, 5 pages

Report number: ccc2020whitepaper_16

Enabling Machine Learning-Ready HPC Ensembles with Merlin

Authors: J. Luc Peterson, Ben Bay, Joe Koning, Peter Robinson, Jessica Semler, Jeremy White, Rushil Anirudh, Kevin Athey, Peer-Timo Bremer, Francesco Di Natale, David Fox, Jim A. Gaffney, Sam A. Jacobs, Bhavya Kailkhura, Bogdan Kustowski, Steven Langer, Brian Spears, Jayaraman Thiagarajan, Brian Van Essen, Jae-Seung Yeom

Abstract: With the growing complexity of computational and experimental facilities, many scientific researchers are turning to machine learning (ML) techniques to analyze large scale ensemble data. With complexities such as multi-component workflows, heterogeneous machine architectures, parallel file systems, and batch scheduling, care must be taken to facilitate this analysis in a high performance computin… ▽ More With the growing complexity of computational and experimental facilities, many scientific researchers are turning to machine learning (ML) techniques to analyze large scale ensemble data. With complexities such as multi-component workflows, heterogeneous machine architectures, parallel file systems, and batch scheduling, care must be taken to facilitate this analysis in a high performance computing (HPC) environment. In this paper, we present Merlin, a workflow framework to enable large ML-friendly ensembles of scientific HPC simulations. By augmenting traditional HPC with distributed compute technologies, Merlin aims to lower the barrier for scientific subject matter experts to incorporate ML into their analysis. In addition to its design, we describe some example applications that Merlin has enabled on leadership-class HPC resources, such as the ML-augmented optimization of nuclear fusion experiments and the calibration of infectious disease models to study the progression of and possible mitigation strategies for COVID-19. △ Less

Submitted 1 July, 2021; v1 submitted 5 December, 2019; originally announced December 2019.

Comments: 28 pages, 9 figures; Submitted to FGCS

Report number: LLNL-JRNL-821884

O2- and CO-Rich Atmospheres for Potentially Habitable Environments on TRAPPIST-1 Planets

Authors: Renyu Hu, Luke Peterson, Eric T. Wolf

Abstract: Small exoplanets of nearby M dwarf stars present the possibility to find and characterize habitable worlds within the next decade. TRAPPIST-1, an ultracool M dwarf star, was recently found to have seven Earth-sized planets of predominantly rocky composition. The planets e, f, and g can have a liquid water ocean on their surface given appropriate atmospheres of N2 and CO2. Particularly, climate mod… ▽ More Small exoplanets of nearby M dwarf stars present the possibility to find and characterize habitable worlds within the next decade. TRAPPIST-1, an ultracool M dwarf star, was recently found to have seven Earth-sized planets of predominantly rocky composition. The planets e, f, and g can have a liquid water ocean on their surface given appropriate atmospheres of N2 and CO2. Particularly, climate models have shown that the planets e and f can sustain a global liquid water ocean, for >=0.2 bar CO2 plus 1 bar N2, or >=2 bars CO2, respectively. These atmospheres are irradiated by ultraviolet emission from the star's moderately active chromosphere, and the consequence of this irradiation is unknown. Here we show that chemical reactions driven by the irradiation can produce and maintain more than 0.2 bar O2 and 0.05 bar CO if the CO2 is >=0.1 bar. The abundance of O2 and CO can rise to more than 1 bar under certain boundary conditions. Because of this O2-CO runaway, habitable environments on the TRAPPIST-1 planets entail an O2- and CO-rich atmosphere with coexisting O3. The only process that would prevent the runaway is direct recombination of O2 and CO in the ocean, a reaction that is facilitated biologically. Our results indicate that O2, O3, and CO should be considered together with CO2 as the primary molecules in the search for atmospheric signatures from temperate and rocky planets of TRAPPIST-1 and other M dwarf stars. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: ApJ accepted

Modeling, simulation and validation of supersonic parachute inflation dynamics during Mars landing

Authors: Daniel Z. Huang, Philip Avery, Charbel Farhat, Jason Rabinovitch, Armen Derkevorkian, Lee D Peterson

Abstract: A high fidelity multi-physics Eulerian computational framework is presented for the simulation of supersonic parachute inflation during Mars landing. Unlike previous investigations in this area, the framework takes into account an initial folding pattern of the parachute, the flow compressibility effect on the fabric material porosity, and the interactions between supersonic fluid flows and the su… ▽ More A high fidelity multi-physics Eulerian computational framework is presented for the simulation of supersonic parachute inflation during Mars landing. Unlike previous investigations in this area, the framework takes into account an initial folding pattern of the parachute, the flow compressibility effect on the fabric material porosity, and the interactions between supersonic fluid flows and the suspension lines. Several adaptive mesh refinement (AMR)-enabled, large edge simulation (LES)-based, simulations of a full-size disk-gap-band (DGB) parachute inflating in the low-density, low-pressure, carbon dioxide (CO2) Martian atmosphere are reported. The comparison of the drag histories and the first peak forces between the simulation results and experimental data collected during the NASA Curiosity Rover's Mars atmospheric entry shows reasonable agreements. Furthermore, a rudimentary material failure analysis is performed to provide an estimate of the safety factor for the parachute decelerator system. The proposed framework demonstrates the potential of using Computational Fluid Dynamics (CFD) and Fluid-Structure Interaction (FSI)-based simulation tools for future supersonic parachute design. △ Less

Submitted 3 December, 2019; originally announced December 2019.

Comments: 24 pages, 12 figures

Parallelizing Training of Deep Generative Models on Massive Scientific Datasets

Authors: Sam Ade Jacobs, Brian Van Essen, David Hysom, Jae-Seung Yeom, Tim Moon, Rushil Anirudh, Jayaraman J. Thiagaranjan, Shusen Liu, Peer-Timo Bremer, Jim Gaffney, Tom Benson, Peter Robinson, Luc Peterson, Brian Spears

Abstract: Training deep neural networks on large scientific data is a challenging task that requires enormous compute power, especially if no pre-trained models exist to initialize the process. We present a novel tournament method to train traditional as well as generative adversarial networks built on LBANN, a scalable deep learning framework optimized for HPC systems. LBANN combines multiple levels of par… ▽ More Training deep neural networks on large scientific data is a challenging task that requires enormous compute power, especially if no pre-trained models exist to initialize the process. We present a novel tournament method to train traditional as well as generative adversarial networks built on LBANN, a scalable deep learning framework optimized for HPC systems. LBANN combines multiple levels of parallelism and exploits some of the worlds largest supercomputers. We demonstrate our framework by creating a complex predictive model based on multi-variate data from high-energy-density physics containing hundreds of millions of images and hundreds of millions of scalar values derived from tens of millions of simulations of inertial confinement fusion. Our approach combines an HPC workflow and extends LBANN with optimized data ingestion and the new tournament-style training algorithm to produce a scalable neural network architecture using a CORAL-class supercomputer. Experimental results show that 64 trainers (1024 GPUs) achieve a speedup of 70.2 over a single trainer (16 GPUs) baseline, and an effective 109% parallel efficiency. △ Less

Submitted 5 October, 2019; originally announced October 2019.

Scalable Topological Data Analysis and Visualization for Evaluating Data-Driven Models in Scientific Applications

Authors: Shusen Liu, Di Wang, Dan Maljovec, Rushil Anirudh, Jayaraman J. Thiagarajan, Sam Ade Jacobs, Brian C. Van Essen, David Hysom, Jae-Seung Yeom, Jim Gaffney, Luc Peterson, Peter B. Robinson, Harsh Bhatia, Valerio Pascucci, Brian K. Spears, Peer-Timo Bremer

Abstract: With the rapid adoption of machine learning techniques for large-scale applications in science and engineering comes the convergence of two grand challenges in visualization. First, the utilization of black box models (e.g., deep neural networks) calls for advanced techniques in exploring and interpreting model behaviors. Second, the rapid growth in computing has produced enormous datasets that re… ▽ More With the rapid adoption of machine learning techniques for large-scale applications in science and engineering comes the convergence of two grand challenges in visualization. First, the utilization of black box models (e.g., deep neural networks) calls for advanced techniques in exploring and interpreting model behaviors. Second, the rapid growth in computing has produced enormous datasets that require techniques that can handle millions or more samples. Although some solutions to these interpretability challenges have been proposed, they typically do not scale beyond thousands of samples, nor do they provide the high-level intuition scientists are looking for. Here, we present the first scalable solution to explore and analyze high-dimensional functions often encountered in the scientific data analysis pipeline. By combining a new streaming neighborhood graph construction, the corresponding topology computation, and a novel data aggregation scheme, namely topology aware datacubes, we enable interactive exploration of both the topological and the geometric aspect of high-dimensional data. Following two use cases from high-energy-density (HED) physics and computational biology, we demonstrate how these capabilities have led to crucial new insights in both applications. △ Less

Submitted 18 July, 2019; originally announced July 2019.

Mean-dependent nonstationary spatial models

Authors: Geoffrey Colin Lee Peterson, Joseph Guinness, Adam Terando, Brian J. Reich

Abstract: Nonstationarity is a major challenge in analyzing spatial data. For example, daily precipitation measurements may have increased variability and decreased spatial smoothness in areas with high mean rainfall. Common nonstationary covariance models introduce parameters specific to each location, giving a highly-parameterized model which is difficult to fit. We develop a nonstationary spatial model t… ▽ More Nonstationarity is a major challenge in analyzing spatial data. For example, daily precipitation measurements may have increased variability and decreased spatial smoothness in areas with high mean rainfall. Common nonstationary covariance models introduce parameters specific to each location, giving a highly-parameterized model which is difficult to fit. We develop a nonstationary spatial model that uses the mean to determine the covariance in a region, resulting in a far simpler, albeit more specialized, model. We explore inferential and predictive properties of the model under various simulated data situations. We show that this model in certain circumstances improves predictions compared to a standard stationary spatial model. We further propose a computationally efficient approximation that has comparable predictive accuracy. We also develop a test for nonstationary data and show it reliably identifies nonstationarity. We apply these methods to daily precipitation in Puerto Rico. △ Less

Submitted 29 May, 2019; originally announced May 2019.

Transfer learning to model inertial confinement fusion experiments

Authors: K. D. Humbird, J. L. Peterson, R. G. McClarren

Abstract: Inertial confinement fusion (ICF) experiments are designed using computer simulations that are approximations of reality, and therefore must be calibrated to accurately predict experimental observations. In this work, we propose a novel nonlinear technique for calibrating from simulations to experiments, or from low fidelity simulations to high fidelity simulations, via "transfer learning". Transf… ▽ More Inertial confinement fusion (ICF) experiments are designed using computer simulations that are approximations of reality, and therefore must be calibrated to accurately predict experimental observations. In this work, we propose a novel nonlinear technique for calibrating from simulations to experiments, or from low fidelity simulations to high fidelity simulations, via "transfer learning". Transfer learning is a commonly used technique in the machine learning community, in which models trained on one task are partially retrained to solve a separate, but related task, for which there is a limited quantity of data. We introduce the idea of hierarchical transfer learning, in which neural networks trained on low fidelity models are calibrated to high fidelity models, then to experimental data. This technique essentially bootstraps the calibration process, enabling the creation of models which predict high fidelity simulations or experiments with minimal computational cost. We apply this technique to a database of ICF simulations and experiments carried out at the Omega laser facility. Transfer learning with deep neural networks enables the creation of models that are more predictive of Omega experiments than simulations alone. The calibrated models accurately predict future Omega experiments, and are used to search for new, optimal implosion designs. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Predicting the time-evolution of multi-physics systems with sequence-to-sequence models

Authors: K. D. Humbird, J. L. Peterson, R. G. McClarren

Abstract: In this work, sequence-to-sequence (seq2seq) models, originally developed for language translation, are used to predict the temporal evolution of complex, multi-physics computer simulations. The predictive performance of seq2seq models is compared to state transition models for datasets generated with multi-physics codes with varying levels of complexity - from simple 1D diffusion calculations to… ▽ More In this work, sequence-to-sequence (seq2seq) models, originally developed for language translation, are used to predict the temporal evolution of complex, multi-physics computer simulations. The predictive performance of seq2seq models is compared to state transition models for datasets generated with multi-physics codes with varying levels of complexity - from simple 1D diffusion calculations to simulations of inertial confinement fusion implosions. Seq2seq models demonstrate the ability to accurately emulate complex systems, enabling the rapid estimation of the evolution of quantities of interest in computationally expensive simulations. △ Less

Submitted 14 November, 2018; originally announced November 2018.

Resampling Forgery Detection Using Deep Learning and A-Contrario Analysis

Authors: Arjuna Flenner, Lawrence Peterson, Jason Bunk, Tajuddin Manhar Mohammed, Lakshmanan Nataraj, B. S. Manjunath

Abstract: The amount of digital imagery recorded has recently grown exponentially, and with the advancement of software, such as Photoshop or Gimp, it has become easier to manipulate images. However, most images on the internet have not been manipulated and any automated manipulation detection algorithm must carefully control the false alarm rate. In this paper we discuss a method to automatically detect lo… ▽ More The amount of digital imagery recorded has recently grown exponentially, and with the advancement of software, such as Photoshop or Gimp, it has become easier to manipulate images. However, most images on the internet have not been manipulated and any automated manipulation detection algorithm must carefully control the false alarm rate. In this paper we discuss a method to automatically detect local resampling using deep learning while controlling the false alarm rate using a-contrario analysis. The automated procedure consists of three primary steps. First, resampling features are calculated for image blocks. A deep learning classifier is then used to generate a heatmap that indicates if the image block has been resampled. We expect some of these blocks to be falsely identified as resampled. We use a-contrario hypothesis testing to both identify if the patterns of the manipulated blocks indicate if the image has been tampered with and to localize the manipulation. We demonstrate that this strategy is effective in indicating if an image has been manipulated and localizing the manipulations. △ Less

Submitted 1 March, 2018; originally announced March 2018.

Comments: arXiv admin note: text overlap with arXiv:1802.03154

Conformal boundary operators, T-curvatures, and conformal fractional Laplacians of odd order

Authors: A. Rod Gover, Lawrence J. Peterson

Abstract: We construct continuously parametrised families of conformally invariant boundary operators on densities. These may also be viewed as conformally covariant boundary operators on functions and generalise to higher orders the first-order conformal Robin operator and an analogous third-order operator of Chang-Qing. Our families include operators of critical order on odd-dimensional boundaries. Combin… ▽ More We construct continuously parametrised families of conformally invariant boundary operators on densities. These may also be viewed as conformally covariant boundary operators on functions and generalise to higher orders the first-order conformal Robin operator and an analogous third-order operator of Chang-Qing. Our families include operators of critical order on odd-dimensional boundaries. Combined with the (conformal Laplacian power) GJMS operators, a suitable selection of the boundary operators yields formally self-adjoint elliptic conformal boundary problems. Working on a conformal manifold with boundary, we show that the operators yield odd-order conformally invariant fractional Laplacian pseudo-differential operators. To do this, we use higher-order conformally invariant Dirichlet-to-Neumann constructions. We also find and construct new curvature quantities associated to our new operator families. These have links to the Branson Q-curvature and include higher-order generalisations of the mean curvature and the T-curvature of Chang-Qing. In the case of the standard conformal hemisphere, the boundary operator construction is particularly simple; the resulting operators provide an elementary construction of families of symmetry breaking intertwinors between the spherical principal series representations of the conformal group of the equator, as studied by Juhl and others. We use our constructions to shed light on some conjectures of Juhl. △ Less

Submitted 22 February, 2018; originally announced February 2018.

Comments: 58 pages

MSC Class: 53A30 (Primary) 53A55; 34B05; 35S05; 53B25; 53B50 (Secondary)

Journal ref: Pacific J. Math. 311 (2021) 277-328

Boosting Image Forgery Detection using Resampling Features and Copy-move analysis

Authors: Tajuddin Manhar Mohammed, Jason Bunk, Lakshmanan Nataraj, Jawadul H. Bappy, Arjuna Flenner, B. S. Manjunath, Shivkumar Chandrasekaran, Amit K. Roy-Chowdhury, Lawrence Peterson

Abstract: Realistic image forgeries involve a combination of splicing, resampling, cloning, region removal and other methods. While resampling detection algorithms are effective in detecting splicing and resampling, copy-move detection algorithms excel in detecting cloning and region removal. In this paper, we combine these complementary approaches in a way that boosts the overall accuracy of image manipula… ▽ More Realistic image forgeries involve a combination of splicing, resampling, cloning, region removal and other methods. While resampling detection algorithms are effective in detecting splicing and resampling, copy-move detection algorithms excel in detecting cloning and region removal. In this paper, we combine these complementary approaches in a way that boosts the overall accuracy of image manipulation detection. We use the copy-move detection method as a pre-filtering step and pass those images that are classified as untampered to a deep learning based resampling detection framework. Experimental results on various datasets including the 2017 NIST Nimble Challenge Evaluation dataset comprising nearly 10,000 pristine and tampered images shows that there is a consistent increase of 8%-10% in detection rates, when copy-move algorithm is combined with different resampling detection algorithms. △ Less

Submitted 19 February, 2018; v1 submitted 9 February, 2018; originally announced February 2018.

Deep neural network initialization with decision trees

Authors: K. D. Humbird, J. L. Peterson, R. G. McClarren

Abstract: In this work a novel, automated process for constructing and initializing deep feed-forward neural networks based on decision trees is presented. The proposed algorithm maps a collection of decision trees trained on the data into a collection of initialized neural networks, with the structures of the networks determined by the structures of the trees. The tree-informed initialization acts as a war… ▽ More In this work a novel, automated process for constructing and initializing deep feed-forward neural networks based on decision trees is presented. The proposed algorithm maps a collection of decision trees trained on the data into a collection of initialized neural networks, with the structures of the networks determined by the structures of the trees. The tree-informed initialization acts as a warm-start to the neural network training process, resulting in efficiently trained, accurate networks. These models, referred to as "deep jointly-informed neural networks" (DJINN), demonstrate high predictive performance for a variety of regression and classification datasets, and display comparable performance to Bayesian hyper-parameter optimization at a lower computational cost. By combining the user-friendly features of decision tree models with the flexibility and scalability of deep neural networks, DJINN is an attractive algorithm for training predictive models on a wide range of complex datasets. △ Less

Submitted 2 July, 2018; v1 submitted 3 July, 2017; originally announced July 2017.

Detection and Localization of Image Forgeries using Resampling Features and Deep Learning

Authors: Jason Bunk, Jawadul H. Bappy, Tajuddin Manhar Mohammed, Lakshmanan Nataraj, Arjuna Flenner, B. S. Manjunath, Shivkumar Chandrasekaran, Amit K. Roy-Chowdhury, Lawrence Peterson

Abstract: Resampling is an important signature of manipulated images. In this paper, we propose two methods to detect and localize image manipulations based on a combination of resampling features and deep learning. In the first method, the Radon transform of resampling features are computed on overlapping image patches. Deep learning classifiers and a Gaussian conditional random field model are then used t… ▽ More Resampling is an important signature of manipulated images. In this paper, we propose two methods to detect and localize image manipulations based on a combination of resampling features and deep learning. In the first method, the Radon transform of resampling features are computed on overlapping image patches. Deep learning classifiers and a Gaussian conditional random field model are then used to create a heatmap. Tampered regions are located using a Random Walker segmentation method. In the second method, resampling features computed on overlapping image patches are passed through a Long short-term memory (LSTM) based network for classification and localization. We compare the performance of detection/localization of both these methods. Our experimental results show that both techniques are effective in detecting and localizing digital image forgeries. △ Less

Submitted 3 July, 2017; originally announced July 2017.

An Enhanced Nonlinear Critical Gradient for Electron Turbulent Transport due to Reversed Magnetic Shear

Authors: J. L. Peterson, G. W. Hammett, D. R. Mikkelsen, H. Y. Yuh, J. Candy, W. Guttenfelder, S. M. Kaye, B. LeBlanc

Abstract: The first nonlinear gyrokinetic simulations of electron internal transport barriers (e-ITBs) in the National Spherical Torus Experiment show that reversed magnetic shear can suppress thermal transport by increasing the nonlinear critical gradient for electron-temperature-gradient-driven turbulence to three times its linear critical value. An interesting feature of this turbulence is nonlinearly dr… ▽ More The first nonlinear gyrokinetic simulations of electron internal transport barriers (e-ITBs) in the National Spherical Torus Experiment show that reversed magnetic shear can suppress thermal transport by increasing the nonlinear critical gradient for electron-temperature-gradient-driven turbulence to three times its linear critical value. An interesting feature of this turbulence is nonlinearly driven off-midplane radial streamers. This work reinforces the experimental observation that magnetic shear is likely an effective way of triggering and sustaining e-ITBs in magnetic fusion devices. △ Less

Submitted 11 May, 2011; originally announced May 2011.

Comments: 4 pages, 5 figures

Report number: PPPL--4621

Future Directions of Research in Geometry: A Summary of the Panel Discussion at the 2007 Midwest Geometry Conference

Authors: Lawrence J. Peterson

Abstract: The 2007 Midwest Geometry Conference included a panel discussion devoted to open problems and the general direction of future research in fields related to the main themes of the conference. This paper summarizes the comments made during the panel discussion. The 2007 Midwest Geometry Conference included a panel discussion devoted to open problems and the general direction of future research in fields related to the main themes of the conference. This paper summarizes the comments made during the panel discussion. △ Less

Submitted 16 August, 2007; originally announced August 2007.

Comments: This is a contribution to the Proceedings of the 2007 Midwest Geometry Conference in honor of Thomas P. Branson, published in SIGMA (Symmetry, Integrability and Geometry: Methods and Applications) at http://www.emis.de/journals/SIGMA/

Journal ref: SIGMA 3 (2007), 081, 7 pages

The ambient obstruction tensor and the conformal deformation complex

Authors: A. Rod Gover, Lawrence J. Peterson

Abstract: The (Fefferman-Graham) ambient obstruction tensor is a conformally invariant symmetric trace-free 2-tensor on even-dimensional Riemannian and pseudo-Riemannian manifolds. The conformal deformation complex is a differential complex related to infinitesimal deformations of conformal structure. We construct a conformally invariant differential operator on algebraic Weyl tensors that gives special c… ▽ More The (Fefferman-Graham) ambient obstruction tensor is a conformally invariant symmetric trace-free 2-tensor on even-dimensional Riemannian and pseudo-Riemannian manifolds. The conformal deformation complex is a differential complex related to infinitesimal deformations of conformal structure. We construct a conformally invariant differential operator on algebraic Weyl tensors that gives special curved analogues of certain operators related to the deformation complex and that, upon application to the Weyl curvature, yields the obstruction tensor. This new definition of the obstruction tensor leads to simple direct proofs that the obstruction tensor is divergence free and vanishes for a class of metrics which (strictly) includes the class of conformally Einstein metrics. Our main constructions are based on the ambient metric of Fefferman-Graham and its relation to the conformal tractor connection. We prove that the obstruction tensor is an obstruction to the formal problem of obtaining a Ricci-flat ambient metric and that it may equivalently be viewed as the obstruction to finding an ambient metric with curvature harmonic for a certain (ambient) form Laplacian. This leads to a new ambient formula for the obstruction in terms of a power of this form Laplacian acting on the ambient curvature. This result leads us to construct Laplacian type operators that generalise the conformal Laplacians of Graham-Jenne-Mason-Sparling. As background to these issues, we give an explicit construction of the deformation complex in dimensions n>3, construct two related (detour) complexes, and establish essential properties of the operators in these. △ Less

Submitted 19 May, 2005; v1 submitted 17 August, 2004; originally announced August 2004.

Comments: 40 pages. No figures. Record of changes: V1, 18 August 2004; V2, 9 November 2004, adding a reference and correcting some typographical errors; V3, 19 May 2005, shortening article, updating references. Version V3 accepted by Pacific Journal of Mathematics

MSC Class: 53A55 (Primary) 22E70; 53A30; 58J10 (Secondary)

Journal ref: Pacific J. Math. 226 (2006) no. 2, 309-351.

Statistical randomization test for QCD intermittency in a single-event distribution

Authors: Leif E. Peterson

Abstract: A randomization test was developed to determine the statistical significance of QCD intermittency in single-event distributions. A total of 96 simulated intermittent distributions based on standard normal Gaussian distributions of size N=500, 1000, 1500, 2000, 4000, 8000, 16000, and 32000 containing induced holes and spikes were tested for intermittency. Non-intermittent null distributions were… ▽ More A randomization test was developed to determine the statistical significance of QCD intermittency in single-event distributions. A total of 96 simulated intermittent distributions based on standard normal Gaussian distributions of size N=500, 1000, 1500, 2000, 4000, 8000, 16000, and 32000 containing induced holes and spikes were tested for intermittency. Non-intermittent null distributions were also simulated as part of the test. A log-linear model was developed to simultaneously test the significance of fit coefficients for the y-intercept and slope contribution to ln(F_2) vs. ln(M) from both the intermittent and null distributions. Statistical power was also assessed for each fit coefficient to reflect the proportion of times out of 1000 tests each coefficient was statistically significant, given the induced effect size and sample size of the Gaussians. Results indicate that the slope of ln(F_2) vs. ln(M) for intermittent distributions increased with decreasing sample size, due to artificially-induced holes occurring in sparse histograms. For intermittent Gaussians with 4000 variates, there was approximately 70% power to detect a slope difference of 0.02 between intermittent and null distributions. For sample sizes of 8000 and greater, there was more than 70% power to detect a slope difference of 0.01. The randomization test performed satisfactorily since the power of the test for intermittency decreased with decreasing sample size. Power was near-zero when the test was applied to null distributions. The randomization test can be used to establish the statistical significance of intermittency in empirical single-event Gaussian distributions. △ Less

Submitted 5 April, 2004; originally announced April 2004.

Comments: 23 pages, 6 figures, 3 table

Conformally invariant powers of the Laplacian, Q-curvature, and tractor calculus

Authors: A. Rod Gover, Lawrence J. Peterson

Abstract: We describe an elementary algorithm for expressing, as explicit formulae in tractor calculus, the conformally invariant GJMS operators due to C.R. Graham et alia. These differential operators have leading part a power of the Laplacian. Conformal tractor calculus is the natural induced bundle calculus associated to the conformal Cartan connection. Applications discussed include standard formulae… ▽ More We describe an elementary algorithm for expressing, as explicit formulae in tractor calculus, the conformally invariant GJMS operators due to C.R. Graham et alia. These differential operators have leading part a power of the Laplacian. Conformal tractor calculus is the natural induced bundle calculus associated to the conformal Cartan connection. Applications discussed include standard formulae for these operators in terms of the Levi-Civita connection and its curvature and a direct definition and formula for T. Branson's so-called Q-curvature (which integrates to a global conformal invariant) as well as generalisations of the operators and the Q-curvature. Among examples, the operators of order 4, 6, and 8 and the related Q-curvatures are treated explicitly. The algorithm exploits the ambient metric construction of Fefferman and Graham and includes a procedure for converting the ambient curvature and its covariant derivatives into tractor calculus expressions. This is partly based on "Standard tractors and the conformal ambient metric construction" (A. Cap and A.R. Gover, math.DG/0207016), where the relationship of the normal standard tractor bundle to the ambient construction is described. △ Less

Submitted 1 February, 2003; v1 submitted 15 January, 2002; originally announced January 2002.

Comments: 42 pages. No figures. Record of changes: V1, 15 January 2002: Original posting. V2, 17 January 2002: Changing comment fields. Leaving abstract and text of article unchanged. V3, 1 February 2003: Minor changes and typographical corrections throughout article

MSC Class: 53B50 (Primary) 53A30; 53A55; 58J70 (Secondary)

Journal ref: Comm. Math. Phys. 235 (2003), no. 2, 339-378.

HEXTE Studies of Sco X-1 Spectra: Detections of Hard X-Ray Tails Beyond 200 keV

Authors: F. D'Amico, W. A. Heindl, R. E. Rothschild, L. E. Peterson, D. E. Gruber, M. Pelling, J. A. Tomsick

Abstract: Using the HEXTE experiment on-board the RXTE satellite, we performed a search for hard X-ray tails in Sco X-1 spectra. We found strong evidence for the presence of such a non-thermal component on several occasions. Using the PCA/RXTE we were able to track the position of the source along the Z diagram, and we observed that the presence of the hard X-ray tail is not confined to a particular regio… ▽ More Using the HEXTE experiment on-board the RXTE satellite, we performed a search for hard X-ray tails in Sco X-1 spectra. We found strong evidence for the presence of such a non-thermal component on several occasions. Using the PCA/RXTE we were able to track the position of the source along the Z diagram, and we observed that the presence of the hard X-ray tail is not confined to a particular region. However, we found a correlation between the power law index of the non-thermal component and the position of the source in the Z diagram, suggesting that the hard X-ray spectrum (i.e., E > 50 keV) becomes flatter as the mass accretion rate increases. We were also able to study the temporal variation of the appearance/absence of the hard X-ray component. With our derived luminosities, we were also able to test the idea that X-ray luminosities can be used to distinguish between X-ray binary systems containing neutron stars and black holes. △ Less

Submitted 22 January, 2001; originally announced January 2001.

Comments: 6 pages, 2 figures; Presented at the 33rd COSPAR Scientific Assembly (Warsaw, Poland, 2000 July); Accepted for publication in Advances in Space Research

The Spectrum of Diffuse Cosmic Hard X-Rays Measured with HEAO-1

Authors: D. E. Gruber, J. L. Matteson, L. E. Peterson, G. V. Jung

Abstract: The spectrum of the diffuse isotropic component of cosmic X-rays over the 13-180 keV range was determined by the UCSD/MIT Hard X-Ray and Gamma-Ray instrument (HEAO A4) on the High Energy Astronomical Observatory-1 (HEAO-1). The instrument consists of a complex of actively shielded and collimated scintillation counters, including the Low Energy Detector set from which the data reported here were… ▽ More The spectrum of the diffuse isotropic component of cosmic X-rays over the 13-180 keV range was determined by the UCSD/MIT Hard X-Ray and Gamma-Ray instrument (HEAO A4) on the High Energy Astronomical Observatory-1 (HEAO-1). The instrument consists of a complex of actively shielded and collimated scintillation counters, including the Low Energy Detector set from which the data reported here were obtained. These data join smoothly with the spectrum at lower energies reported by the GSFC HEAO A2 instrument and with that measured to 400 keV by the HEAO A4 Medium Energy Detectors. The HEAO data set also joins the recent results from COMPTEL on the Compton Gamma-Ray Observatory in the 1-10 MeV range, which failed to confirm the existence of an "MeV bump" in this range. Although the spectrum over the entire range 3 keV < E < 100 GeV can be fit by a simple empirical analytic expression, the origin is likely due to a number of distinct source components. The prevailing idea for the origin is that the hard X-ray spectrum is due to X-rays from various AGN components, particularly Seyfert galaxies extending to cosmological distances, and that the low energy gamma-rays may be due to emission from type 1a Supernovae, also integrated to cosmological distances. The higher energy gamma-ray spectrum defined by EGRET, also on the CGRO, may be due to unresolved gamma-ray emitting blazars. Models of production by these source components, extrapolated to the present epoch, must reproduce the observationally derived spectrum. △ Less

Submitted 31 March, 1999; originally announced March 1999.

Comments: LaTeX, 16pp, 3figs

Report number: SP-98-25

The X-ray Spectrum of the Plerionic System PSR B1509-58/MSH 15-52

Authors: D. Marsden, P. R. Blanco, D. E. Gruber, W. A. Heindl, M. R. Pelling, L. E. Peterson, R. E. Rothschild, A. H. Rots, K. Jahoda, D. J. Macomb

Abstract: We present the results of observations of the PSR B1509$-$58/MSH 15$-$52 system in X-rays ($2-250$ keV) by the Rossi X-ray Timing Explorer. The spectra of the peak of the pulsed component (radio phase $0.17-0.53$) is fit by a power law of photon index $1.36\pm0.01$, with no evidence of a high energy spectral break seen up to $\sim200$ keV. For the off-pulse spectral component, the spectrum from… ▽ More We present the results of observations of the PSR B1509$-$58/MSH 15$-$52 system in X-rays ($2-250$ keV) by the Rossi X-ray Timing Explorer. The spectra of the peak of the pulsed component (radio phase $0.17-0.53$) is fit by a power law of photon index $1.36\pm0.01$, with no evidence of a high energy spectral break seen up to $\sim200$ keV. For the off-pulse spectral component, the spectrum from $2-250$ keV is fit by a power law of photon index $2.215\pm0.005$. An iron emission line at 6.7 keV with an equivalent width of 129 eV improves the fit, but only at a marginal significance. Thermal bremsstrahlung and Raymond-Smith models produce much worse fits to the unpulsed data. The lack of a high energy spectral break in the pulsed emission implies an efficiency of $\geq 3%$ in the conversion of pulsar spindown energy to pulsed X-rays in the system. △ Less

Submitted 10 October, 1997; originally announced October 1997.

Comments: Accepted by ApJ Letters, 18 pages, 4 postscript figures