Showing 1–50 of 63 results for author: Leake, J

Compatibility of Real-Rooted Polynomials with Mixed Signs

Authors: Jonathan Leake, Nick Ryder

Abstract: We characterize compatible families of real-rooted polynomials, allowing both positive and negative leading coefficients. Our characterization naturally generalizes the same-sign characterization used by Chudnovsky and Seymour in their famous 2007 paper proving the real-rootedness of independence polynomials of claw-free graphs, thus fully settling a question left open in their paper. Our methods… ▽ More We characterize compatible families of real-rooted polynomials, allowing both positive and negative leading coefficients. Our characterization naturally generalizes the same-sign characterization used by Chudnovsky and Seymour in their famous 2007 paper proving the real-rootedness of independence polynomials of claw-free graphs, thus fully settling a question left open in their paper. Our methods are generally speaking elementary, utilizing mainly linear algebra and the established theory of interlacing polynomials, with a bit of invariant theory. △ Less

Submitted 5 August, 2024; v1 submitted 23 July, 2024; originally announced July 2024.

MSC Class: 26C10; 05A20

Capacity bounds on integral flows and the Kostant partition function

Authors: Jonathan Leake, Alejandro H. Morales

Abstract: The type $A$ Kostant partition function is an important combinatorial object with various applications: it counts integer flows on the complete directed graph, computes Hilbert series of spaces of diagonal harmonics, and can be used to compute weight and tensor product multiplicities of representations. In this paper we study asymptotics of the Kostant partition function, improving on various prev… ▽ More The type $A$ Kostant partition function is an important combinatorial object with various applications: it counts integer flows on the complete directed graph, computes Hilbert series of spaces of diagonal harmonics, and can be used to compute weight and tensor product multiplicities of representations. In this paper we study asymptotics of the Kostant partition function, improving on various previously known lower bounds and settling conjectures of O'Neill and Yip. Our methods build upon recent results and techniques of Brändén-Leake-Pak, who used Lorentzian polynomials and Gurvits' capacity method to bound the number of lattice points of transportation and flow polytopes. Finally, we also give new two-sided bounds using the Lidskii formulas from subdivisions of flow polytopes. △ Less

Submitted 24 June, 2024; v1 submitted 11 June, 2024; originally announced June 2024.

Comments: 34 pages, 3 figures, and an appendix, v2. fixed typos and added a footnote about an announced improvement to Theorem 1.1 by Balashov, Bulavenko, Dranowski, and Molybog that we learned about after v1. was posted

MSC Class: 05A16; 05A20; 52B05; 52A39 (Primary) 52B20; 32A08; 52B55; 28D20 (Secondary)

The Onset of Magnetic Reconnection in Dynamically Evolving Current Sheets

Authors: James Leake, Lars Daldorff, James Klimchuk

Abstract: We present the first results of three-dimensional (3D) numerical magnetohydrodynamic (MHD) simulations of the onset of magnetic reconnection via the tearing instability in dynamically thinning current sheets in the solar corona. In all our simulations, the onset of the non-linear tearing instability, which leads to the break-up of the thinning current sheet, does not occur until after the instabil… ▽ More We present the first results of three-dimensional (3D) numerical magnetohydrodynamic (MHD) simulations of the onset of magnetic reconnection via the tearing instability in dynamically thinning current sheets in the solar corona. In all our simulations, the onset of the non-linear tearing instability, which leads to the break-up of the thinning current sheet, does not occur until after the instability growth time becomes faster than the dynamic thinning time. Furthermore, as in previous 3D MHD simulations of static current sheets in the corona, for some parameters, the amount of magnetic shear is a fundamental switch-on parameter, which has consequences for coronal heating models. These results open up the possibility of using observable quantities of coronal current sheets to predict when they will break-up and release magnetic energy to power various energetic phenomena and/or heat the atmosphere. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Solar Eruptions Triggered by Flux Emergence Below or Near a Coronal Flux Rope

Authors: T. Török, M. G. Linton, J. E. Leake, Z. Mikić, R. Lionello, V. S. Titov, C. Downs

Abstract: Observations have shown a clear association of filament/prominence eruptions with the emergence of magnetic flux in or near filament channels. Magnetohydrodynamic (MHD) simulations have been employed to systematically study the conditions under which such eruptions occur. These simulations to date have modeled filament channels as two-dimensional (2D) flux ropes or 3D uniformly sheared arcades. He… ▽ More Observations have shown a clear association of filament/prominence eruptions with the emergence of magnetic flux in or near filament channels. Magnetohydrodynamic (MHD) simulations have been employed to systematically study the conditions under which such eruptions occur. These simulations to date have modeled filament channels as two-dimensional (2D) flux ropes or 3D uniformly sheared arcades. Here we present MHD simulations of flux emergence into a more realistic configuration consisting of a bipolar active region containing a line-tied 3D flux rope. We use the coronal flux-rope model of Titov et al. (2014) as the initial condition and drive our simulations by imposing boundary conditions extracted from a flux-emergence simulation by Leake et al. (2013). We identify three mechanisms that determine the evolution of the system: (i) reconnection displacing foot points of field lines overlying the coronal flux rope, (ii) changes of the ambient field due to the intrusion of new flux at the boundary, and (iii) interaction of the (axial) electric currents in the pre-existing and newly emerging flux systems. The relative contributions and effects of these mechanisms depend on the properties of the pre-existing and emerging flux systems. Here we focus on the location and orientation of the emerging flux relative to the coronal flux rope. Varying these parameters, we investigate under which conditions an eruption of the latter is triggered. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 18 pages, 8 figures, accepted for publication by The Astrophysical Journal

From Trees to Polynomials and Back Again: New Capacity Bounds with Applications to TSP

Authors: Leonid Gurvits, Nathan Klein, Jonathan Leake

Abstract: We give simply exponential lower bounds on the probabilities of a given strongly Rayleigh distribution, depending only on its expectation. This resolves a weak version of a problem left open by Karlin-Klein-Oveis Gharan in their recent breakthrough work on metric TSP, and this resolution leads to a minor improvement of their approximation factor for metric TSP. Our results also allow for a more st… ▽ More We give simply exponential lower bounds on the probabilities of a given strongly Rayleigh distribution, depending only on its expectation. This resolves a weak version of a problem left open by Karlin-Klein-Oveis Gharan in their recent breakthrough work on metric TSP, and this resolution leads to a minor improvement of their approximation factor for metric TSP. Our results also allow for a more streamlined analysis of the algorithm. To achieve these new bounds, we build upon the work of Gurvits-Leake on the use of the productization technique for bounding the capacity of a real stable polynomial. This technique allows one to reduce certain inequalities for real stable polynomials to products of affine linear forms, which have an underlying matrix structure. In this paper, we push this technique further by characterizing the worst-case polynomials via bipartitioned forests. This rigid combinatorial structure yields a clean induction argument, which implies our stronger bounds. In general, we believe the results of this paper will lead to further improvement and simplification of the analysis of various combinatorial and probabilistic bounds and algorithms. △ Less

Submitted 9 May, 2024; v1 submitted 15 November, 2023; originally announced November 2023.

Comments: ICALP 2024

Simulating the photospheric to coronal plasma using magnetohydrodyanamic characteristics I: data-driven boundary conditions

Authors: Lucas A. Tarr, N. Dylan Kee, Mark G. Linton, Peter W. Schuck, James E. Leake

Abstract: We develop a general description of how information propagates through a magnetohydrodynamic (MHD) system based on the method of characteristics and use that to formulate numerical boundary conditions (BCs) that are intrinsically consistent with the MHD equations. Our formulation includes two major advances for simulations of the Sun. First, we derive data-driven BCs that optimally match the state… ▽ More We develop a general description of how information propagates through a magnetohydrodynamic (MHD) system based on the method of characteristics and use that to formulate numerical boundary conditions (BCs) that are intrinsically consistent with the MHD equations. Our formulation includes two major advances for simulations of the Sun. First, we derive data-driven BCs that optimally match the state of the plasma inferred from a time series of observations of a boundary (e.g., the solar photosphere). Second, our method directly handles random noise and systematic bias in the observations, and finds a solution for the boundary evolution that is strictly consistent with MHD and maximally consistent with the observations. We validate the method against a Ground Truth (GT) simulation of an expanding spheromak. The data-driven simulation can reproduce the GT simulation above the photosphere with high fidelity when driven at high cadence. Errors progressively increase for lower driving cadence until a threshold cadence is reached and the driven simulation can no longer accurately reproduce the GT simulation. However, our characteristic formulation of the BCs still requires adherence of the boundary evolution to the MHD equations even when the driven solution departs from the true solution in the driving layer. That increasing departure clearly indicates when additional information at the boundary is needed to fully specify the correct evolution of the system. The method functions even when no information about the evolution of some variables on the lower boundary is available, albeit with a further decrease in fidelity. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: Submitted for publication to the Astrophysical Journal; 77 pages, 12 figures

The Thickness of Electric Current Sheets and Implications for Coronal Heating

Authors: James A. Klimchuk, James E. Leake, Lars K. S. Daldorff, Craig D. Johnston

Abstract: The thickness of current sheets is extremely important, especially as it relates to the onset of fast magnetic reconnection. Onset determines how much magnetic free energy can build up in a field before it is explosively released. This has implications for many phenomena on the Sun and throughout the universe, including the heating of the solar corona. Significant effort has been devoted to the qu… ▽ More The thickness of current sheets is extremely important, especially as it relates to the onset of fast magnetic reconnection. Onset determines how much magnetic free energy can build up in a field before it is explosively released. This has implications for many phenomena on the Sun and throughout the universe, including the heating of the solar corona. Significant effort has been devoted to the question of whether equilibrium current sheets in realistic geometries have finite or zero thickness. Using a simple force balance analysis, we show why current sheets without a guide field (2D) and with a guide field that is invariant in the guide field direction (2.5D) cannot be in equilibrium if they have both finite thickness and finite length. We then estimate the conditions under which the tension of a curved line-tied guide field can facilitate equilibrium in 3D sheets that are finite in all dimensions. Finally, we argue that some quasi-statically evolving current sheets undergoing slow stressing (e.g., when the coronal magnetic field is subjected to photospheric boundary driving) may reach a critical shear, at which point they lose equilibrium, spontaneously collapse, and reconnect. The critical shear is generally consistent with the heating requirements of solar active regions. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 15 pages, 5 figures, accepted by Frontiers in Physics

Lorentzian polynomials on cones

Authors: Petter Brändén, Jonathan Leake

Abstract: Inspired by the theory of hyperbolic polynomials and Hodge theory, we develop the theory of Lorentzian polynomials on cones. This notion captures the Hodge-Riemann relations of degree zero and one. Motivated by fundamental properties of volume polynomials of Chow rings of simplicial fans, we define a class of multivariate polynomials which we call hereditary polynomials. We give a complete and eas… ▽ More Inspired by the theory of hyperbolic polynomials and Hodge theory, we develop the theory of Lorentzian polynomials on cones. This notion captures the Hodge-Riemann relations of degree zero and one. Motivated by fundamental properties of volume polynomials of Chow rings of simplicial fans, we define a class of multivariate polynomials which we call hereditary polynomials. We give a complete and easily checkable characterization of hereditary Lorentzian polynomials. This characterization is used to give elementary and simple proofs of the Heron-Rota-Welsh conjecture for the characteristic polynomial of a matroid, and the Alexandrov-Fenchel inequalities for convex bodies. We then characterize Chow rings of simplicial fans which satisfy the Hodge-Riemann relations of degree zero and one, and we prove that this property only depends on the support of the fan. Several different characterizations of Lorentzian polynomials on cones are provided. △ Less

Submitted 25 April, 2023; originally announced April 2023.

The Multiview Observatory for Solar Terrestrial Science (MOST)

Authors: N. Gopalswamy, S. Christe, S. F. Fung, Q. Gong, J. R. Gruesbeck, L. K. Jian, S. G. Kanekal, C. Kay, T. A. Kucera, J. E. Leake, L. Li, P. Makela, P. Nikulla, N. L. Reginald, A. Shih, S. K. Tadikonda, N. Viall, L. B. Wilson III, S. Yashiro, L. Golub, E. DeLuca, K. Reeves, A. C. Sterling, A. R. Winebarger, C. DeForest , et al. (32 additional authors not shown)

Abstract: We report on a study of the Multiview Observatory for Solar Terrestrial Science (MOST) mission that will provide comprehensive imagery and time series data needed to understand the magnetic connection between the solar interior and the solar atmosphere/inner heliosphere. MOST will build upon the successes of SOHO and STEREO missions with new views of the Sun and enhanced instrument capabilities. T… ▽ More We report on a study of the Multiview Observatory for Solar Terrestrial Science (MOST) mission that will provide comprehensive imagery and time series data needed to understand the magnetic connection between the solar interior and the solar atmosphere/inner heliosphere. MOST will build upon the successes of SOHO and STEREO missions with new views of the Sun and enhanced instrument capabilities. This article is based on a study conducted at NASA Goddard Space Flight Center that determined the required instrument refinement, spacecraft accommodation, launch configuration, and flight dynamics for mission success. MOST is envisioned as the next generation great observatory positioned to obtain three-dimensional information of large-scale heliospheric structures such as coronal mass ejections, stream interaction regions, and the solar wind itself. The MOST mission consists of 2 pairs of spacecraft located in the vicinity of Sun-Earth Lagrange points L4 (MOST1, MOST3) and L5 (MOST2 and MOST4). The spacecraft stationed at L4 (MOST1) and L5 (MOST2) will each carry seven remote-sensing and three in-situ instrument suites, including a novel radio package known as the Faraday Effect Tracker of Coronal and Heliospheric structures (FETCH). MOST3 and MOST4 will carry only the FETCH instruments and are positioned at variable locations along the Earth orbit up to 20° ahead of L4 and 20° behind L5, respectively. FETCH will have polarized radio transmitters and receivers on all four spacecraft to measure the magnetic content of solar wind structures propagating from the Sun to Earth using the Faraday rotation technique. The MOST mission will be able to sample the magnetized plasma throughout the Sun-Earth connected space during the mission lifetime over a solar cycle. △ Less

Submitted 10 December, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Comments: 42 pages, 19 figures, 8 tables, to appear in J. Atmospheric and Solar Terrestrial Physics

Deterministic Approximation Algorithms for Volumes of Spectrahedra

Authors: Mahmut Levent Doğan, Jonathan Leake, Mohan Ravichandran

Abstract: We give a method for computing asymptotic formulas and approximations for the volumes of spectrahedra, based on the maximum-entropy principle from statistical physics. The method gives an approximate volume formula based on a single convex optimization problem of minimizing $-\log \det P$ over the spectrahedron. Spectrahedra can be described as affine slices of the convex cone of positive semi-def… ▽ More We give a method for computing asymptotic formulas and approximations for the volumes of spectrahedra, based on the maximum-entropy principle from statistical physics. The method gives an approximate volume formula based on a single convex optimization problem of minimizing $-\log \det P$ over the spectrahedron. Spectrahedra can be described as affine slices of the convex cone of positive semi-definite (PSD) matrices, and the method yields efficient deterministic approximation algorithms and asymptotic formulas whenever the number of affine constraints is sufficiently dominated by the dimension of the PSD cone. Our approach is inspired by the work of Barvinok and Hartigan who used an analogous framework for approximately computing volumes of polytopes. Spectrahedra, however, possess a remarkable feature not shared by polytopes, a new fact that we also prove: central sections of the set of density matrices (the quantum version of the simplex) all have asymptotically the same volume. This allows for very general approximation algorithms, which apply to large classes of naturally occurring spectrahedra. We give two main applications of this method. First, we apply this method to what we call the "multi-way Birkhoff spectrahedron" and obtain an explicit asymptotic formula for its volume. This spectrahedron is the set of quantum states with maximal entanglement (i.e., the quantum states having univariant quantum marginals equal to the identity matrix) and is the quantum analog of the multi-way Birkhoff polytope. Second, we apply this method to explicitly compute the asymptotic volume of central sections of the set of density matrices. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Variable-wavelength quick scanning nano-focused X-ray microscopy for in situ strain and tilt mapping

Authors: Marie-ingrid Richard, Thomas W Cornelius, Florian Lauraux, Jean-Baptiste Molin, Christoph Kirchlechner, Steven J Leake, Jérôme Carnis, Tobias U Schülli, Ludovic Thilly, Olivier Thomas

Abstract: Compression of micro-pillars is followed in situ by a quick nano-focused X-ray scanning microscopy technique combined with three-dimensional reciprocal space mapping. Compared to other attempts using 2 X-ray nanobeams, it avoids any motion or vibration that would lead to a destruction of the sample. The technique consists of scanning both the energy of the incident nano-focused X-ray beam and the… ▽ More Compression of micro-pillars is followed in situ by a quick nano-focused X-ray scanning microscopy technique combined with three-dimensional reciprocal space mapping. Compared to other attempts using 2 X-ray nanobeams, it avoids any motion or vibration that would lead to a destruction of the sample. The technique consists of scanning both the energy of the incident nano-focused X-ray beam and the in-plane translations of the focusing optics along the X-ray beam. Here, we demonstrate the approach by imaging the strain and lattice orientation of Si micro-pillars and their pedestals during in situ compression. Varying the energy of the incident beam instead of rocking the sample and mapping the focusing optics instead of moving the sample supplies a vibration-free measurement of the reciprocal space maps without removal of the mechanical load. The maps of strain and lattice orientation are in good agreement with the ones recorded by ordinary rocking-curve scans. Variable-wavelength quick scanning X-ray microscopy opens the route for in situ strain and tilt mapping towards more diverse and complex materials environments, especially where sample manipulation is difficult. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Journal ref: Small, Wiley-VCH Verlag, 2020

The Role of Reconnection in the Onset of Solar Eruptions

Authors: James Leake, Mark Linton, Spiro Antiochos

Abstract: Solar eruptive events such as coronal mass ejections and eruptive flares are frequently associated with the emergence of magnetic flux from the convection zone into the corona. We use three dimensional magnetohydrodynamic numerical simulations to study the interaction of coronal magnetic fields with emerging flux and determine the conditions that lead to eruptive activity. A simple parameter study… ▽ More Solar eruptive events such as coronal mass ejections and eruptive flares are frequently associated with the emergence of magnetic flux from the convection zone into the corona. We use three dimensional magnetohydrodynamic numerical simulations to study the interaction of coronal magnetic fields with emerging flux and determine the conditions that lead to eruptive activity. A simple parameter study is performed, varying the relative angle between emerging magnetic flux and a pre-existing coronal dipole field. We find that in all cases, the emergence results in a sheared magnetic arcade that transitions to a twisted coronal flux rope via low-lying magnetic reconnection. This structure, however, is constrained by its own outer field, and so is non-eruptive in the absence of reconnection with the overlying coronal field. The amount of this overlying reconnection is determined by the relative angle between the emerged and pre-existing fields. The reconnection between emerging and pre-existing fields is necessary to generate sufficient expansion of the emerging structure so that flare-like reconnection below the coronal flux rope becomes strong enough to trigger its release. Our results imply that the relative angle is the key parameter in determining whether the resultant active regions exhibit eruptive behavior, and is thus a potentially useful candidate for predicting eruptions in newly emerging active regions. More generally, our results demonstrate that the detailed interaction between the convection zone/photosphere and the corona must be calculated self-consistently in order to model solar eruptions accurately. △ Less

Submitted 23 May, 2022; originally announced May 2022.

Imaging the Breathing of a Platinum Nanoparticle in Electrochemical Environment

Authors: Clément Atlan, Corentin Chatelier, Maxime Dupraz, Isaac Martens, Arnaud Viola, Ni Li, Lu Gao, Steven J. Leake, Tobias U. Schülli, Joël Eymery, Frédéric Maillard, Marie-Ingrid Richard

Abstract: Surface strain is widely used in gas phase catalysis and electrocatalysis to control the binding energies of adsorbates on metal surfaces. However, $in$ $situ$ or $operando$ strain measurements are experimentally challenging, especially on nanomaterials. Here, we take advantage of the 4$^{th}$ generation Extremely Brilliant Source at the European Synchrotron Radiation Facility (ESRF-EBS, Grenoble,… ▽ More Surface strain is widely used in gas phase catalysis and electrocatalysis to control the binding energies of adsorbates on metal surfaces. However, $in$ $situ$ or $operando$ strain measurements are experimentally challenging, especially on nanomaterials. Here, we take advantage of the 4$^{th}$ generation Extremely Brilliant Source at the European Synchrotron Radiation Facility (ESRF-EBS, Grenoble, France) to quantify the distribution of strain inside a Pt nanoparticle, and to determine its morphology in an electrochemical environment. Our results show for the first time evidence of heterogeneous and potential-dependent strain distribution between highly-coordinated ({100} and {111} facets) and under-coordinated atoms (edges and corners) as well as evidence of strain propagation from the surface to the bulk of the nanoparticle. These results provide dynamic structural insights to better simulate and design efficient nanocatalysts for energy storage and conversion applications. △ Less

Submitted 2 February, 2023; v1 submitted 14 March, 2022; originally announced March 2022.

Impact of 3D Structure on Magnetic Reconnection

Authors: Lars K. S. Daldorff, James E. Leake, James A. Klimchuk

Abstract: Results from 2.5D and 3D studies of the onset and development of the tearing instability are presented, using high fidelity resistive MHD simulations. A limited parameter study of the strength of the reconnecting field (or shear angle) was performed. An initially simple 1D equilibrium was used, consisting of a modified force-free current sheet, with periodic boundary conditions in all directions.… ▽ More Results from 2.5D and 3D studies of the onset and development of the tearing instability are presented, using high fidelity resistive MHD simulations. A limited parameter study of the strength of the reconnecting field (or shear angle) was performed. An initially simple 1D equilibrium was used, consisting of a modified force-free current sheet, with periodic boundary conditions in all directions. In all cases, the linear and non-linear evolution led to a primary current sheet between two large flux ropes. The global reconnection rate during this later stage was analyzed in all simulations. It was found that in 2.5D the primary current sheet fragmented due to plasmoids, and as expected, the global reconnection rate, calculated using multiple methods, increases with the strength of the reconnecting field due to a stronger Alfvén speed. In 3D, the presence of interacting oblique modes of the tearing instability complicates the simple 2.5D picture, entangling the magnetic field of the inflow and introducing a negative effect on the reconnection rate. The two competing effects of stronger Alfvén speed and entangling, which both increase with the strength of the reconnecting field, resulted in a decrease in the reconnection rate with increasing reconnecting field. For all simulations, the 3D rates were less than in 2.5D, but suggest that as one goes to weak reconnecting field (or strong guide field), the system becomes more 2.5D like and the 2.5D and 3D rates converge. These results have relevance to situations like nano-flare heating and flare current sheets in the corona. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 15 pages, 10 figures

Lorentzian polynomials on cones and the Heron-Rota-Welsh conjecture

Authors: Petter Brändén, Jonathan Leake

Abstract: We give a short proof of the log-concavity of the coefficients of the reduced characteristic polynomial of a matroid. The proof uses an extension of the theory of Lorentzian polynomials to convex cones, and reproves the Hodge-Riemann relations of degree one for the Chow ring of a matroid. We give a short proof of the log-concavity of the coefficients of the reduced characteristic polynomial of a matroid. The proof uses an extension of the theory of Lorentzian polynomials to convex cones, and reproves the Hodge-Riemann relations of degree one for the Chow ring of a matroid. △ Less

Submitted 11 October, 2021; v1 submitted 1 October, 2021; originally announced October 2021.

Optimization and Sampling Under Continuous Symmetry: Examples and Lie Theory

Authors: Jonathan Leake, Nisheeth K. Vishnoi

Abstract: In the last few years, the notion of symmetry has provided a powerful and essential lens to view several optimization or sampling problems that arise in areas such as theoretical computer science, statistics, machine learning, quantum inference, and privacy. Here, we present two examples of nonconvex problems in optimization and sampling where continuous symmetries play -- implicitly or explicitly… ▽ More In the last few years, the notion of symmetry has provided a powerful and essential lens to view several optimization or sampling problems that arise in areas such as theoretical computer science, statistics, machine learning, quantum inference, and privacy. Here, we present two examples of nonconvex problems in optimization and sampling where continuous symmetries play -- implicitly or explicitly -- a key role in the development of efficient algorithms. These examples rely on deep and hidden connections between nonconvex symmetric manifolds and convex polytopes, and are heavily generalizable. To formulate and understand these generalizations, we then present an introduction to Lie theory -- an indispensable mathematical toolkit for capturing and working with continuous symmetries. We first present the basics of Lie groups, Lie algebras, and the adjoint actions associated with them, and we also mention the classification theorem for Lie algebras. Subsequently, we present Kostant's convexity theorem and show how it allows us to reduce linear optimization problems over orbits of Lie groups to linear optimization problems over polytopes. Finally, we present the Harish-Chandra and the Harish-Chandra--Itzykson--Zuber (HCIZ) formulas, which convert partition functions (integrals) over Lie groups into sums over the corresponding (discrete) Weyl groups, enabling efficient sampling algorithms. △ Less

Submitted 2 September, 2021; originally announced September 2021.

Comments: This article is to supplement the talks by the authors at the Bootcamp in the semester on Geometric Methods for Optimization and Sampling at the Simons Institute for the Theory of Computing

Bragg Coherent Imaging of nanoprecipitates: role of superstructure reflections

Authors: Maxime Dupraz, Steven J. Leake, Marie-Ingrid Richard

Abstract: Coherent precipitation of ordered phases is responsible for providing exceptional high temperature mechanical properties in a wide range of compositionally complex alloys (CCAs). Ordered phases are also essential to enhance the magnetic or catalytic properties of alloyed nanoparticles. The present work aims at demonstrating the relevance of Bragg coherent diffraction imaging (BCDI) to study bulk a… ▽ More Coherent precipitation of ordered phases is responsible for providing exceptional high temperature mechanical properties in a wide range of compositionally complex alloys (CCAs). Ordered phases are also essential to enhance the magnetic or catalytic properties of alloyed nanoparticles. The present work aims at demonstrating the relevance of Bragg coherent diffraction imaging (BCDI) to study bulk and thin film samples or isolated nanoparticles containing coherent nanoprecipitates / ordered phases. Crystals of a few tens of nanometres are modelled with realistic interatomic potentials and relaxed after introduction of coherent ordered nanoprecipitates. Diffraction patterns from fundamental and superstructure reflections are calculated in the kinematic approximation and used as input to retrieve the strain fields using algorithmic inversion. We first tackle the case of single nanoprecipitates and show that the strain field distribution from the ordered phase is retrieved very accurately. Then, we investigate the influence of the order parameter S on the strain field retrieved from the superstructure reflections and evidence that a very accurate strain distribution can be retrieved for partially ordered phases with large and inhomogeneous strains. In a subsequent section, we evaluate the relevance of BCDI for the study of systems containing many precipitates and demonstrate that the technique is relevant for such systems. Finally, we discuss the experimental feasibility of using BCDI to image ordered phases, in the light of the new possibilities offered by the 4 th generation synchrotron sources. △ Less

Submitted 24 February, 2021; originally announced February 2021.

Comments: 33 pages, 12 figures

Journal ref: September 2020 Journal of Applied Crystallography 53(5):1353-1369

The Rise and Emergence of Untwisted Toroidal Flux Ropes on the Sun

Authors: Kalman J. Knizhnik, James E. Leake, Mark G. Linton, Sally Dacie

Abstract: Magnetic flux ropes (MFRs) rising buoyantly through the Sun's convection zone are thought to be subject to viscous forces preventing them from rising coherently. Numerous studies have suggested that MFRs require a minimum twist in order to remain coherent during their rise. Furthermore, even MFRs that get to the photosphere may be unable to successfully emerge into the corona unless they are at le… ▽ More Magnetic flux ropes (MFRs) rising buoyantly through the Sun's convection zone are thought to be subject to viscous forces preventing them from rising coherently. Numerous studies have suggested that MFRs require a minimum twist in order to remain coherent during their rise. Furthermore, even MFRs that get to the photosphere may be unable to successfully emerge into the corona unless they are at least moderately twisted, since the magnetic pressure gradient needs to overcome the weight of the photospheric plasma. To date, however, no lower limit has been placed on the critical minimum twist required for an MFR to rise coherently through the convection zone or emerge through the photosphere. In this paper, we simulate an untwisted toroidal MFR which is able to rise from the convection zone and emerge through the photosphere as an active region that resembles those observed on the Sun. We show that untwisted MFRs can remain coherent during their rise and then pile-up near the photosphere, triggering the undular instability, allowing the MFR to emerge through the photosphere. We propose that the toroidal geometry of our MFR is critical for its coherent rise. Upon emerging, a pair of lobes rises into the corona which interact and reconnect, resulting in a localized high speed jet. The resulting photospheric magnetogram displays the characteristic salt-and-pepper structure often seen in observations. Our major result is that MFRs need not be twisted to rise coherently through the convection zone and emerge through the photosphere. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: Accepted to ApJ

Sampling Matrices from Harish-Chandra-Itzykson-Zuber Densities with Applications to Quantum Inference and Differential Privacy

Authors: Jonathan Leake, Colin S. McSwiggen, Nisheeth K. Vishnoi

Abstract: Given two $n \times n$ Hermitian matrices $Y$ and $Λ$, the Harish-Chandra-Itzykson-Zuber (HCIZ) distribution on the unitary group $\text{U}(n)$ is $e^{\text{tr}(UΛU^*Y)}dμ(U)$, where $μ$ is the Haar measure on $\text{U}(n)$. The density $e^{\text{tr}(UΛU^*Y)}$ is known as the HCIZ density. Random unitary matrices distributed according to the HCIZ density are important in various settings in physic… ▽ More Given two $n \times n$ Hermitian matrices $Y$ and $Λ$, the Harish-Chandra-Itzykson-Zuber (HCIZ) distribution on the unitary group $\text{U}(n)$ is $e^{\text{tr}(UΛU^*Y)}dμ(U)$, where $μ$ is the Haar measure on $\text{U}(n)$. The density $e^{\text{tr}(UΛU^*Y)}$ is known as the HCIZ density. Random unitary matrices distributed according to the HCIZ density are important in various settings in physics and random matrix theory. However, the basic question of efficient sampling from the HCIZ distribution has remained open. We present two efficient algorithms to sample matrices from distributions that are close to the HCIZ distribution. The first algorithm outputs samples that are $ξ$-close in total variation distance and requires polynomially many arithmetic operations in $\log 1/ξ$ and the number of bits needed to encode $Y$ and $Λ$. The second algorithm comes with a stronger guarantee that the samples are $ξ$-close in infinity divergence, but the number of arithmetic operations depends polynomially on $1/ξ$, the number of bits needed to encode $Y$ and $Λ$, and the differences of the largest and the smallest eigenvalues of $Y$ and $Λ$. HCIZ densities can also be viewed as exponential densities on $\text{U}(n)$-orbits, and these densities have been studied in statistics, machine learning, and theoretical computer science. Thus our results have the following applications: 1) an efficient algorithm to sample from complex versions of matrix Langevin distributions studied in statistics, 2) an efficient algorithm to sample from continuous max-entropy distributions on unitary orbits, which implies an efficient algorithm to sample a pure quantum state from the entropy-maximizing ensemble representing a given density matrix, and 3) an efficient algorithm for differentially private rank-$k$ approximation, with improved utility bounds for $k>1$. △ Less

Submitted 6 April, 2021; v1 submitted 10 November, 2020; originally announced November 2020.

Comments: Full version of a paper appearing in STOC 2021

On the Computability of Continuous Maximum Entropy Distributions: Adjoint Orbits of Lie Groups

Authors: Jonathan Leake, Nisheeth K. Vishnoi

Abstract: Given a point $A$ in the convex hull of a given adjoint orbit $\mathcal{O}(F)$ of a compact Lie group $G$, we give a polynomial time algorithm to compute the probability density supported on $\mathcal{O}(F)$ whose expectation is $A$ and that minimizes the Kullback-Leibler divergence to the $G$-invariant measure on $\mathcal{O}(F)$. This significantly extends the recent work of the authors (STOC 20… ▽ More Given a point $A$ in the convex hull of a given adjoint orbit $\mathcal{O}(F)$ of a compact Lie group $G$, we give a polynomial time algorithm to compute the probability density supported on $\mathcal{O}(F)$ whose expectation is $A$ and that minimizes the Kullback-Leibler divergence to the $G$-invariant measure on $\mathcal{O}(F)$. This significantly extends the recent work of the authors (STOC 2020) who presented such a result for the manifold of rank $k$-projections which is a specific adjoint orbit of the unitary group $\mathrm{U}(n)$. Our result relies on the ellipsoid method-based framework proposed in prior work; however, to apply it to the general setting of compact Lie groups, we need tools from Lie theory. For instance, properties of the adjoint representation are used to find the defining equalities of the minimal affine space containing the convex hull of $\mathcal{O}(F)$, and to establish a bound on the optimal dual solution. Also, the Harish-Chandra integral formula is used to obtain an evaluation oracle for the dual objective function. While the Harish-Chandra integral formula allows us to write certain integrals over the adjoint orbit of a Lie group as a sum of a small number of determinants, it is only defined for elements of a chosen Cartan subalgebra of the Lie algebra $\mathfrak{g}$ of $G.$ We show how it can be applied to our setting with the help of Kostant's convexity theorem. Further, the convex hull of an adjoint orbit is a type of orbitope, and the orbitopes studied in this paper are known to be spectrahedral. Thus our main result can be viewed as extending the maximum entropy framework to a class of spectrahedra. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Decoding the Pre-Eruptive Magnetic Field Configurations of Coronal Mass Ejections

Authors: S. Patsourakos, A. Vourlidas, T. Török, B. Kliem, S. K. Antiochos, V. Archontis, G. Aulanier, X. Cheng, G. Chintzoglou, M. K. Georgoulis, L. M. Green, J. E. Leake, R. Moore, A. Nindos, P. Syntelis, S. L. Yardley, V. Yurchyshyn, J. Zhang

Abstract: A clear understanding of the nature of the pre-eruptive magnetic field configurations of Coronal Mass Ejections (CMEs) is required for understanding and eventually predicting solar eruptions. Only two, but seemingly disparate, magnetic configurations are considered viable; namely, sheared magnetic arcades (SMA) and magnetic flux ropes (MFR). They can form via three physical mechanisms (flux emerge… ▽ More A clear understanding of the nature of the pre-eruptive magnetic field configurations of Coronal Mass Ejections (CMEs) is required for understanding and eventually predicting solar eruptions. Only two, but seemingly disparate, magnetic configurations are considered viable; namely, sheared magnetic arcades (SMA) and magnetic flux ropes (MFR). They can form via three physical mechanisms (flux emergence, flux cancellation, helicity condensation) . Whether the CME culprit is an SMA or an MFR, however, has been strongly debated for thirty years. We formed an International Space Science Institute (ISSI) team to address and resolve this issue and report the outcome here. We review the status of the field across modeling and observations, identify the open and closed issues, compile lists of SMA and MFR observables to be tested against observations and outline research activities to close the gaps in our current understanding. We propose that the combination of multi-viewpoint multi-thermal coronal observations and multi-height vector magnetic field measurements is the optimal approach for resolving the issue conclusively. We demonstrate the approach using MHD simulations and synthetic coronal images. Our key conclusion is that the differentiation of pre-eruptive configurations in terms of SMAs and MFRs seems artificial. Both observations and modeling can be made consistent if the pre-eruptive configuration exists in a hybrid state that is continuously evolving from an SMA to an MFR. Thus, the 'dominant' nature of a given configuration will largely depend on its evolutionary stage (SMA-like early-on, MFR-like near the eruption). △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: Space Science Reviews, accepted for publication

Lower bounds for contingency tables via Lorentzian polynomials

Authors: Petter Brändén, Jonathan Leake, Igor Pak

Abstract: We present a new lower bound on the number of contingency tables, improving upon and extending previous lower bounds by Barvinok and Gurvits. As an application, we obtain new lower bounds on the volumes of flow and transportation polytopes. Our proofs are based on recent results on Lorentzian polynomials. We present a new lower bound on the number of contingency tables, improving upon and extending previous lower bounds by Barvinok and Gurvits. As an application, we obtain new lower bounds on the volumes of flow and transportation polytopes. Our proofs are based on recent results on Lorentzian polynomials. △ Less

Submitted 13 February, 2021; v1 submitted 13 August, 2020; originally announced August 2020.

Comments: 28 pages

Capacity Lower Bounds via Productization

Authors: Leonid Gurvits, Jonathan Leake

Abstract: We give a sharp lower bound on the capacity of a real stable polynomial, depending only on the value of its gradient at $x = 1$. This result implies a sharp improvement to a similar inequality proved by Linial-Samorodnitsky-Wigderson in 2000, which was crucial to the analysis of their permanent approximation algorithm. Such inequalities have played an important role in the recent work on operator… ▽ More We give a sharp lower bound on the capacity of a real stable polynomial, depending only on the value of its gradient at $x = 1$. This result implies a sharp improvement to a similar inequality proved by Linial-Samorodnitsky-Wigderson in 2000, which was crucial to the analysis of their permanent approximation algorithm. Such inequalities have played an important role in the recent work on operator scaling and its generalizations and applications, and in fact we use our bound to construct a new scaling algorithm for real stable polynomials. In addition, we give a strong improvement on previous lower bounds of the capacity of a non-homogeneous real stable polynomial, depending only on the value of its gradient at $x = 1$. Crucially, this new bound is independent of the degree of the polynomial, and has singly exponential dependence on the number of variables. This compares favorably to the bounds used recently in the fantastic work of Karlin-Klein-Oveis Gharan to give an improved approximation factor for metric TSP, where this dependence is doubly exponential. Such bounds were conjectured to exist by the authors, and thus our new bound should imply further improvement to the approximation factor for metric TSP. The new technique we develop to prove this bound is productization, which says that any real stable polynomial can be approximated at any point in the positive orthant by a product of linear forms. Beyond the results of this paper, our main hope is that this new technique will allow us to avoid "frightening technicalities", in the words of Laurent and Schrijver, that often accompany combinatorial lower bounds. △ Less

Submitted 21 September, 2022; v1 submitted 16 July, 2020; originally announced July 2020.

Comments: STOC 2021; this version has a new and improved bound which should imply improvements to the metric TSP results of Karlin-Klein-Oveis Gharan

On the computability of continuous maximum entropy distributions with applications

Authors: Jonathan Leake, Nisheeth K. Vishnoi

Abstract: We initiate a study of the following problem: Given a continuous domain $Ω$ along with its convex hull $\mathcal{K}$, a point $A \in \mathcal{K}$ and a prior measure $μ$ on $Ω$, find the probability density over $Ω$ whose marginal is $A$ and that minimizes the KL-divergence to $μ$. This framework gives rise to several extremal distributions that arise in mathematics, quantum mechanics, statistics,… ▽ More We initiate a study of the following problem: Given a continuous domain $Ω$ along with its convex hull $\mathcal{K}$, a point $A \in \mathcal{K}$ and a prior measure $μ$ on $Ω$, find the probability density over $Ω$ whose marginal is $A$ and that minimizes the KL-divergence to $μ$. This framework gives rise to several extremal distributions that arise in mathematics, quantum mechanics, statistics, and theoretical computer science. Our technical contributions include a polynomial bound on the norm of the optimizer of the dual problem that holds in a very general setting and relies on a "balance" property of the measure $μ$ on $Ω$, and exact algorithms for evaluating the dual and its gradient for several interesting settings of $Ω$ and $μ$. Together, along with the ellipsoid method, these results imply polynomial-time algorithms to compute such KL-divergence minimizing distributions in several cases. Applications of our results include: 1) an optimization characterization of the Goemans-Williamson measure that is used to round a positive semidefinite matrix to a vector, 2) the computability of the entropic barrier for polytopes studied by Bubeck and Eldan, and 3) a polynomial-time algorithm to compute the barycentric quantum entropy of a density matrix that was proposed as an alternative to von Neumann entropy in the 1970s: this corresponds to the case when $Ω$ is the set of rank one projections matrices and $μ$ corresponds to the Haar measure on the unit sphere. Our techniques generalize to the setting of Hermitian rank $k$ projections using the Harish-Chandra-Itzykson-Zuber formula, and are applicable even beyond, to adjoint orbits of compact Lie groups. △ Less

Submitted 15 April, 2020; originally announced April 2020.

Comments: 50 pages, STOC 2020

The Onset of 3D Magnetic Reconnection and Heating in the Solar Corona

Authors: James Leake, Lars Daldorff, James Klimchuk

Abstract: Magnetic reconnection, a fundamentally important process in many aspects of astrophysics, is believed to be initiated by the tearing instability of an electric current sheet, a region where magnetic field abruptly changes direction and electric currents build up. Recent studies have suggested that the amount of magnetic shear in these structures is a critical parameter for the switch-on nature of… ▽ More Magnetic reconnection, a fundamentally important process in many aspects of astrophysics, is believed to be initiated by the tearing instability of an electric current sheet, a region where magnetic field abruptly changes direction and electric currents build up. Recent studies have suggested that the amount of magnetic shear in these structures is a critical parameter for the switch-on nature of magnetic reconnection in the solar atmosphere, at fluid spatial scales much larger than kinetic scales. We present results of simulations of reconnection in 3D current sheets with conditions appropriate to the solar corona. Using high-fidelity simulations, we follow the evolution of the linear and non-linear 3D tearing instability, leading to reconnection. We find that, depending on the parameter space, magnetic shear can play a vital role in the onset of significant energy release and heating via non-linear tearing. Two regimes in our study exist, dependent on whether the current sheet is longer or shorter than the wavelength of the fastest growing parallel mode (in the corresponding infinite system), thus determining whether sub-harmonics are present in the actual system. In one regime, where the fastest growing parallel mode has sub-harmonics, the non-linear interaction of these sub-harmonics and the coalescence of 3D plasmoids dominates the non-linear evolution, with magnetic shear playing only a weak role in the amount of energy released. In the second regime, where the fastest growing parallel mode has no-sub-harmonics, then only strongly sheared current sheets, where oblique mode are strong enough to compete with the dominant parallel mode, show any significant energy release. We expect both regimes to exist on the Sun, and so our results have important consequences for the the question of reconnection onset in different solar physics applications. △ Less

Submitted 25 February, 2020; v1 submitted 9 January, 2020; originally announced January 2020.

Comments: 11 pages9 figures

Evidence of Charge Density Wave transverse pinning by x-ray micro-diffraction

Authors: E. Bellec, I. Gonzalez-Vallejo, V. L. R. Jacques, A. A. Sinchenko, A. P. Orlov, P. Monceau, S. J. Leake, D. Le Bolloc'h

Abstract: Incommensurate charge density waves (CDW) have the extraordinary ability to display non-Ohmic behavior when submitted to an external field. The mechanism leading to this non trivial dynamics is still not well understood, although recent experimental studies tend to prove that it is due to solitonic transport. Solitons could come from the relaxation of the strained CDW within an elastic-to-plastic… ▽ More Incommensurate charge density waves (CDW) have the extraordinary ability to display non-Ohmic behavior when submitted to an external field. The mechanism leading to this non trivial dynamics is still not well understood, although recent experimental studies tend to prove that it is due to solitonic transport. Solitons could come from the relaxation of the strained CDW within an elastic-to-plastic transition. However, the nucleation process and the transport of these charged topological objects have never been observed at the local scale until now. In this letter, we use in-situ scanning x-ray micro-diffraction with micrometer resolution of a NbSe$_3$ sample designed to have sliding and non-sliding areas. Direct imaging of the charge density wave deformation is obtained using an analytical approach based on the phase gradient to disentangle the transverse from the longitudinal components over a large surface of a hundred microns size. We show that the CDW dissociates itself from the host lattice in the sliding regime and displays a large transverse deformation, ten times larger than the longitudinal one and strongly dependent on the amplitude and the direction of the applied currents. This deformation continuously extends across the macroscopic sample dimensions, over a distance 10 000 times greater than the CDW wavelength despite the presence of strong defects while remaining strongly pinned by the lateral surfaces. This 2D quantitative study highlights the prominent role of shear effect that should play a significant role in the nucleation of solitons. △ Less

Submitted 17 February, 2020; v1 submitted 9 May, 2019; originally announced May 2019.

Journal ref: Phys. Rev. B 101, 125122 (2020)

Reconstructing Extreme Space Weather from Planet Hosting Stars

Authors: V. S. Airapetian, V. Adibekyan, M. Ansdell, D. Alexander, T. Bastian, S. Boro Saikia, A. S. Brun, O. Cohen, M. Cuntz, W. Danchi, J. Davenport, J. DeNolfo, R. DeVore, C. F. Dong, J. J. Drake, K. France, F. Fraschetti, K. Herbst, K. Garcia-Sage, M. Gillon, A. Glocer, J. L. Grenfell, G. Gronoff, N. Gopalswamy, M. Guedel , et al. (58 additional authors not shown)

Abstract: The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditi… ▽ More The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: White Paper submitted to the Astronomy and Astrophysics Decadal Survey (Astro2020), 8 pages, 1 figure

On the Further Structure of the Finite Free Convolutions

Authors: Jonathan Leake, Nick Ryder

Abstract: Since the celebrated resolution of Kadison-Singer (via the Paving Conjecture) by Marcus, Spielman, and Srivastava, much study has been devoted to further understanding and generalizing the techniques of their proof. Specifically, their barrier method was crucial to achieving the required polynomial root bounds on the finite free convolution. But unfortunately this method required individual analys… ▽ More Since the celebrated resolution of Kadison-Singer (via the Paving Conjecture) by Marcus, Spielman, and Srivastava, much study has been devoted to further understanding and generalizing the techniques of their proof. Specifically, their barrier method was crucial to achieving the required polynomial root bounds on the finite free convolution. But unfortunately this method required individual analysis for each usage, and the existence of a larger encapsulating framework is an important open question. In this paper, we make steps toward such a framework by generalizing their root bound to all differential operators. We further conjecture a large class of root bounds, the resolution of which would require for more robust techniques. We further give an important counterexample to a very natural multivariate version of their bound, which if true would have implied tight bounds for the Paving Conjecture. △ Less

Submitted 13 November, 2018; originally announced November 2018.

Comments: 15 pages, 1 figure

Time variations of the non-potential and volume-threading magnetic helicities

Authors: Luis Linan, Étienne Pariat, Kostas Moraitis, Gherardo Valori, James E. Leake

Abstract: Relative magnetic helicity is a gauge invariant quantity suitable for the study of the magnetic helicity content of heliospheric plasmas. Relative magnetic helicity can be decomposed uniquely into two gauge invariant quantities, the magnetic helicity of the non-potential component of the field, and a complementary volume-threading helicity. Recent analysis of numerical experiments simulating the g… ▽ More Relative magnetic helicity is a gauge invariant quantity suitable for the study of the magnetic helicity content of heliospheric plasmas. Relative magnetic helicity can be decomposed uniquely into two gauge invariant quantities, the magnetic helicity of the non-potential component of the field, and a complementary volume-threading helicity. Recent analysis of numerical experiments simulating the generation of solar eruptions have shown that the ratio of the non-potential helicity to the total relative helicity is a clear marker of the eruptivity of the magnetic system, and that the high value of that quantity could be a sufficient condition for the onset of the instability generating the eruptions. The present study introduces the first analytical examination of the time variations of these non-potential and volume-threading helicities. The validity of the analytical formulas derived are confirmed with analysis of three-dimensional (3D) magnetohydrodynamics (MHD) simulations of solar coronal dynamics. Both the analytical investigation, and the numerical application show that, unlike magnetic helicity, the non-potential and the volume-threading helicities are not conserved quantities, even in the ideal MHD regime. A term corresponding to the transformation between the non-potential and volume-threading helicities frequently dominates their dynamics. This finding has an important consequence for their estimation in the solar corona: unlike with relative helicity, their volume coronal evolution cannot be ascertained by the flux of these quantities through the volume's boundaries. Only techniques extrapolating the 3D coronal field will enable both the proper study of the non-potential and volume-threading helicities, and the observational analysis of helicity-based solar-eruptivity proxies. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 24 pages, 14 figures, accepted for publication in the Astrophysical Journal

Sequential eruptions triggered by flux emergence - observations and modeling

Authors: Sally Dacie, Tibor Torok, Pascal Demoulin, Mark Linton, Cooper Downs, Lidia van Driel-Gesztelyi, David Long, James Leake

Abstract: We describe and analyze observations by the Solar Dynamics Observatory of the emergence of a small, bipolar active region within an area of unipolar magnetic flux that was surrounded by a circular, quiescent filament. Within only eight hours of the start of the emergence, a partial splitting of the filament and two consecutive coronal mass ejections took place. We argue that all three dynamic even… ▽ More We describe and analyze observations by the Solar Dynamics Observatory of the emergence of a small, bipolar active region within an area of unipolar magnetic flux that was surrounded by a circular, quiescent filament. Within only eight hours of the start of the emergence, a partial splitting of the filament and two consecutive coronal mass ejections took place. We argue that all three dynamic events occurred as a result of particular magnetic-reconnection episodes between the emerging bipole and the pre-existing coronal magnetic field. In order to substantiate our interpretation, we consider three-dimensional magnetohydrodynamic simulations that model the emergence of magnetic flux in the vicinity of a large-scale coronal flux rope. The simulations qualitatively reproduce most of the reconnection episodes suggested by the observations; as well as the filament-splitting, the first eruption, and the formation of sheared/twisted fields that may have played a role in the second eruption. Our results suggest that the position of emerging flux with respect to the background magnetic configuration is a crucial factor for the resulting evolution, while previous results suggest that parameters such as the orientation or the amount of emerging flux are important as well. This poses a challenge for predicting the onset of eruptions that are triggered by flux emergence, and it calls for a detailed survey of the relevant parameter space by means of numerical simulations. △ Less

Submitted 29 June, 2018; originally announced July 2018.

Comments: 17 pages, 8 figures, accepted for publication by the Astrophysical Journal

Counting Matchings via Capacity Preserving Operators

Authors: Leonid Gurvits, Jonathan Leake

Abstract: The notion of the capacity of a polynomial was introduced by Gurvits around 2005, originally to give drastically simplified proofs of the Van der Waerden lower bound for permanents of doubly stochastic matrices and Schrijver's inequality for perfect matchings of regular bipartite graphs. Since this seminal work, the notion of capacity has been utilized to bound various combinatorial quantities and… ▽ More The notion of the capacity of a polynomial was introduced by Gurvits around 2005, originally to give drastically simplified proofs of the Van der Waerden lower bound for permanents of doubly stochastic matrices and Schrijver's inequality for perfect matchings of regular bipartite graphs. Since this seminal work, the notion of capacity has been utilized to bound various combinatorial quantities and to give polynomial-time algorithms to approximate such quantities (e.g., the number of bases of a matroid). These types of results are often proven by giving bounds on how much a particular differential operator can change the capacity of a given polynomial. In this paper, we unify the theory surrounding such capacity preserving operators by giving tight capacity preservation bounds for all nondegenerate real stability preservers. We then use this theory to give a new proof of a recent result of Csikvári, which settled Friedland's lower matching conjecture. △ Less

Submitted 16 October, 2020; v1 submitted 12 April, 2018; originally announced April 2018.

Exploring Extreme Space Weather Factors of Exoplanetary Habitability

Authors: V. S. Airapetian, V. Adibekyan, M. Ansdell, O. Cohen, M. Cuntz, W. Danchi, C. F. Dong, J. J. Drake, A. Fahrenbach, K. France, K. Garcia-Sage, A. Glocer, J. L. Grenfell, G. Gronoff, H. Hartnett, W. Henning, N. R. Hinkel, A. G. Jensen, M. Jin, P. Kalas, S. R. Kane, K. Kobayashi, R. Kopparapu, J. Leake, M. López-Puertas , et al. (24 additional authors not shown)

Abstract: It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems.… ▽ More It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years. △ Less

Submitted 9 March, 2018; originally announced March 2018.

Comments: 6 pages, the white paper submitted to the US National Academy of Sciences call on Exoplanet Science Strategy

Domain Wall Orientations in Ferroelectric Superlattices Probed with Synchrotron X-Ray Diffraction

Authors: Marios Hadjimichael, Edoardo Zatterin, Stéphanie Fernandez-Peña, Steven J. Leake, Pavlo Zubko

Abstract: Ferroelectric domains in PbTiO$_3$/SrTiO$_3$ superlattices were studied using synchrotron X-ray diffraction. Macroscopic measurements revealed a change in the domain wall orientation from $\left\lbrace 100 \right\rbrace $ to $\left\lbrace 110 \right\rbrace $ crystallographic planes with increasing temperature. The temperature range of this reorientation depends on the ferroelectric layer thickness… ▽ More Ferroelectric domains in PbTiO$_3$/SrTiO$_3$ superlattices were studied using synchrotron X-ray diffraction. Macroscopic measurements revealed a change in the domain wall orientation from $\left\lbrace 100 \right\rbrace $ to $\left\lbrace 110 \right\rbrace $ crystallographic planes with increasing temperature. The temperature range of this reorientation depends on the ferroelectric layer thickness and domain period. Using a nanofocused beam, local changes in domain wall orientation within the buried ferroelectric layers were imaged, both in structurally uniform regions of the sample and near defect sites and argon ion etched patterns. Domain walls were found to exhibit preferential alignment with the straight edges of the etched patterns as well as with structural features associated with defect sites. The distribution of out-of-plane lattice parameters was mapped around one such feature, showing that it is accompanied by inhomogeneous strain and large strain gradients. △ Less

Submitted 5 January, 2018; originally announced January 2018.

A Hydrodynamic Model of Alfvénic Wave Heating in a Coronal Loop and its Chromospheric Footpoints

Authors: Jeffrey W. Reep, Alexander J. B. Russell, Lucas A. Tarr, James E. Leake

Abstract: Alfvénic waves have been proposed as an important energy transport mechanism in coronal loops, capable of delivering energy to both the corona and chromosphere and giving rise to many observed features, of flaring and quiescent regions. In previous work, we established that resistive dissipation of waves (ambipolar diffusion) can drive strong chromospheric heating and evaporation, capable of produ… ▽ More Alfvénic waves have been proposed as an important energy transport mechanism in coronal loops, capable of delivering energy to both the corona and chromosphere and giving rise to many observed features, of flaring and quiescent regions. In previous work, we established that resistive dissipation of waves (ambipolar diffusion) can drive strong chromospheric heating and evaporation, capable of producing flaring signatures. However, that model was based on a simplified assumption that the waves propagate instantly to the chromosphere, an assumption which the current work removes. Via a ray tracing method, we have implemented traveling waves in a field-aligned hydrodynamic simulation that dissipate locally as they propagate along the field line. We compare this method to and validate against the magnetohydrodynamics code Lare3D. We then examine the importance of travel times to the dynamics of the loop evolution, finding that (1) the ionization level of the plasma plays a critical role in determining the location and rate at which waves dissipate; (2) long duration waves effectively bore a hole into the chromosphere, allowing subsequent waves to penetrate deeper than previously expected, unlike an electron beam whose energy deposition rises in height as evaporation reduces the mean-free paths of the electrons; (3) the dissipation of these waves drives a pressure front that propagates to deeper depths, unlike energy deposition by an electron beam. △ Less

Submitted 17 December, 2017; originally announced December 2017.

Comments: Accepted to ApJ

Connecting the q-Multiplicative Convolution and the Finite Difference Convolution

Authors: Jonathan Leake, Nick Ryder

Abstract: In a recent paper, Brändén, Krasikov, and Shapiro consider root location preservation properties of finite difference operators. To this end, the authors describe a natural polynomial convolution operator and conjecture that it preserves root mesh properties. We prove this conjecture using two methods. The first develops a novel connection between the additive (Walsh) and multiplicative (Grace-Sze… ▽ More In a recent paper, Brändén, Krasikov, and Shapiro consider root location preservation properties of finite difference operators. To this end, the authors describe a natural polynomial convolution operator and conjecture that it preserves root mesh properties. We prove this conjecture using two methods. The first develops a novel connection between the additive (Walsh) and multiplicative (Grace-Szegö) convolutions, which can be generically used to transfer results from multiplicative to additive. We then use this to transfer an analogous result, due to Lamprecht, which demonstrates logarithmic root mesh preservation properties of a certain $q$-multiplicative convolution operator. The second method proves the result directly using a modification of Lamprecht's proof of the logarithmic root mesh result. We present his original argument in a streamlined fashion and then make the appropriate alterations to apply it to the additive case. △ Less

Submitted 7 December, 2017; originally announced December 2017.

Electric-Current Neutralization, Magnetic Shear, and Eruptive Activity in Solar Active Regions

Authors: Yang Liu, Xudong Sun, Tibor Török, Viacheslav S. Titov, James E. Leake

Abstract: The physical conditions that determine whether or not solar active regions (ARs) produce strong flares and coronal mass ejections (CMEs) are not yet well understood. Here we investigate the association between electric-current neutralization, magnetic shear along polarity inversion lines (PILs), and eruptive activity in four ARs; two emerging and two well-developed ones. We find that the CME-produ… ▽ More The physical conditions that determine whether or not solar active regions (ARs) produce strong flares and coronal mass ejections (CMEs) are not yet well understood. Here we investigate the association between electric-current neutralization, magnetic shear along polarity inversion lines (PILs), and eruptive activity in four ARs; two emerging and two well-developed ones. We find that the CME-producing ARs are characterized by a strongly non-neutralized total current, while the total current in the ARs that did not produce CMEs is almost perfectly neutralized. The difference in the PIL-shear between these two groups is much less pronounced, which suggests that the degree of current-neutralization may serve as a better proxy for assessing the ability of ARs to produce CMEs. △ Less

Submitted 15 August, 2017; originally announced August 2017.

Comments: ApJL in press

A Representation Theoretic Explanation of the Borcea-Brändén Characterization

Authors: Jonathan Leake

Abstract: In 2009, Borcea and Brändén characterize all linear operators on multivariate polynomials which preserve the property of being non-vanishing (stable) on products of prescribed open circular regions. We give a representation theoretic interpretation of their findings, which generalizes and simplifies their result and leads to a conceptual unification of many related results in polynomial stability… ▽ More In 2009, Borcea and Brändén characterize all linear operators on multivariate polynomials which preserve the property of being non-vanishing (stable) on products of prescribed open circular regions. We give a representation theoretic interpretation of their findings, which generalizes and simplifies their result and leads to a conceptual unification of many related results in polynomial stability theory. At the heart of this unification is a generalized Grace's theorem which addresses polynomials whose roots are all contained in some real interval or ray. This generalization allows us to extend the Borcea-Brändén result to characterize a certain subclass of the linear operators which preserve such polynomials. △ Less

Submitted 6 August, 2021; v1 submitted 19 June, 2017; originally announced June 2017.

Testing predictors of eruptivity using parametric flux emergence simulations

Authors: Chloé Guennou, Etienne Pariat, Nicole Vilmer, James E. Leake

Abstract: Solar flares and coronal mass ejections (CMEs) are among the most energetic events in the solar system, impacting the near-Earth environment. Flare productivity is empirically known to be correlated with the size and complexity of active regions. Several indicators, based on magnetic-field data from active regions, have been tested for flare forecasting in recent years. None of these indicators, o… ▽ More Solar flares and coronal mass ejections (CMEs) are among the most energetic events in the solar system, impacting the near-Earth environment. Flare productivity is empirically known to be correlated with the size and complexity of active regions. Several indicators, based on magnetic-field data from active regions, have been tested for flare forecasting in recent years. None of these indicators, or combinations thereof, have yet demonstrated an unambiguous eruption or flare criterion. Furthermore, numerical simulations have been only barely used to test the predictability of these parameters. In this context, we used the 3D parametric MHD numerical simulations of the self-consistent formation of the flux emergence of a twisted flux tube, inducing the formation of stable and unstable magnetic flux ropes of Leake (2013, 2014). We use these numerical simulations to investigate the eruptive signatures observable in various magnetic scalar parameters and provide highlights on data analysis processing. Time series of 2D photospheric-like magnetograms are used from parametric simulations of stable and unstable flux emergence, to compute a list of about 100 different indicators. This list includes parameters previously used for operational forecasting, physical parameters used for the first time, as well as new quantities specifically developed for this purpose. Our results indicate that only parameters measuring the total non-potentiality of active regions associated with magnetic inversion line properties, such as the Falconer parameters $L_{ss}$, $WL_{ss}$, $L_{sg}$ and $WL_{sg}$, as well as the new current integral $WL_{sc}$ and length $L_{sc}$ parameters, present a significant ability to distinguish the eruptive cases of the model from the non-eruptive cases, possibly indicating that they are promising flare and eruption predictors. △ Less

Submitted 15 June, 2017; originally announced June 2017.

Comments: 46 pages, 16 figures, accepted for publication in Space Weather and Space Climate on June, 8th

Magnetic Flux Emergence and Decay Rates for Preceder and Follower Sunspots Observed with HMI

Authors: A. A. Norton, E. H. Jones, M. G. Linton, J. E. Leake

Abstract: We quantify the emergence and decay rates of preceder (p) and follower (f) sunspots within ten active regions from 2010-2014 using Space-weather Helioseismic Magnetic Imager Active Region Patch data. The sunspots are small- to mid-sized regions and contain a signed flux within a single polarity sunspot of $(1.1-6.5)\times 10^{21}$ Mx. The net unsigned flux within the regions, including plage, rang… ▽ More We quantify the emergence and decay rates of preceder (p) and follower (f) sunspots within ten active regions from 2010-2014 using Space-weather Helioseismic Magnetic Imager Active Region Patch data. The sunspots are small- to mid-sized regions and contain a signed flux within a single polarity sunspot of $(1.1-6.5)\times 10^{21}$ Mx. The net unsigned flux within the regions, including plage, ranges from $(5.1-20)\times 10^{21}$ Mx. Rates are calculated with and without intensity contours to differentiate between sunspot formation and flux emergence. Signed flux emergence rates, calculated with intensity contours, for the p (f) spots average $6.8$ (4.9) $\times 10^{19}$ Mx h$^{-1}$, while decay rates are $-1.9 (-3.4)\times 10^{19}$ Mx h$^{-1}$. The mean, signed flux emergence rate of the regions, including plage, is $7.1 \times 10^{19}$ Mx h$^{-1}$ for a mean peak flux of $5.9 \times 10^{21}$ Mx. Using a synthesis of these results and others reported previously, there is a clear trend for larger flux regions to emerge faster than smaller ones. Observed emergence rates ($dφ/dt$, Mx h$^{-1}$) scale with total signed peak flux, $\tildeφ_{max}$, as a power law with an exponent of 0.36, i.e., $dφ/dt = A \tildeφ_{max}^{0.36}$. The observed rates may assist in constraining the boundary and initial conditions in simulations which already demonstrate increased rates for flux tubes with higher buoyancy and twist, or in the presence of a strong upflow. Overall, the observed emergence rates are smaller than those in simulations, which may indicate a slower rise of the flux in the interior than captured in simulations. △ Less

Submitted 4 May, 2017; originally announced May 2017.

Comments: Accepted by ApJ, 17 pages, 6 figures

Relative magnetic helicity as a diagnostic of solar eruptivity

Authors: E. Pariat, J. E. Leake, G. Valori, M. G. Linton, F. P. Zuccarello, K. Dalmasse

Abstract: The discovery of clear criteria that can deterministically describe the eruptive state of a solar active region would lead to major improvements on space weather predictions. Using series of numerical simulations of the emergence of a magnetic flux rope in a magnetized coronal, leading either to eruptions or to stable configurations, we test several global scalar quantities for the ability to disc… ▽ More The discovery of clear criteria that can deterministically describe the eruptive state of a solar active region would lead to major improvements on space weather predictions. Using series of numerical simulations of the emergence of a magnetic flux rope in a magnetized coronal, leading either to eruptions or to stable configurations, we test several global scalar quantities for the ability to discriminate between the eruptive and the non-eruptive simulations. From the magnetic field generated by the three-dimensional magnetohydrodynamical simulations, we compute and analyse the evolution of the magnetic flux, of the magnetic energy and its decomposition into potential and free energies, and of the relative magnetic helicity and its decomposition. Unlike the magnetic flux and magnetic energies, magnetic helicities are able to markedly distinguish the eruptive from the non-eruptive simulations. We find that the ratio of the magnetic helicity of the current-carrying magnetic field to the total relative helicity presents the highest values for the eruptive simulations, in the pre-eruptive phase only. We observe that the eruptive simulations do not possess the highest value of total magnetic helicity. In the framework of our numerical study, the magnetic energies and the total relative helicity do not correspond to good eruptivity proxies. Our study highlights that the ratio of magnetic helicities diagnoses very clearly the eruptive potential of our parametric simulations. Our study shows that magnetic-helicity-based quantities may be very efficient for the prediction of solar eruptions. △ Less

Submitted 30 March, 2017; originally announced March 2017.

Journal ref: A&A 601, A125 (2017)

Testing the Accuracy of Data Driven MHD Simulations of Active Region Evolution

Authors: James E. Leake, Mark G. Linton, Peter W. Schuck

Abstract: Models for the evolution of the solar coronal magnetic field are vital for understanding solar activity, yet the best measurements of magnetic field lie at the photosphere, necessitating the development of coronal models which are \textit{"data-driven"} at the photosphere. We present an investigation to determine the feasibility and accuracy of such methods. Our validation framework uses a simulat… ▽ More Models for the evolution of the solar coronal magnetic field are vital for understanding solar activity, yet the best measurements of magnetic field lie at the photosphere, necessitating the development of coronal models which are \textit{"data-driven"} at the photosphere. We present an investigation to determine the feasibility and accuracy of such methods. Our validation framework uses a simulation of active region (AR) formation, modeling the emergence of magnetic flux from the convection zone to the corona, as a ground-truth dataset, to supply both the photospheric information, and to perform the validation of the data-driven method. We focus our investigation on how the data-driven model accuracy depends on the temporal frequency of the driving data. The Helioseismic and Magnetic Imager on NASA's Solar {\jl Dynamics} Observatory produces full-disk vector magnetic field measurements at a 12 minute cadence. Using our framework we show that ARs that emerge over 25 hours can be modeled by the data-driving method with only $\sim$1\% error in the free magnetic energy, assuming the photospheric information is specified every 12 minutes. However, for rapidly evolving features, under-sampling of the dynamics at this cadence leads to a strobe effect, generating large electric currents and incorrect coronal morphology and energies. We derive a sampling condition for the driving cadence based on the evolution of these small-scale features, and show that higher-cadence driving can lead to acceptable errors. Future work will investigate the source of errors associated with deriving plasma variables from the photospheric magnetograms as well as other sources of errors such as reduced resolution and instrument bias and noise. △ Less

Submitted 22 February, 2017; originally announced February 2017.

Science Objective: Understanding Energy Transport by Alfvénic Waves in Solar Flares

Authors: Jeffrey W. Reep, Harry P. Warren, James E. Leake, Lucas A. Tarr, Alexander J. B. Russell, Graham S. Kerr, Hugh S. Hudson

Abstract: Solar flares are driven by the release of magnetic energy from reconnection events in the solar corona, whereafter energy is transported to the chromosphere, heating the plasma and causing the characteristic radiative losses. In the collisional thick-target model, electrons accelerated to energies exceeding 10 keV traverse the corona and impact the chromosphere, where they deposit their energy thr… ▽ More Solar flares are driven by the release of magnetic energy from reconnection events in the solar corona, whereafter energy is transported to the chromosphere, heating the plasma and causing the characteristic radiative losses. In the collisional thick-target model, electrons accelerated to energies exceeding 10 keV traverse the corona and impact the chromosphere, where they deposit their energy through collisions with the much denser plasma in the lower atmosphere. While there are undoubtedly high energy non-thermal electrons accelerated in flares, it is unclear whether these electron beams are the sole mechanism of energy transport, or whether they only dominate in certain phases of the flare's evolution. Alfvénic waves are generated during the post-reconnection relaxation of magnetic field lines, so it is important to examine their role in energy transport. △ Less

Submitted 6 February, 2017; originally announced February 2017.

Comments: Next Generation Solar Physics Mission White paper, 4 pages, 2 figures

Magnetoacoustic Waves in a Stratified Atmosphere with a Magnetic Null Point

Authors: Lucas A. Tarr, Mark Linton, James Leake

Abstract: We perform nonlinear MHD simulations to study the propagation of magnetoacoustic waves from the photosphere to the low corona. We focus on a 2D system with a gravitationally stratified atmosphere and three photospheric concentrations of magnetic flux that give rise to a magnetic null point with a magnetic dome topology. We find that a single wavepacket introduced at the lower boundary splits into… ▽ More We perform nonlinear MHD simulations to study the propagation of magnetoacoustic waves from the photosphere to the low corona. We focus on a 2D system with a gravitationally stratified atmosphere and three photospheric concentrations of magnetic flux that give rise to a magnetic null point with a magnetic dome topology. We find that a single wavepacket introduced at the lower boundary splits into multiple secondary wavepackets. A portion of the packet refracts towards the null due to the varying Alfvén speed. Waves incident on the equipartition contour surrounding the null, where the sound and Alfvén speeds coincide, partially transmit, reflect, and mode convert between branches of the local dispersion relation. Approximately $15.5\%$ of the wavepacket's initial energy ($E_{input}$) converges on the null, mostly as a fast magnetoacoustic wave. Conversion is very efficient: $70\%$ of the energy incident on the null is converted to slow modes propagating away from the null, $7\%$ leaves as a fast wave, and the remaining $23\%$ (0.036$E_{input}$) is locally dissipated. The acoustic energy leaving the null is strongly concentrated along field lines near each of the null's four separatrices. The portion of the wavepacket that refracts towards the null, and the amount of current accumulation, depends on the vertical and horizontal wavenumbers and the centroid position of the wavepacket as it crosses the photosphere. Regions that refract towards or away from the null do not simply coincide with regions of open versus closed magnetic field or the local field orientation. We also modeled wavepacket propagation using a WKB method and found that it agrees qualitatively, though not quantitatively, with the results of the numerical simulation. △ Less

Submitted 21 December, 2016; v1 submitted 18 November, 2016; originally announced November 2016.

Comments: 33 pages, 11 figures, submitted to ApJ

The Contribution of Coronal Jets To The Solar Wind

Authors: R. Lionello, T. TÖrÖk, V. S. Titov, J. E. Leake, Z. MikiĆ, J. A. Linker, M. G. Linton

Abstract: Transient collimated plasma eruptions in the solar corona, commonly known as coronal (or X-ray) jets, are among the most interesting manifestations of solar activity. It has been suggested that these events contribute to the mass and energy content of the corona and solar wind, but the extent of these contributions remains uncertain. We have recently modeled the formation and evolution of coronal… ▽ More Transient collimated plasma eruptions in the solar corona, commonly known as coronal (or X-ray) jets, are among the most interesting manifestations of solar activity. It has been suggested that these events contribute to the mass and energy content of the corona and solar wind, but the extent of these contributions remains uncertain. We have recently modeled the formation and evolution of coronal jets using a three-dimensional (3D) magnetohydrodynamic (MHD) code with thermodynamics in a large spherical domain that includes the solar wind. Our model is coupled to 3D MHD flux-emergence simulations, i.e, we use boundary conditions provided by such simulations to drive a time-dependent coronal evolution. The model includes parametric coronal heating, radiative losses, and thermal conduction, which enables us to simulate the dynamics and plasma properties of coronal jets in a more realistic manner than done so far. Here we employ these simulations to calculate the amount of mass and energy transported by coronal jets into the outer corona and inner heliosphere. Based on observed jet-occurrence rates, we then estimate the total contribution of coronal jets to the mass and energy content of the solar wind to (0.4-3.0) % and (0.3-1.0) %, respectively. Our results are largely consistent with the few previous rough estimates obtained from observations, supporting the conjecture that coronal jets provide only a small amount of mass and energy to the solar wind. We emphasize, however, that more advanced observations and simulations are needed to substantiate this conjecture. △ Less

Submitted 10 October, 2016; originally announced October 2016.

Generalizations of the Matching Polynomial to the Multivariate Independence Polynomial

Authors: Jonathan Leake, Nick Ryder

Abstract: We generalize two main theorems of matching polynomials of undirected simple graphs, namely, real-rootedness and the Heilmann-Lieb root bound. Viewing the matching polynomial of a graph $G$ as the independence polynomial of the line graph of $G$, we determine conditions for the extension of these theorems to the independence polynomial of any graph. In particular, we show that a stability-like pro… ▽ More We generalize two main theorems of matching polynomials of undirected simple graphs, namely, real-rootedness and the Heilmann-Lieb root bound. Viewing the matching polynomial of a graph $G$ as the independence polynomial of the line graph of $G$, we determine conditions for the extension of these theorems to the independence polynomial of any graph. In particular, we show that a stability-like property of the multivariate independence polynomial characterizes claw-freeness. Finally, we give and extend multivariate versions of Godsil's theorems on the divisibility of matching polynomials of trees related to $G$. △ Less

Submitted 19 February, 2019; v1 submitted 3 October, 2016; originally announced October 2016.

Comments: This research was supported by NSF Grant CCF-155375, 20 pages

Mixed Determinants and the Kadison-Singer problem

Authors: Jonathan Leake, Mohan Ravichandran

Abstract: We adapt the arguments of Marcus, Spielman and Srivastava in their proof of the Kadison-Singer problem to prove improved paving estimates. Working with Anderson's paving formulation of Kadison-Singer instead of Weaver's vector balancing version, we show that the machinery of interlacing polyomials due to Marcus, Spielman and Srivastava works in this setting as well. The relevant expected character… ▽ More We adapt the arguments of Marcus, Spielman and Srivastava in their proof of the Kadison-Singer problem to prove improved paving estimates. Working with Anderson's paving formulation of Kadison-Singer instead of Weaver's vector balancing version, we show that the machinery of interlacing polyomials due to Marcus, Spielman and Srivastava works in this setting as well. The relevant expected characteristic polynomials turn out to be related to the so called "mixed determinants" that have been carefully studied by Borcea and Branden. This technique allows us to show that any projection with diagonal entries strictly less than $\frac{1}{4}$ can be two paved, matching recent results of Bownik, Casazza, Marcus and Speegle, though our estimates are asymptotically weaker. We also show that any projection with diagonal entries at most $\frac{1}{2}$ can be four paved, yielding improvements over currently known estimates. We also relate the problem of finding optimal paving estimates to bounding the root intervals of a natural one parameter deformation of the characteristic polynomial of a matrix that turns out to have some remarkable combinatorial properties. △ Less

Submitted 1 November, 2018; v1 submitted 14 September, 2016; originally announced September 2016.

Comments: 26 pages, no figures

Sunquake Generation by Coronal Magnetic Restructuring

Authors: Alexander J. B. Russell, Michaela K. Mooney, James E. Leake, Hugh S. Hudson

Abstract: Sunquakes are the surface signatures of acoustic waves in the Sun's interior that are produced by some but not all flares and coronal mass ejections (CMEs). This paper explores a mechanism for sunquake generation by the magnetic field changes that occur during flares and CMEs, using MHD simulations with a semiempirical FAL-C atmosphere to demonstrate the generation of acoustic waves in the interio… ▽ More Sunquakes are the surface signatures of acoustic waves in the Sun's interior that are produced by some but not all flares and coronal mass ejections (CMEs). This paper explores a mechanism for sunquake generation by the magnetic field changes that occur during flares and CMEs, using MHD simulations with a semiempirical FAL-C atmosphere to demonstrate the generation of acoustic waves in the interior in response to changing magnetic tilt in the corona. We find that Alfven-sound resonance combined with the ponderomotive force produces acoustic waves in the interior with sufficient energy to match sunquake observations when the magnetic field angle changes by the order of 10 degrees in a region where the coronal field strength is a few hundred gauss or more. The most energetic sunquakes are produced when the coronal field is strong, while the variation of magnetic field strength with height and the time scale of the tilt change are of secondary importance. △ Less

Submitted 22 August, 2016; v1 submitted 26 February, 2016; originally announced February 2016.

Comments: 6 pages, 3 figures; accepted to ApJ

Modeling Jets in the Corona and Solar Wind

Authors: T. Torok, R. Lionello, V. S. Titov, J. E. Leake, Z. Mikic, J. A. Linker, M. G. Linton

Abstract: Coronal jets are transient, collimated eruptions that occur in regions of predominantly open magnetic field in the solar corona. Our understanding of these events has greatly evolved in recent years but several open questions, such as the contribution of coronal jets to the solar wind, remain. Here we present an overview of the observations and numerical modeling of coronal jets, followed by a bri… ▽ More Coronal jets are transient, collimated eruptions that occur in regions of predominantly open magnetic field in the solar corona. Our understanding of these events has greatly evolved in recent years but several open questions, such as the contribution of coronal jets to the solar wind, remain. Here we present an overview of the observations and numerical modeling of coronal jets, followed by a brief description of "next-generation" simulations that include an advanced description of the energy transfer in the corona ("thermodynamic MHD"), large spherical computational domains, and the solar wind. These new models will allow us to address some of the open questions. △ Less

Submitted 30 November, 2015; originally announced November 2015.

Comments: 12 pages, 4 figures, to appear in ASP proceedings

Toward A Self Consistent MHD Model of Chromospheres and Winds From Late Type Evolved Stars

Authors: V. S. Airapetian, J. E. Leake, K. G. Carpenter

Abstract: We present the first magnetohydrodynamic model of the stellar chromospheric heating and acceleration of the outer atmospheres of cool evolved stars, using alpha Tau as a case study. We used a 1.5D MHD code with a generalized Ohm's law that accounts for the effects of partial ionization in the stellar atmosphere to study Alfven wave dissipation and wave reflection. We have demonstrated that due to… ▽ More We present the first magnetohydrodynamic model of the stellar chromospheric heating and acceleration of the outer atmospheres of cool evolved stars, using alpha Tau as a case study. We used a 1.5D MHD code with a generalized Ohm's law that accounts for the effects of partial ionization in the stellar atmosphere to study Alfven wave dissipation and wave reflection. We have demonstrated that due to inclusion of the effects of ion-neutral collisions in magnetized weakly ionized chromospheric plasma on resistivity and the appropriate grid resolution, the numerical resistivity becomes 1-2 orders of magnitude smaller than the physical resistivity. The motions introduced by non-linear transverse Alfven waves can explain non-thermally broadened and non-Gaussian profiles of optically thin UV lines forming in the stellar chromosphere of alpha Tau and other late-type giant and supergiant stars. The calculated heating rates in the stellar chromosphere due to resistive (Joule) dissipation of electric currents, induced by upward propagating non-linear Alfven waves, are consistent with observational constraints on the net radiative losses in UV lines and the continuum from alpha Tau. At the top of the chromosphere, Alfven waves experience significant reflection, producing downward propagating transverse waves that interact with upward propagating waves and produce velocity shear in the chromosphere. Our simulations also suggest that momentum deposition by non-linear Alfven waves becomes significant in the outer chromosphere at 1 stellar radius from the photosphere. The calculated terminal velocity and the mass loss rate are consistent with the observationally derived wind properties in alpha Tau. △ Less

Submitted 12 September, 2014; originally announced September 2014.

Comments: To be appeared in 18th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun Proceedings of Lowell Observatory (9-13 June 2014) Edited by G. van Belle & H. Harris

Simulations of Emerging Magnetic Flux. II: The formation of Unstable Coronal Flux Ropes and the Initiation of CMEs

Authors: J. E. Leake, M. G. Linton, S. K. Antiochos

Abstract: We present results from 3D magnetohydrodynamic (MHD) simulations of the emergence of a twisted convection zone flux tube into a pre-existing coronal dipole field. As in previous simulations, following the partial emergence of the sub-surface flux into the corona, a combination of vortical motions and internal magnetic reconnection forms a coronal flux rope. Then, in the simulations presented here,… ▽ More We present results from 3D magnetohydrodynamic (MHD) simulations of the emergence of a twisted convection zone flux tube into a pre-existing coronal dipole field. As in previous simulations, following the partial emergence of the sub-surface flux into the corona, a combination of vortical motions and internal magnetic reconnection forms a coronal flux rope. Then, in the simulations presented here, external reconnection between the emerging field and the pre-existing dipole coronal field allows further expansion of the coronal flux rope into the corona. After sufficient expansion, internal reconnection occurs beneath the coronal flux rope axis, and the flux rope erupts up to the top boundary of the simulation domain ($\sim$ 36 Mm above the surface). We find that the presence of a pre-existing field, orientated in a direction to facilitate reconnection with the emerging field, is vital to the fast rise of the coronal flux rope. The simulations shown in this paper are able to self-consistently create many of the surface and coronal signatures used by coronal mass ejection (CME) models. These signatures include: surface shearing and rotational motions; quadrupolar geometry above the surface; central sheared arcades reconnecting with oppositely orientated overlying dipole fields; the formation of coronal flux ropes underlying potential coronal field; and internal reconnection which resembles the classical flare reconnection scenario. This suggests that proposed mechanisms for the initiation of a CME, such as "magnetic breakout", are operating during the emergence of new active regions. △ Less

Submitted 26 March, 2014; v1 submitted 7 February, 2014; originally announced February 2014.