Showing 1–50 of 70 results for author: Johnson, M A

Evaluation of Provenance Serialisations for Astronomical Provenance

Authors: Michael A. C. Johnson, Marcus Paradies, Hans-Rainer Klöckner, Albina Muzafarova, Kristen Lackeos, David J. Champion, Marta Dembska, Sirko Schindler

Abstract: Provenance data from astronomical pipelines are instrumental in establishing trust and reproducibility in the data processing and products. In addition, astronomers can query their provenance to answer questions routed in areas such as anomaly detection, recommendation, and prediction. The next generation of astronomical survey telescopes such as the Vera Rubin Observatory or Square Kilometre Arra… ▽ More Provenance data from astronomical pipelines are instrumental in establishing trust and reproducibility in the data processing and products. In addition, astronomers can query their provenance to answer questions routed in areas such as anomaly detection, recommendation, and prediction. The next generation of astronomical survey telescopes such as the Vera Rubin Observatory or Square Kilometre Array, are capable of producing peta to exabyte scale data, thereby amplifying the importance of even small improvements to the efficiency of provenance storage or querying. In order to determine how astronomers should store and query their provenance data, this paper reports on a comparison between the turtle and JSON provenance serialisations. The triple store Apache Jena Fuseki and the graph database system Neo4j were selected as representative database management systems (DBMS) for turtle and JSON, respectively. Simulated provenance data was uploaded to and queried over each DBMS and the metrics measured for comparison were the accuracy and timing of the queries as well as the data upload times. It was found that both serialisations are competent for this purpose, and both have similar query accuracy. The turtle provenance was found to be more efficient at storing and uploading the data. Regarding queries, for small datasets ($<$5MB) and simple information retrieval queries, the turtle serialisation was also found to be more efficient. However, queries for JSON serialised provenance were found to be more efficient for more complex queries which involved matching patterns across the DBMS, this effect scaled with the size of the queried provenance. △ Less

Submitted 19 July, 2024; originally announced July 2024.

Comments: 9 pages, 8 figures, to be published in the 16th International Workshop on Theory and Practice of Provenance

Pipeline Provenance for Analysis, Evaluation, Trust or Reproducibility

Authors: Michael A. C. Johnson, Hans-Rainer Klöckner, Albina Muzafarova, Kristen Lackeos, David J. Champion, Marta Dembska, Sirko Schindler, Marcus Paradies

Abstract: Data volumes and rates of research infrastructures will continue to increase in the upcoming years and impact how we interact with their final data products. Little of the processed data can be directly investigated and most of it will be automatically processed with as little user interaction as possible. Capturing all necessary information of such processing ensures reproducibility of the final… ▽ More Data volumes and rates of research infrastructures will continue to increase in the upcoming years and impact how we interact with their final data products. Little of the processed data can be directly investigated and most of it will be automatically processed with as little user interaction as possible. Capturing all necessary information of such processing ensures reproducibility of the final results and generates trust in the entire process. We present PRAETOR, a software suite that enables automated generation, modelling, and analysis of provenance information of Python pipelines. Furthermore, the evaluation of the pipeline performance, based upon a user defined quality matrix in the provenance, enables the first step of machine learning processes, where such information can be fed into dedicated optimisation procedures. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 4 pages, 3 figures

Journal ref: Research Notes of the AAS 8.4 (2024): 100

Orbital Stability of Smooth Solitary Waves for the Novikov Equation

Authors: Brett Ehrman, Mathew A. Johnson, Stéphane Lafortune

Abstract: We study the orbital stability of smooth solitary wave solutions of the Novikov equation, which is a Camassa-Holm type equation with cubic nonlinearities. These solitary waves are shown to exist as a one-parameter family (up to spatial translations) parameterized by their asymptotic endstate, and are encoded as critical points of a particular action functional. As an important step in our analysis… ▽ More We study the orbital stability of smooth solitary wave solutions of the Novikov equation, which is a Camassa-Holm type equation with cubic nonlinearities. These solitary waves are shown to exist as a one-parameter family (up to spatial translations) parameterized by their asymptotic endstate, and are encoded as critical points of a particular action functional. As an important step in our analysis we must study the spectrum the Hessian of this action functional, which turns out to be a nonlocal integro-differential operator acting on $L^2(\mathbb{R})$. We provide a combination of analytical and numerical evidence that the necessary spectral hypotheses always holds for the Novikov equation. Together with a detailed study of the associated Vakhitov-Kolokolov condition, our analysis indicates that all smooth solitary wave solutions of the Novikov equation are nonlinearly orbitally stable. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: 34 pages, 4 figures

Rapid cryogenic characterisation of 1024 integrated silicon quantum dots

Authors: Edward J. Thomas, Virginia N. Ciriano-Tejel, David F. Wise, Domenic Prete, Mathieu de Kruijf, David J. Ibberson, Grayson M. Noah, Alberto Gomez-Saiz, M. Fernando Gonzalez-Zalba, Mark A. I. Johnson, John J. L. Morton

Abstract: Quantum computers are nearing the thousand qubit mark, with the current focus on scaling to improve computational performance. As quantum processors grow in complexity, new challenges arise such as the management of device variability and the interface with supporting electronics. Spin qubits in silicon quantum dots are poised to address these challenges with their proven control fidelities and po… ▽ More Quantum computers are nearing the thousand qubit mark, with the current focus on scaling to improve computational performance. As quantum processors grow in complexity, new challenges arise such as the management of device variability and the interface with supporting electronics. Spin qubits in silicon quantum dots are poised to address these challenges with their proven control fidelities and potential for compatibility with large-scale integration. Here, we demonstrate the integration of 1024 silicon quantum dots with on-chip digital and analogue electronics, all operating below 1 K. A high-frequency analogue multiplexer provides fast access to all devices with minimal electrical connections, enabling characteristic data across the quantum dot array to be acquired in just 5 minutes. We achieve this by leveraging radio-frequency reflectometry with state-of-the-art signal integrity, reaching a minimum integration time of 160 ps. Key quantum dot parameters are extracted by fast automated machine learning routines to assess quantum dot yield and understand the impact of device design. We find correlations between quantum dot parameters and room temperature transistor behaviour that may be used as a proxy for in-line process monitoring. Our results show how rapid large-scale studies of silicon quantum devices can be performed at lower temperatures and measurement rates orders of magnitude faster than current probing techniques, and form a platform for the future on-chip addressing of large scale qubit arrays. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Comments: Main text: 14 pages, 8 figures, 1 table Supplementary: 8 pages, 6 figures

Orbital Stability of Periodic Traveling Waves in the $b$-Camassa-Holm Equation

Authors: Brett Ehrman, Mathew A. Johnson

Abstract: In this paper, we identify criteria that guarantees the nonlinear orbital stability of a given periodic traveling wave solution within the b-family Camassa-Holm equation. These periodic waves exist as 3-parameter families (up to spatial translations) of smooth traveling wave solutions, and their stability criteria are expressed in terms of Jacobians of the conserved quantities with respect to thes… ▽ More In this paper, we identify criteria that guarantees the nonlinear orbital stability of a given periodic traveling wave solution within the b-family Camassa-Holm equation. These periodic waves exist as 3-parameter families (up to spatial translations) of smooth traveling wave solutions, and their stability criteria are expressed in terms of Jacobians of the conserved quantities with respect to these parameters. The stability criteria utilizes a general Hamiltonian structure which exists for every $b>1$, and hence applies outside of the completely integrable cases ($b=2$ and $b=3$). △ Less

Submitted 19 February, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 26 pages, 1 figure. This update fixes some typos, and adds discussion regarding numerical calculations and verification of our stability conditions

Tomography of entangling two-qubit logic operations in exchange-coupled donor electron spin qubits

Authors: Holly G. Stemp, Serwan Asaad, Mark R. van Blankenstein, Arjen Vaartjes, Mark A. I. Johnson, Mateusz T. Mądzik, Amber J. A. Heskes, Hannes R. Firgau, Rocky Y. Su, Chih Hwan Yang, Arne Laucht, Corey I. Ostrove, Kenneth M. Rudinger, Kevin Young, Robin Blume-Kohout, Fay E. Hudson, Andrew S. Dzurak, Kohei M. Itoh, Alexander M. Jakob, Brett C. Johnson, David N. Jamieson, Andrea Morello

Abstract: Scalable quantum processors require high-fidelity universal quantum logic operations in a manufacturable physical platform. Donors in silicon provide atomic size, excellent quantum coherence and compatibility with standard semiconductor processing, but no entanglement between donor-bound electron spins has been demonstrated to date. Here we present the experimental demonstration and tomography of… ▽ More Scalable quantum processors require high-fidelity universal quantum logic operations in a manufacturable physical platform. Donors in silicon provide atomic size, excellent quantum coherence and compatibility with standard semiconductor processing, but no entanglement between donor-bound electron spins has been demonstrated to date. Here we present the experimental demonstration and tomography of universal 1- and 2-qubit gates in a system of two weakly exchange-coupled electrons, bound to single phosphorus donors introduced in silicon by ion implantation. We surprisingly observe that the exchange interaction has no effect on the qubit coherence. We quantify the fidelity of the quantum operations using gate set tomography (GST), and we use the universal gate set to create entangled Bell states of the electrons spins, with fidelity ~ 93%, and concurrence 0.91 +/- 0.08. These results form the necessary basis for scaling up donor-based quantum computers. △ Less

Submitted 2 March, 2024; v1 submitted 27 September, 2023; originally announced September 2023.

Error channels in quantum nondemolition measurements on spin systems

Authors: Benjamin Joecker, Holly G. Stemp, Irene Fernández de Fuentes, Mark A. I. Johnson, Andrea Morello

Abstract: Quantum nondemolition (QND) measurements are a precious resource for quantum information processing. Repetitive QND measurements can boost the fidelity of qubit preparation and measurement, even when the underlying single-shot measurements are of low fidelity. However, this fidelity boost is limited by the degree in which the physical system allows for a truly QND process -- slight deviations from… ▽ More Quantum nondemolition (QND) measurements are a precious resource for quantum information processing. Repetitive QND measurements can boost the fidelity of qubit preparation and measurement, even when the underlying single-shot measurements are of low fidelity. However, this fidelity boost is limited by the degree in which the physical system allows for a truly QND process -- slight deviations from ideal QND measurement result in bit flip errors (`quantum jumps') if the measurement is repeated too often. Here, we develop a theoretical framework to understand and quantify the resulting error arising from deviation from perfect QND measurement in model spin qubit systems. We first develop our model on the ubiquitous example of exchange-coupled electron spins qubits tunnel-coupled to a charge reservoir. We then extend it to electron-nuclear spin systems, to illustrate the crucial similarities and differences between the two limits. Applied to the well-understood platform of a donor nuclear spin in silicon, the model shows excellent agreement with experiments. For added generality, we conclude the work by considering the effect of anisotropic spin couplings. △ Less

Submitted 27 November, 2023; v1 submitted 26 July, 2023; originally announced July 2023.

Journal ref: Phys.Rev.B 109 (2024) 085302

Nonlinear Subharmonic Dynamics of Spectrally Stable Lugiato-Lefever Periodic Waves

Authors: Mariana Haragus, Mathew A. Johnson, Wesley R. Perkins, Björn de Rijk

Abstract: We study the nonlinear dynamics of perturbed, spectrally stable $T$-periodic stationary solutions of the Lugiato-Lefever equation (LLE), a damped nonlinear Schrödinger equation with forcing that arises in nonlinear optics. It is known that for each $N\in\mathbb{N}$, such a $T$-periodic wave train is (orbitally) asymptotically stable against $NT$-periodic, i.e. subharmonic, perturbations. Unfortuna… ▽ More We study the nonlinear dynamics of perturbed, spectrally stable $T$-periodic stationary solutions of the Lugiato-Lefever equation (LLE), a damped nonlinear Schrödinger equation with forcing that arises in nonlinear optics. It is known that for each $N\in\mathbb{N}$, such a $T$-periodic wave train is (orbitally) asymptotically stable against $NT$-periodic, i.e. subharmonic, perturbations. Unfortunately, in such results both the allowable size of initial perturbations as well as the exponential decay rates of perturbations depend on $N$ and, in fact, tend to zero as $N\to\infty$, leading to a lack of uniformity in the period of the perturbation. In recent work, the authors performed a delicate decomposition of the associated linearized solution operator and obtained linear estimates which are uniform in $N$. The dynamical description suggested by this uniform linear theory indicates that the corresponding nonlinear iteration can only be closed if one allows for a spatio-temporal phase modulation of the underlying wave. However, such a modulated perturbation is readily seen to satisfy a quasilinear equation, yielding an inherent loss of regularity. We regain regularity by transferring a nonlinear damping estimate, which has recently been obtained for the LLE in the case of localized perturbations to the case of subharmonic perturbations. Thus, we obtain a nonlinear, subharmonic stability result for periodic stationary solutions of the LLE that is uniform in $N$. This in turn yields an improved nonuniform subharmonic stability result providing an $N$-independent ball of initial perturbations which eventually exhibit exponential decay at an $N$-dependent rate. Finally, we argue that our results connect in the limit $N \to \infty$ to previously established stability results against localized perturbations, thereby unifying existing theories. △ Less

Submitted 23 July, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 31 pages. Minor update

Navigating the 16-dimensional Hilbert space of a high-spin donor qudit with electric and magnetic fields

Authors: Irene Fernández de Fuentes, Tim Botzem, Mark A. I. Johnson, Arjen Vaartjes, Serwan Asaad, Vincent Mourik, Fay E. Hudson, Kohei M. Itoh, Brett C. Johnson, Alexander M. Jakob, Jeffrey C. McCallum, David N. Jamieson, Andrew S. Dzurak, Andrea Morello

Abstract: Efficient scaling and flexible control are key aspects of useful quantum computing hardware. Spins in semiconductors combine quantum information processing with electrons, holes or nuclei, control with electric or magnetic fields, and scalable coupling via exchange or dipole interaction. However, accessing large Hilbert space dimensions has remained challenging, due to the short-distance nature of… ▽ More Efficient scaling and flexible control are key aspects of useful quantum computing hardware. Spins in semiconductors combine quantum information processing with electrons, holes or nuclei, control with electric or magnetic fields, and scalable coupling via exchange or dipole interaction. However, accessing large Hilbert space dimensions has remained challenging, due to the short-distance nature of the interactions. Here, we present an atom-based semiconductor platform where a 16-dimensional Hilbert space is built by the combined electron-nuclear states of a single antimony donor in silicon. We demonstrate the ability to navigate this large Hilbert space using both electric and magnetic fields, with gate fidelity exceeding 99.8% on the nuclear spin, and unveil fine details of the system Hamiltonian and its susceptibility to control and noise fields. These results establish high-spin donors as a rich platform for practical quantum information and to explore quantum foundations. △ Less

Submitted 14 June, 2023; v1 submitted 12 June, 2023; originally announced June 2023.

Comments: 31 pages and 19 figures including Supplementary Materials

Journal ref: Nat Commun 15, 1380 (2024)

Modulational Instability in the Ostrovsky Equation and Related Models

Authors: Bhavna, Mathew A. Johnson, Ashish Kumar Pandey

Abstract: We study the modulational instability of small-amplitude periodic traveling wave solutions in a generalized Ostrovsky equation. Specifically, we investigate the invertibility of the associated linearized operator in the vicinity of the origin and derive a modulational instability index that depends on the dispersion and nonlinearity. We then show that the small-amplitude periodic traveling waves o… ▽ More We study the modulational instability of small-amplitude periodic traveling wave solutions in a generalized Ostrovsky equation. Specifically, we investigate the invertibility of the associated linearized operator in the vicinity of the origin and derive a modulational instability index that depends on the dispersion and nonlinearity. We then show that the small-amplitude periodic traveling waves of the Ostrovsky equation exhibit modulational instability if the wavenumber is greater than a critical value which agrees with previously found numerical results in \cite{Whitfield2014Rotation-inducedWaves} both qualitatively and quantitatively. We also study the effects of surface tension on modulational instability using the index. △ Less

Submitted 14 May, 2023; originally announced May 2023.

Parma: Confidential Containers via Attested Execution Policies

Authors: Matthew A. Johnson, Stavros Volos, Ken Gordon, Sean T. Allen, Christoph M. Wintersteiger, Sylvan Clebsch, John Starks, Manuel Costa

Abstract: Container-based technologies empower cloud tenants to develop highly portable software and deploy services in the cloud at a rapid pace. Cloud privacy, meanwhile, is important as a large number of container deployments operate on privacy-sensitive data, but challenging due to the increasing frequency and sophistication of attacks. State-of-the-art confidential container-based designs leverage proc… ▽ More Container-based technologies empower cloud tenants to develop highly portable software and deploy services in the cloud at a rapid pace. Cloud privacy, meanwhile, is important as a large number of container deployments operate on privacy-sensitive data, but challenging due to the increasing frequency and sophistication of attacks. State-of-the-art confidential container-based designs leverage process-based trusted execution environments (TEEs), but face security and compatibility issues that limits their practical deployment. We propose Parma, an architecture that provides lift-and-shift deployment of unmodified containers while providing strong security protection against a powerful attacker who controls the untrusted host and hypervisor. Parma leverages VM-level isolation to execute a container group within a unique VM-based TEE. Besides container integrity and user data confidentiality and integrity, Parma also offers container attestation and execution integrity based on an attested execution policy. Parma execution policies provide an inductive proof over all future states of the container group. This proof, which is established during initialization, forms a root of trust that can be used for secure operations within the container group without requiring any modifications of the containerized workflow itself (aside from the inclusion of the execution policy.) We evaluate Parma on AMD SEV-SNP processors by running a diverse set of workloads demonstrating that workflows exhibit 0-26% additional overhead in performance over running outside the enclave, with a mean 13% overhead on SPEC2017, while requiring no modifications to their program code. Adding execution policies introduces less than 1% additional overhead. Furthermore, we have deployed Parma as the underlying technology driving Confidential Containers on Azure Container Instances. △ Less

Submitted 7 March, 2023; v1 submitted 8 February, 2023; originally announced February 2023.

Comments: 12 pages, 6 figures, 2 tables

Design and Performance of a Novel Low Energy Multi-Species Beamline for the ALPHA Antihydrogen Experiment

Authors: C. J. Baker, W. Bertsche, A. Capra, C. L. Cesar, M. Charlton, A. J. Christensen, R. Collister, A. Cridland Mathad, S. Eriksson, A. Evans, N. Evetts, S. Fabbri, J. Fajans, T. Friesen, M. C. Fujiwara, D. R. Gill, P. Grandemange, P. Granum, J. S. Hangst, M. E. Hayden, D. Hodgkinson, C. A. Isaac, M. A. Johnson, J. M. Jones, S. A. Jones , et al. (25 additional authors not shown)

Abstract: The ALPHA Collaboration, based at the CERN Antiproton Decelerator, has recently implemented a novel beamline for low-energy ($\lesssim$ 100 eV) positron and antiproton transport between cylindrical Penning traps that have strong axial magnetic fields. Here, we describe how a combination of semianalytical and numerical calculations were used to optimise the layout and design of this beamline. Using… ▽ More The ALPHA Collaboration, based at the CERN Antiproton Decelerator, has recently implemented a novel beamline for low-energy ($\lesssim$ 100 eV) positron and antiproton transport between cylindrical Penning traps that have strong axial magnetic fields. Here, we describe how a combination of semianalytical and numerical calculations were used to optimise the layout and design of this beamline. Using experimental measurements taken during the initial commissioning of the instrument, we evaluate its performance and validate the models used for its development. By combining data from a range of sources, we show that the beamline has a high transfer efficiency, and estimate that the percentage of particles captured in the experiments from each bunch is (78 $\pm$ 3)% for up to $10^{5}$ antiprotons, and (71 $\pm$ 5)% for bunches of up to $10^{7}$ positrons. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: 15 pages, 15 figures

Transverse phase space tomography in the CLARA accelerator test facility using image compression and machine learning

Authors: Andrzej Wolski, Mark A. Johnson, Matthew King, Boris L. Militsyn, Peter H. Williams

Abstract: We describe a novel technique, based on image compression and machine learning, for transverse phase space tomography in two degrees of freedom in an accelerator beamline. The technique has been used in the CLARA accelerator test facility at Daresbury Laboratory: results from the machine learning method are compared with those from a conventional tomography algorithm (algebraic reconstruction), ap… ▽ More We describe a novel technique, based on image compression and machine learning, for transverse phase space tomography in two degrees of freedom in an accelerator beamline. The technique has been used in the CLARA accelerator test facility at Daresbury Laboratory: results from the machine learning method are compared with those from a conventional tomography algorithm (algebraic reconstruction), applied to the same data. The use of machine learning allows reconstruction of the 4D phase space distribution of the beam to be carried out much more rapidly than using conventional tomography algorithms, and also enables the use of image compression to reduce significantly the size of the data sets involved in the analysis. Results from the machine learning technique are at least as good as those from the algebraic reconstruction tomography in characterising the beam behaviour, in terms of the variation of the beam size in response to variation of the quadrupole strengths. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Beating the thermal limit of qubit initialization with a Bayesian Maxwell's demon

Authors: Mark A. I. Johnson, Mateusz T. Mądzik, Fay E. Hudson, Kohei M. Itoh, Alexander M. Jakob, David N. Jamieson, Andrew Dzurak, Andrea Morello

Abstract: Fault-tolerant quantum computing requires initializing the quantum register in a well-defined fiducial state. In solid-state systems, this is typically achieved through thermalization to a cold reservoir, such that the initialization fidelity is fundamentally limited by temperature. Here, we present a method of preparing a fiducial quantum state with a confidence beyond the thermal limit. It is ba… ▽ More Fault-tolerant quantum computing requires initializing the quantum register in a well-defined fiducial state. In solid-state systems, this is typically achieved through thermalization to a cold reservoir, such that the initialization fidelity is fundamentally limited by temperature. Here, we present a method of preparing a fiducial quantum state with a confidence beyond the thermal limit. It is based on real time monitoring of the qubit through a negative-result measurement -- the equivalent of a `Maxwell's demon' that triggers the experiment only upon the appearance of a qubit in the lowest-energy state. We experimentally apply it to initialize an electron spin qubit in silicon, achieving a ground-state initialization fidelity of 98.9(4)%, corresponding to a 20$\times$ reduction in initialization error compared to the unmonitored system. A fidelity approaching 99.9% could be achieved with realistic improvements in the bandwidth of the amplifier chain or by slowing down the rate of electron tunneling from the reservoir. We use a nuclear spin ancilla, measured in quantum nondemolition mode, to prove the value of the electron initialization fidelity far beyond the intrinsic fidelity of the electron readout. However, the method itself does not require an ancilla for its execution, saving the need for additional resources. The quantitative analysis of the initialization fidelity reveals that a simple picture of spin-dependent electron tunneling does not correctly describe the data. Our digital `Maxwell's demon' can be applied to a wide range of quantum systems, with minimal demands on control and detection hardware. △ Less

Submitted 1 November, 2022; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 15 pages, 7 figures

Journal ref: Physical Review X 12, (2022)

Astronomical Pipeline Provenance: A Use Case Evaluation

Authors: Michael A. C. Johnson, Marcus Paradies, Marta Dembska, Kristen Lackeos, Hans-Rainer Klöckner, David J. Champion, Sirko Schindler

Abstract: In this decade astronomy is undergoing a paradigm shift to handle data from next generation observatories such as the Square Kilometre Array (SKA) or the Vera C. Rubin Observatory (LSST). Producing real time data streams of up to 10 TB/s and data products of the order of 600 Pbytes/year, the SKA will be the biggest civil data producing machine of the world that demands novel solutions on how these… ▽ More In this decade astronomy is undergoing a paradigm shift to handle data from next generation observatories such as the Square Kilometre Array (SKA) or the Vera C. Rubin Observatory (LSST). Producing real time data streams of up to 10 TB/s and data products of the order of 600 Pbytes/year, the SKA will be the biggest civil data producing machine of the world that demands novel solutions on how these data volumes can be stored and analysed. Through the use of complex, automated pipelines the provenance of this real time data processing is key to establish confidence within the system, its final data products, and ultimately its scientific results. The intention of this paper is to lay the foundation for making an automated provenance generation tool for astronomical/data-processing pipelines. We therefore present a use case analysis, specific to the astronomical needs which addresses the issues of trust and reproducibility as well as other ulterior use cases which are of interest to astronomers. This analysis is subsequently used as the basis to discuss the requirements, challenges, and opportunities involved in designing both the tool and the associated provenance model. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: 9 pages, 1 table

ACM Class: H.3.3; J.2

Journal ref: 13th International Workshop on Theory and Practice of Provenance (2021) 54

Subharmonic Dynamics of Wave Trains in the Korteweg-de Vries / Kuramoto-Sivashinsky Equation

Authors: Mathew A. Johnson, Wesley R. Perkins

Abstract: We study the stability and nonlinear local dynamics of spectrally stable periodic wave trains of the Korteweg-de Vries / Kuramoto-Sivashinsky equation when subjected to classes of periodic perturbations. It is known that for each $N\in\mathbb{N}$, such a $T$-periodic wave train is asymptotically stable to $NT$-periodic, i.e., subharmonic, perturbations, in the sense that initially nearby data will… ▽ More We study the stability and nonlinear local dynamics of spectrally stable periodic wave trains of the Korteweg-de Vries / Kuramoto-Sivashinsky equation when subjected to classes of periodic perturbations. It is known that for each $N\in\mathbb{N}$, such a $T$-periodic wave train is asymptotically stable to $NT$-periodic, i.e., subharmonic, perturbations, in the sense that initially nearby data will converge asymptotically to a small Galilean boost of the underlying wave, with exponential rates of decay. However, both the allowable size of initial perturbations and the exponential rates of decay depend on $N$ and, in fact, tend to zero as $N\to\infty$, leading to a lack of uniformity in such subharmonic stability results. Our goal here is to build upon a recent methodology introduced by the authors in the reaction-diffusion setting and achieve a subharmonic stability result which is uniform in $N$. This work is motivated by the dynamics of such wave trains when subjected to perturbations which are localized (i.e., integrable on the line). △ Less

Submitted 17 September, 2021; originally announced September 2021.

Comments: 26 pages, 1 figure

Multi-wavelength Optical and NIR Variability Analysis of the Blazar PKS 0027-426

Authors: E. Guise, S. F. Hönig, T. Almeyda, K. Horne, M. Kishimoto, M. Aguena, S. Allam, F. Andrade-Oliveira, J. Asorey, M. Banerji, E. Bertin, B. Boulderstone, D. Brooks, D. L. Burke, A. Carnero Rosell, D. Carollo, M. Carrasco Kind, J. Carretero, M. Costanzi, L. N. da Costa, T. M. Davis, J. De Vicente, P. Doel, S. Everett, I. Ferrero , et al. (40 additional authors not shown)

Abstract: We present multi-wavelength spectral and temporal variability analysis of PKS 0027-426 using optical griz observations from DES (Dark Energy Survey) between 2013-2018 and VOILETTE (VEILS Optical Light curves of Extragalactic TransienT Events) between 2018-2019 and near infrared (NIR) JKs observations from VEILS (VISTAExtragalactic Infrared Legacy Survey) between 2017-2019. Multiple methods of cros… ▽ More We present multi-wavelength spectral and temporal variability analysis of PKS 0027-426 using optical griz observations from DES (Dark Energy Survey) between 2013-2018 and VOILETTE (VEILS Optical Light curves of Extragalactic TransienT Events) between 2018-2019 and near infrared (NIR) JKs observations from VEILS (VISTAExtragalactic Infrared Legacy Survey) between 2017-2019. Multiple methods of cross-correlation of each combination of light curve provides measurements of possible lags between optical-optical, optical-NIR, and NIR-NIR emission, for each observation season and for the entire observational period. Inter-band time lag measurements consistently suggest either simultaneous emission or delays between emission regions on timescales smaller than the cadences of observations. The colour-magnitude relation between each combination of filters was also studied to determine the spectral behaviour of PKS 0027-426. Our results demonstrate complex colour behaviour that changes between bluer when brighter (BWB), stable when brighter (SWB) and redder when brighter (RWB) trends over different timescales and using different combinations of optical filters. Additional analysis of the optical spectra is performed to provide further understanding of this complex spectral behaviour. △ Less

Submitted 25 November, 2021; v1 submitted 30 August, 2021; originally announced August 2021.

Comments: Accepted for publication in MNRAS. 33 pages, 28 figures, 11 tables

Nonlinear Modulational Dynamics of Spectrally Stable Lugiato-Lefever Periodic Waves

Authors: Mariana Haragus, Mathew A. Johnson, Wesley R. Perkins, Björn de Rijk

Abstract: We consider the nonlinear stability of spectrally stable periodic waves in the Lugiato-Lefever equation (LLE), a damped nonlinear Schrödinger equation with forcing that arises in nonlinear optics. So far, nonlinear stability of such solutions has only been established against co-periodic perturbations by exploiting the existence of a spectral gap. In this paper, we consider perturbations which are… ▽ More We consider the nonlinear stability of spectrally stable periodic waves in the Lugiato-Lefever equation (LLE), a damped nonlinear Schrödinger equation with forcing that arises in nonlinear optics. So far, nonlinear stability of such solutions has only been established against co-periodic perturbations by exploiting the existence of a spectral gap. In this paper, we consider perturbations which are localized, i.e., integrable on the line. Such localized perturbations naturally yield the absence of a spectral gap, so we must rely on a substantially different method with origins in the stability analysis of periodic waves in reaction-diffusion systems. The relevant linear estimates have been obtained in recent work by the first three authors through a delicate decomposition of the associated linearized solution operator. Since its most critical part just decays diffusively, the nonlinear iteration can only be closed if one allows for a spatio-temporal phase modulation. However, the modulated perturbation satisfies a quasilinear equation yielding an apparent loss of regularity. To overcome this obstacle, we incorporate tame estimates on the unmodulated perturbation, which satisfies a semilinear equation in which no derivatives are lost, yet where decay is too slow to close an independent iteration scheme. We obtain nonlinear stability of periodic steady waves in the LLE against localized perturbations with precisely the same decay rates as predicted by the linear theory. △ Less

Submitted 10 February, 2022; v1 submitted 3 June, 2021; originally announced June 2021.

Comments: 32 pages, updated exposition

Solitary waves in a Whitham equation with small surface tension

Authors: Mathew A. Johnson, Tien Truong, Miles H. Wheeler

Abstract: Using a nonlocal version of the center manifold theorem and a normal form reduction, we prove the existence of small-amplitude generalized solitary-wave solutions and modulated solitary-wave solutions to the steady gravity-capillary Whitham equation with weak surface tension. Through the application of the center manifold theorem, the nonlocal equation for the solitary wave profiles is reduced to… ▽ More Using a nonlocal version of the center manifold theorem and a normal form reduction, we prove the existence of small-amplitude generalized solitary-wave solutions and modulated solitary-wave solutions to the steady gravity-capillary Whitham equation with weak surface tension. Through the application of the center manifold theorem, the nonlocal equation for the solitary wave profiles is reduced to a four-dimensional system of ODEs inheriting reversibility. Along particular parameter curves, relating directly to the classical gravity-capillary water wave problem, the associated linear operator is seen to undergo either a reversible $0^{2+}(i k_0)$ bifurcation or a reversible $(i s)^2$ bifurcation. Through a normal form transformation, the reduced system of ODEs along each relevant parameter curve is seen to be well approximated by a truncated system retaining only second-order or third-order terms. These truncated systems relate directly to systems obtained in the study of the full gravity-capillary water wave equation and, as such, the existence of generalized and modulated solitary waves for the truncated systems is guaranteed by classical works, and they are readily seen to persist as solutions of the gravity-capillary Whitham equation due to reversibility. Consequently, this work illuminates further connections between the gravity-capillary Whitham equation and the full two-dimensional gravity-capillary water wave problem. △ Less

Submitted 23 September, 2021; v1 submitted 3 March, 2021; originally announced March 2021.

Comments: 35 pages, 4 figures. Some mathematical typos fixed

Subharmonic Dynamics of Wave Trains in Reaction Diffusion Systems

Authors: Mathew A. Johnson, Wesley R. Perkins

Abstract: We investigate the stability and nonlinear local dynamics of spectrally stable wave trains in reaction-diffusion systems. For each $N\in\mathbb{N}$, such $T$-periodic traveling waves are easily seen to be nonlinearly asymptotically stable (with asymptotic phase) with exponential rates of decay when subject to $NT$-periodic, i.e., subharmonic, perturbations. However, both the allowable size of pert… ▽ More We investigate the stability and nonlinear local dynamics of spectrally stable wave trains in reaction-diffusion systems. For each $N\in\mathbb{N}$, such $T$-periodic traveling waves are easily seen to be nonlinearly asymptotically stable (with asymptotic phase) with exponential rates of decay when subject to $NT$-periodic, i.e., subharmonic, perturbations. However, both the allowable size of perturbations and the exponential rates of decay depend on $N$, and, in particular, they tend to zero as $N\to\infty$, leading to a lack of uniformity in such subharmonic stability results. In this work, we build on recent work by the authors and introduce a methodology that allows us to achieve a stability result for subharmonic perturbations which is uniform in $N$. Our work is motivated by the dynamics of such waves when subject to perturbations which are localized (i.e. integrable on the line), which has recently received considerable attention by many authors. △ Less

Submitted 11 March, 2021; v1 submitted 25 September, 2020; originally announced September 2020.

Comments: 24 pages. Updated exposition and incorporated referee comments

Linear Modulational and Subharmonic Dynamics of Spectrally Stable Lugiato-Lefever Periodic Waves

Authors: Mariana Haragus, Mathew A. Johnson, Wesley R. Perkins

Abstract: We study the linear dynamics of spectrally stable $T$-periodic stationary solutions of the Lugiato-Lefever equation (LLE), a damped nonlinear Schrödinger equation with forcing that arises in nonlinear optics. Such $T$-periodic solutions are nonlinearly stable to $NT$-periodic, i.e. subharmonic, perturbations for each $N\in\mathbb{N}$ with exponential decay rates of perturbations of the form… ▽ More We study the linear dynamics of spectrally stable $T$-periodic stationary solutions of the Lugiato-Lefever equation (LLE), a damped nonlinear Schrödinger equation with forcing that arises in nonlinear optics. Such $T$-periodic solutions are nonlinearly stable to $NT$-periodic, i.e. subharmonic, perturbations for each $N\in\mathbb{N}$ with exponential decay rates of perturbations of the form $e^{-δ_N t}$. However, both the exponential rates of decay $δ_N$ and the allowable size of the initial perturbations tend to $0$ as $N\to\infty$, so that this result is non-uniform in $N$ and, in fact, empty in the limit $N=\infty$. The primary goal of this paper is to introduce a methodology, in the context of the LLE, by which a uniform stability result for subharmonic perturbations may be achieved, at least at the linear level. The obtained uniform decay rates are shown to agree precisely with the polynomial decay rates of localized, i.e. integrable on the real line, perturbations of such spectrally stable periodic solutions of the LLE. This work both unifies and expands on several existing works in the literature concerning the stability and dynamics of such waves, and sets forth a general methodology for studying such problems in other contexts. △ Less

Submitted 14 January, 2021; v1 submitted 7 July, 2020; originally announced July 2020.

Comments: 36 pages, 2 figures. Minor typos fixed, some exposition updated

Measuring out-of-time-ordered correlation functions without reversing time evolution

Authors: Philip Daniel Blocher, Serwan Asaad, Vincent Mourik, Mark A. I. Johnson, Andrea Morello, Klaus Mølmer

Abstract: Out-of-time-ordered correlation functions (OTOCs) play a crucial role in the study of thermalization, entanglement, and quantum chaos, as they quantify the scrambling of quantum information due to complex interactions. As a consequence of their out-of-time-ordered nature, OTOCs are difficult to measure experimentally. Here we propose an OTOC measurement protocol that does not rely on the reversal… ▽ More Out-of-time-ordered correlation functions (OTOCs) play a crucial role in the study of thermalization, entanglement, and quantum chaos, as they quantify the scrambling of quantum information due to complex interactions. As a consequence of their out-of-time-ordered nature, OTOCs are difficult to measure experimentally. Here we propose an OTOC measurement protocol that does not rely on the reversal of time evolution and is easy to implement in a range of experimental settings. The protocol accounts for both pure and mixed initial states, and is applicable to systems that interact with environmental degrees of freedom. We demonstrate the application of our protocol by the characterization of scrambling in a periodically-driven spin that exhibits quantum chaos. △ Less

Submitted 26 October, 2022; v1 submitted 9 March, 2020; originally announced March 2020.

Comments: 12 pages, 6 figures. Updated figure labels and author affiliation. Made a few minor edits to the main text

Journal ref: Phys. Rev. A 106, 042429 (2022)

Stability of Traveling wave solutions of Nonlinear Dispersive equations of NLS type

Authors: Katelyn Plaisier Leisman, Jared C Bronski, Mathew A Johnson, Robert Marangell

Abstract: In this paper we present a rigorous modulational stability theory for periodic traveling wave solutions to equations of nonlinear Schrödinger (NLS) type. We first argue that, for Hamiltonian dispersive equations with a non-singular symplectic form and $d$ conserved quantities (in addition to the Hamiltonian), one expects that generically ${\mathcal L}$, the linearization around a periodic travelin… ▽ More In this paper we present a rigorous modulational stability theory for periodic traveling wave solutions to equations of nonlinear Schrödinger (NLS) type. We first argue that, for Hamiltonian dispersive equations with a non-singular symplectic form and $d$ conserved quantities (in addition to the Hamiltonian), one expects that generically ${\mathcal L}$, the linearization around a periodic traveling wave, will have a $2d$ dimensional generalized kernel, with a particular Jordan structure: The kernel of ${\mathcal L}$ is expected to be $d$ dimensional, the first generalized kernel is expected to be $d$ dimensional, and there are expected to be no higher generalized kernels. The breakup of this $2d$ dimensional kernel under perturbations arising from a change in boundary conditions dictates the modulational stability or instability of the underlying periodic traveling wave. This general picture is worked out in detail for the case of equations of NLS type. We give explicit genericity conditions that guarantee that the Jordan form is the generic one: these take the form of non-vanishing determinants of certain matrices whose entries can be expressed in terms of a finite number of moments of the traveling wave solution. Assuming that these genericity conditions are met we give a normal form for the small eigenvalues that result from the break-up of the generalized kernel, in the form of the eigenvalues of a quadratic matrix pencil. We compare these results to direct numerical simulation in a number of cases of interest: the cubic and quintic NLS equations, for focusing and defocusing nonlinearities, subject to both longitudinal and transverse perturbations. The longitudinal stability of traveling waves of the cubic NLS has been previously studied using the integrability: in this case our results agree with those in the literature. All of the remaining cases appear to be new. △ Less

Submitted 5 November, 2019; v1 submitted 11 October, 2019; originally announced October 2019.

Comments: 59 pages, 25 subfigures of 10 figures

MSC Class: 37K45 (Primary) 35B35; 35B10; 37K05 (Secondary)

Coherent electrical control of a single high-spin nucleus in silicon

Authors: Serwan Asaad, Vincent Mourik, Benjamin Joecker, Mark A. I. Johnson, Andrew D. Baczewski, Hannes R. Firgau, Mateusz T. Mądzik, Vivien Schmitt, Jarryd J. Pla, Fay E. Hudson, Kohei M. Itoh, Jeffrey C. McCallum, Andrew S. Dzurak, Arne Laucht, Andrea Morello

Abstract: Nuclear spins are highly coherent quantum objects. In large ensembles, their control and detection via magnetic resonance is widely exploited, e.g. in chemistry, medicine, materials science and mining. Nuclear spins also featured in early ideas and demonstrations of quantum information processing. Scaling up these ideas requires controlling individual nuclei, which can be detected when coupled to… ▽ More Nuclear spins are highly coherent quantum objects. In large ensembles, their control and detection via magnetic resonance is widely exploited, e.g. in chemistry, medicine, materials science and mining. Nuclear spins also featured in early ideas and demonstrations of quantum information processing. Scaling up these ideas requires controlling individual nuclei, which can be detected when coupled to an electron. However, the need to address the nuclei via oscillating magnetic fields complicates their integration in multi-spin nanoscale devices, because the field cannot be localized or screened. Control via electric fields would resolve this problem, but previous methods relied upon transducing electric signals into magnetic fields via the electron-nuclear hyperfine interaction, which severely affects the nuclear coherence. Here we demonstrate the coherent quantum control of a single antimony (spin-7/2) nucleus, using localized electric fields produced within a silicon nanoelectronic device. The method exploits an idea first proposed in 1961 but never realized experimentally with a single nucleus. Our results are quantitatively supported by a microscopic theoretical model that reveals how the purely electrical modulation of the nuclear electric quadrupole interaction, in the presence of lattice strain, results in coherent nuclear spin transitions. The spin dephasing time, 0.1 seconds, surpasses by orders of magnitude those obtained via methods that require a coupled electron spin for electrical drive. These results show that high-spin quadrupolar nuclei could be deployed as chaotic models, strain sensors and hybrid spin-mechanical quantum systems using all-electrical controls. Integrating electrically controllable nuclei with quantum dots could pave the way to scalable nuclear- and electron-spin-based quantum computers in silicon that operate without the need for oscillating magnetic fields. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: Main text and figures followed by methods, extended data, and supplementary information

Journal ref: Nature 579, 205 (2020)

Modulational Instability of Viscous Fluid Conduit Periodic Waves

Authors: Mathew A. Johnson, Wesley R. Perkins

Abstract: In this paper, we are interested in studying the modulational dynamics of interfacial waves rising buoyantly along a conduit of a viscous liquid. Formally, the behavior of modulated periodic waves on large space and time scales may be described through the use of Whitham modulation theory. The application of Whitham theory, however, is based on formal asymptotic (WKB) methods, thus removing a laye… ▽ More In this paper, we are interested in studying the modulational dynamics of interfacial waves rising buoyantly along a conduit of a viscous liquid. Formally, the behavior of modulated periodic waves on large space and time scales may be described through the use of Whitham modulation theory. The application of Whitham theory, however, is based on formal asymptotic (WKB) methods, thus removing a layer of rigor that would otherwise support their predictions. In this study, we aim at rigorously verifying the predictions of the Whitham theory, as it pertains to the modulational stability of periodic waves, in the context of the so-called conduit equation, a nonlinear dispersive PDE governing the evolution of the circular interface separating a light, viscous fluid rising buoyantly through a heavy, more viscous, miscible fluid at small Reynolds numbers. In particular, using rigorous spectral perturbation theory, we connect the predictions of Whitham theory to the rigorous spectral (in particular, modulational) stability of the underlying wave trains. This makes rigorous recent formal results on the conduit equation obtained by Maiden and Hoefer. △ Less

Submitted 4 November, 2019; v1 submitted 11 April, 2019; originally announced April 2019.

Comments: 31 pages, 2 figures. Minor typos corrected

On the bifurcation diagram of the capillary-gravity Whitham equation

Authors: Mats Ehrnström, Mathew A. Johnson, Ola I. H. Maehlen, Filippo Remonato

Abstract: We study the bifurcation of periodic travelling waves of the capillary-gravity Whitham equation. This is a nonlinear pseudo-differential equation that combines the canonical shallow water nonlinearity with the exact (unidirectional) dispersion for finite-depth capillary-gravity waves. Starting from the line of zero solutions, we give a complete description of all small periodic solutions, unimodal… ▽ More We study the bifurcation of periodic travelling waves of the capillary-gravity Whitham equation. This is a nonlinear pseudo-differential equation that combines the canonical shallow water nonlinearity with the exact (unidirectional) dispersion for finite-depth capillary-gravity waves. Starting from the line of zero solutions, we give a complete description of all small periodic solutions, unimodal as well bimodal, using simple and double bifurcation via Lyapunov--Schmidt reductions. Included in this study is the resonant case when one wavenumber divides another. Some bifurcation formulas are studied, enabling us, in almost all cases, to continue the unimodal bifurcation curves into global curves. By characterizing the range of the surface tension parameter for which the integral kernel corresponding to the linear dispersion operator is completely monotone (and therefore positive and convex; the threshold value for this to happen turns out to be \(T = \frac{4}{π^2}\), not the critical Bond number \(\frac{1}{3}\)), we are able to say something about the nodal properties of solutions, even in the presence of surface tension. Finally, we present a few general results for the equation and discuss, in detail, the complete bifurcation diagram as far as it is known from analytical and numerical evidence. Interestingly, we find, analytically, secondary bifurcation curves connecting different branches of solutions; and, numerically, that all supercritical waves preserve their basic nodal structure, converging asymptotically in \(L^2(\SM)\) (but not in \(L^\infty\)) towards one of the two constant solution curves. △ Less

Submitted 11 January, 2019; originally announced January 2019.

MSC Class: 35Q35; 37K50; 76N10

Electron spin relaxation of single phosphorus donors in metal-oxide-semiconductor nanoscale devices

Authors: Stefanie B. Tenberg, Serwan Asaad, Mateusz T. Mądzik, Mark A. I. Johnson, Benjamin Joecker, Arne Laucht, Fay E. Hudson, Kohei M. Itoh, A. Malwin Jakob, Brett C. Johnson, David N. Jamieson, Jeffrey C. McCallum, Andrew S. Dzurak, Robert Joynt, Andrea Morello

Abstract: We analyze the electron spin relaxation rate $1/T_1$ of individual ion-implanted $^{31}$P donors, in a large set of metal-oxide-semiconductor (MOS) silicon nanoscale devices, with the aim of identifying spin relaxation mechanisms peculiar to the environment of the spins. The measurements are conducted at low temperatures ($T\approx 100$~mK), as a function of external magnetic field $B_0$ and donor… ▽ More We analyze the electron spin relaxation rate $1/T_1$ of individual ion-implanted $^{31}$P donors, in a large set of metal-oxide-semiconductor (MOS) silicon nanoscale devices, with the aim of identifying spin relaxation mechanisms peculiar to the environment of the spins. The measurements are conducted at low temperatures ($T\approx 100$~mK), as a function of external magnetic field $B_0$ and donor electrochemical potential $μ_{\rm D}$. We observe a magnetic field dependence of the form $1/T_1\propto B_0^5$ for $B_0\gtrsim 3\,$ T, corresponding to the phonon-induced relaxation typical of donors in the bulk. However, the relaxation rate varies by up to two orders of magnitude between different devices. We attribute these differences to variations in lattice strain at the location of the donor. For $B_0\lesssim 3\,$T, the relaxation rate changes to $1/T_1\propto B_0$ for two devices. This is consistent with relaxation induced by evanescent-wave Johnson noise created by the metal structures fabricated above the donors. At such low fields, where $T_1>1\,$s, we also observe and quantify the spurious increase of $1/T_1$ when the electrochemical potential of the spin excited state $|\uparrow\rangle$ comes in proximity to empty states in the charge reservoir, leading to spin-dependent tunneling that resets the spin to $|\downarrow\rangle$. These results give precious insights into the microscopic phenomena that affect spin relaxation in MOS nanoscale devices, and provide strategies for engineering spin qubits with improved spin lifetimes. △ Less

Submitted 26 March, 2019; v1 submitted 17 December, 2018; originally announced December 2018.

Comments: 14 pages, 8 figures. Version 2 has extensive text revisions, typo corrections, and two extra authors

Journal ref: Phys. Rev. B 99, 205306 (2019)

The Plane's The Thing: The Case for Wide-Fast-Deep Coverage of the Galactic Plane and Bulge

Authors: Jay Strader, Elias Aydi, Christopher Britt, Adam Burgasser, Laura Chomiuk, Will Clarkson, Brian D. Fields, Poshak Gandhi, Leo Girardi, John Gizis, Jacob Hogan, Michael A. C. Johnson, James Lauroesch, Michael Liu, Tom Maccarone, Peregrine McGehee, Dante Minniti, Koji Mukai, C. Tanner Murphey, Alexandre Roman-Lopez, Simone Scaringi, Jennifer Sobeck, Kirill Sokolovsky, Xilu Wang

Abstract: We argue that the exclusion of the Galactic Plane and Bulge from the uniform wide-fast-deep (WFD) LSST survey cadence is fundamentally inconsistent with two of the main science drivers of LSST: Mapping the Milky Way and Exploring the Transient Optical Sky. We outline the philosophical basis for this claim and then describe a number of important science goals that can only be addressed by WFD-like… ▽ More We argue that the exclusion of the Galactic Plane and Bulge from the uniform wide-fast-deep (WFD) LSST survey cadence is fundamentally inconsistent with two of the main science drivers of LSST: Mapping the Milky Way and Exploring the Transient Optical Sky. We outline the philosophical basis for this claim and then describe a number of important science goals that can only be addressed by WFD-like coverage of the Plane and Bulge. △ Less

Submitted 29 November, 2018; originally announced November 2018.

Comments: white paper on LSST cadence optimization, 9 pages

On the Dynamics of Traveling Fronts Arising in Nanoscale Pattern Formation

Authors: Mathew A. Johnson, Gregory D. Lyng, Connor Smith

Abstract: We study the stability and dynamics of traveling-front solutions of a modified Kuramoto--Sivashinsky equation arising in the modeling of nanoscale ripple patterns that form when a nominally flat solid surface is bombarded with a broad ion beam at an oblique angle of incidence. Structurally, the linearized operators associated with these fronts have unstable essential spectrum---corresponding to in… ▽ More We study the stability and dynamics of traveling-front solutions of a modified Kuramoto--Sivashinsky equation arising in the modeling of nanoscale ripple patterns that form when a nominally flat solid surface is bombarded with a broad ion beam at an oblique angle of incidence. Structurally, the linearized operators associated with these fronts have unstable essential spectrum---corresponding to instability of the spatially asymptotic states---and stable point spectrum---corresponding to stability of the transition profile of the front. We show that these waves are linearly orbitally asymptotically stable in appropriate exponentially weighted spaces. While the technical device of exponential weights allows us to accommodate the unstable essential spectrum of individual waves in our linear analysis, it does not shed light on the long-time pattern formation that is observed experimentally and in numerical simulations. To begin to address this issue, we consider a periodic array of unstable front and back solutions. While not an exact solution of the governing equation, this periodic pattern mimics experimentally observed phenomena. Our numerical experiments suggest that the convecting instabilities associated with each individual wave are damped as they pass through transition layers and that this stabilization mechanism underlies the pattern formation seen in experiments. △ Less

Submitted 2 July, 2019; v1 submitted 4 October, 2018; originally announced October 2018.

Comments: 31 pages, 14 Figures. Minor update in exposition

Prospecting Period Measurements with LSST - Low Mass X-ray Binaries as a Test Case

Authors: Michael A. C. Johnson, Poshak Gandhi, Adriane P. Chapman, Luc Moreau, Philip A. Charles, William I. Clarkson, Adam B. Hill

Abstract: The Large Synoptic Survey Telescope (LSST) will provide for unbiased sampling of variability properties of objects with $r$ mag $<$ 24. This should allow for those objects whose variations reveal their orbital periods ($P_{orb}$), such as low mass X-ray binaries (LMXBs) and related objects, to be examined in much greater detail and with uniform systematic sampling. However, the baseline LSST obser… ▽ More The Large Synoptic Survey Telescope (LSST) will provide for unbiased sampling of variability properties of objects with $r$ mag $<$ 24. This should allow for those objects whose variations reveal their orbital periods ($P_{orb}$), such as low mass X-ray binaries (LMXBs) and related objects, to be examined in much greater detail and with uniform systematic sampling. However, the baseline LSST observing strategy has temporal sampling that is not optimised for such work in the Galaxy. Here we assess four candidate observing strategies for measurement of $P_{orb}$ in the range 10 minutes to 50 days. We simulate multi-filter quiescent LMXB lightcurves including ellipsoidal modulation and stochastic flaring, and then sample these using LSST's operations simulator (OpSim) over the (mag, $P_{orb}$) parameter space, and over five sightlines sampling a range of possible reddening values. The percentage of simulated parameter space with correctly returned periods ranges from $\sim$23 %, for the current baseline strategy, to $\sim$70 % for the two simulated specialist strategies. Convolving these results with a $P_{orb}$ distribution, a modelled Galactic spatial distribution and reddening maps, we conservatively estimate that the most recent version of the LSST baseline strategy will allow $P_{orb}$ determination for $\sim$18 % of the Milky Way's LMXB population, whereas strategies that do not reduce observations of the Galactic Plane can improve this dramatically to $\sim$32 %. This increase would allow characterisation of the full binary population by breaking degeneracies between suggested $P_{orb}$ distributions in the literature. Our results can be used in the ongoing assessment of the effectiveness of various potential cadencing strategies. △ Less

Submitted 8 January, 2019; v1 submitted 24 September, 2018; originally announced September 2018.

Comments: Replacement after addressing minor corrections from the referee - mainly improvements in clarification

Generalized Solitary Waves in the Gravity-Capillary Whitham Equation

Authors: Mathew A. Johnson, J. Douglas Wright

Abstract: We study the existence of traveling wave solutions to a unidirectional shallow water model which incorporates the full linear dispersion relation for both gravitational and capillary restoring forces. Using functional analytic techniques, we show that for small surface tension (corresponding to Bond numbers between $0$ and ${1}/{3}$) there exists small amplitude solitary waves that decay to asympt… ▽ More We study the existence of traveling wave solutions to a unidirectional shallow water model which incorporates the full linear dispersion relation for both gravitational and capillary restoring forces. Using functional analytic techniques, we show that for small surface tension (corresponding to Bond numbers between $0$ and ${1}/{3}$) there exists small amplitude solitary waves that decay to asymptotically small periodic waves at spatial infinity. The size of the oscillations in the far field are shown to be small beyond all algebraic orders in the amplitude of the wave. We also present numerical evidence, based on the recent analytical work of Hur \& Johnson, that the asymptotic end states are modulationally stable for all Bond numbers between $0$ and $1/3$. △ Less

Submitted 30 July, 2018; originally announced July 2018.

MSC Class: 76B45

Gaia DR2 Distances and Peculiar Velocities for Galactic Black Hole Transients

Authors: Poshak Gandhi, Anjali Rao, Michael A. C. Johnson, John A. Paice, Thomas J. Maccarone

Abstract: We report on a first census of Galactic black hole X-ray binary (BHXRB) properties with the second data release (DR2) of {\em Gaia}, focusing on dynamically confirmed and strong candidate black hole transients. DR2 provides five-parameter astrometric solutions including position, parallax and proper motion for 11 of a sample of 24 systems. Distance estimates are tested with parallax inversion as w… ▽ More We report on a first census of Galactic black hole X-ray binary (BHXRB) properties with the second data release (DR2) of {\em Gaia}, focusing on dynamically confirmed and strong candidate black hole transients. DR2 provides five-parameter astrometric solutions including position, parallax and proper motion for 11 of a sample of 24 systems. Distance estimates are tested with parallax inversion as well as Bayesian inference. We derive an empirically motivated characteristic scale length of $L$=2.17$\pm$0.12 kpc for this BHXRB population to infer distances based upon an exponentially decreasing space density prior. Geometric DR2 parallaxes provide new, independent distance estimates, but the faintness of this population in quiescence results in relatively large fractional distance uncertainties. Despite this, DR2 estimates generally agree with literature distances. The most discrepant case is BW Cir, for which detailed studies of the donor star have suggested a distant location at >~25 kpc. A large DR2 measured parallax and relatively high proper motion instead prefer significantly smaller distances, suggesting that the source may instead be amongst the nearest of XRBs. However, both distances create problems for interpretation of the source, and follow-up data are required to resolve its true nature. DR2 also provides a first distance estimate to one source, MAXI J1820+070, and novel proper motion estimates for 7 sources. Peculiar velocities relative to Galactic rotation exceed $\sim$ 50 km s$^{-1}$ for the bulk of the sample, with a median system kinetic energy of peculiar motion of $\sim$ 5 $\times$ 10$^{47}$ erg. BW Cir could be a new high-velocity BHXRB if its astrometry is confirmed. A putative anti-correlation between peculiar velocity and black hole mass is found, as expected in mass-dependent BH kick formation channels, but this trend remains weak in the DR2 data. △ Less

Submitted 5 February, 2019; v1 submitted 30 April, 2018; originally announced April 2018.

Comments: MNRAS in press

Spectral Stability of Inviscid Roll Waves

Authors: Mathew A. Johnson, Pascal Noble, L. Miguel Rodrigues, Zhao Yang, Kevin Zumbrun

Abstract: We carry out a systematic analytical and numerical study of spectral stability of discontinuous roll wave solutions of the inviscid Saint Venant equations, based on a periodic Evans-Lopatinski determinant analogous to the periodic Evans function of Gardner in the (smooth) viscous case, obtaining a complete spectral stability diagram useful in hydraulic engineering and related applications. In part… ▽ More We carry out a systematic analytical and numerical study of spectral stability of discontinuous roll wave solutions of the inviscid Saint Venant equations, based on a periodic Evans-Lopatinski determinant analogous to the periodic Evans function of Gardner in the (smooth) viscous case, obtaining a complete spectral stability diagram useful in hydraulic engineering and related applications. In particular, we obtain an explicit low-frequency stability boundary, which, moreover, matches closely with its (numerically-determined) counterpart in the viscous case. This is seen to be related to but not implied by the associated formal first-order Whitham modulation equations. △ Less

Submitted 9 March, 2018; originally announced March 2018.

Comments: 45 pages, 11 figures

Numerical Bifurcation and Spectral Stability of Wavetrains in Bidirectional Whitham Models

Authors: Kyle M. Claassen, Mathew A. Johnson

Abstract: We consider several different bidirectional Whitham equations that have recently appeared in the literature. Each of these models combine the full two-way dispersion relation from the incompressible Euler equations with a canonical shallow water nonlinearity, providing nonlocal model equations that may be expected to exhibit some of the interesting high-frequency phenomena present in the Euler equ… ▽ More We consider several different bidirectional Whitham equations that have recently appeared in the literature. Each of these models combine the full two-way dispersion relation from the incompressible Euler equations with a canonical shallow water nonlinearity, providing nonlocal model equations that may be expected to exhibit some of the interesting high-frequency phenomena present in the Euler equations that standard "long-wave" theories fail to capture. Of particular interest here is the existence and stability of periodic traveling wave solutions in such models. Using numerical bifurcation techniques we construct global bifurcation diagrams for each system and compare the global structure of branches, together with the possibility of bifurcation branches terminating in a "highest" singular (peaked/cusped) wave. We also numerically approximate the stability spectrum along these bifurcation branches and compare the stability predictions of these models. Our results confirm a number of analytical results concerning the stability of asymptotically small waves in these models and provide new insights into the existence and stability of large amplitude waves. △ Less

Submitted 9 April, 2018; v1 submitted 26 October, 2017; originally announced October 2017.

Comments: 36 pages, 19 figures

AMPNet: Asynchronous Model-Parallel Training for Dynamic Neural Networks

Authors: Alexander L. Gaunt, Matthew A. Johnson, Maik Riechert, Daniel Tarlow, Ryota Tomioka, Dimitrios Vytiniotis, Sam Webster

Abstract: New types of machine learning hardware in development and entering the market hold the promise of revolutionizing deep learning in a manner as profound as GPUs. However, existing software frameworks and training algorithms for deep learning have yet to evolve to fully leverage the capability of the new wave of silicon. We already see the limitations of existing algorithms for models that exploit s… ▽ More New types of machine learning hardware in development and entering the market hold the promise of revolutionizing deep learning in a manner as profound as GPUs. However, existing software frameworks and training algorithms for deep learning have yet to evolve to fully leverage the capability of the new wave of silicon. We already see the limitations of existing algorithms for models that exploit structured input via complex and instance-dependent control flow, which prohibits minibatching. We present an asynchronous model-parallel (AMP) training algorithm that is specifically motivated by training on networks of interconnected devices. Through an implementation on multi-core CPUs, we show that AMP training converges to the same accuracy as conventional synchronous training algorithms in a similar number of epochs, but utilizes the available hardware more efficiently even for small minibatch sizes, resulting in significantly shorter overall training times. Our framework opens the door for scaling up a new class of deep learning models that cannot be efficiently trained today. △ Less

Submitted 22 June, 2017; v1 submitted 27 May, 2017; originally announced May 2017.

Comments: 17 pages, 13 figures

Existence of a highest wave in a fully dispersive two-way shallow water model

Authors: Mats Ehrnström, Mathew A. Johnson, Kyle M. Claassen

Abstract: We consider the existence of periodic traveling waves in a bidirectional Whitham equation, combining the full two-way dispersion relation from the incompressible Euler equations with a canonical shallow water nonlinearity. Of particular interest is the existence of a highest, cusped, traveling wave solution, which we obtain as a limiting case at the end of the main bifurcation branch of $2π$-perio… ▽ More We consider the existence of periodic traveling waves in a bidirectional Whitham equation, combining the full two-way dispersion relation from the incompressible Euler equations with a canonical shallow water nonlinearity. Of particular interest is the existence of a highest, cusped, traveling wave solution, which we obtain as a limiting case at the end of the main bifurcation branch of $2π$-periodic traveling wave solutions continuing from the zero state. Unlike the unidirectional Whitham equation, containing only one branch of the full Euler dispersion relation, where such a highest wave behaves like $|x|^{1/2}$ near its crest, the cusped waves obtained here behave like $|x\log|x||$. Although the linear operator involved in this equation can be easily represented in terms of an integral operator, it maps continuous functions out of the Hölder and Lipschitz scales of function spaces by introducing logarithmic singularities. Since the nonlinearity is also of higher order than that of the unidirectional Whitham equation, several parts of our proofs and results deviate from those of the corresponding unidirectional equation, with the analysis of the logarithmic singularity being the most subtle component. This paper is part of a longer research programme for understanding the interplay between nonlinearities and dispersion in the formation of large-amplitude waves and their singularities. △ Less

Submitted 12 April, 2018; v1 submitted 8 October, 2016; originally announced October 2016.

Comments: 34 pages, 5 figures

Oscillation estimates of eigenfunctions via the combinatorics of noncrossing partitions

Authors: Vera Mikyoung Hur, Mathew A. Johnson, Jeremy L. Martin

Abstract: We study oscillations in the eigenfunctions for a fractional Schrödinger operator on the real line. An argument in the spirit of Courant's nodal domain theorem applies to an associated local problem in the upper half plane and provides a bound on the number of nodal domains for the extensions of the eigenfunctions. Using the combinatorial properties of noncrossing partitions, we turn the nodal dom… ▽ More We study oscillations in the eigenfunctions for a fractional Schrödinger operator on the real line. An argument in the spirit of Courant's nodal domain theorem applies to an associated local problem in the upper half plane and provides a bound on the number of nodal domains for the extensions of the eigenfunctions. Using the combinatorial properties of noncrossing partitions, we turn the nodal domain bound into an estimate for the number of sign changes in the eigenfunctions. We discuss applications in the periodic setting and the Steklov problem on planar domains. △ Less

Submitted 5 September, 2017; v1 submitted 7 September, 2016; originally announced September 2016.

Comments: 20 pages

MSC Class: 35P05; 05A18

Journal ref: Discrete Analysis 2017:13, 20 pp

Near-Field studies of anisotropic variations and temperature induced structural changes in a supported single lipid bilayer

Authors: Merrell A. Johnson, Ricardo S. Decca

Abstract: Temperature controlled Polarization Modulation Near-Field Scanning Optical Microscopy (PM-NSOM) measurements of a single supported $L_{β^{\prime}}$ 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid bilayer are presented. The effective retardance ($ΔS = \frac{2 π(n_e-n_o)t}λ$), where $t$ is the thickness of the bilayer and $λ$ is the wavelength of light used and the direction of the projecti… ▽ More Temperature controlled Polarization Modulation Near-Field Scanning Optical Microscopy (PM-NSOM) measurements of a single supported $L_{β^{\prime}}$ 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid bilayer are presented. The effective retardance ($ΔS = \frac{2 π(n_e-n_o)t}λ$), where $t$ is the thickness of the bilayer and $λ$ is the wavelength of light used and the direction of the projection of the acyl chains ($θ$) were measured simultaneously. From $ΔS$ the birefringence ($n_e-n_o$) of the bilayer was determined. A change of $ΔS$ $\sim$ (3.8 $\pm$ 0.3)~mrad at the transition temperature $T_m \sim$ 41$^{\circ}$C between the gel $L_{β^{\prime}}$ to liquid disorder $L_α$ was observed in a single planar bilayer. This agrees well with previous values of ($n_e-n_o$) in the $L_{β^{\prime}}$ phase and translates to an assumed $\langle φ\rangle $ $\sim$ 32$^{\circ}$ when T $<T_m$ and 0$^{\circ}$ when T $>T_m$. Evidence of supper heating and supper cooling are presented, along with a discussion of the behavior that occurs around $T_m$. △ Less

Submitted 23 August, 2016; originally announced August 2016.

Comments: 8 pages, 7 figures

Journal ref: Phys. Rev. E 101, 032416 (2020)

Nondegeneracy and Stability of Antiperiodic Bound States for Fractional Nonlinear Schrödinger Equations

Authors: Kyle M. Claassen, Mathew A. Johnson

Abstract: We consider the existence and stability of real-valued, spatially antiperiodic standing wave solutions to a family of nonlinear Schrödinger equations with fractional dispersion and power-law nonlinearity. As a key technical result, we demonstrate that the associated linearized operator is nondegenerate when restricted to antiperiodic perturbations, i.e. that its kernel is generated by the translat… ▽ More We consider the existence and stability of real-valued, spatially antiperiodic standing wave solutions to a family of nonlinear Schrödinger equations with fractional dispersion and power-law nonlinearity. As a key technical result, we demonstrate that the associated linearized operator is nondegenerate when restricted to antiperiodic perturbations, i.e. that its kernel is generated by the translational and gauge symmetries of the governing evolution equation. In the process, we provide a characterization of the antiperiodic ground state eigenfunctions for linear fractional Schrödinger operators on $\mathbb{R}$ with real-valued, periodic potentials as well as a Sturm-Liouville type oscillation theory for the higher antiperiodic eigenfunctions. △ Less

Submitted 20 November, 2018; v1 submitted 20 June, 2016; originally announced June 2016.

Comments: 46 pages, 2 figures

Winning versus losing during gambling and its neural correlates

Authors: Pierre Sacré, Matthew S. D. Kerr, Sandya Subramanian, Kevin Kahn, Jorge Gonzalez-Martinez, Matthew A. Johnson, John T. Gale, Sridevi V. Sarma

Abstract: Humans often make decisions which maximize an internal utility function. For example, humans often maximize their expected reward when gambling and this is considered as a "rational" decision. However, humans tend to change their betting strategies depending on how they "feel". If someone has experienced a losing streak, they may "feel" that they are more likely to win on the next hand even though… ▽ More Humans often make decisions which maximize an internal utility function. For example, humans often maximize their expected reward when gambling and this is considered as a "rational" decision. However, humans tend to change their betting strategies depending on how they "feel". If someone has experienced a losing streak, they may "feel" that they are more likely to win on the next hand even though the odds of the game have not changed. That is, their decisions are driven by their emotional state. In this paper, we investigate how the human brain responds to wins and losses during gambling. Using a combination of local field potential recordings in human subjects performing a financial decision-making task, spectral analyses, and non-parametric cluster statistics, we investigated whether neural responses in different cognitive and limbic brain areas differ between wins and losses after decisions are made. In eleven subjects, the neural activity modulated significantly between win and loss trials in one brain region: the anterior insula ($p=0.01$). In particular, gamma activity (30-70 Hz) increased in the anterior insula when subjects just realized that they won. Modulation of metabolic activity in the anterior insula has been observed previously in functional magnetic resonance imaging studies during decision making and when emotions are elicited. However, our study is able to characterize temporal dynamics of electrical activity in this brain region at the millisecond resolution while decisions are made and after outcomes are revealed. △ Less

Submitted 10 February, 2016; originally announced February 2016.

Note on the stability of viscous roll-waves

Authors: Blake Barker, Mathew A. Johnson, Pascal Noble, L. Miguel Rodrigues, Kevin Zumbrun

Abstract: In this note, we announce a complete classification of stability of periodic roll-wave solutions of the viscous shallow-water equations, from their onset at Froude number $F\approx 2$ up to the infinite-Froude limit. For intermediate Froude numbers, we obtain numerically a particularly simple power-law relation between $F$ and the boundaries of the region of stable periods, that appears potentiall… ▽ More In this note, we announce a complete classification of stability of periodic roll-wave solutions of the viscous shallow-water equations, from their onset at Froude number $F\approx 2$ up to the infinite-Froude limit. For intermediate Froude numbers, we obtain numerically a particularly simple power-law relation between $F$ and the boundaries of the region of stable periods, that appears potentially useful in hydraulic engineering applications. In the asymptotic regime $F\to 2$ (onset), we provide an analytic expression of the stability boundaries whereas in the limit $F\to\infty$, we show that roll-waves are always unstable. △ Less

Submitted 7 June, 2016; v1 submitted 5 October, 2015; originally announced October 2015.

Comments: 6 pages, 2 figures

Modulational Instability in the Whitham Equation with Surface Tension and Vorticity

Authors: Vera Mikyoung Hur, Mathew A. Johnson

Abstract: We study modulational stability and instability in the Whitham equation, combining the dispersion relation of water waves and a nonlinearity of the shallow water equations, and modified to permit the effects of surface tension and constant vorticity. When the surface tension coefficient is large, we show that a periodic traveling wave of sufficiently small amplitude is unstable to long wavelength… ▽ More We study modulational stability and instability in the Whitham equation, combining the dispersion relation of water waves and a nonlinearity of the shallow water equations, and modified to permit the effects of surface tension and constant vorticity. When the surface tension coefficient is large, we show that a periodic traveling wave of sufficiently small amplitude is unstable to long wavelength perturbations if the wave number is greater than a critical value, and stable otherwise, similarly to the Benjamin-Feir instability of gravity waves. In the case of weak surface tension, we find intervals of stable and unstable wave numbers, whose boundaries are associated with the extremum of the group velocity, the resonance between the first and second harmonics, the resonance between long and short waves, and a resonance between dispersion and the nonlinearity. For each constant vorticity we show that a periodic traveling wave of sufficiently small amplitude is unstable if the wave number is greater than a critical value, and stable otherwise. Moreover it can be made stable for a sufficiently large vorticity. The results agree with those based upon numerical computations or formal multiple-scale expansions for the physical problem. △ Less

Submitted 26 August, 2015; v1 submitted 24 June, 2015; originally announced June 2015.

Comments: 17 pages, 4 figures. References updated

Stability of viscous St. Venant roll-waves: from onset to infinite-Froude number limit

Authors: Blake Barker, Mathew A. Johnson, Pascal Noble, L. Miguel Rodrigues, Kevin Zumbrun

Abstract: We study the spectral stability of roll-wave solutions of the viscous St. Venant equations modeling inclined shallow-water flow, both at onset in the small-Froude number or "weakly unstable" limit $F\to 2^+$ and for general values of the Froude number $F$, including the limit $F\to +\infty$. In the former, $F\to 2^+$, limit, the shallow water equations are formally approximated by a Korteweg de Vr… ▽ More We study the spectral stability of roll-wave solutions of the viscous St. Venant equations modeling inclined shallow-water flow, both at onset in the small-Froude number or "weakly unstable" limit $F\to 2^+$ and for general values of the Froude number $F$, including the limit $F\to +\infty$. In the former, $F\to 2^+$, limit, the shallow water equations are formally approximated by a Korteweg de Vries/Kuramoto-Sivashinsky (KdV-KS) equation that is a singular perturbation of the standard Korteweg de Vries (KdV) equation modeling horizontal shallow water flow. Our main analytical result is to rigorously validate this formal limit, showing that stability as $F\to 2^+$ is equivalent to stability of the corresponding KdV-KS waves in the KdV limit. Together with recent results obtained for KdV-KS by Johnson--Noble--Rodrigues--Zumbrun and Barker, this gives not only the first rigorous verification of stability for any single viscous St. Venant roll wave, but a complete classification of stability in the weakly unstable limit. In the remainder of the paper, we investigate numerically and analytically the evolution of the stability diagram as Froude number increases to infinity. Notably, we find transition at around $F=2.3$ from weakly unstable to different, large-$F$ behavior, with stability determined by simple power law relations. The latter stability criteria are potentially useful in hydraulic engineering applications, for which typically $2.5\leq F\leq 6.0$. △ Less

Submitted 31 March, 2015; v1 submitted 3 March, 2015; originally announced March 2015.

Comments: 38 pages. Important upgrade made in the slope condition, and added Numerical Observation 2 that was inadvertently left out in previous draft

Modulational Instability in Equations of KdV Type

Authors: Jared C. Bronski, Vera Mikyoung Hur, Mathew A. Johnson

Abstract: It is a matter of experience that nonlinear waves in dispersive media, propagating primarily in one direction, may appear periodic in small space and time scales, but their characteristics --- amplitude, phase, wave number, etc. --- slowly vary in large space and time scales. In the 1970's, Whitham developed an asymptotic (WKB) method to study the effects of small "modulations" on nonlinear period… ▽ More It is a matter of experience that nonlinear waves in dispersive media, propagating primarily in one direction, may appear periodic in small space and time scales, but their characteristics --- amplitude, phase, wave number, etc. --- slowly vary in large space and time scales. In the 1970's, Whitham developed an asymptotic (WKB) method to study the effects of small "modulations" on nonlinear periodic wave trains. Since then, there has been a great deal of work aiming at rigorously justifying the predictions from Whitham's formal theory. We discuss recent advances in the mathematical understanding of the dynamics, in particular, the instability of slowly modulated wave trains for nonlinear dispersive equations of KdV type. △ Less

Submitted 12 January, 2015; originally announced January 2015.

Comments: 40 pages. To appear in upcoming title in Lecture Notes in Physics

Modulational Instability in the Whitham Equation for Water Waves

Authors: Vera Mikyoung Hur, Mathew A. Johnson

Abstract: We show that periodic traveling waves with sufficiently small amplitudes of the Whitham equation, which incorporates the dispersion relation of surface water waves and the nonlinearity of the shallow water equations,are spectrally unstable to long wavelengths perturbations if the wave number is greater than a critical value, bearing out the Benjamin-Feir instability of Stokes waves; they are spect… ▽ More We show that periodic traveling waves with sufficiently small amplitudes of the Whitham equation, which incorporates the dispersion relation of surface water waves and the nonlinearity of the shallow water equations,are spectrally unstable to long wavelengths perturbations if the wave number is greater than a critical value, bearing out the Benjamin-Feir instability of Stokes waves; they are spectrally stable to square integrable perturbations otherwise. The proof involves a spectral perturbation of the associated linearized operator with respect to the Floquet exponent and the small amplitude parameter. We extend the results to related, nonlinear dispersive equations. △ Less

Submitted 14 May, 2014; v1 submitted 5 December, 2013; originally announced December 2013.

Comments: 20 pages, 2 figures. Replacement incorporates referees' comments

Stability of Periodic Traveling Waves for Nonlinear Dispersive Equations

Authors: Vera Mikyoung Hur, Mathew A. Johnson

Abstract: We consider the stability and instability of periodic travling waves for Korteweg-de Vries type equations with fractional dispersion and other nonlinear dispersive equations. We establish that a constrained minimizer for the related variational problem is nonlinearly stable to period preserving perturbations. We also discuss when the associated linearized equation admits exponentially growing solu… ▽ More We consider the stability and instability of periodic travling waves for Korteweg-de Vries type equations with fractional dispersion and other nonlinear dispersive equations. We establish that a constrained minimizer for the related variational problem is nonlinearly stable to period preserving perturbations. We also discuss when the associated linearized equation admits exponentially growing solutions. The proof utilizes the variational structure of the equation. △ Less

Submitted 12 January, 2015; v1 submitted 19 March, 2013; originally announced March 2013.

Comments: 27 pages. Missing assumption added to stability theorem. arXiv admin note: text overlap with arXiv:1303.2605

Behavior of periodic solutions of viscous conservation laws under localized and nonlocalized perturbations

Authors: Mathew A. Johnson, Pascal Noble, L. Miguel Rodrigues, Kevin Zumbrun

Abstract: We establish nonlinear stability and asymptotic behavior of traveling periodic waves of viscous conservation laws under localized perturbations or nonlocalized perturbations asymptotic to constant shifts in phase, showing that long-time behavior is governed by an associated second-order formal Whitham modulation system. A key point is to identify the way in which initial perturbations translate to… ▽ More We establish nonlinear stability and asymptotic behavior of traveling periodic waves of viscous conservation laws under localized perturbations or nonlocalized perturbations asymptotic to constant shifts in phase, showing that long-time behavior is governed by an associated second-order formal Whitham modulation system. A key point is to identify the way in which initial perturbations translate to initial data for this formal system, a task accomplished by detailed estimates on the linearized solution operator about the background wave. Notably, our approach gives both a common theoretical treatment and a complete classification in terms of "phase-coupling" or "-decoupling" of general systems of conservation or balance laws, encompassing cases that had previously been studied separately or not at all. At the same time, our refined description of solutions gives the new result of nonlinear asymptotic stability with respect to localized perturbations in the phase-decoupled case, further distinguishing behavior in the different cases. An interesting technical aspect of our analysis is that for systems of conservation laws the Whitham modulation description is of system rather than scalar form, as a consequence of which renormalization methods such as have been used to treat the reaction-diffusion case in general do not seem to apply. △ Less

Submitted 9 November, 2012; originally announced November 2012.

Comments: 63 pages

Stability of Small Periodic Waves in Fractional KdV Type Equations

Authors: Mathew A. Johnson

Abstract: We consider the effects of varying dispersion and nonlinearity on the stability of periodic traveling wave solutions of nonlinear PDE of KdV-type, including generalized KdV and Benjamin-Ono equations. In this investigation, we consider the spectral stability of such solutions that arise as small perturbations of an equilibrium state. A key feature of our analysis is the development of a nonlocal F… ▽ More We consider the effects of varying dispersion and nonlinearity on the stability of periodic traveling wave solutions of nonlinear PDE of KdV-type, including generalized KdV and Benjamin-Ono equations. In this investigation, we consider the spectral stability of such solutions that arise as small perturbations of an equilibrium state. A key feature of our analysis is the development of a nonlocal Floquet-like theory that is suitable to analyze the $L^2(\RM)$ spectrum of the associated linearized operators. Using spectral perturbation theory then, we derive a relationship between the power of the nonlinearity and the symbol of the fractional dispersive operator that determines the spectral stability and instability to arbitrary small localized perturbations. △ Less

Submitted 19 March, 2013; v1 submitted 8 October, 2012; originally announced October 2012.

Comments: An error in Lemma 2.1 has been corrected, and the remaining analysis and results have been modified accordingly

An Instability Index Theory for Quadratic Pencils and Applications

Authors: Jared Bronski, Mathew A. Johnson, Todd Kapitula

Abstract: Primarily motivated by the stability analysis of nonlinear waves in second-order in time Hamiltonian systems, in this paper we develop an instability index theory for quadratic operator pencils acting on a Hilbert space. In an extension of the known theory for linear pencils, explicit connections are made between the number of eigenvalues of a given quadratic operator pencil with positive real par… ▽ More Primarily motivated by the stability analysis of nonlinear waves in second-order in time Hamiltonian systems, in this paper we develop an instability index theory for quadratic operator pencils acting on a Hilbert space. In an extension of the known theory for linear pencils, explicit connections are made between the number of eigenvalues of a given quadratic operator pencil with positive real parts to spectral information about the individual operators comprising the coefficients of the spectral parameter in the pencil. As an application, we apply the general theory developed here to yield spectral and nonlinear stability/instability results for abstract second-order in time wave equations. More specifically, we consider the problem of the existence and stability of spatially periodic waves for the "good" Boussinesq equation. In the analysis our instability index theory provides an explicit, and somewhat surprising, connection between the stability of a given periodic traveling wave solution of the "good" Boussinesq equation and the stability of the same periodic profile, but with different wavespeed, in the nonlinear dynamics of a related generalized Korteweg-de Vries equation. △ Less

Submitted 16 July, 2012; originally announced July 2012.

Comments: 25 pages, 3 figures

Nonlinear modulational stability of periodic traveling-wave solutions of the generalized Kuramoto-Sivashinsky equation

Authors: Blake Barker, Mathew A. Johnson, Pascal Noble, L. Miguel Rodrigues, Kevin Zumbrun

Abstract: In this paper we consider the spectral and nonlinear stability of periodic traveling wave solutions of a generalized Kuramoto-Sivashinsky equation. In particular, we resolve the long-standing question of nonlinear modulational stability by demonstrating that spectrally stable waves are nonlinearly stable when subject to small localized (integrable) perturbations. Our analysis is based upon detaile… ▽ More In this paper we consider the spectral and nonlinear stability of periodic traveling wave solutions of a generalized Kuramoto-Sivashinsky equation. In particular, we resolve the long-standing question of nonlinear modulational stability by demonstrating that spectrally stable waves are nonlinearly stable when subject to small localized (integrable) perturbations. Our analysis is based upon detailed estimates of the linearized solution operator, which are complicated by the fact that the (necessarily essential) spectrum of the associated linearization intersects the imaginary axis at the origin. We carry out a numerical Evans function study of the spectral problem and find bands of spectrally stable periodic traveling waves, in close agreement with previous numerical studies of Frisch-She-Thual, Bar-Nepomnyashchy, Chang-Demekhin-Kopelevich, and others carried out by other techniques. We also compare predictions of the associated Whitham modulation equations, which formally describe the dynamics of weak large scale perturbations of a periodic wave train, with numerical time evolution studies, demonstrating their effectiveness at a practical level. For the reader's convenience, we include in an appendix the corresponding treatment of the Swift-Hohenberg equation, a nonconservative counterpart of the generalized Kuramoto-Sivashinsky equation for which the nonlinear stability analysis is considerably simpler, together with numerical Evans function analyses extending spectral stability analyses of Mielke and Schneider. △ Less

Submitted 16 March, 2012; originally announced March 2012.

Comments: 78 pages, 11 figures