Showing 1–50 of 182 results for author: Schroder, B

Parallel-in-time solution of hyperbolic PDE systems via characteristic-variable block preconditioning

Authors: H. De Sterck, R. D. Falgout, O. A. Krzysik, J. B. Schroder

Abstract: We consider the parallel-in-time solution of hyperbolic partial differential equation (PDE) systems in one spatial dimension, both linear and nonlinear. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner appli… ▽ More We consider the parallel-in-time solution of hyperbolic partial differential equation (PDE) systems in one spatial dimension, both linear and nonlinear. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that inter-variable coupling for characteristic variables is weak relative to intra-variable coupling, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small. For an $\ell$-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of $\ell$ scalar linear(ized)-advection-like problems, each being associated with a different characteristic wave-speed in the underlying linear(ized) PDE. We approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media, and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Reconfigurable Intelligent Surfaces for 6G Mobile Networks: An Industry R&D Perspective

Authors: Maik Sode, Michael Ponschab, Lucas N. Ribeiro, Sven Haesloop, Ehsan Tohidi, Michael Peter, Sławomir Stańczak, Bilal H. Mohamed, Wilhelm Keusgen, Heinz Mellein, Eslam Yassin, Bernd Schroeder

Abstract: The reconfigurable intelligent surface (RIS) technology is a potential solution to enhance network capacity and coverage without significant investment into additional infrastructure in 6G networks. This work highlights the interest of the mobile communication industry on RIS, and discusses the development of liquid crystal-based RIS for improved energy efficiency and coverage in the millimeter-wa… ▽ More The reconfigurable intelligent surface (RIS) technology is a potential solution to enhance network capacity and coverage without significant investment into additional infrastructure in 6G networks. This work highlights the interest of the mobile communication industry on RIS, and discusses the development of liquid crystal-based RIS for improved energy efficiency and coverage in the millimeter-wave band. Furthermore, the paper discusses perspectives and insights from an industry R\&D point of view, addressing relevant use cases, technical requirements, implementation challenges, and practical considerations for RIS deployment optimization in the context of 6G networks. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Scaling Properties of Liquid Dynamics Predicted from a Single Configuration: Pseudoisomorphs for Harmonic-Bonded Molecules

Authors: Zahraa Sheydaafar, Jeppe C. Dyre, Thomas B. Schrøder

Abstract: Isomorphs are curves in the thermodynamic phase diagram of invariant excess entropy, structure, and dynamics, while pseudoisomorphs are curves of invariant structure and dynamics, but not of the excess entropy. The latter curves have been shown to exist in molecular models with flexible bonds [Olsen et al., J. Chem. Phys. 145, 241103 (2016)]. We here present three force-based methods to trace out… ▽ More Isomorphs are curves in the thermodynamic phase diagram of invariant excess entropy, structure, and dynamics, while pseudoisomorphs are curves of invariant structure and dynamics, but not of the excess entropy. The latter curves have been shown to exist in molecular models with flexible bonds [Olsen et al., J. Chem. Phys. 145, 241103 (2016)]. We here present three force-based methods to trace out pseudoisomorphs based on a single configuration and test them on the asymmetric dumbbell and 10-bead Lennard-Jones chain models with bonds modeled as harmonic springs. The three methods are based on requiring that particle forces, center-of-mass forces, and torques, respectively, are invariant in reduced units. For each of the two investigated models we identify a method that works well for tracing out pseudoisomorphs. △ Less

Submitted 16 June, 2024; originally announced June 2024.

Resonant emittance mixing of flat beams in plasma accelerators

Authors: Severin Diederichs, Carlo Benedetti, Angel Ferran Pousa, Alexander Sinn, Jens Osterhoff, Carl B. Schroeder, Maxence Thévenet

Abstract: Plasma accelerators sustain large field gradients and could enable future compact linear colliders. To achieve the required high luminosity, linear colliders rely on flat beams to avoid potentially deleterious beamstrahlung effects. Here, we show that flat beams in plasma accelerators can be subject to beam quality degradation due to emittance mixing caused by transverse coupling in the wakefields… ▽ More Plasma accelerators sustain large field gradients and could enable future compact linear colliders. To achieve the required high luminosity, linear colliders rely on flat beams to avoid potentially deleterious beamstrahlung effects. Here, we show that flat beams in plasma accelerators can be subject to beam quality degradation due to emittance mixing caused by transverse coupling in the wakefields. When there is a resonance between the betatron oscillations in the horizontal and vertical planes for the beam particles in a coupled wakefield, the transverse emittances fully exchange, leading to a round beam. Depending on the mechanism causing the resonance, the use of laser drivers, flat particle beam drivers, or hollow plasma channels can avoid the resonance and mitigate the emittance deterioration. △ Less

Submitted 9 March, 2024; originally announced March 2024.

Generalized Optimal AMG Convergence Theory for Nonsymmetric and Indefinite Problems

Authors: Ahsan Ali, James Brannick, Karsten Kahl, Oliver A. Krzysik, Jacob B. Schroder, Ben S. Southworth

Abstract: Algebraic multigrid (AMG) is known to be an effective solver for many sparse symmetric positive definite (SPD) linear systems. For SPD systems, the convergence theory of AMG is well-understood in terms of the $A$-norm but in a nonsymmetric setting such an energy norm is non-existent. For this reason, convergence of AMG for nonsymmetric systems of equations remains an open area of research. Existin… ▽ More Algebraic multigrid (AMG) is known to be an effective solver for many sparse symmetric positive definite (SPD) linear systems. For SPD systems, the convergence theory of AMG is well-understood in terms of the $A$-norm but in a nonsymmetric setting such an energy norm is non-existent. For this reason, convergence of AMG for nonsymmetric systems of equations remains an open area of research. Existing nonsymmetric AMG algorithms in this setting mostly rely on heuristics motivated by SPD convergence theory. In the SPD setting, the classical form of optimal AMG interpolation provides a useful insight in determining the two grid convergence rate of the method. In this work, we discuss a generalization of the optimal AMG convergence theory targeting nonsymmetric problems by constructing a $2\times 2$ block symmetric indefinite system so that the Petrov-Galerkin AMG process for the nonsymmetric matrix $A$ can be recast as a Galerkin AMG process for a symmetric indefinite system. We show that using this generalization of the optimal interpolation theory, one can obtain the same identity for the two-grid convergence rate as that derived in the SPD setting for optimal interpolation. We also provide supporting numerical results for the convergence result and nonsymmetric advection-diffusion problems. △ Less

Submitted 20 January, 2024; originally announced January 2024.

Comments: 10 pages, 2 figures, submitted as a student paper for the 18th Copper Mountain Conference on Iterative Methods

Parallel-in-time solution of scalar nonlinear conservation laws

Authors: H. De Sterck, R. D. Falgout, O. A. Krzysik, J. B. Schroder

Abstract: We consider the parallel-in-time solution of scalar nonlinear conservation laws in one spatial dimension. The equations are discretized in space with a conservative finite-volume method using weighted essentially non-oscillatory (WENO) reconstructions, and in time with high-order explicit Runge-Kutta methods. The solution of the global, discretized space-time problem is sought via a nonlinear iter… ▽ More We consider the parallel-in-time solution of scalar nonlinear conservation laws in one spatial dimension. The equations are discretized in space with a conservative finite-volume method using weighted essentially non-oscillatory (WENO) reconstructions, and in time with high-order explicit Runge-Kutta methods. The solution of the global, discretized space-time problem is sought via a nonlinear iteration that uses a novel linearization strategy in cases of non-differentiable equations. Under certain choices of discretization and algorithmic parameters, the nonlinear iteration coincides with Newton's method, although, more generally, it is a preconditioned residual correction scheme. At each nonlinear iteration, the linearized problem takes the form of a certain discretization of a linear conservation law over the space-time domain in question. An approximate parallel-in-time solution of the linearized problem is computed with a single multigrid reduction-in-time (MGRIT) iteration. The MGRIT iteration employs a novel coarse-grid operator that is a modified conservative semi-Lagrangian discretization and generalizes those we have developed previously for non-conservative scalar linear hyperbolic problems. Numerical tests are performed for the inviscid Burgers and Buckley--Leverett equations. For many test problems, the solver converges in just a handful of iterations with convergence rate independent of mesh resolution, including problems with (interacting) shocks and rarefactions. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Reconstruction of the Ranks of the Nonextremal Cards and of Ordered Sets with a Minmax Pair of Pseudo-Similar Points

Authors: Bernd S. W. Schröder

Abstract: For every ordered set, we reconstruct the deck obtained by removal of the elements of rank r that are neither minimal nor maximal. Consequently, we also reconstruct the deck obtained by removal of the extremal, that is, minimal or maximal, elements. Finally, we reconstruct the ordered sets with a minmax pair of pseudo-similar points. For every ordered set, we reconstruct the deck obtained by removal of the elements of rank r that are neither minimal nor maximal. Consequently, we also reconstruct the deck obtained by removal of the extremal, that is, minimal or maximal, elements. Finally, we reconstruct the ordered sets with a minmax pair of pseudo-similar points. △ Less

Submitted 29 January, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

A high-intensity laser-based positron source

Authors: S. S. Bulanov, C. Benedetti, D. Terzani, C. B. Schroeder, E. Esarey, T. Blackburn, M. Marklund

Abstract: Plasma based acceleration is considered a promising concept for the next generation of linear electron-positron colliders. Despite the great progress achieved over last twenty years in laser technology, laser and beam driven particle acceleration, and special target availability, positron acceleration remains significantly underdeveloped if compared to electron acceleration. This is due to both th… ▽ More Plasma based acceleration is considered a promising concept for the next generation of linear electron-positron colliders. Despite the great progress achieved over last twenty years in laser technology, laser and beam driven particle acceleration, and special target availability, positron acceleration remains significantly underdeveloped if compared to electron acceleration. This is due to both the specifics of the plasma-based acceleration, and the lack of adequate positron sources tailored for the subsequent plasma based acceleration. Here a positron source based on the collision of a high energy electron beam with a high intensity laser pulse is proposed. The source relies on the subsequent multi-photon Compton and Breit-Wheeleer processes to generate an electron-positron pair out of a high energy photon emitted by an electron. Due to the strong dependence of the Breit-Wheeler process rate on photon energy and field strength, positrons are created with low divergence in a small volume around the peak of the laser pulse. The resulting low emittance in the submicron range potentially makes such positron source interesting for collider applications. △ Less

Submitted 17 November, 2023; originally announced November 2023.

Comments: 28 pages, 10 figures, 2 tables

Emittance-preserving acceleration of high-quality positron beams using warm plasma filaments

Authors: Severin Diederichs, Carlo Benedetti, Eric Esarey, Alexander Sinn, Jens Osterhoff, Carl B. Schroeder, Maxence Thévenet

Abstract: Preserving the quality of positron beams in plasma-based accelerators, where wakefields are generated in electron filaments, is challenging. These wakefields are characterized by transversely non-linear focusing fields and non-uniform accelerating fields. However, a nonzero plasma temperature linearizes the transverse wakefield within the central region of the electron filament. In this study, we… ▽ More Preserving the quality of positron beams in plasma-based accelerators, where wakefields are generated in electron filaments, is challenging. These wakefields are characterized by transversely non-linear focusing fields and non-uniform accelerating fields. However, a nonzero plasma temperature linearizes the transverse wakefield within the central region of the electron filament. In this study, we employ 3D particle-in-cell simulations with mesh refinement to demonstrate that beams with emittances on the order of tens of nanometers are contained within the linearized region of the transverse wakefield. This enables emittance preservation to one percent, while positron beams with the same charge and micrometer emittances, which sample the non-linear part of the transverse wakefield, experience a relative emittance growth of ten percent. Additionally, we observe a significant reduction in the growth rate of the slice energy spread for the tens of nanometers emittance beams in comparison to the micrometer emittance beams. The utilization of warm plasmas in conjunction with low-emittance beams opens up new avenues for enhancing the beam quality across various plasma-based positron acceleration approaches. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: To be submitted as a proceedings for the 6th European Advanced Accelerator Concepts workshop

A compact, all-optical positron production and collection scheme

Authors: Davide Terzani, Carlo Benedetti, Stepan S. Bulanov, Carl B. Schroeder, Eric Esarey

Abstract: In this paper we discuss a compact, laser-plasma-based scheme for the generation of positron beams suitable to be implemented in an all-optical setup. A laser-plasma-accelerated electron beam hits a solid target producing electron-positron pairs via bremsstrahlung. The back of the target serves as a plasma mirror to in-couple a laser pulse into a plasma stage located right after the mirror where t… ▽ More In this paper we discuss a compact, laser-plasma-based scheme for the generation of positron beams suitable to be implemented in an all-optical setup. A laser-plasma-accelerated electron beam hits a solid target producing electron-positron pairs via bremsstrahlung. The back of the target serves as a plasma mirror to in-couple a laser pulse into a plasma stage located right after the mirror where the laser drives a plasma wave (or wakefield). By properly choosing the delay between the laser and the electron beam the positrons produced in the target can be trapped in the wakefield, where they are focused and accelerated during the transport, resulting in a collimated beam. This approach minimizes the ballistic propagation time and enhances the trapping efficiency. The system can be used as an injector of positron beams and has potential applications in the development of a future, compact, plasma-based electron-positron linear collider. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Space-Time Block Preconditioning for Incompressible Resistive Magnetohydrodynamics

Authors: Federico Danieli, Ben S. Southworth, Jacob B. Schroder

Abstract: This work develops a novel all-at-once space-time preconditioning approach for resistive magnetohydrodynamics (MHD), with a focus on model problems targeting fusion reactor design. We consider parallel-in-time due to the long time domains required to capture the physics of interest, as well as the complexity of the underlying system and thereby computational cost of long-time integration. To ameli… ▽ More This work develops a novel all-at-once space-time preconditioning approach for resistive magnetohydrodynamics (MHD), with a focus on model problems targeting fusion reactor design. We consider parallel-in-time due to the long time domains required to capture the physics of interest, as well as the complexity of the underlying system and thereby computational cost of long-time integration. To ameliorate this cost by using many processors, we thus develop a novel approach to solving the whole space-time system that is parallelizable in both space and time. We develop a space-time block preconditioning for resistive MHD, following the space-time block preconditioning concept first introduced by Danieli et al. in 2022 for incompressible flow, where an effective preconditioner for classic sequential time-stepping is extended to the space-time setting. The starting point for our derivation is the continuous Schur complement preconditioner by Cyr et al. in 2021, which we proceed to generalise in order to produce, to our knowledge, the first space-time block preconditioning approach for the challenging equations governing incompressible resistive MHD. The numerical results are promising for the model problems of island coalescence and tearing mode, with the overhead computational cost associated with space-time preconditioning versus sequential time-stepping being modest and primarily in the range of 2x-5x, which is low for parallel-in-time schemes in general. Additionally, the scaling results for inner (linear) and outer (nonlinear) iterations are flat in the case of fixed time-step size and only grow very slowly in the case of time-step refinement. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 25 pages, 4 figures, 3 tables

MSC Class: 65F08; 65Y05; 76W05; 65M60; 65M22

Small Representations of Permutation Groups by Embedding the Domain in an Ordered Set

Authors: Bernd Schröder

Abstract: We present a representation for permutation groups as the automorphism group of an ordered set $U$ such that the automorphism group's action on a subset $T\subseteq U$ is the permutation group itself. For many imprimitive permutation groups, the representation uses fewer points than the smallest representations to date. We present a representation for permutation groups as the automorphism group of an ordered set $U$ such that the automorphism group's action on a subset $T\subseteq U$ is the permutation group itself. For many imprimitive permutation groups, the representation uses fewer points than the smallest representations to date. △ Less

Submitted 15 September, 2023; v1 submitted 20 August, 2023; originally announced August 2023.

Towards the Automorphism Conjecture I: Combinatorial Control

Authors: Bernd S. W. Schröder

Abstract: This paper exploits adjacencies between the orbits of an ordered set P and a consequence of the classification of finite simple groups to, in many cases, exponentially bound the number of automorphisms. Results clearly identify the structures which currently prevent the proof of such an exponential bound, or which indeed inflate the number of automorphisms beyond such a bound. This is a first step… ▽ More This paper exploits adjacencies between the orbits of an ordered set P and a consequence of the classification of finite simple groups to, in many cases, exponentially bound the number of automorphisms. Results clearly identify the structures which currently prevent the proof of such an exponential bound, or which indeed inflate the number of automorphisms beyond such a bound. This is a first step towards a possible resolution of the Automorphism Conjecture for ordered sets. △ Less

Submitted 11 September, 2023; v1 submitted 10 August, 2023; originally announced August 2023.

Constrained Local Approximate Ideal Restriction for Advection-Diffusion Problems

Authors: Ahsan Ali, James Brannick, Karsten Kahl, Oliver A. Krzysik, Jacob B. Schroder, Ben S. Southworth

Abstract: This paper focuses on developing a reduction-based algebraic multigrid method that is suitable for solving general (non)symmetric linear systems and is naturally robust from pure advection to pure diffusion. Initial motivation comes from a new reduction-based algebraic multigrid (AMG) approach, $\ell$AIR (local approximate ideal restriction), that was developed for solving advection-dominated prob… ▽ More This paper focuses on developing a reduction-based algebraic multigrid method that is suitable for solving general (non)symmetric linear systems and is naturally robust from pure advection to pure diffusion. Initial motivation comes from a new reduction-based algebraic multigrid (AMG) approach, $\ell$AIR (local approximate ideal restriction), that was developed for solving advection-dominated problems. Though this new solver is very effective in the advection dominated regime, its performance degrades in cases where diffusion becomes dominant. This is consistent with the fact that in general, reduction-based AMG methods tend to suffer from growth in complexity and/or convergence rates as the problem size is increased, especially for diffusion dominated problems in two or three dimensions. Motivated by the success of $\ell$AIR in the advective regime, our aim in this paper is to generalize the AIR framework with the goal of improving the performance of the solver in diffusion dominated regimes. To do so, we propose a novel way to combine mode constraints as used commonly in energy minimization AMG methods with the local approximation of ideal operators used in $\ell$AIR. The resulting constrained $\ell$AIR (C$\ell$AIR) algorithm is able to achieve fast scalable convergence on advective and diffusive problems. In addition, it is able to achieve standard low complexity hierarchies in the diffusive regime through aggressive coarsening, something that has been previously difficult for reduction-based methods. △ Less

Submitted 14 May, 2024; v1 submitted 1 July, 2023; originally announced July 2023.

Comments: Revised form published in the SIAM Journal on Scientific computing (SPECIAL SECTION Copper Mountain 2023), 22 pages, 7 Figures

MSC Class: 65N55 (Primary); 65N22; 65F08; 65F10 (Secondary)

Journal ref: SIAM J. SCI. COMPUT. 2024 Vol. 0, No. 0, pp. S96-S122

Evolution of Relativistic Pair Beams: Implications for Laboratory and TeV Astrophysics

Authors: Marvin Beck, Oindrila Ghosh, Florian Grüner, Martin Pohl, Carl B. Schroeder, Günter Sigl, Ryan D. Stark, Benno Zeitler

Abstract: Missing cascades from TeV blazar beams indicate that collective plasma effects may play a significant role in their energy loss. It is possible to mimic the evolution of such highly energetic pair beams in laboratory experiments using modern accelerators. The fate of the beam is governed by two different processes, energy loss through the unstable mode and energetic broadening of the pair beam thr… ▽ More Missing cascades from TeV blazar beams indicate that collective plasma effects may play a significant role in their energy loss. It is possible to mimic the evolution of such highly energetic pair beams in laboratory experiments using modern accelerators. The fate of the beam is governed by two different processes, energy loss through the unstable mode and energetic broadening of the pair beam through diffusion in momentum space. We chalk out this evolution using a Fokker-Planck approach in which the drift and the diffusion terms respectively describe these phenomena in a compact form. We present particle-in-cell simulations to trace the complete evolution of the unstable beam-plasma system for a generic narrow Gaussian pair beam for which the growth rate is reactive. We show that the instability leads to an energetic broadening of the pair beam, slowing down the instability growth in the linear phase, in line with the analytical and numerical solutions of the Fokker-Planck equation. Whereas in a laboratory experiment the change in the momentum distribution is an easily measured observable as a feedback of the instability, the consequence of diffusive broadening in an astrophysical scenario can be translated to an increase in the opening angle of the pair beam. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: 20 pages, 10 figures

The Automorphism Conjecture for Ordered Sets of Width $\leq 11$ (Version 2)

Authors: Bernd S. W. Schröder

Abstract: We introduce a recursive method to deconstruct the automorphism group of an ordered set. By connecting this method with deep results for permutation groups, we prove the Automorphism Conjecture for ordered sets of width less than or equal to $11$. Subsequent investigations show that the method presented here could lead to a resolution of the Automorphism Conjecture. We introduce a recursive method to deconstruct the automorphism group of an ordered set. By connecting this method with deep results for permutation groups, we prove the Automorphism Conjecture for ordered sets of width less than or equal to $11$. Subsequent investigations show that the method presented here could lead to a resolution of the Automorphism Conjecture. △ Less

Submitted 22 May, 2023; originally announced May 2023.

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

Strong-Field QED Experiments using the BELLA PW Laser Dual Beamlines

Authors: M. Turner, S. S. Bulanov, C. Benedetti, A. J. Gonsalves, W. P. Leemans, K. Nakamura, J. van Tilborg, C. B. Schroeder, C. G. R. Geddes, E. Esarey

Abstract: The Petawatt (PW) laser facility of the Berkeley Lab Laser Accelerator (BELLA) Center has recently commissioned its second laser pulse transport line. This new beamline can be operated in parallel with the first beamline and enables strong-field quantum electrodynamics (SF-QED) experiments at BELLA. In this paper, we present an overview of the upgraded BELLA PW facility with a SF-QED experimental… ▽ More The Petawatt (PW) laser facility of the Berkeley Lab Laser Accelerator (BELLA) Center has recently commissioned its second laser pulse transport line. This new beamline can be operated in parallel with the first beamline and enables strong-field quantum electrodynamics (SF-QED) experiments at BELLA. In this paper, we present an overview of the upgraded BELLA PW facility with a SF-QED experimental layout in which intense laser pulses collide with GeV-class laser-wakefield-accelerated electron beams. We present simulation results showing that experiments will allow the study of laser-particle interactions from the classical to the SF-QED regime with a nonlinear quantum parameter of up to $χかい\sim$2. In addition, we show that experiments will enable the study and production of GeV-class, mrad-divergence positron beams via the Breit-Wheeler process. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 5 figures

The Automorphism Conjecture for Ordered Sets of Width $\leq 11$

Authors: Bernd Schröder

Abstract: We prove the automorphism conjecture for ordered sets of width less than or equal to 11. The proof supports the meta conjecture that a large number of automorphisms is achievable only as some type of product of independent automorphisms on highly symmetric subsets. We prove the automorphism conjecture for ordered sets of width less than or equal to 11. The proof supports the meta conjecture that a large number of automorphisms is achievable only as some type of product of independent automorphisms on highly symmetric subsets. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 57 pages, 6 figures

MSC Class: 06A07

Efficient multigrid reduction-in-time for method-of-lines discretizations of linear advection

Authors: H. De Sterck, R. D. Falgout, O. A. Krzysik, J. B. Schroder

Abstract: Parallel-in-time methods for partial differential equations (PDEs) have been the subject of intense development over recent decades, particularly for diffusion-dominated problems. It has been widely reported in the literature, however, that many of these methods perform quite poorly for advection-dominated problems. Here we analyze the particular iterative parallel-in-time algorithm of multigrid r… ▽ More Parallel-in-time methods for partial differential equations (PDEs) have been the subject of intense development over recent decades, particularly for diffusion-dominated problems. It has been widely reported in the literature, however, that many of these methods perform quite poorly for advection-dominated problems. Here we analyze the particular iterative parallel-in-time algorithm of multigrid reduction-in-time (MGRIT) for discretizations of constant-wave-speed linear advection problems. We focus on common method-of-lines discretizations that employ upwind finite differences in space and Runge-Kutta methods in time. Using a convergence framework we developed in previous work, we prove for a subclass of these discretizations that, if using the standard approach of rediscretizing the fine-grid problem on the coarse grid, robust MGRIT convergence with respect to CFL number and coarsening factor is not possible. This poor convergence and non-robustness is caused, at least in part, by an inadequate coarse-grid correction for smooth Fourier modes known as characteristic components.We propose an alternative coarse-grid that provides a better correction of these modes. This coarse-grid operator is related to previous work and uses a semi-Lagrangian discretization combined with an implicitly treated truncation error correction. Theory and numerical experiments show the coarse-grid operator yields fast MGRIT convergence for many of the method-of-lines discretizations considered, including for both implicit and explicit discretizations of high order. Parallel results demonstrate substantial speed-up over sequential time-stepping. △ Less

Submitted 20 March, 2023; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: v2 updates v1 with a small amount of rewriting and with include parallel timing results

High-dimensional neural network potentials for accurate vibrational frequencies: The formic acid dimer benchmark

Authors: Dilshana Shanavas Rasheeda, Alberto Martín Santa Daría, Benjamin Schröder, Edit Mátyus, Jörg Behler

Abstract: In recent years, machine learning potentials (MLP) for atomistic simulations have attracted a lot of attention in chemistry and materials science. Many new approaches have been developed with the primary aim to transfer the accuracy of electronic structure calculations to large condensed systems containing thousands of atoms. In spite of these advances, the reliability of modern MLPs in reproducin… ▽ More In recent years, machine learning potentials (MLP) for atomistic simulations have attracted a lot of attention in chemistry and materials science. Many new approaches have been developed with the primary aim to transfer the accuracy of electronic structure calculations to large condensed systems containing thousands of atoms. In spite of these advances, the reliability of modern MLPs in reproducing the subtle details of the multi-dimensional potential-energy surface is still difficult to assess for such systems. On the other hand, moderately sized systems enabling the application of tools for thorough and systematic quality-control are nowadays rarely investigated. In this work we use benchmark-quality harmonic and anharmonic vibrational frequencies as a sensitive probe for the validation of high-dimensional neural network potentials. For the case of the formic acid dimer, a frequently studied model system for which stringent spectroscopic data became recently available, we show that high-quality frequencies can be obtained from state-of-the-art calculations in excellent agreement with coupled cluster theory and experimental data. △ Less

Submitted 8 September, 2022; v1 submitted 7 September, 2022; originally announced September 2022.

Multigrid Reduction in Time for Chaotic Dynamical Systems

Authors: David A. Vargas, Robert D. Falgout, Stefanie Günther, Jacob B. Schroder

Abstract: As CPU clock speeds have stagnated and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes can be a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel in Time schemes, such as the… ▽ More As CPU clock speeds have stagnated and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes can be a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel in Time schemes, such as the Multigrid Reduction in Time (MGRIT) method, remedy this by parallelizing across time-steps, and have shown promising results for parabolic problems. However, chaotic problems have proved more difficult, since chaotic initial value problems (IVPs) are inherently ill-conditioned. MGRIT relies on a hierarchy of successively coarser time-grids to iteratively correct the solution on the finest time-grid, but due to the nature of chaotic systems, small inaccuracies on the coarser levels can be greatly magnified and lead to poor coarse-grid corrections. Here we introduce a modified MGRIT algorithm based on an existing quadratically converging nonlinear extension to the multigrid Full Approximation Scheme (FAS), as well as a novel time-coarsening scheme. Together, these approaches better capture long-term chaotic behavior on coarse-grids and greatly improve convergence of MGRIT for chaotic IVPs. Further, we introduce a novel low memory variant of the algorithm for solving chaotic PDEs with MGRIT which not only solves the IVP, but also provides estimates for the unstable Lyapunov vectors of the system. We provide supporting numerical results for the Lorenz system and demonstrate parallel speedup for the chaotic Kuramoto- Sivashinsky partial differential equation over a significantly longer time-domain than in previous works. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: 24 pages, 10 figures. arXiv admin note: text overlap with arXiv:2201.10441

Report number: LLNL-JRNL-838414 MSC Class: 65M22; 65M55

Parallel Energy-Minimization Prolongation for Algebraic Multigrid

Authors: Carlo Janna, Andrea Franceschini, Jacob B. Schroder, Luke Olson

Abstract: Algebraic multigrid (AMG) is one of the most widely used solution techniques for linear systems of equations arising from discretized partial differential equations. The popularity of AMG stems from its potential to solve linear systems in almost linear time, that is with an O(n) complexity, where n is the problem size. This capability is crucial at the present, where the increasing availability o… ▽ More Algebraic multigrid (AMG) is one of the most widely used solution techniques for linear systems of equations arising from discretized partial differential equations. The popularity of AMG stems from its potential to solve linear systems in almost linear time, that is with an O(n) complexity, where n is the problem size. This capability is crucial at the present, where the increasing availability of massive HPC platforms pushes for the solution of very large problems. The key for a rapidly converging AMG method is a good interplay between the smoother and the coarse-grid correction, which in turn requires the use of an effective prolongation. From a theoretical viewpoint, the prolongation must accurately represent near kernel components and, at the same time, be bounded in the energy norm. For challenging problems, however, ensuring both these requirements is not easy and is exactly the goal of this work. We propose a constrained minimization procedure aimed at reducing prolongation energy while preserving the near kernel components in the span of interpolation. The proposed algorithm is based on previous energy minimization approaches utilizing a preconditioned restricted conjugate gradients method, but has new features and a specific focus on parallel performance and implementation. It is shown that the resulting solver, when used for large real-world problems from various application fields, exhibits excellent convergence rates and scalability and outperforms at least some more traditional AMG approaches. △ Less

Submitted 5 August, 2022; originally announced August 2022.

Comments: 24 pages, 4 figures

MSC Class: 65N55; 68W10; 65F08

Improving the Reliability of Next Generation SSDs using WOM-v Codes

Authors: Shehbaz Jaffer, Kaveh Mahdaviani, Bianca Schroeder

Abstract: High density Solid State Drives, such as QLC drives, offer increased storage capacity, but a magnitude lower Program and Erase (P/E) cycles, limiting their endurance and hence usability. We present the design and implementation of non-binary, Voltage-Based Write-Once-Memory (WOM-v) Codes to improve the lifetime of QLC drives. First, we develop a FEMU based simulator test-bed to evaluate the gains… ▽ More High density Solid State Drives, such as QLC drives, offer increased storage capacity, but a magnitude lower Program and Erase (P/E) cycles, limiting their endurance and hence usability. We present the design and implementation of non-binary, Voltage-Based Write-Once-Memory (WOM-v) Codes to improve the lifetime of QLC drives. First, we develop a FEMU based simulator test-bed to evaluate the gains of WOM-v codes on real world workloads. Second, we propose and implement two optimizations, an efficient garbage collection mechanism and an encoding optimization to drastically improve WOM-v code endurance without compromising performance. A careful evaluation, including microbenchmarks and trace-driven evaluation, demonstrates that WOM-v codes can reduce Erase cycles for QLC drives by 4.4x-11.1x for real world workloads with minimal performance overheads resulting in improved QLC SSD lifetime. △ Less

Submitted 23 July, 2022; originally announced August 2022.

Comments: 15 pages, 13 Figures, Published at USENIX FAST'22

Journal ref: 20th USENIX Conference on File and Storage Technologies (FAST) 2022

Self-stabilizing positron acceleration in a plasma column

Authors: S. Diederichs, C. Benedetti, M. Thévenet, E. Esarey, J. Osterhoff, C. B. Schroeder

Abstract: Plasma accelerators sustain extreme field gradients, and potentially enable future compact linear colliders. Although tremendous progress has been achieved in accelerating electron beams in a plasma accelerator, positron acceleration with collider-relevant parameters is challenging. A recently proposed positron acceleration scheme relying on the wake generated by an electron drive beam in a plasma… ▽ More Plasma accelerators sustain extreme field gradients, and potentially enable future compact linear colliders. Although tremendous progress has been achieved in accelerating electron beams in a plasma accelerator, positron acceleration with collider-relevant parameters is challenging. A recently proposed positron acceleration scheme relying on the wake generated by an electron drive beam in a plasma column has been shown to be able to accelerate positron witness beams with low emittance and low energy spread. However, since this scheme relies on cylindrical symmetry, it is possibly prone to transverse instabilities that could lead, ultimately, to beam break-up. In this article, we show that the witness beam itself is subject to various damping mechanisms and, therefore, this positron acceleration scheme is inherently stable towards misalignment of the drive and witness beams. This enables stable, high-quality plasma-based positron acceleration. △ Less

Submitted 23 June, 2022; originally announced June 2022.

A toy model for viscous liquid dynamics

Authors: Filip Samuelsen, Lorenzo Costigliola, Thomas B. Schrøder

Abstract: A simple model for viscous liquid dynamics is introduced. Consider the surface of the union of hyper-spheres centered at random positions inside a hypercube with periodic boundary conditions. It is argued and demonstrated by numerical simulations that at high dimensions geodetic flows on this surface is a good model for viscous liquid dynamics. It is shown that this simple model exhibits viscous d… ▽ More A simple model for viscous liquid dynamics is introduced. Consider the surface of the union of hyper-spheres centered at random positions inside a hypercube with periodic boundary conditions. It is argued and demonstrated by numerical simulations that at high dimensions geodetic flows on this surface is a good model for viscous liquid dynamics. It is shown that this simple model exhibits viscous dynamics for densities above the percolation threshold in $8$, $12$ and $16$ dimensions. Thus the slowing down of the dynamics, measured by the mean-squared displacement, extends to several orders of magnitude similarly to what is observed in other models for viscous dynamics. Furthermore, the shape of the mean-squared displacement is to a very good approximation the same as for the standard model in simulations of viscous liquids: the Kob-Andersen binary Lennard Jones mixture. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 4 pages, 5 Figures

Online Charge Measurement for Petawatt Laser-Driven Ion Acceleration

Authors: Laura D. Geulig, Lieselotte Obst-Huebl, Kei Nakamura, Jianhui Bin, Qing Ji, Sven Steinke, Antoine Snijders, Jian-Hua Mao, Eleanor Blakely, Anthony J. Gonsalves, Stepan S. Bulanov, Jeroen van Tilborg, Carl B. Schroeder, Cameron G. R. Geddes, Eric Esarey, Markus Roth, Thomas Schenkel

Abstract: Laser-driven ion beams have gained considerable attention for their potential use in multidisciplinary research and technology. Pre-clinical studies into their radiobiological effectiveness have established the prospect of using laser-driven ion beams for radiotherapy. In particular, research into the beneficial effects of ultra-high instantaneous dose rates is enabled by the high ion bunch charge… ▽ More Laser-driven ion beams have gained considerable attention for their potential use in multidisciplinary research and technology. Pre-clinical studies into their radiobiological effectiveness have established the prospect of using laser-driven ion beams for radiotherapy. In particular, research into the beneficial effects of ultra-high instantaneous dose rates is enabled by the high ion bunch charge and uniquely short bunch lengths present for laser-driven ion beams. Such studies require reliable, online dosimetry methods to monitor the bunch charge for every laser shot to ensure that the prescribed dose is accurately applied to the biological sample. In this paper we present the first successful use of an Integrating Current Transformer (ICT) for laser-driven ion accelerators. This is a non-invasive diagnostic to measure the charge of the accelerated ion bunch. It enables online dose measurements in radiobiological experiments and facilitates ion beam tuning, in particular, optimization of the laser ion source and alignment of the proton transport beamline. We present the ICT implementation and the correlation with other diagnostics such as radiochromic films, a Thomson parabola spectrometer and a scintillator. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Towards Generalizable Semantic Product Search by Text Similarity Pre-training on Search Click Logs

Authors: Zheng Liu, Wei Zhang, Yan Chen, Weiyi Sun, Tianchuan Du, Benjamin Schroeder

Abstract: Recently, semantic search has been successfully applied to e-commerce product search and the learned semantic space(s) for query and product encoding are expected to generalize to unseen queries or products. Yet, whether generalization can conveniently emerge has not been thoroughly studied in the domain thus far. In this paper, we examine several general-domain and domain-specific pre-trained Rob… ▽ More Recently, semantic search has been successfully applied to e-commerce product search and the learned semantic space(s) for query and product encoding are expected to generalize to unseen queries or products. Yet, whether generalization can conveniently emerge has not been thoroughly studied in the domain thus far. In this paper, we examine several general-domain and domain-specific pre-trained Roberta variants and discover that general-domain fine-tuning does not help generalization, which aligns with the discovery of prior art. Proper domain-specific fine-tuning with clickstream data can lead to better model generalization, based on a bucketed analysis of a publicly available manual annotated query-product pair da △ Less

Submitted 28 April, 2022; v1 submitted 11 April, 2022; originally announced April 2022.

Comments: 9 pages. ACL 2022 ECNLP Workshop

ACM Class: F.2.2; I.2.7

Journal ref: ECNLP 2022

Whitepaper submitted to Snowmass21: Advanced accelerator linear collider demonstration facility at intermediate energy

Authors: C. Benedetti, S. S. Bulanov, E. Esarey, C. G. R. Geddes A. J. Gonsalves, P. M. Jacobs, S. Knapen, B. Nachman, K. Nakamura, S. Pagan Griso, C. B. Schroeder, D. Terzani, J. van Tilborg, M. Turner, W. -M. Yao, R. Bernstein, V. Shiltsev, S. J. Gessner, M. J. Hogan, T. Nelson, C. Jing, I. Low, X. Lu, R. Yoshida, C. Lee, P. Meade , et al. (8 additional authors not shown)

Abstract: It is widely accepted that the next lepton collider beyond a Higgs factory would require center-of-mass energy of the order of up to 15 TeV. Since, given reasonable space and cost restrictions, conventional accelerator technology reaches its limits near this energy, high-gradient advanced acceleration concepts are attractive. Advanced and novel accelerators (ANAs) are leading candidates due to the… ▽ More It is widely accepted that the next lepton collider beyond a Higgs factory would require center-of-mass energy of the order of up to 15 TeV. Since, given reasonable space and cost restrictions, conventional accelerator technology reaches its limits near this energy, high-gradient advanced acceleration concepts are attractive. Advanced and novel accelerators (ANAs) are leading candidates due to their ability to produce acceleration gradients on the order of 1--100~GV/m, leading to compact acceleration structures. Over the last 10-15 years significant progress has been achieved in accelerating electron beams by ANAs. For example, the demonstration of several-GeV electron beams from laser-powered capillary discharge waveguides, as well as the proof-of-principle coupling of two accelerating structures powered by different laser pulses, has increased interest in ANAs as a viable technology to be considered for a compact, TeV-class, lepton linear collider. However, intermediate facilities are required to test the technology and demonstrate key subsystems. A 20-100 GeV center-of-mass energy ANA-based lepton collider can be a possible candidate for an intermediate facility. Apart from being a test beam facility for accelerator and detector studies, this collider will provide opportunities to study muon and proton beam acceleration, investigate charged particle interactions with extreme electromagnetic fields (relevant for beam delivery system designs and to study the physics at the interaction point), as well as precision Quantum Chromodynamics and Beyond the Standard Model physics measurements. Possible applications of this collider include the studies of $γがんまγがんま$ and $e$-ion collider designs. △ Less

Submitted 15 April, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: contribution to Snowmass 2021

Journal ref: INST 19 T01010 (2024)

Snowmass Whitepaper AF6: Plasma-Based Particle Sources

Authors: M. Fuchs, B. A. Shadwick, N. Vafaei-Najafabadi, A. G. R. Thomas, G. Andonian, M. Büscher, A. Lehrach, O. Apsimon, G. Xia, D. Filippetto, C. B. Schroeder, M. C. Downer

Abstract: High-brightness beams generated by particle sources based on advanced accelerator concepts have the potential to become an essential part of future accelerator technology. High-gradient accelerators can generate and rapidly accelerate particle beams to relativistic energies while minimizing irreversible detrimental effects to the beam brightness that occur at low beam energies. Due to the high acc… ▽ More High-brightness beams generated by particle sources based on advanced accelerator concepts have the potential to become an essential part of future accelerator technology. High-gradient accelerators can generate and rapidly accelerate particle beams to relativistic energies while minimizing irreversible detrimental effects to the beam brightness that occur at low beam energies. Due to the high accelerating gradients, these novel accelerators are also significantly more compact than conventional technology. The beam parameters of these particle sources are largely determined by the injection and subsequent acceleration processes. While there has been significant progress crucial parameters that are required for a future collider or more near-term applications, including X-ray free-electron lasers (XFELs), such as a sufficiently small energy spread and small emittance for bunches with a high charge and at high pulse repetition rate. Major research and development efforts are required to realize these approaches for a front-end injector for a future collider in order to address these limitations. In particular, this includes methods to control and manipulate the phase-space and spin degrees-of-freedom of ultrashort LWFA electron bunches with high accuracy, methods that increase the laser-to-electron beam efficiency and increased repetition rate. This also includes the development of high-resolution diagnostics, such as full 6D phase-space measurements, beam polarimetry and high-fidelity simulation tools. A further increase in beam luminosity can be achieve through emittance damping. For future colliders, the damping rings might be replaced by a substantially more compact plasma-based approach. Here, plasma wigglers are used to achieve similar damping performance but over a two orders of magnitude reduced length. △ Less

Submitted 17 March, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: contribution to Snowmass 2021

Linear colliders based on laser-plasma accelerators

Authors: C. Benedetti, S. S. Bulanov, E. Esarey, C. G. R. Geddes, A. J. Gonsalves, A. Huebl, R. Lehe, K. Nakamura, C. B. Schroeder, D. Terzani, J. van Tilborg, M. Turner, J. -L. Vay, T. Zhou, F. Albert, J. Bromage, E. M. Campbell, D. H. Froula, J. P. Palastro, J. Zuegel, D. Bruhwiler, N. M. Cook, B. Cros, M. C. Downer, M. Fuchs , et al. (18 additional authors not shown)

Abstract: White paper to the Proceedings of the U.S. Particle Physics Community Planning Exercise (Snowmass 2021): Linear colliders based on laser-plasma accelerators White paper to the Proceedings of the U.S. Particle Physics Community Planning Exercise (Snowmass 2021): Linear colliders based on laser-plasma accelerators △ Less

Submitted 4 July, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021, Accelerator Frontier

Beam dynamics challenges in linear colliders based on laser-plasma accelerators

Authors: C. B. Schroeder, C. Benedetti, S. S. Bulanov, D. Terzani, E. Esarey, C. G. R. Geddes

Abstract: In this paper we discuss design considerations and beam dynamics challenges associated with laser-driven plasma-based accelerators as applied to multi-TeV-scale linear colliders. Plasma accelerators provide ultra-high gradients and ultra-short bunches, offering the potential for compact linacs and reduced power requirements. We show that stable, efficient acceleration with beam quality preservatio… ▽ More In this paper we discuss design considerations and beam dynamics challenges associated with laser-driven plasma-based accelerators as applied to multi-TeV-scale linear colliders. Plasma accelerators provide ultra-high gradients and ultra-short bunches, offering the potential for compact linacs and reduced power requirements. We show that stable, efficient acceleration with beam quality preservation is possible in the nonlinear bubble regime of laser-plasma accelerators using beam shaping. Ion motion, naturally occuring for dense beams (i.e., low emittance and high energy) severely damps transverse beam instabilities. Coulomb scattering by the background ions is considered and it is shown that the strong focusing in the plasma strongly suppresses scattering-induced emittance growth. Betatron radiation emission from the transverse motion of the beam in the plasma will result in beam power loss and energy spread growth; however for sub-100 nm emittances, the beam power loss and energy spread growth will be sub-percent for multi-TeV-class plasma linacs. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Toward Parallel in Time for Chaotic Dynamical Systems

Authors: David A. Vargas, Robert D. Falgout, Stefanie Günther, Jacob B. Schroder

Abstract: As CPU clock speeds have stagnated, and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes are a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel in Time schemes, such as the Mu… ▽ More As CPU clock speeds have stagnated, and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes are a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel in Time schemes, such as the Multigrid Reduction in Time (MGRIT) method, remedy this by parallelizing across time-steps, and have shown promising results for parabolic problems. However, chaotic problems have proved more difficult, since chaotic initial value problems are inherently ill-conditioned. MGRIT relies on a hierarchy of successively coarser time-grids to iteratively correct the solution on the finest time-grid, but due to the nature of chaotic systems, subtle inaccuracies on the coarser levels can lead to poor coarse-grid corrections. Here we propose a modification to nonlinear FAS multigrid, as well as a novel time-coarsening scheme, which together better capture long term behavior on coarse grids and greatly improve convergence of MGRIT for chaotic initial value problems. We provide supporting numerical results for the Lorenz system model problem. △ Less

Submitted 25 January, 2022; originally announced January 2022.

Report number: LLNL-CONF-830992 MSC Class: 65M22; 65M55

Weighted Relaxation for Multigrid Reduction in Time

Authors: Masumi Sugiyama, Jacob B. Schroder, Ben S. Southworth, Stephanie Friedhoff

Abstract: Based on current trends in computer architectures, faster compute speeds must come from increased parallelism rather than increased clock speeds, which are currently stagnate. This situation has created the well-known bottleneck for sequential time-integration, where each individual time-value (i.e., time-step) is computed sequentially. One approach to alleviate this and achieve parallelism in tim… ▽ More Based on current trends in computer architectures, faster compute speeds must come from increased parallelism rather than increased clock speeds, which are currently stagnate. This situation has created the well-known bottleneck for sequential time-integration, where each individual time-value (i.e., time-step) is computed sequentially. One approach to alleviate this and achieve parallelism in time is with multigrid. In this work, we consider multigrid-reduction-in-time (MGRIT), a multilevel method applied to the time dimension that computes multiple time-steps in parallel. Like all multigrid methods, MGRIT relies on the complementary relationship between relaxation on a fine-grid and a correction from the coarse grid to solve the problem. All current MGRIT implementations are based on unweighted-Jacobi relaxation; here we introduce the concept of weighted relaxation to MGRIT. We derive new convergence bounds for weighted relaxation, and use this analysis to guide the selection of relaxation weights. Numerical results then demonstrate that non-unitary relaxation weights consistently yield faster convergence rates and lower iteration counts for MGRIT when compared with unweighted relaxation. In most cases, weighted relaxation yields a 10%-20% saving in iterations. For A-stable integration schemes, results also illustrate that under-relaxation can restore convergence in some cases where unweighted relaxation is not convergent. △ Less

Submitted 5 July, 2021; originally announced July 2021.

Comments: 32 pages, 38 figures

Report number: LA-UR-21-26114 MSC Class: 65F10; 65M55

Scaling properties of liquid dynamics predicted from a single configuration: Small rigid molecules

Authors: Zahraa Sheydaafar, Jeppe C. Dyre, Thomas B. Schrøder

Abstract: Isomorphs are curves in the thermodynamic phase diagram along which structure and dynamics are invariant to a good approximation. There are two main ways to trace out isomorphs, the configurational-adiabat method and the direct-isomorph-check method. Recently a new method based on the scaling properties of forces was introduced and shown to work very well for atomic systems [T. B. Schroder, Phys.… ▽ More Isomorphs are curves in the thermodynamic phase diagram along which structure and dynamics are invariant to a good approximation. There are two main ways to trace out isomorphs, the configurational-adiabat method and the direct-isomorph-check method. Recently a new method based on the scaling properties of forces was introduced and shown to work very well for atomic systems [T. B. Schroder, Phys. Rev. Lett. 129, 245501 (2022)]. A unique feature of this method is that it only requires a single equilibrium configuration for tracing out an isomorph. We here test generalizations of this method to molecular systems and compare to simulations of three simple molecular models: the asymmetric dumbbell model of two Lennard-Jones spheres, the symmetric inverse-power-law dumbbell model, and the Lewis-Wahnström o-terphenyl model. We introduce and test two force-based and one torque-based methods, all of which require just a single configuration for tracing out an isomorph. Overall, the method based on requiring invariant center-of-mass reduced forces works best. △ Less

Submitted 8 April, 2023; v1 submitted 28 May, 2021; originally announced May 2021.

Comments: 20 pages, 12 figures

Journal ref: J. Phys. Chem. B 127, 2837 (2023)

Predicting scaling properties from a single fluid configuration

Authors: Thomas B. Schrøder

Abstract: Time-dependent dynamical properties of a fluid can not be estimated from a single configuration without performing a simulation. Here we show, however, that the scaling properties of both structure and dynamics can be predicted from a single configuration. The new method is demonstrated to work very well for equilibrium dynamics of the Kob-Andersen Binary Lennard-Jones mixture. Furthermore, the me… ▽ More Time-dependent dynamical properties of a fluid can not be estimated from a single configuration without performing a simulation. Here we show, however, that the scaling properties of both structure and dynamics can be predicted from a single configuration. The new method is demonstrated to work very well for equilibrium dynamics of the Kob-Andersen Binary Lennard-Jones mixture. Furthermore, the method is applied to isobaric cooling where the liquid falls out of equilibrium and forms a glass, demonstrating that the method requires neither equilibrium nor constant volume conditions to work, in contrast to existing methods. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Comments: 11 pages, 5 figures

Time-periodic steady-state solution of fluid-structure interaction and cardiac flow problems through multigrid-reduction-in-time

Authors: Andreas Hessenthaler, Robert D. Falgout, Jacob B. Schroder, Adelaide de Vecchi, David Nordsletten, Oliver Röhrle

Abstract: In this paper, a time-periodic MGRIT algorithm is proposed as a means to reduce the time-to-solution of numerical algorithms by exploiting the time periodicity inherent to many applications in science and engineering. The time-periodic MGRIT algorithm is applied to a variety of linear and nonlinear single- and multiphysics problems that are periodic-in-time. It is demonstrated that the proposed pa… ▽ More In this paper, a time-periodic MGRIT algorithm is proposed as a means to reduce the time-to-solution of numerical algorithms by exploiting the time periodicity inherent to many applications in science and engineering. The time-periodic MGRIT algorithm is applied to a variety of linear and nonlinear single- and multiphysics problems that are periodic-in-time. It is demonstrated that the proposed parallel-in-time algorithm can obtain the same time-periodic steady-state solution as sequential time-stepping. It is shown that the required number of MGRIT iterations can be estimated a priori and that the new MGRIT variant can significantly and consistently reduce the time-to-solution compared to sequential time-stepping, irrespective of the number of dimensions, linear or nonlinear PDE models, single-physics or coupled problems and the employed computing resources. The numerical experiments demonstrate that the time-periodic MGRIT algorithm enables a greater level of parallelism yielding faster turnaround, and thus, facilitating more complex and more realistic problems to be solved. △ Less

Submitted 14 November, 2021; v1 submitted 1 May, 2021; originally announced May 2021.

Comments: 29 pages; 7 pages Supplementary Materials

Report number: LLNL-JRNL-820515

Effect of Nozzle Curvature on Supersonic Gas Jets Used in Laser-Plasma Acceleration

Authors: Ocean Zhou, Hai-En Tsai, Tobias M. Ostermayr, Liona Fan-Chiang, Jeroen van Tilborg, Carl B. Schroeder, Eric Esarey, Cameron G. R. Geddes

Abstract: Supersonic gas jets produced by converging-diverging (C-D) nozzles are commonly used as targets for laser-plasma acceleration (LPA) experiments. A major point of interest for these targets is the gas density at the region of interaction where the laser ionizes the gas plume to create a plasma, providing the acceleration structure. Tuning the density profiles at this interaction region is crucial t… ▽ More Supersonic gas jets produced by converging-diverging (C-D) nozzles are commonly used as targets for laser-plasma acceleration (LPA) experiments. A major point of interest for these targets is the gas density at the region of interaction where the laser ionizes the gas plume to create a plasma, providing the acceleration structure. Tuning the density profiles at this interaction region is crucial to LPA optimization. A "flat-top" density profile is desired at this line of interaction to control laser propagation and high energy electron acceleration, while a short high-density profile is often preferred for acceleration of lower-energy tightly-focused laser-plasma interactions. A particular design parameter of interest is the curvature of the nozzle's diverging section. We examine three nozzle designs with different curvatures: the concave "bell", straight conical and convex "trumpet" nozzles. We demonstrate that, at mm-scale distances from the nozzle exit, the trumpet and straight nozzles, if optimized, produce "flat-top" density profiles whereas the bell nozzle creates focused regions of gas with higher densities. An optimization procedure for the trumpet nozzle is derived and compared to the straight nozzle optimization process. We find that the trumpet nozzle, by providing an extra parameter of control through its curvature, is more versatile for creating flat-top profiles and its optimization procedure is more refined compared to the straight nozzle and the straight nozzle optimization process. We present results for different nozzle designs from computational fluid dynamics (CFD) simulations performed with the program ANSYS Fluent and verify them experimentally using neutral density interferometry. △ Less

Submitted 15 September, 2021; v1 submitted 27 April, 2021; originally announced April 2021.

Comments: 14 pages, 16 figures, 5 tables

Journal ref: Physics of Plasmas Vol 28, Issue 9, Pg. 093107 (2021)

Does mesoscopic elasticity control viscous slowing down in glassforming liquids?

Authors: Geert Kapteijns, David Richard, Eran Bouchbinder, Thomas B. Schrøder, Jeppe C. Dyre, Edan Lerner

Abstract: The dramatic slowing down of relaxation dynamics of liquids approaching the glass transition remains a highly debated problem, where the crux of the puzzle resides in the elusive increase of the activation barrier $ΔでるたE(T)$ with decreasing temperature $T$. A class of theoretical frameworks -- known as elastic models -- attribute this temperature dependence to the variations of the liquid's macroscop… ▽ More The dramatic slowing down of relaxation dynamics of liquids approaching the glass transition remains a highly debated problem, where the crux of the puzzle resides in the elusive increase of the activation barrier $ΔでるたE(T)$ with decreasing temperature $T$. A class of theoretical frameworks -- known as elastic models -- attribute this temperature dependence to the variations of the liquid's macroscopic elasticity, quantified by the high-frequency shear modulus $G_\infty(T)$. While elastic models find some support in a number of experimental studies, these models do not take into account the spatial structures, length scales, and heterogeneity associated with structural relaxation in supercooled liquids. Here, we propose that viscous slowing down is controlled by a mesoscopic elastic stiffness $κかっぱ(T)$, defined as the characteristic stiffness of response fields to local dipole forces in the liquid's underlying inherent states. First, we show that $κかっぱ(T)$ -- which is intimately related to the energy and length scales characterizing quasilocalized, nonphononic excitations in glasses -- increases more strongly with decreasing $T$ than the macroscopic inherent state shear modulus $G(T)$ in several computer liquids. Second, we show that the simple relation $ΔでるたE(T)\proptoκかっぱ(T)$ holds remarkably well for some computer liquids, implying a direct connection between the liquid's underlying mesoscopic elasticity and enthalpic energy barriers. On the other hand, we show that for other computer liquids, the above relation fails. Finally, we provide strong evidence that what distinguishes computer liquids in which the $ΔでるたE(T) \propto κかっぱ(T)$ relation holds, from those in which it does not, is that the latter feature highly granular potential energy landscapes, where many sub-basins separated by low activation barriers exist. [Rest of abstract abridged] △ Less

Submitted 17 August, 2021; v1 submitted 21 March, 2021; originally announced March 2021.

Comments: 15 pages, 12 figures. Accepted manuscript

Journal ref: J. Chem. Phys. 155, 074502 (2021)

Radial Density Profile and Stability of Capillary Discharge Plasma Waveguides of Lengths up to 40 Centimeters

Authors: M. Turner, A. J. Gonsalves, S. S. Bulvanov, C. Benedetti, N. A. Bobrova, V. A. Gasilov, P. V. Sasorov, G. Korn, K. Nakamura, J. van Tilborg, C. G. Geddes, C. B. Schroeder, E. Esarey

Abstract: We measured the parameter reproducibility and radial electron density profile of capillary discharge waveguides with diameters of 650 um to 2 mm and lengths of 9 to 40 cm. To our knowledge, 40 cm is the longest discharge capillary plasma waveguide to date. This length is important for >= 10 GeV electron energy gain in a single laser driven plasma wakefield acceleration (LPA) stage. Evaluation of w… ▽ More We measured the parameter reproducibility and radial electron density profile of capillary discharge waveguides with diameters of 650 um to 2 mm and lengths of 9 to 40 cm. To our knowledge, 40 cm is the longest discharge capillary plasma waveguide to date. This length is important for >= 10 GeV electron energy gain in a single laser driven plasma wakefield acceleration (LPA) stage. Evaluation of waveguide parameter variations showed that their focusing strength was stable and reproducible to <0.2% and their average on-axis plasma electron density to <1%. These variations explain only a small fraction of LPA electron bunch variations observed in experiments to date. Measurements of laser pulse centroid oscillations revealed that the radial channel profile rises faster than parabolic and are in excellent agreement with magneto-hydro-dynamic simulation results. We show that the effects of non-parabolic contributions on Gaussian pulse propagation were negligible when the pulse was approximately matched to the channel. However, they affected pulse propagation for a non-matched configuration in which the waveguide was used as a plasma telescope to change the focused laser pulse spot size. △ Less

Submitted 14 January, 2021; v1 submitted 10 December, 2020; originally announced December 2020.

Prediction-Based Power Oversubscription in Cloud Platforms

Authors: Alok Kumbhare, Reza Azimi, Ioannis Manousakis, Anand Bonde, Felipe Frujeri, Nithish Mahalingam, Pulkit Misra, Seyyed Ahmad Javadi, Bianca Schroeder, Marcus Fontoura, Ricardo Bianchini

Abstract: Datacenter designers rely on conservative estimates of IT equipment power draw to provision resources. This leaves resources underutilized and requires more datacenters to be built. Prior work has used power capping to shave the rare power peaks and add more servers to the datacenter, thereby oversubscribing its resources and lowering capital costs. This works well when the workloads and their ser… ▽ More Datacenter designers rely on conservative estimates of IT equipment power draw to provision resources. This leaves resources underutilized and requires more datacenters to be built. Prior work has used power capping to shave the rare power peaks and add more servers to the datacenter, thereby oversubscribing its resources and lowering capital costs. This works well when the workloads and their server placements are known. Unfortunately, these factors are unknown in public clouds, forcing providers to limit the oversubscription so that performance is never impacted. In this paper, we argue that providers can use predictions of workload performance criticality and virtual machine (VM) resource utilization to increase oversubscription. This poses many challenges, such as identifying the performance-critical workloads from black-box VMs, creating support for criticality-aware power management, and increasing oversubscription while limiting the impact of capping. We address these challenges for the hardware and software infrastructures of Microsoft Azure. The results show that we enable a 2x increase in oversubscription with minimum impact to critical workloads. △ Less

Submitted 29 October, 2020; originally announced October 2020.

Clocking plasmon nanofocusing by THz near-field streaking

Authors: Lara Wimmer, Benjamin Schröder, Murat Sivis, Georg Herink, Claus Ropers

Abstract: We apply terahertz (THz) near-field streaking in a nanofocusing geometry to investigate plasmon polariton propagation on the shaft of a conical nanotip. By evaluating the delay between a streaking spectrogram for plasmon-induced photoemission with a measurement for direct apex excitation, we obtain an average plasmon group velocity, which is in agreement with numerical simulations. Combining plasm… ▽ More We apply terahertz (THz) near-field streaking in a nanofocusing geometry to investigate plasmon polariton propagation on the shaft of a conical nanotip. By evaluating the delay between a streaking spectrogram for plasmon-induced photoemission with a measurement for direct apex excitation, we obtain an average plasmon group velocity, which is in agreement with numerical simulations. Combining plasmon-induced photoemission with THz near-field streaking facilitates extensive control over localized photoelectron sources for time-resolved imaging and diffraction. △ Less

Submitted 25 August, 2020; originally announced August 2020.

Comments: 4 pages, 3 figures

Journal ref: Appl. Phys. Lett. 111, 131102 (2017)

Powering Electronic Devices from Salt Gradients in AA Battery-Sized Stacks of Hydrogel-Infused Paper

Authors: Anirvan Guha, Trevor J. Kalkus, Thomas B. H. Schroeder, Oliver G. Willis, Michael Mayer

Abstract: Strongly electric fish use gradients of ions within their bodies to generate stunning external electrical discharges; the most powerful of these organisms, the Atlantic torpedo ray, can produce pulses of over 1 kW from its electric organs. Despite extensive study of this phenomenon in nature, the development of artificial power generation schemes based on ion gradients for portable, wearable, or i… ▽ More Strongly electric fish use gradients of ions within their bodies to generate stunning external electrical discharges; the most powerful of these organisms, the Atlantic torpedo ray, can produce pulses of over 1 kW from its electric organs. Despite extensive study of this phenomenon in nature, the development of artificial power generation schemes based on ion gradients for portable, wearable, or implantable human use has remained out of reach. Previously, inspired by the electric eel, we developed an artificial electric organ that generated electricity from ion gradients within stacked hydrogels and, like the eel, was optimized to deliver large voltages that exceeded 100 V. Due to its high internal resistance, the current of this power source was, however, too low to power standard electronics. Here we introduce an artificial electric organ that takes inspiration from the unique morphologies of torpedo rays for maximal current output. This power source uses a hybrid material of hydrogel-infused paper to create, organize, and reconfigure stacks of thin, arbitrarily large gel films both in series and in parallel. The resulting increase in electrical power by almost two orders of magnitude compared to the original eel-inspired design makes it possible to power electronic devices and establishes that biology's mechanism of generating significant electrical power can now be realized from benign and soft materials in a portable size. △ Less

Submitted 28 May, 2020; originally announced May 2020.

Comments: 34 pages, 10 figures, 1 table

Structured Query-Based Image Retrieval Using Scene Graphs

Authors: Brigit Schroeder, Subarna Tripathi

Abstract: A structured query can capture the complexity of object interactions (e.g. 'woman rides motorcycle') unlike single objects (e.g. 'woman' or 'motorcycle'). Retrieval using structured queries therefore is much more useful than single object retrieval, but a much more challenging problem. In this paper we present a method which uses scene graph embeddings as the basis for an approach to image retriev… ▽ More A structured query can capture the complexity of object interactions (e.g. 'woman rides motorcycle') unlike single objects (e.g. 'woman' or 'motorcycle'). Retrieval using structured queries therefore is much more useful than single object retrieval, but a much more challenging problem. In this paper we present a method which uses scene graph embeddings as the basis for an approach to image retrieval. We examine how visual relationships, derived from scene graphs, can be used as structured queries. The visual relationships are directed subgraphs of the scene graph with a subject and object as nodes connected by a predicate relationship. Notably, we are able to achieve high recall even on low to medium frequency objects found in the long-tailed COCO-Stuff dataset, and find that adding a visual relationship-inspired loss boosts our recall by 10% in the best case. △ Less

Submitted 13 May, 2020; originally announced May 2020.

Comments: Accepted to Diagram Image Retrieval and Analysis (DIRA) Workshop at CVPR 2020

Multilevel Initialization for Layer-Parallel Deep Neural Network Training

Authors: Eric C. Cyr, Stefanie Günther, Jacob B. Schroder

Abstract: This paper investigates multilevel initialization strategies for training very deep neural networks with a layer-parallel multigrid solver. The scheme is based on the continuous interpretation of the training problem as a problem of optimal control, in which neural networks are represented as discretizations of time-dependent ordinary differential equations. A key goal is to develop a method able… ▽ More This paper investigates multilevel initialization strategies for training very deep neural networks with a layer-parallel multigrid solver. The scheme is based on the continuous interpretation of the training problem as a problem of optimal control, in which neural networks are represented as discretizations of time-dependent ordinary differential equations. A key goal is to develop a method able to intelligently initialize the network parameters for the very deep networks enabled by scalable layer-parallel training. To do this, we apply a refinement strategy across the time domain, that is equivalent to refining in the layer dimension. The resulting refinements create deep networks, with good initializations for the network parameters coming from the coarser trained networks. We investigate the effectiveness of such multilevel "nested iteration" strategies for network training, showing supporting numerical evidence of reduced run time for equivalent accuracy. In addition, we study whether the initialization strategies provide a regularizing effect on the overall training process and reduce sensitivity to hyperparameters and randomness in initial network parameters. △ Less

Submitted 18 December, 2019; originally announced December 2019.

Using Image Priors to Improve Scene Understanding

Authors: Brigit Schroeder, Hanlin Tang, Alexandre Alahi

Abstract: Semantic segmentation algorithms that can robustly segment objects across multiple camera viewpoints are crucial for assuring navigation and safety in emerging applications such as autonomous driving. Existing algorithms treat each image in isolation, but autonomous vehicles often revisit the same locations or maintain information from the immediate past. We propose a simple yet effective method f… ▽ More Semantic segmentation algorithms that can robustly segment objects across multiple camera viewpoints are crucial for assuring navigation and safety in emerging applications such as autonomous driving. Existing algorithms treat each image in isolation, but autonomous vehicles often revisit the same locations or maintain information from the immediate past. We propose a simple yet effective method for leveraging these image priors to improve semantic segmentation of images from sequential driving datasets. We examine several methods to fuse these temporal scene priors, and introduce a prior fusion network that is able to learn how to transfer this information. The prior fusion model improves the accuracy over the non-prior baseline from 69.1% to 73.3% for dynamic classes, and from 88.2% to 89.1% for static classes. Compared to models such as FCN-8, our prior method achieves the same accuracy with 5 times fewer parameters. We used a simple encoder decoder backbone, but this general prior fusion method could be applied to more complex semantic segmentation backbones. We also discuss how structured representations of scenes in the form of a scene graph could be leveraged as priors to further improve scene understanding. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: Accepted to Women in Computer Vision (WiCV) Workshop at CVPR 2019

Triplet-Aware Scene Graph Embeddings

Authors: Brigit Schroeder, Subarna Tripathi, Hanlin Tang

Abstract: Scene graphs have become an important form of structured knowledge for tasks such as for image generation, visual relation detection, visual question answering, and image retrieval. While visualizing and interpreting word embeddings is well understood, scene graph embeddings have not been fully explored. In this work, we train scene graph embeddings in a layout generation task with different forms… ▽ More Scene graphs have become an important form of structured knowledge for tasks such as for image generation, visual relation detection, visual question answering, and image retrieval. While visualizing and interpreting word embeddings is well understood, scene graph embeddings have not been fully explored. In this work, we train scene graph embeddings in a layout generation task with different forms of supervision, specifically introducing triplet super-vision and data augmentation. We see a significant performance increase in both metrics that measure the goodness of layout prediction, mean intersection-over-union (mIoU)(52.3% vs. 49.2%) and relation score (61.7% vs. 54.1%),after the addition of triplet supervision and data augmentation. To understand how these different methods affect the scene graph representation, we apply several new visualization and evaluation methods to explore the evolution of the scene graph embedding. We find that triplet supervision significantly improves the embedding separability, which is highly correlated with the performance of the layout prediction model. △ Less

Submitted 19 September, 2019; originally announced September 2019.

Comments: Accepted to Scene Graph Representation Learning workshop at ICCV 2019

Multiple-colliding laser pulses as a basis for studying high-field high-energy physics

Authors: J. Magnusson, A. Gonoskov, M. Marklund, T. Zh. Esirkepov, J. K. Koga, K. Kondo, M. Kando, S. V. Bulanov, G. Korn, C. G. R. Geddes, C. B. Schroeder, E. Esarey, S. S. Bulanov

Abstract: Apart from maximizing the strength of optical electromagnetic fields achievable at high-intensity laser facilities, the collision of several phase-matched laser pulses has been theoretically identified as a trigger of and way to study various phenomena. These range from the basic processes of strong-field quantum electrodynamics to the extraordinary dynamics of the generated electron-positron plas… ▽ More Apart from maximizing the strength of optical electromagnetic fields achievable at high-intensity laser facilities, the collision of several phase-matched laser pulses has been theoretically identified as a trigger of and way to study various phenomena. These range from the basic processes of strong-field quantum electrodynamics to the extraordinary dynamics of the generated electron-positron plasmas. This has paved the way for several experimental proposals aimed at both fundamental studies of matter at extreme conditions and the creation of particle and radiation sources. Because of the unprecedented capabilities of such sources they have the potential to open up new opportunities for experimental studies in nuclear and quark-gluon physics. We here perform a systematic analysis of different regimes and opportunities achievable with the concept of multiple-colliding laser pulses (MCLP), for both current and upcoming laser facilities. We reveal that several distinct regimes could be within reach of multi-PW laser facilities. △ Less

Submitted 12 June, 2019; originally announced June 2019.

Comments: 15 pages, 12 figures

Journal ref: Phys. Rev. A 100, 063404 (2019)

Solid-like mean-square displacement in glass-forming liquids

Authors: Thomas B. Schrøder, Jeppe C. Dyre

Abstract: It was recently shown that the real part of the frequency-dependent fluidity for several glass-forming liquids of different chemistry conforms to the prediction of the random barrier model (RBM) devised for ac electrical conduction in disordered solids [S. P. Bierwirth \textit{et al.}, Phys. Rev. Lett. {\bf 119}, 248001 (2017)]. Inspired by these results we introduce a crystallization-resistant mo… ▽ More It was recently shown that the real part of the frequency-dependent fluidity for several glass-forming liquids of different chemistry conforms to the prediction of the random barrier model (RBM) devised for ac electrical conduction in disordered solids [S. P. Bierwirth \textit{et al.}, Phys. Rev. Lett. {\bf 119}, 248001 (2017)]. Inspired by these results we introduce a crystallization-resistant modification of the Kob-Andersen binary Lennard-Jones mixture for which the results of extensive graphics-processing unit (GPU)-based molecular-dynamics simulations are presented. We find that the low-temperature mean-square displacement is fitted well by the RBM prediction, which involves no shape parameters. This finding highlights the challenge of explaining why a simple model based on hopping of non-interacting particles in a fixed random energy landscape can reproduce the complex and highly cooperative dynamics of glass-forming liquids. △ Less

Submitted 7 March, 2020; v1 submitted 27 May, 2019; originally announced May 2019.

Journal ref: J. Chem. Phys. 152, 141101 (2020)

Ion Acceleration in Laser Generated Mega Tesla Magnetic Vortex

Authors: Jaehong Park, Stepan S. Bulanov, Jianhui Bin, Qing Ji, Sven Steinke, Jean-Luc Vay, Cameron G. R. Geddes, Carl B. Schroeder, Wim P. Leemans, Thomas Schenkel, Eric Esarey

Abstract: Magnetic Vortex Acceleration (MVA) from near critical density targets is one of the promising schemes of laser-driven ion acceleration. 3D particle-in-cell simulations are used to explore a more extensive laser-target parameter space than previously reported on in the literature as well as to study the laser pulse coupling to the target, the structure of the fields, and the properties of the accel… ▽ More Magnetic Vortex Acceleration (MVA) from near critical density targets is one of the promising schemes of laser-driven ion acceleration. 3D particle-in-cell simulations are used to explore a more extensive laser-target parameter space than previously reported on in the literature as well as to study the laser pulse coupling to the target, the structure of the fields, and the properties of the accelerated ion beam in the MVA scheme. The efficiency of acceleration depends on the coupling of the laser energy to the self-generated channel in the target. The accelerated proton beams demonstrate high level of collimation with achromatic angular divergence, and carry a significant amount of charge. For PW-class lasers, this acceleration regime provides favorable scaling of maximum ion energy with laser power for optimized interaction parameters. The mega Tesla-level magnetic fields generated by the laser-driven co-axial plasma structure in the target are prerequisite for accelerating protons to the energy of several hundred MeV. △ Less

Submitted 5 April, 2019; originally announced April 2019.

Comments: submitted to Physics of Plasmas

The Role of Energy Minimization in Algebraic Multigrid Interpolation

Authors: James Brannick, Scott P. MacLachlan, Jacob B. Schroder, Ben S. Southworth

Abstract: Algebraic multigrid (AMG) methods are powerful solvers with linear or near-linear computational complexity for certain classes of linear systems, Ax=b. Broadening the scope of problems that AMG can effectively solve requires the development of improved interpolation operators. Such development is often based on AMG convergence theory. However, convergence theory in AMG tends to have a disconnect w… ▽ More Algebraic multigrid (AMG) methods are powerful solvers with linear or near-linear computational complexity for certain classes of linear systems, Ax=b. Broadening the scope of problems that AMG can effectively solve requires the development of improved interpolation operators. Such development is often based on AMG convergence theory. However, convergence theory in AMG tends to have a disconnect with AMG in practice due to the practical constraints of (i) maintaining matrix sparsity in transfer and coarse-grid operators, and (ii) retaining linear complexity in the setup and solve phase. This paper presents a review of fundamental results in AMG convergence theory, followed by a discussion on how these results can be used to motivate interpolation operators in practice. A general weighted energy minimization functional is then proposed to form interpolation operators, and a novel `diagonal' preconditioner for Sylvester- or Lyapunov-type equations developed simultaneously. Although results based on the weighted energy minimization typically underperform compared to a fully constrained energy minimization, numerical results provide new insight into the role of energy minimization and constraint vectors in AMG interpolation. △ Less

Submitted 13 February, 2019; originally announced February 2019.