Showing 1–50 of 65 results for author: McKenzie, D

Normalizing Basis Functions: Approximate Stationary Models for Large Spatial Data

Authors: Antony Sikorski, Daniel McKenzie, Douglas Nychka

Abstract: In geostatistics, traditional spatial models often rely on the Gaussian Process (GP) to fit stationary covariances to data. It is well known that this approach becomes computationally infeasible when dealing with large data volumes, necessitating the use of approximate methods. A powerful class of methods approximate the GP as a sum of basis functions with random coefficients. Although this techni… ▽ More In geostatistics, traditional spatial models often rely on the Gaussian Process (GP) to fit stationary covariances to data. It is well known that this approach becomes computationally infeasible when dealing with large data volumes, necessitating the use of approximate methods. A powerful class of methods approximate the GP as a sum of basis functions with random coefficients. Although this technique offers computational efficiency, it does not inherently guarantee a stationary covariance. To mitigate this issue, the basis functions can be "normalized" to maintain a constant marginal variance, avoiding unwanted artifacts and edge effects. This allows for the fitting of nearly stationary models to large, potentially non-stationary datasets, providing a rigorous base to extend to more complex problems. Unfortunately, the process of normalizing these basis functions is computationally demanding. To address this, we introduce two fast and accurate algorithms to the normalization step, allowing for efficient prediction on fine grids. The practical value of these algorithms is showcased in the context of a spatial analysis on a large dataset, where significant computational speedups are achieved. While implementation and testing are done specifically within the LatticeKrig framework, these algorithms can be adapted to other basis function methods operating on regular grids. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Determining the nanoflare heating frequency of an X-ray Bright Point observed by MaGIXS

Authors: Biswajit Mondal, P. S. Athiray, Amy R. Winebarger, Sabrina L. Savage, Ken Kobayashi, Stephen Bradshaw, Will Barnes, Patrick R. Champey, Peter Cheimets, Jaroslav Dudik, Leon Golub, Helen E. Mason, David E. McKenzie, Christopher S. Moore, Chad Madsen, Katharine K. Reeves, Paola Testa, Genevieve D. Vigil, Harry P. Warren, Robert W. Walsh, Giulio Del Zanna

Abstract: Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multi-million kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission… ▽ More Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multi-million kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission from an X-ray bright point (XBP) observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with nearest-available observations from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO) and X-Ray Telescope (XRT) onboard Hinode observatory. The length and magnetic field strength of the coronal loops are derived from the linear-force-free extrapolation of the observed photospheric magnetogram by Helioseismic and Magnetic Imager (HMI) onboard SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the loop length and magnetic field strength. The simulation results are then compared and matched against the measured intensity from AIA, XRT, and MaGIXS. Our model results indicate the observed emissions from the XBP under study could be well matched by a distribution of nanoflares with average delay times 1500 s to 3000 s, which suggest that the heating is dominated by high-frequency events. Further, we demonstrate the high sensitivity of MaGIXS and XRT to diagnose the heating frequency using this method, while AIA passbands are found to be the least sensitive. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: Accepted for publication in the Astrophysical Journal (ApJ)

A crystallographic approach to symmetry-breaking in fluid layers

Authors: John F. Rudge, Dan McKenzie

Abstract: Symmetry-breaking bifurcations, where a flow state with a certain symmetry undergoes a transition to state with a different symmetry, are ubiquitous in fluid mechanics. Much can be understood about the nature of these transitions from symmetry alone, using the theory of groups and their representations. Here we show how the extensive databases on groups in crystallography can be exploited to yield… ▽ More Symmetry-breaking bifurcations, where a flow state with a certain symmetry undergoes a transition to state with a different symmetry, are ubiquitous in fluid mechanics. Much can be understood about the nature of these transitions from symmetry alone, using the theory of groups and their representations. Here we show how the extensive databases on groups in crystallography can be exploited to yield insights into fluid-dynamical problems. In particular, we demonstrate the application of the crystallographic layer groups to problems in fluid layers, using thermal convection as an example. Crystallographic notation provides a concise and unambiguous description of the symmetries involved, and we advocate its broader use by the fluid dynamics community. △ Less

Submitted 9 April, 2024; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: 29 pages, 10 figures, 3 supplementary tables. Minor edits from initial submission including new figure 1

Fermat Distances: Metric Approximation, Spectral Convergence, and Clustering Algorithms

Authors: Nicolás García Trillos, Anna Little, Daniel McKenzie, James M. Murphy

Abstract: We analyze the convergence properties of Fermat distances, a family of density-driven metrics defined on Riemannian manifolds with an associated probability measure. Fermat distances may be defined either on discrete samples from the underlying measure, in which case they are random, or in the continuum setting, in which they are induced by geodesics under a density-distorted Riemannian metric. We… ▽ More We analyze the convergence properties of Fermat distances, a family of density-driven metrics defined on Riemannian manifolds with an associated probability measure. Fermat distances may be defined either on discrete samples from the underlying measure, in which case they are random, or in the continuum setting, in which they are induced by geodesics under a density-distorted Riemannian metric. We prove that discrete, sample-based Fermat distances converge to their continuum analogues in small neighborhoods with a precise rate that depends on the intrinsic dimensionality of the data and the parameter governing the extent of density weighting in Fermat distances. This is done by leveraging novel geometric and statistical arguments in percolation theory that allow for non-uniform densities and curved domains. Our results are then used to prove that discrete graph Laplacians based on discrete, sample-driven Fermat distances converge to corresponding continuum operators. In particular, we show the discrete eigenvalues and eigenvectors converge to their continuum analogues at a dimension-dependent rate, which allows us to interpret the efficacy of discrete spectral clustering using Fermat distances in terms of the resulting continuum limit. The perspective afforded by our discrete-to-continuum Fermat distance analysis leads to new clustering algorithms for data and related insights into efficient computations associated to density-driven spectral clustering. Our theoretical analysis is supported with numerical simulations and experiments on synthetic and real image data. △ Less

Submitted 7 July, 2023; originally announced July 2023.

Fundamentals of impulsive energy release in the corona

Authors: Albert Y. Shih, Lindsay Glesener, Säm Krucker, Silvina Guidoni, Steven Christe, Katharine K. Reeves, Szymon Gburek, Amir Caspi, Meriem Alaoui, Joel Allred, Marina Battaglia, Wayne Baumgartner, Brian Dennis, James Drake, Keith Goetz, Leon Golub, Iain Hannah, Laura Hayes, Gordon Holman, Andrew Inglis, Jack Ireland, Graham Kerr, James Klimchuk, David McKenzie, Christopher S. Moore , et al. (8 additional authors not shown)

Abstract: It is essential that there be coordinated and co-optimized observations in X-rays, gamma-rays, and EUV during the peak of solar cycle 26 (~2036) to significantly advance our understanding of impulsive energy release in the corona. The open questions include: What are the physical origins of space-weather events? How are particles accelerated at the Sun? How is impulsively released energy transport… ▽ More It is essential that there be coordinated and co-optimized observations in X-rays, gamma-rays, and EUV during the peak of solar cycle 26 (~2036) to significantly advance our understanding of impulsive energy release in the corona. The open questions include: What are the physical origins of space-weather events? How are particles accelerated at the Sun? How is impulsively released energy transported throughout the solar atmosphere? How is the solar corona heated? Many of the processes involved in triggering, driving, and sustaining solar eruptive events -- including magnetic reconnection, particle acceleration, plasma heating, and energy transport in magnetized plasmas -- also play important roles in phenomena throughout the Universe. This set of observations can be achieved through a single flagship mission or, with foreplanning, through a combination of major missions (e.g., the previously proposed FIERCE mission concept). △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 5 pages, 1 figure

Journal ref: Bulletin of the AAS, Vol. 55, Issue 3, Whitepaper #364 (5pp); 2023 July 31

It begins with a boundary: A geometric view on probabilistically robust learning

Authors: Leon Bungert, Nicolás García Trillos, Matt Jacobs, Daniel McKenzie, Đorđe Nikolić, Qingsong Wang

Abstract: Although deep neural networks have achieved super-human performance on many classification tasks, they often exhibit a worrying lack of robustness towards adversarially generated examples. Thus, considerable effort has been invested into reformulating Empirical Risk Minimization (ERM) into an adversarially robust framework. Recently, attention has shifted towards approaches which interpolate betwe… ▽ More Although deep neural networks have achieved super-human performance on many classification tasks, they often exhibit a worrying lack of robustness towards adversarially generated examples. Thus, considerable effort has been invested into reformulating Empirical Risk Minimization (ERM) into an adversarially robust framework. Recently, attention has shifted towards approaches which interpolate between the robustness offered by adversarial training and the higher clean accuracy and faster training times of ERM. In this paper, we take a fresh and geometric view on one such method -- Probabilistically Robust Learning (PRL) (Robey et al., ICML, 2022). We propose a geometric framework for understanding PRL, which allows us to identify a subtle flaw in its original formulation and to introduce a family of probabilistic nonlocal perimeter functionals to address this. We prove existence of solutions using novel relaxation methods and study properties as well as local limits of the introduced perimeters. △ Less

Submitted 30 May, 2023; originally announced May 2023.

Magnetic Energy Powers the Corona: How We Can Understand its 3D Storage & Release

Authors: Amir Caspi, Daniel B. Seaton, Roberto Casini, Cooper Downs, Sarah E. Gibson, Holly Gilbert, Lindsay Glesener, Silvina E. Guidoni, J. Marcus Hughes, David McKenzie, Joseph Plowman, Katharine K. Reeves, Pascal Saint-Hilaire, Albert Y. Shih, Matthew J. West

Abstract: The coronal magnetic field is the prime driver behind many as-yet unsolved mysteries: solar eruptions, coronal heating, and the solar wind, to name a few. It is, however, still poorly observed and understood. We highlight key questions related to magnetic energy storage, release, and transport in the solar corona, and their relationship to these important problems. We advocate for new and multi-po… ▽ More The coronal magnetic field is the prime driver behind many as-yet unsolved mysteries: solar eruptions, coronal heating, and the solar wind, to name a few. It is, however, still poorly observed and understood. We highlight key questions related to magnetic energy storage, release, and transport in the solar corona, and their relationship to these important problems. We advocate for new and multi-point co-optimized measurements, sensitive to magnetic field and other plasma parameters, spanning from optical to $γがんま$-ray wavelengths, to bring closure to these long-standing and fundamental questions. We discuss how our approach can fully describe the 3D magnetic field, embedded plasma, particle energization, and their joint evolution to achieve these objectives. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 16 pages, 3 figures

Journal ref: Bulletin of the AAS, Vol. 55, Issue 3, Whitepaper #049 (16pp); 2023 July 31

Improving Multi-Dimensional Data Formats, Access, and Assimilation Tools for the Twenty-First Century

Authors: Daniel B. Seaton, Amir Caspi, Roberto Casini, Cooper Downs, Sarah E. Gibson, Holly Gilbert, Lindsay Glesener, Silvina E. Guidoni, J. Marcus Hughes, David McKenzie, Joseph Plowman, Katharine K. Reeves, Pascal Saint-Hilaire, Albert Y. Shih, Matthew J. West

Abstract: Heliophysics image data largely relies on a forty-year-old ecosystem built on the venerable Flexible Image Transport System (FITS) data standard. While many in situ measurements use newer standards, they are difficult to integrate with multiple data streams required to develop global understanding. Additionally, most data users still engage with data in much the same way as they did decades ago. H… ▽ More Heliophysics image data largely relies on a forty-year-old ecosystem built on the venerable Flexible Image Transport System (FITS) data standard. While many in situ measurements use newer standards, they are difficult to integrate with multiple data streams required to develop global understanding. Additionally, most data users still engage with data in much the same way as they did decades ago. However, contemporary missions and models require much more complex support for 3D multi-parameter data, robust data assimilation strategies, and integration of multiple individual data streams required to derive complete physical characterizations of the Sun and Heliospheric plasma environment. In this white paper we highlight some of the 21$^\mathsf{st}$ century challenges for data frameworks in heliophysics, consider an illustrative case study, and make recommendations for important steps the field can take to modernize its data products and data usage models. Our specific recommendations include: (1) Investing in data assimilation capability to drive advanced data-constrained models, (2) Investing in new strategies for integrating data across multiple instruments to realize measurements that cannot be produced from single observations, (3) Rethinking old data use paradigms to improve user access, develop deep understanding, and decrease barrier to entry for new datasets, and (4) Investing in research on data formats better suited for multi-dimensional data and cloud-based computing. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 9 pages, 3 figures

Journal ref: Bulletin of the AAS, Vol. 55, Issue 3, Whitepaper #361 (9pp); 2023 July 31

COMPLETE: A flagship mission for complete understanding of 3D coronal magnetic energy release

Authors: Amir Caspi, Daniel B. Seaton, Roberto Casini, Cooper Downs, Sarah E. Gibson, Holly Gilbert, Lindsay Glesener, Silvina E. Guidoni, J. Marcus Hughes, David McKenzie, Joseph Plowman, Katharine K. Reeves, Pascal Saint-Hilaire, Albert Y. Shih, Matthew J. West

Abstract: COMPLETE is a flagship mission concept combining broadband spectroscopic imaging and comprehensive magnetography from multiple viewpoints around the Sun to enable tomographic reconstruction of 3D coronal magnetic fields and associated dynamic plasma properties, which provide direct diagnostics of energy release. COMPLETE re-imagines the paradigm for solar remote-sensing observations through purpos… ▽ More COMPLETE is a flagship mission concept combining broadband spectroscopic imaging and comprehensive magnetography from multiple viewpoints around the Sun to enable tomographic reconstruction of 3D coronal magnetic fields and associated dynamic plasma properties, which provide direct diagnostics of energy release. COMPLETE re-imagines the paradigm for solar remote-sensing observations through purposefully co-optimized detectors distributed on multiple spacecraft that operate as a single observatory, linked by a comprehensive data/model assimilation strategy to unify individual observations into a single physical framework. We describe COMPLETE's science goals, instruments, and mission implementation. With targeted investment by NASA, COMPLETE is feasible for launch in 2032 to observe around the maximum of Solar Cycle 26. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 10 pages, 6 figures, 1 table

Journal ref: Bulletin of the AAS, Vol. 55, Issue 3, Whitepaper #048 (10pp); 2023 July 31

Evidence for the Operation of the Hanle and Magneto-Optical Effects in the Scattering Polarization Signals Observed by CLASP2 Across the Mg II h and k Lines

Authors: Ryohko Ishikawa, Javier Trujillo Bueno, Ernest Alsina Ballester, Luca Belluzzi, Tanausú del Pino Alemán, David E. McKenzie, Frédéric Auchère, Ken Kobayashi, Takenori J. Okamoto, Laurel A. Rachmeler, Donguk Song

Abstract: Radiative transfer investigations of the solar Mg II h and k resonance lines around 280~nm showed that, while their circular polarization (Stokes V) signals arise from the Zeeman effect, the linear polarization profiles (Stokes Q and U) are dominated by the scattering of anisotropic radiation and the Hanle and magneto-optical (MO) effects. Using the unprecedented observations of the Mg II and Mn I… ▽ More Radiative transfer investigations of the solar Mg II h and k resonance lines around 280~nm showed that, while their circular polarization (Stokes V) signals arise from the Zeeman effect, the linear polarization profiles (Stokes Q and U) are dominated by the scattering of anisotropic radiation and the Hanle and magneto-optical (MO) effects. Using the unprecedented observations of the Mg II and Mn I resonance lines obtained by the Chromospheric LAyer Spectro-Polarimeter (CLASP2), here we investigate how the linear polarization signals at different wavelengths (i.e., at the center, and at the near and far wings of the k line) vary with the longitudinal component of the magnetic field ($B_{L}$) at their approximate height of formation. The $B_{L}$ is estimated from the V signals in the aforementioned spectral lines. Particular attention is given to the following quantities that are expected to be influenced by the presence of magnetic fields through the Hanle and MO effects: the sign of the U signals, the total linear polarization amplitude ($LP$) and its direction ($χかい$) with respect to a reference direction. We find that at the center and near wings of the $k$ line, the behavior of these quantities is significantly different in the observed quiet and plage regions, and that both $LP$ and $χかい$ seem to depend on $B_{L}$. These observational results are indicative of the operation of the Hanle effect △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 26 pages, 18 figures, accepted for publication in the Astrophysical Journal

Learning to Solve Integer Linear Programs with Davis-Yin Splitting

Authors: Daniel McKenzie, Samy Wu Fung, Howard Heaton

Abstract: In many applications, a combinatorial problem must be repeatedly solved with similar, but distinct parameters. Yet, the parameters $w$ are not directly observed; only contextual data $d$ that correlates with $w$ is available. It is tempting to use a neural network to predict $w$ given $d$. However, training such a model requires reconciling the discrete nature of combinatorial optimization with th… ▽ More In many applications, a combinatorial problem must be repeatedly solved with similar, but distinct parameters. Yet, the parameters $w$ are not directly observed; only contextual data $d$ that correlates with $w$ is available. It is tempting to use a neural network to predict $w$ given $d$. However, training such a model requires reconciling the discrete nature of combinatorial optimization with the gradient-based frameworks used to train neural networks. When the problem in question is an Integer Linear Program (ILP), one approach to overcome this training issue is to consider a continuous relaxation of the combinatorial problem. While existing methods utilizing this approach have shown to be highly effective on small problems, they do not always scale well to large problems. In this work, we draw on ideas from modern convex optimization to design a network and training scheme which scales effortlessly to problems with thousands of variables. Our experiments verify the computational advantage our proposed method enjoys on two representative problems, namely the shortest path problem and the knapsack problem. △ Less

Submitted 21 March, 2024; v1 submitted 30 January, 2023; originally announced January 2023.

Tomography of a solar plage with the Tenerife Inversion Code

Authors: Hao Li, T. del Pino Alemán, J. Trujillo Bueno, R. Ishikawa, E. Alsina Ballester, David E. McKenzie, Frédéric Auchère, Ken Kobayashi, Takenori J. Okamoto, Laurel A. Rachmeler, Donguk Song

Abstract: We apply the Tenerife Inversion Code (TIC) to the plage spectropolarimetric observations obtained by the Chromospheric LAyer SpectroPolarimeter (CLASP2). These unprecedented data consist of full Stokes profiles in the spectral region around the Mg II h and k lines for a single slit position, with around two thirds of the 200 arcsec slit crossing a plage region and the rest crossing an enhanced net… ▽ More We apply the Tenerife Inversion Code (TIC) to the plage spectropolarimetric observations obtained by the Chromospheric LAyer SpectroPolarimeter (CLASP2). These unprecedented data consist of full Stokes profiles in the spectral region around the Mg II h and k lines for a single slit position, with around two thirds of the 200 arcsec slit crossing a plage region and the rest crossing an enhanced network. A former analysis of these data had allowed us to infer the longitudinal component of the magnetic field by applying the weak field approximation (WFA) to the circular polarization profiles, and to assign the inferred magnetic fields to different layers of the solar atmosphere based on the results of previous theoretical radiative transfer investigations. In this work, we apply the recently developed TIC to the same data. We obtain the stratified model atmosphere that fits the intensity and circular polarization profiles at each position along the spectrograph slit and we compare our results for the longitudinal component of the magnetic field with the previously obtained WFA results, highlighting the generally good agreement in spite of the fact that the WFA is known to produce an underestimation when applied to the outer lobes of the Mg II h and k circular polarization profiles. Finally, we use the inverted model atmospheres to give a rough estimation of the energy that could be carried by Alfvèn waves propagating along the chromosphere in the plage and network regions, showing that it is sufficient to compensate the estimated energy losses in the chromosphere of solar active regions. △ Less

Submitted 31 January, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: Accepted for publication in The Astrophysical Journal

The First Flight of the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS)

Authors: Sabrina L. Savage, Amy R. Winebarger, Ken Kobayashi, P. S. Athiray, Dyana Beabout, Leon Golub, Robert W. Walsh, Brent Beabout, Stephen Bradshaw, Alexander R. Bruccoleri, Patrick R. Champey, Peter Cheimets, Jonathan Cirtain, Edward DeLuca, Giulio Del Zanna, Anthony Guillory, Harlan Haight, Ralf K. Heilmann, Edward Hertz, William Hogue, Jeffery Kegley, Jeffery Kolodziejczak, Chad Madsen, Helen Mason, David E. McKenzie , et al. (12 additional authors not shown)

Abstract: The Marshall Grazing Incidence X-ray Spectrometer (MaGIXS) sounding rocket experiment launched on July 30, 2021 from the White Sands Missile Range in New Mexico. MaGIXS is a unique solar observing telescope developed to capture X-ray spectral images, in the 6 - 24 Angstrom wavelength range, of coronal active regions. Its novel design takes advantage of recent technological advances related to fabr… ▽ More The Marshall Grazing Incidence X-ray Spectrometer (MaGIXS) sounding rocket experiment launched on July 30, 2021 from the White Sands Missile Range in New Mexico. MaGIXS is a unique solar observing telescope developed to capture X-ray spectral images, in the 6 - 24 Angstrom wavelength range, of coronal active regions. Its novel design takes advantage of recent technological advances related to fabricating and optimizing X-ray optical systems as well as breakthroughs in inversion methodologies necessary to create spectrally pure maps from overlapping spectral images. MaGIXS is the first instrument of its kind to provide spatially resolved soft X-ray spectra across a wide field of view. The plasma diagnostics available in this spectral regime make this instrument a powerful tool for probing solar coronal heating. This paper presents details from the first MaGIXS flight, the captured observations, the data processing and inversion techniques, and the first science results. △ Less

Submitted 7 December, 2022; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: 20 pages, 18 figures

Adapting Zeroth Order Algorithms for Comparison-Based Optimization

Authors: Isha Slavin, Daniel McKenzie

Abstract: Comparison-Based Optimization (CBO) is an optimization paradigm that assumes only very limited access to the objective function f(x). Despite the growing relevance of CBO to real-world applications, this field has received little attention as compared to the adjacent field of Zeroth-Order Optimization (ZOO). In this work we propose a relatively simple method for converting ZOO algorithms to CBO al… ▽ More Comparison-Based Optimization (CBO) is an optimization paradigm that assumes only very limited access to the objective function f(x). Despite the growing relevance of CBO to real-world applications, this field has received little attention as compared to the adjacent field of Zeroth-Order Optimization (ZOO). In this work we propose a relatively simple method for converting ZOO algorithms to CBO algorithms, thus greatly enlarging the pool of known algorithms for CBO. Via PyCUTEst, we benchmarked these algorithms against a suite of unconstrained problems. We then used hyperparameter tuning to determine optimal values of the parameters of certain algorithms, and utilized visualization tools such as heat maps and line graphs for purposes of interpretation. All our code is available at https://github.com/ishaslavin/Comparison_Based_Optimization. △ Less

Submitted 23 March, 2023; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: Pending review at SIURO

Theory of magnon-polaritons in quantum Ising materials

Authors: R. D. Mckenzie, M. Libersky, D. M. Silevitch, T. F. Rosenbaum

Abstract: We present a theory of magnon-polaritons in quantum Ising materials, and develop a formalism describing the coupling between light and matter as an Ising system is tuned through its quantum critical point. The theory is applied to Ising materials having multilevel single-site Hamiltonians, in which multiple magnon modes are present, such as the insulating Ising magnet LiHoF$_4$ . We find that the… ▽ More We present a theory of magnon-polaritons in quantum Ising materials, and develop a formalism describing the coupling between light and matter as an Ising system is tuned through its quantum critical point. The theory is applied to Ising materials having multilevel single-site Hamiltonians, in which multiple magnon modes are present, such as the insulating Ising magnet LiHoF$_4$ . We find that the magnon-photon coupling strengths may be tuned by the applied transverse field, with the coupling between the soft mode present in the quantum Ising material and a photonic resonator mode diverging at the quantum critical point of the material. A fixed system of spins will not exhibit the diamagnetic response expected when light is coupled to mobile spins or atoms. Without the diamagnetic response, one expects a divergent magnon-photon coupling strength to lead to a superradiant quantum phase transition. However, this neglects the effects of damping and decoherence present in any real system. We show that damping and decoherence may block the superradiant quantum phase transition, and lead to weak coupling between the soft magnon mode and the resonator mode. The results of the theory are applied to experimental data on the model system LiHoF$_4$ in a microwave resonator. △ Less

Submitted 18 September, 2022; originally announced September 2022.

Comments: 22 pages, 6 figures

Quiet Sun Center to Limb Variation of the Linear Polarization Observed by CLASP2 Across the Mg II h & k Lines

Authors: L. A. Rachmeler, J. Trujillo Bueno, D. E. McKenzie, R. Ishikawa, F. Auchere, K. Kobayashi, R. Kano, T. J. Okamoto, C. W. Bethge, D. Song, E. Alsina Ballester, L. Belluzzi, T. del Pino Aleman, M. Yoshida, T. Shimizu, A. Winebarger, A. R. Kobelski, G. D. Vigil, B. De Pontieu, Y. Suematsu, N. Narukage, M. Kubo, T. Sakao, H. Hara, J. Stepan , et al. (3 additional authors not shown)

Abstract: The CLASP2 (Chromospheric LAyer SpectroPolarimeter 2) sounding rocket mission was launched on 2019 April 11. CLASP2 measured the four Stokes parameters of the Mg II h & k spectral region around 2800 Angstroms along a 200 arcsecond slit at three locations on the solar disk, achieving the first spatially and spectrally resolved observations of the solar polarization in this near ultraviolet region.… ▽ More The CLASP2 (Chromospheric LAyer SpectroPolarimeter 2) sounding rocket mission was launched on 2019 April 11. CLASP2 measured the four Stokes parameters of the Mg II h & k spectral region around 2800 Angstroms along a 200 arcsecond slit at three locations on the solar disk, achieving the first spatially and spectrally resolved observations of the solar polarization in this near ultraviolet region. The focus of the work presented here is the center-to-limb variation of the linear polarization across these resonance lines, which is produced by the scattering of anisotropic radiation in the solar atmosphere. The linear polarization signals of the Mg II h & k lines are sensitive to the magnetic field from the low to the upper chromosphere through the Hanle and magneto-optical effects. We compare the observations to theoretical predictions from radiative transfer calculations in unmagnetized semi-empirical models, arguing that magnetic fields and horizontal inhomogeneities are needed to explain the observed polarization signals and spatial variations. This comparison is an important step in both validating and refining our understanding of the physical origin of these polarization signatures, and also in paving the way toward future space telescopes for probing the magnetic fields of the solar upper atmosphere via ultraviolet spectropolarimetry. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Comments: 14 pages, 5 figures, to be published in the Astrophysical Journal (ApJ), initially submitted May 2022, revised submission July 2022

Journal ref: ApJ 936 67 (2022)

Parallel plasma loops and the energization of the solar corona

Authors: Hardi Peter, Lakshmi Pradeep Chitta, Feng Chen, David I. Pontin, Amy R. Winebarger, Leon Golub, Sabrina L. Savage, Laurel A. Rachmeler, Ken Kobayashi, David H. Brooks, Jonathan W. Cirtain, Bart De Pontieu, David E. McKenzie, Richard J. Morton, Paola Testa, Sanjiv K. Tiwari, Robert W. Walsh, Harry P. Warren

Abstract: The outer atmosphere of the Sun is composed of plasma heated to temperatures well in excess of the visible surface. We investigate short cool and warm (<1 MK) loops seen in the core of an active region to address the role of field-line braiding in energising these structures. We report observations from the High-resolution Coronal imager (Hi-C) that have been acquired in a coordinated campaign wit… ▽ More The outer atmosphere of the Sun is composed of plasma heated to temperatures well in excess of the visible surface. We investigate short cool and warm (<1 MK) loops seen in the core of an active region to address the role of field-line braiding in energising these structures. We report observations from the High-resolution Coronal imager (Hi-C) that have been acquired in a coordinated campaign with the Interface Region Imaging Spectrograph (IRIS). In the core of the active region, the 172 A band of Hi-C and the 1400 A channel of IRIS show plasma loops at different temperatures that run in parallel. There is a small but detectable spatial offset of less than 1 arcsec between the loops seen in the two bands. Most importantly, we do not see observational signatures that these loops might be twisted around each other. Considering the scenario of magnetic braiding, our observations of parallel loops imply that the stresses put into the magnetic field have to relax while the braiding is applied: the magnetic field never reaches a highly braided state on these length-scales comparable to the separation of the loops. This supports recent numerical 3D models of loop braiding in which the effective dissipation is sufficiently large that it keeps the magnetic field from getting highly twisted within a loop. △ Less

Submitted 31 May, 2022; originally announced May 2022.

Comments: Accepted for publication in the Astrophysical Journal, 24 pages, 18 figures

From the simplex to the sphere: Faster constrained optimization using the Hadamard parametrization

Authors: Qiuwei Li, Daniel McKenzie, Wotao Yin

Abstract: The standard simplex in R^n, also known as the probability simplex, is the set of nonnegative vectors whose entries sum up to 1. They frequently appear as constraints in optimization problems that arise in machine learning, statistics, data science, operations research, and beyond. We convert the standard simplex to the unit sphere and thus transform the corresponding constrained optimization prob… ▽ More The standard simplex in R^n, also known as the probability simplex, is the set of nonnegative vectors whose entries sum up to 1. They frequently appear as constraints in optimization problems that arise in machine learning, statistics, data science, operations research, and beyond. We convert the standard simplex to the unit sphere and thus transform the corresponding constrained optimization problem into an optimization problem on a simple, smooth manifold. We show that KKT points and strict-saddle points of the minimization problem on the standard simplex all correspond to those of the transformed problem, and vice versa. So, solving one problem is equivalent to solving the other problem. Then, we propose several simple, efficient, and projection-free algorithms using the manifold structure. The equivalence and the proposed algorithm can be extended to optimization problems with unit simplex, weighted probability simplex, or `1-norm sphere constraints. Numerical experiments between the new algorithms and existing ones show the advantages of the new approach △ Less

Submitted 30 August, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: Version 2

Curvature-Aware Derivative-Free Optimization

Authors: Bumsu Kim, HanQin Cai, Daniel McKenzie, Wotao Yin

Abstract: The paper discusses derivative-free optimization (DFO), which involves minimizing a function without access to gradients or directional derivatives, only function evaluations. Classical DFO methods, which mimic gradient-based methods, such as Nelder-Mead and direct search have limited scalability for high-dimensional problems. Zeroth-order methods have been gaining popularity due to the demands of… ▽ More The paper discusses derivative-free optimization (DFO), which involves minimizing a function without access to gradients or directional derivatives, only function evaluations. Classical DFO methods, which mimic gradient-based methods, such as Nelder-Mead and direct search have limited scalability for high-dimensional problems. Zeroth-order methods have been gaining popularity due to the demands of large-scale machine learning applications, and the paper focuses on the selection of the step size $αあるふぁ_k$ in these methods. The proposed approach, called Curvature-Aware Random Search (CARS), uses first- and second-order finite difference approximations to compute a candidate $αあるふぁ_{+}$. We prove that for strongly convex objective functions, CARS converges linearly provided that the search direction is drawn from a distribution satisfying very mild conditions. We also present a Cubic Regularized variant of CARS, named CARS-CR, which converges in a rate of $\mathcal{O}(k^{-1})$ without the assumption of strong convexity. Numerical experiments show that CARS and CARS-CR match or exceed the state-of-the-arts on benchmark problem sets. △ Less

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

Comments: 31 pages, 9 figures

MSC Class: 49M15; 65K05; 68Q25; 90C56

Thermodynamic Evolution of Solar Flare Supra-arcade Downflows

Authors: Z. F. Li, X. Cheng, M. D. Ding, Katharine K. Reeves, DeOndre Kittrell, Mark Weber, David E. McKenzie

Abstract: Solar flares are rapid energy release phenomena that appear as bright ribbons in the chromosphere and high-temperature loops in the corona, respectively. Supra-arcade Downflows (SADs) are plasma voids that first come out above the flare loops and then move quickly towards the flare loop top during the decay phase of the flare. In our work, we study 20 SADs appearing in three flares. By differentia… ▽ More Solar flares are rapid energy release phenomena that appear as bright ribbons in the chromosphere and high-temperature loops in the corona, respectively. Supra-arcade Downflows (SADs) are plasma voids that first come out above the flare loops and then move quickly towards the flare loop top during the decay phase of the flare. In our work, we study 20 SADs appearing in three flares. By differential emission measure (DEM) analysis, we calculate the DEM weighted average temperature and emission measure (EM) of the front region and the main body of SADs. It is found that the temperatures of the SAD front and body tend to increase during the course of SADs flowing downwards. The relationship between the pressure and temperature fits well with the adiabatic equation for both the SAD front and body, suggesting that the heating of SADs is mainly caused by adiabatic compression. Moreover, we also estimate the velocities of SADs via the Fourier Local Correlation Tracking (FLCT) method and find that increase of the temperature of the SAD front presents a correlation with the decrease of the SAD kinetic energy, while the SAD body does not, implying that the viscous process may also heat the SAD front in spite of a limited role. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: 12 pages, 10 figures; published by ApJ

Operator Splitting for Learning to Predict Equilibria in Convex Games

Authors: Daniel McKenzie, Howard Heaton, Qiuwei Li, Samy Wu Fung, Stanley Osher, Wotao Yin

Abstract: Systems of competing agents can often be modeled as games. Assuming rationality, the most likely outcomes are given by an equilibrium (e.g. a Nash equilibrium). In many practical settings, games are influenced by context, i.e. additional data beyond the control of any agent (e.g. weather for traffic and fiscal policy for market economies). Often the exact game mechanics are unknown, yet vast amoun… ▽ More Systems of competing agents can often be modeled as games. Assuming rationality, the most likely outcomes are given by an equilibrium (e.g. a Nash equilibrium). In many practical settings, games are influenced by context, i.e. additional data beyond the control of any agent (e.g. weather for traffic and fiscal policy for market economies). Often the exact game mechanics are unknown, yet vast amounts of historical data consisting of (context, equilibrium) pairs are available, raising the possibility of learning a solver which predicts the equilibria given only the context. We introduce Nash Fixed Point Networks (N-FPNs), a class of neural networks that naturally output equilibria. Crucially, N- FPNs employ a constraint decoupling scheme to handle complicated agent action sets while avoiding expensive projections. Empirically, we find N-FPNs are compatible with the recently developed Jacobian-Free Backpropagation technique for training implicit networks, making them significantly faster and easier to train than prior models. Our experiments show N-FPNs are capable of scaling to problems orders of magnitude larger than existing learned game solvers. △ Less

Submitted 11 June, 2024; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: To appear in SIMODS

JFB: Jacobian-Free Backpropagation for Implicit Networks

Authors: Samy Wu Fung, Howard Heaton, Qiuwei Li, Daniel McKenzie, Stanley Osher, Wotao Yin

Abstract: A promising trend in deep learning replaces traditional feedforward networks with implicit networks. Unlike traditional networks, implicit networks solve a fixed point equation to compute inferences. Solving for the fixed point varies in complexity, depending on provided data and an error tolerance. Importantly, implicit networks may be trained with fixed memory costs in stark contrast to feedforw… ▽ More A promising trend in deep learning replaces traditional feedforward networks with implicit networks. Unlike traditional networks, implicit networks solve a fixed point equation to compute inferences. Solving for the fixed point varies in complexity, depending on provided data and an error tolerance. Importantly, implicit networks may be trained with fixed memory costs in stark contrast to feedforward networks, whose memory requirements scale linearly with depth. However, there is no free lunch -- backpropagation through implicit networks often requires solving a costly Jacobian-based equation arising from the implicit function theorem. We propose Jacobian-Free Backpropagation (JFB), a fixed-memory approach that circumvents the need to solve Jacobian-based equations. JFB makes implicit networks faster to train and significantly easier to implement, without sacrificing test accuracy. Our experiments show implicit networks trained with JFB are competitive with feedforward networks and prior implicit networks given the same number of parameters. △ Less

Submitted 24 December, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Mapping Solar Magnetic Fields from the Photosphere to the Base of the Corona

Authors: Ryohko Ishikawa, Javier Trujillo Bueno, Tanausu del Pino Aleman, Takenori J. Okamoto, David E. McKenzie, Frederic Auchere, Ryouhei Kano, Donguk Song, Masaki Yoshida, Laurel A. Rachmeler, Ken Kobayashi, Hirohisa Hara, Masahito Kubo, Noriyuki Narukage, Taro Sakao, Toshifumi Shimizu, Yoshinori Suematsu, Christian Bethge, Bart De Pontieu, Alberto Sainz Dalda, Genevieve D. Vigil, Amy Winebarger, Ernest Alsina Ballester, Luca Belluzzi, Jiri Stepan , et al. (3 additional authors not shown)

Abstract: Routine ultraviolet imaging of the Sun's upper atmosphere shows the spectacular manifestation of solar activity; yet we remain blind to its main driver, the magnetic field. Here we report unprecedented spectropolarimetric observations of an active region plage and its surrounding enhanced network, showing circular polarization in ultraviolet (Mg II $h$ & $k$ and Mn I) and visible (Fe I) lines. We… ▽ More Routine ultraviolet imaging of the Sun's upper atmosphere shows the spectacular manifestation of solar activity; yet we remain blind to its main driver, the magnetic field. Here we report unprecedented spectropolarimetric observations of an active region plage and its surrounding enhanced network, showing circular polarization in ultraviolet (Mg II $h$ & $k$ and Mn I) and visible (Fe I) lines. We infer the longitudinal magnetic field from the photosphere to the very upper chromosphere. At the top of the plage chromosphere the field strengths reach more than 300 gauss, strongly correlated with the Mg II $k$ line core intensity and the electron pressure. This unique mapping shows how the magnetic field couples the different atmospheric layers and reveals the magnetic origin of the heating in the plage chromosphere. △ Less

Submitted 2 March, 2021; originally announced March 2021.

Comments: 50 pages, 11 figures, 1 table, published in Science Advances

Journal ref: Science Advances, Vol. 7, no. 8, eabe8406, 2021

A Zeroth-Order Block Coordinate Descent Algorithm for Huge-Scale Black-Box Optimization

Authors: HanQin Cai, Yuchen Lou, Daniel McKenzie, Wotao Yin

Abstract: We consider the zeroth-order optimization problem in the huge-scale setting, where the dimension of the problem is so large that performing even basic vector operations on the decision variables is infeasible. In this paper, we propose a novel algorithm, coined ZO-BCD, that exhibits favorable overall query complexity and has a much smaller per-iteration computational complexity. In addition, we di… ▽ More We consider the zeroth-order optimization problem in the huge-scale setting, where the dimension of the problem is so large that performing even basic vector operations on the decision variables is infeasible. In this paper, we propose a novel algorithm, coined ZO-BCD, that exhibits favorable overall query complexity and has a much smaller per-iteration computational complexity. In addition, we discuss how the memory footprint of ZO-BCD can be reduced even further by the clever use of circulant measurement matrices. As an application of our new method, we propose the idea of crafting adversarial attacks on neural network based classifiers in a wavelet domain, which can result in problem dimensions of over 1.7 million. In particular, we show that crafting adversarial examples to audio classifiers in a wavelet domain can achieve the state-of-the-art attack success rate of 97.9%. △ Less

Submitted 11 June, 2021; v1 submitted 21 February, 2021; originally announced February 2021.

Comments: Accepted to ICML 2021

Journal ref: Proceedings of the 38th International Conference on Machine Learning, PMLR 139:1193-1203, 2021

Direct Observation of Collective Electronuclear Modes About a Quantum Critical Point

Authors: M. Libersky, R. D. McKenzie, D. M. Silevitch, P. C. E. Stamp, T. F. Rosenbaum

Abstract: We directly measure the low energy excitation modes of the quantum Ising magnet LiHoF$_4$ using microwave spectroscopy. Instead of a single electronic mode, we find a set of collective electronuclear modes, in which the spin-$1/2$ Ising electronic spins hybridize with the bath of spin-$7/2$ Ho nuclear spins. The lowest-lying electronuclear mode softens at the approach to the quantum critical point… ▽ More We directly measure the low energy excitation modes of the quantum Ising magnet LiHoF$_4$ using microwave spectroscopy. Instead of a single electronic mode, we find a set of collective electronuclear modes, in which the spin-$1/2$ Ising electronic spins hybridize with the bath of spin-$7/2$ Ho nuclear spins. The lowest-lying electronuclear mode softens at the approach to the quantum critical point, even in the presence of disorder. This softening is rapidly quenched by a longitudinal magnetic field. Similar electronuclear structures should exist in other spin-based quantum Ising systems. △ Less

Submitted 13 October, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 12 pages, 5 figures. Accepted for publication in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. vol 127 article 207202 (2021)

Balancing Geometry and Density: Path Distances on High-Dimensional Data

Authors: Anna Little, Daniel McKenzie, James Murphy

Abstract: New geometric and computational analyses of power-weighted shortest-path distances (PWSPDs) are presented. By illuminating the way these metrics balance density and geometry in the underlying data, we clarify their key parameters and discuss how they may be chosen in practice. Comparisons are made with related data-driven metrics, which illustrate the broader role of density in kernel-based unsupe… ▽ More New geometric and computational analyses of power-weighted shortest-path distances (PWSPDs) are presented. By illuminating the way these metrics balance density and geometry in the underlying data, we clarify their key parameters and discuss how they may be chosen in practice. Comparisons are made with related data-driven metrics, which illustrate the broader role of density in kernel-based unsupervised and semi-supervised machine learning. Computationally, we relate PWSPDs on complete weighted graphs to their analogues on weighted nearest neighbor graphs, providing high probability guarantees on their equivalence that are near-optimal. Connections with percolation theory are developed to establish estimates on the bias and variance of PWSPDs in the finite sample setting. The theoretical results are bolstered by illustrative experiments, demonstrating the versatility of PWSPDs for a wide range of data settings. Throughout the paper, our results require only that the underlying data is sampled from a low-dimensional manifold, and depend crucially on the intrinsic dimension of this manifold, rather than its ambient dimension. △ Less

Submitted 7 June, 2021; v1 submitted 16 December, 2020; originally announced December 2020.

MSC Class: 05C85; 05C80 ACM Class: I.5.3

The gray body approximation for radiative heat transfer in evacuated tube solar collectors: Effects of envelope infrared transparency

Authors: Mark A. George, Noboru Takamure, David R. McKenzie

Abstract: A theoretical and experimental analysis is carried out of radiative heat transfer in the coaxial geometry of evacuated tube solar collectors. The gray body approximation implicit in the use of an effective emissivity does not strictly apply to evacuated tube solar collectors due to selective absorber coating and partially transmitting outer glass in the thermal infrared, especially when constructe… ▽ More A theoretical and experimental analysis is carried out of radiative heat transfer in the coaxial geometry of evacuated tube solar collectors. The gray body approximation implicit in the use of an effective emissivity does not strictly apply to evacuated tube solar collectors due to selective absorber coating and partially transmitting outer glass in the thermal infrared, especially when constructed from borosilicate. We develop analytic expressions for the heat transfer through the outer envelope and show the equations no longer follow a simple form where an effective emissivity for the system can be defined. To test all approximations in practice, an experiment is performed using an evacuated solar collector manufactured in the 1980s by the Nitto Kohki company in Japan using the effective emissivity approximation to determine the typical heat transfer characteristics using net radiative heat flows in both directions. This method enabled a good fit to temperature-time data for cooling and heating of the inner tube for temperatures between 10 and 85°C. The results confirm that the effective emittance method can be used in situations with typical glass wall thickness and temperatures. At temperatures greater 100°C, the spectral distribution of the emitted radiation falls significantly within the transmitting region of the outer glass and can no longer be neglected. The work has verified the stability of the vacuum in this type of collector as the tube still functions well, maintaining a low emissivity and good vacuum after approximately 40 years in storage conditions. △ Less

Submitted 25 March, 2022; v1 submitted 24 November, 2020; originally announced November 2020.

Comments: 15 pages, 8 figures

Journal ref: Journal of Applied Physics 131, 125001 (2022)

Who killed Lilly Kane? A case study in applying knowledge graphs to crime fiction

Authors: Mariam Alaverdian, William Gilroy, Veronica Kirgios, Xia Li, Carolina Matuk, Daniel Mckenzie, Tachin Ruangkriengsin, Andrea Bertozzi, Jeffrey Brantingham

Abstract: We present a preliminary study of a knowledge graph created from season one of the television show Veronica Mars, which follows the eponymous young private investigator as she attempts to solve the murder of her best friend Lilly Kane. We discuss various techniques for mining the knowledge graph for clues and potential suspects. We also discuss best practice for collaboratively constructing knowle… ▽ More We present a preliminary study of a knowledge graph created from season one of the television show Veronica Mars, which follows the eponymous young private investigator as she attempts to solve the murder of her best friend Lilly Kane. We discuss various techniques for mining the knowledge graph for clues and potential suspects. We also discuss best practice for collaboratively constructing knowledge graphs from television shows. △ Less

Submitted 23 November, 2020; originally announced November 2020.

A One-bit, Comparison-Based Gradient Estimator

Authors: HanQin Cai, Daniel Mckenzie, Wotao Yin, Zhenliang Zhang

Abstract: We study zeroth-order optimization for convex functions where we further assume that function evaluations are unavailable. Instead, one only has access to a $\textit{comparison oracle}$, which given two points $x$ and $y$ returns a single bit of information indicating which point has larger function value, $f(x)$ or $f(y)$. By treating the gradient as an unknown signal to be recovered, we show how… ▽ More We study zeroth-order optimization for convex functions where we further assume that function evaluations are unavailable. Instead, one only has access to a $\textit{comparison oracle}$, which given two points $x$ and $y$ returns a single bit of information indicating which point has larger function value, $f(x)$ or $f(y)$. By treating the gradient as an unknown signal to be recovered, we show how one can use tools from one-bit compressed sensing to construct a robust and reliable estimator of the normalized gradient. We then propose an algorithm, coined SCOBO, that uses this estimator within a gradient descent scheme. We show that when $f(x)$ has some low dimensional structure that can be exploited, SCOBO outperforms the state-of-the-art in terms of query complexity. Our theoretical claims are verified by extensive numerical experiments. △ Less

Submitted 23 April, 2022; v1 submitted 6 October, 2020; originally announced October 2020.

Journal ref: Applied and Computational Harmonic Analysis, 60 (2022): 242-266

Critical Science Plan for the Daniel K. Inouye Solar Telescope (DKIST)

Authors: Mark P. Rast, Nazaret Bello González, Luis Bellot Rubio, Wenda Cao, Gianna Cauzzi, Edward DeLuca, Bart De Pontieu, Lyndsay Fletcher, Sarah E. Gibson, Philip G. Judge, Yukio Katsukawa, Maria D. Kazachenko, Elena Khomenko, Enrico Landi, Valentin Martínez Pillet, Gordon J. D. Petrie, Jiong Qiu, Laurel A. Rachmeler, Matthias Rempel, Wolfgang Schmidt, Eamon Scullion, Xudong Sun, Brian T. Welsch, Vincenzo Andretta, Patrick Antolin , et al. (62 additional authors not shown)

Abstract: The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With… ▽ More The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With this Critical Science Plan (CSP) we attempt to anticipate some of what those capabilities will enable, providing a snapshot of some of the scientific pursuits that the Daniel K. Inouye Solar Telescope hopes to engage as start-of-operations nears. The work builds on the combined contributions of the DKIST Science Working Group (SWG) and CSP Community members, who generously shared their experiences, plans, knowledge and dreams. Discussion is primarily focused on those issues to which DKIST will uniquely contribute. △ Less

Submitted 20 August, 2020; v1 submitted 18 August, 2020; originally announced August 2020.

The High Inclination Solar Mission

Authors: K. Kobayashi, L. Johnson, H. Thomas, S. McIntosh, D. McKenzie, J. Newmark, A. Heaton, J. Carr, M. Baysinger, Q. Bean, L. Fabisinski, P. Capizzo, K. Clements, S. Sutherlin, J. Garcia, K. Medina, D. Turse

Abstract: The High Inclination Solar Mission (HISM) is a concept for an out-of-the-ecliptic mission for observing the Sun and the heliosphere. The mission profile is largely based on the Solar Polar Imager concept: initially spiraling in to a 0.48 AUえーゆー ecliptic orbit, then increasing the orbital inclination at a rate of $\sim 10$ degrees per year, ultimately reaching a heliographic inclination of $>$75 degree… ▽ More The High Inclination Solar Mission (HISM) is a concept for an out-of-the-ecliptic mission for observing the Sun and the heliosphere. The mission profile is largely based on the Solar Polar Imager concept: initially spiraling in to a 0.48 AUえーゆー ecliptic orbit, then increasing the orbital inclination at a rate of $\sim 10$ degrees per year, ultimately reaching a heliographic inclination of $>$75 degrees. The orbital profile is achieved using solar sails derived from the technology currently being developed for the Solar Cruiser mission, currently under development. HISM remote sensing instruments comprise an imaging spectropolarimeter (Doppler imager / magnetograph) and a visible light coronagraph. The in-situ instruments include a Faraday cup, an ion composition spectrometer, and magnetometers. Plasma wave measurements are made with electrical antennas and high speed magnetometers. The $7,000\,\mathrm{m}^2$ sail used in the mission assessment is a direct extension of the 4-quadrant $1,666\,\mathrm{m}^2$ Solar Cruiser design and employs the same type of high strength composite boom, deployment mechanism, and membrane technology. The sail system modelled is spun (~1 rpm) to assure required boom characteristics with margin. The spacecraft bus features a fine-pointing 3-axis stabilized instrument platform that allows full science observations as soon as the spacecraft reaches a solar distance of 0.48 AUえーゆー. △ Less

Submitted 4 June, 2020; originally announced June 2020.

The drivers of active region outflows into the slow solar wind

Authors: David H. Brooks, Amy R. Winebarger, Sabrina Savage, Harry P. Warren, Bart De Pontieu, Hardi Peter, Jonathan W. Cirtain, Leon Golub, Ken Kobayashi, Scott W. McIntosh, David McKenzie, Richard Morton, Laurel Rachmeler, Paola Testa, Sanjiv Tiwari, Robert Walsh

Abstract: Plasma outflows from the edges of active regions have been suggested as a possible source of the slow solar wind. Spectroscopic measurements show that these outflows have an enhanced elemental composition, which is a distinct signature of the slow wind. Current spectroscopic observations, however, do not have sufficient spatial resolution to distinguish what structures are being measured or to det… ▽ More Plasma outflows from the edges of active regions have been suggested as a possible source of the slow solar wind. Spectroscopic measurements show that these outflows have an enhanced elemental composition, which is a distinct signature of the slow wind. Current spectroscopic observations, however, do not have sufficient spatial resolution to distinguish what structures are being measured or to determine the driver of the outflows. The High-resolution Coronal Imager (Hi-C) flew on a sounding rocket in May, 2018, and observed areas of active region outflow at the highest spatial resolution ever achieved (250 km). Here we use the Hi-C data to disentangle the outflow composition signatures observed with the Hinode satellite during the flight. We show that there are two components to the outflow emission: a substantial contribution from expanded plasma that appears to have been expelled from closed loops in the active region core, and a second contribution from dynamic activity in active region plage, with a composition signature that reflects solar photospheric abundances. The two competing drivers of the outflows may explain the variable composition of the slow solar wind. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: To be published in The Astrophysical Journal. Figures 1, 2, 3 and 12 are reduced resolution to meet size requirements. The original figures will appear in the published version

Zeroth-Order Regularized Optimization (ZORO): Approximately Sparse Gradients and Adaptive Sampling

Authors: HanQin Cai, Daniel Mckenzie, Wotao Yin, Zhenliang Zhang

Abstract: We consider the problem of minimizing a high-dimensional objective function, which may include a regularization term, using (possibly noisy) evaluations of the function. Such optimization is also called derivative-free, zeroth-order, or black-box optimization. We propose a new $\textbf{Z}$eroth-$\textbf{O}$rder $\textbf{R}$egularized $\textbf{O}$ptimization method, dubbed ZORO. When the underlying… ▽ More We consider the problem of minimizing a high-dimensional objective function, which may include a regularization term, using (possibly noisy) evaluations of the function. Such optimization is also called derivative-free, zeroth-order, or black-box optimization. We propose a new $\textbf{Z}$eroth-$\textbf{O}$rder $\textbf{R}$egularized $\textbf{O}$ptimization method, dubbed ZORO. When the underlying gradient is approximately sparse at an iterate, ZORO needs very few objective function evaluations to obtain a new iterate that decreases the objective function. We achieve this with an adaptive, randomized gradient estimator, followed by an inexact proximal-gradient scheme. Under a novel approximately sparse gradient assumption and various different convex settings, we show the (theoretical and empirical) convergence rate of ZORO is only logarithmically dependent on the problem dimension. Numerical experiments show that ZORO outperforms the existing methods with similar assumptions, on both synthetic and real datasets. △ Less

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

Journal ref: SIAM Journal on Optimization 32, no. 2 (2022): 687-714

Is the High-Resolution Coronal Imager Resolving Coronal Strands? Results from AR 12712

Authors: Thomas Williams, Robert W. Walsh, Amy R. Winebarger, David H. Brooks, Jonathan W. Cirtain, Bart Depontieu, Leon Golub, Ken Kobayashi, David E. Mckenzie, Richard J. Morton, Hardi Peter, Laurel A. Rachmeler, Sabrina L. Savage, Paola Testa, Sanjiv K. Tiwari, Harry P. Warren, Benjamin J. Watkinson

Abstract: Following the success of the first mission, the High-Resolution Coronal Imager (Hi-C) was launched for a third time (Hi-C 2.1) on 29th May 2018 from the White Sands Missile Range, NM, USA. On this occasion, 329 seconds of 17.2 nm data of target active region AR 12712 was captured with a cadence of ~4s, and a plate scale of 0.129''/pixel. Using data captured by Hi-C 2.1 and co-aligned observations… ▽ More Following the success of the first mission, the High-Resolution Coronal Imager (Hi-C) was launched for a third time (Hi-C 2.1) on 29th May 2018 from the White Sands Missile Range, NM, USA. On this occasion, 329 seconds of 17.2 nm data of target active region AR 12712 was captured with a cadence of ~4s, and a plate scale of 0.129''/pixel. Using data captured by Hi-C 2.1 and co-aligned observations from SDO/AIA 17.1 nm we investigate the widths of 49 coronal strands. We search for evidence of substructure within the strands that is not detected by AIA, and further consider whether these strands are fully resolved by Hi-C 2.1. With the aid of Multi-Scale Gaussian Normalization (MGN), strands from a region of low-emission that can only be visualized against the contrast of the darker, underlying moss are studied. A comparison is made between these low-emission strands with those from regions of higher emission within the target active region. It is found that Hi-C 2.1 can resolve individual strands as small as ~202km, though more typical strands widths seen are ~513km. For coronal strands within the region of low-emission, the most likely width is significantly narrower than the high-emission strands at ~388km. This places the low-emission coronal strands beneath the resolving capabilities of SDO/AIA, highlighting the need of a permanent solar observatory with the resolving power of Hi-C. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: 19 pages, 10 images, 3 tables

Measuring and modeling the rate of separator reconnection between an emerging and existing active region

Authors: Marika I. McCarthy, Dana W. Longcope, Anna Malanushenko, David E. McKenzie

Abstract: Magnetic reconnection occurs when new flux emerges into the corona and becomes incorporated into the existing coronal field. A new active region (AR) emerging in the vicinity of an existing AR provides a convenient laboratory in which reconnection of this kind can be quantified. We use high time-cadence 171 $Å$ data from SDO/AIA focused on new/old active region pair 11147/11149, to quantify reconn… ▽ More Magnetic reconnection occurs when new flux emerges into the corona and becomes incorporated into the existing coronal field. A new active region (AR) emerging in the vicinity of an existing AR provides a convenient laboratory in which reconnection of this kind can be quantified. We use high time-cadence 171 $Å$ data from SDO/AIA focused on new/old active region pair 11147/11149, to quantify reconnection. We identify new loops as brightenings within a strip of pixels between the regions. This strategy is premised on the assumption that the energy brightening a loop originates in magnetic reconnection. We catalog 301 loops observed in the 48-hour time period beginning with the emergence of AR 11149. The rate at which these loops appear between the two ARs is used to calculate the reconnection rate between them. We then fit these loops with magnetic field, solving for each loop's field strength, geometry, and twist (via its proxy, coronal $αあるふぁ$). We find the rate of newly-brightened flux overestimates the flux which could be undergoing reconnection. This excess can be explained by our finding that the interconnecting region is not at its lowest energy (constant-$αあるふぁ$) state; the extrapolations exhibit loop-to-loop variation in $αあるふぁ$. This flux overestimate may result from the slow emergence of AR 11149, allowing time for Taylor relaxation internal to the domain of the reconnected flux to bring the $αあるふぁ$ distribution towards a single value which provides another mechanism for brightening loops after they are first created. △ Less

Submitted 14 November, 2019; originally announced November 2019.

Comments: Accepted for publication to ApJ

Gravity, topography, and melt generation rates from simple 3D models of mantle convection

Authors: Matthew E. Lees, John F. Rudge, Dan McKenzie

Abstract: Convection in fluid layers at high Rayleigh number (Ra $\sim 10^6$) have a spoke pattern planform. Instabilities in the bottom thermal boundary layer develop into hot rising sheets of fluid, with a component of radial flow towards a central upwelling plume. The sheets form the "spokes" of the pattern, and the plumes the "hubs". Such a pattern of flow is expected to occur beneath plate interiors on… ▽ More Convection in fluid layers at high Rayleigh number (Ra $\sim 10^6$) have a spoke pattern planform. Instabilities in the bottom thermal boundary layer develop into hot rising sheets of fluid, with a component of radial flow towards a central upwelling plume. The sheets form the "spokes" of the pattern, and the plumes the "hubs". Such a pattern of flow is expected to occur beneath plate interiors on Earth, but it remains a challenge to use observations to place constraints on the convective planform of the mantle. Here we present predictions of key surface observables (gravity, topography, and rates of melt generation) from simple 3D numerical models of convection in a fluid layer. These models demonstrate that gravity and topography have only limited sensitivity to the spokes, and mostly reflect the hubs (the rising and sinking plumes). By contrast, patterns of melt generation are more sensitive to short wavelength features in the flow. There is the potential to have melt generation along the spokes, but at a rate which is relatively small compared with that at the hubs. Such melting of spokes can only occur when the lithosphere is sufficiently thin ($\lesssim 80$ km) and mantle water contents are sufficiently high ($\gtrsim 100$ ppm). The distribution of volcanism across the Middle East, Arabia and Africa north of equator suggests that it results from such spoke pattern convection. △ Less

Submitted 5 March, 2020; v1 submitted 7 November, 2019; originally announced November 2019.

Comments: Very minor changes from previous submission

Hi-C 2.1 Observations of Jetlet-like Events at Edges of Solar Magnetic Network Lane

Authors: Navdeep K. Panesar, Alphonse C. Sterling, Ronald L. Moore, Amy R. Winebarger, Sanjiv K. Tiwari, Sabrina L. Savage, Leon Golub, Laurel A. Rachmeler, Ken Kobayashi, David H. Brooks, Jonathan W. Cirtain, Bart De Pontieu, David E. McKenzie, Richard J. Morton, Hardi Peter, Paola Testa, Robert W. Walsh, Harry P. Warren

Abstract: We present high-resolution, high-cadence observations of six, fine-scale, on-disk jet-like events observed by the High-resolution Coronal Imager 2.1 (Hi-C 2.1) during its sounding-rocket flight. We combine the Hi-C 2.1 images with images from SDO/AIA, and IRIS, and investigate each event's magnetic setting with co-aligned line-of-sight magnetograms from SDO/HMI. We find that: (i) all six events ar… ▽ More We present high-resolution, high-cadence observations of six, fine-scale, on-disk jet-like events observed by the High-resolution Coronal Imager 2.1 (Hi-C 2.1) during its sounding-rocket flight. We combine the Hi-C 2.1 images with images from SDO/AIA, and IRIS, and investigate each event's magnetic setting with co-aligned line-of-sight magnetograms from SDO/HMI. We find that: (i) all six events are jetlet-like (having apparent properties of jetlets), (ii) all six are rooted at edges of magnetic network lanes, (iii) four of the jetlet-like events stem from sites of flux cancelation between majority-polarity network flux and merging minority-polarity flux, and (iv) four of the jetlet-like events show brightenings at their bases reminiscent of the base brightenings in coronal jets. The average spire length of the six jetlet-like events (9,000$\pm$3000km) is three times shorter than that for IRIS jetlets (27,000$\pm$8000km). While not ruling out other generation mechanisms, the observations suggest that at least four of these events may be miniature versions of both larger-scale coronal jets that are driven by minifilament eruptions and still-larger-scale solar eruptions that are driven by filament eruptions. Therefore, we propose that our Hi-C events are driven by the eruption of a tiny sheared-field flux rope, and that the flux-rope field is built and triggered to erupt by flux cancelation. △ Less

Submitted 6 November, 2019; originally announced November 2019.

Comments: 10 pages, 5 figures, 1 Table, accepted for publication in ApJ Letters

Fine-scale explosive energy release at sites of prospective magnetic flux cancellation in the core of the solar active region observed by Hi-C 2.1, IRIS and SDO

Authors: Sanjiv K. Tiwari, Navdeep K. Panesar, Ronald L. Moore, Bart De Pontieu, Amy R. Winebarger, Leon Golub, Sabrina L. Savage, Laurel A. Rachmeler, Ken Kobayashi, Paola Testa, Harry P. Warren, David H. Brooks, Jonathan W. Cirtain, David E. McKenzie, Richard J. Morton, Hardi Peter, Robert W. Walsh

Abstract: The second Hi-C flight (Hi-C2.1) provided unprecedentedly-high spatial and temporal resolution ($\sim$250km, 4.4s) coronal EUV images of Fe IX/X emission at 172 Å, of AR 12712 on 29-May-2018, during 18:56:21-19:01:56 UT. Three morphologically-different types (I: dot-like, II: loop-like, III: surge/jet-like) of fine-scale sudden-brightening events (tiny microflares) are seen within and at the ends… ▽ More The second Hi-C flight (Hi-C2.1) provided unprecedentedly-high spatial and temporal resolution ($\sim$250km, 4.4s) coronal EUV images of Fe IX/X emission at 172 Å, of AR 12712 on 29-May-2018, during 18:56:21-19:01:56 UT. Three morphologically-different types (I: dot-like, II: loop-like, III: surge/jet-like) of fine-scale sudden-brightening events (tiny microflares) are seen within and at the ends of an arch filament system in the core of the AR. Although type Is (not reported before) resemble IRIS-bombs (in size, and brightness wrt surroundings), our dot-like events are apparently much hotter, and shorter in span (70s). We complement the 5-minute-duration Hi-C2.1 data with SDO/HMI magnetograms, SDO/AIA EUV images, and IRIS UV spectra and slit-jaw images to examine, at the sites of these events, brightenings and flows in the transition-region and corona and evolution of magnetic flux in the photosphere. Most, if not all, of the events are seated at sites of opposite-polarity magnetic flux convergence (sometimes driven by adjacent flux emergence), implying likely flux cancellation at the microflare's polarity inversion line. In the IRIS spectra and images, we find confirming evidence of field-aligned outflow from brightenings at the ends of loops of the arch filament system. In types I and II the explosion is confined, while in type III the explosion is ejective and drives jet-like outflow. The light-curves from Hi-C, AIA and IRIS peak nearly simultaneously for many of these events and none of the events display a systematic cooling sequence as seen in typical coronal flares, suggesting that these tiny brightening-events have chromospheric/transition-region origin. △ Less

Submitted 4 November, 2019; originally announced November 2019.

Comments: 35 pages, 25 figures, 1 table, accepted for publication in ApJ

Atomic-Scale Patterning of Arsenic in Silicon by Scanning Tunneling Microscopy

Authors: Taylor J. Z. Stock, Oliver Warschkow, Procopios C. Constantinou, Juerong Li, Sarah Fearn, Eleanor Crane, Emily V. S. Hofmann, Alexander Kölker, David R. McKenzie, Steven R. Schofield, Neil J. Curson

Abstract: Over the last two decades, prototype devices for future classical and quantum computing technologies have been fabricated, by using scanning tunneling microscopy and hydrogen resist lithography to position phosphorus atoms in silicon with atomic-scale precision. Despite these successes, phosphine remains the only donor precursor molecule to have been demonstrated as compatible with the hydrogen re… ▽ More Over the last two decades, prototype devices for future classical and quantum computing technologies have been fabricated, by using scanning tunneling microscopy and hydrogen resist lithography to position phosphorus atoms in silicon with atomic-scale precision. Despite these successes, phosphine remains the only donor precursor molecule to have been demonstrated as compatible with the hydrogen resist lithography technique. The potential benefits of atomic-scale placement of alternative dopant species have, until now, remained unexplored. In this work, we demonstrate successful fabrication of atomic-scale structures of arsenic-in-silicon. Using a scanning tunneling microscope tip, we pattern a monolayer hydrogen mask to selectively place arsenic atoms on the Si(001) surface using arsine as the precursor molecule. We fully elucidate the surface chemistry and reaction pathways of arsine on Si(001), revealing significant differences to phosphine. We explain how these differences result in enhanced surface immobilization and in-plane confinement of arsenic compared to phosphorus, and a dose-rate independent arsenic saturation density of $0.24{\pm}0.04$ ML. We demonstrate the successful encapsulation of arsenic delta-layers using silicon molecular beam epitaxy, and find electrical characteristics that are competitive with equivalent structures fabricated with phosphorus. Arsenic delta-layers are also found to offer improvement in out-of-plane confinement compared to similarly prepared phosphorus layers, while still retaining >80% carrier activation and sheet resistances of $<2 kΩおめが/{\square}$. These excellent characteristics of arsenic represent opportunities to enhance existing capabilities of atomic-scale fabrication of dopant structures in silicon, and are particularly important for three-dimensional devices, where vertical control of the position of device components is critical. △ Less

Submitted 15 October, 2019; originally announced October 2019.

The High-Resolution Coronal Imager, Flight 2.1

Authors: Laurel A. Rachmeler, Amy R. Winebarger, Sabrina L. Savage, Leon Golub, Ken Kobayashi, Genevieve D. Vigil, David H. Brooks, Jonathan W. Cirtain, Bart De Pontieu, David E. McKenzie, Richard J. Morton, Hardi Peter, Paola Testa, Sanjiv K. Tiwari, Robert W. Walsh, Harry P. Warren, Caroline Alexander, Darren Ansell, Brent L. Beabout, Dyana L. Beabout, Christian W. Bethge, Patrick R. Champey, Peter N. Cheimets, Mark A. Cooper, Helen K. Creel , et al. (27 additional authors not shown)

Abstract: The third flight of the High-Resolution Coronal Imager (Hi-C 2.1) occurred on May 29, 2018, the Sounding Rocket was launched from White Sands Missile Range in New Mexico. The instrument has been modified from its original configuration (Hi-C 1) to observe the solar corona in a passband that peaks near 172 Angstrom and uses a new, custom-built low-noise camera. The instrument targeted Active Region… ▽ More The third flight of the High-Resolution Coronal Imager (Hi-C 2.1) occurred on May 29, 2018, the Sounding Rocket was launched from White Sands Missile Range in New Mexico. The instrument has been modified from its original configuration (Hi-C 1) to observe the solar corona in a passband that peaks near 172 Angstrom and uses a new, custom-built low-noise camera. The instrument targeted Active Region 12712, and captured 78 images at a cadence of 4.4 sec (18:56:22 - 19:01:57 UT; 5 min and 35 sec observing time). The image spatial resolution varies due to quasi-periodic motion blur from the rocket; sharp images contain resolved features of at least 0.47 arcsec. There are coordinated observations from multiple ground- and space-based telescopes providing an unprecedented opportunity to observe the mass and energy coupling between the chromosphere and the corona. Details of the instrument and the data set are presented in this paper. △ Less

Submitted 12 September, 2019; originally announced September 2019.

Comments: 26 pages, 15 figures, submitted to Solar Physics

Journal ref: Solar Physics, Vol 294, Article number 174, 2019

Power Weighted Shortest Paths for Clustering Euclidean Data

Authors: Daniel Mckenzie, Steven Damelin

Abstract: We study the use of power weighted shortest path distance functions for clustering high dimensional Euclidean data, under the assumption that the data is drawn from a collection of disjoint low dimensional manifolds. We argue, theoretically and experimentally, that this leads to higher clustering accuracy. We also present a fast algorithm for computing these distances. We study the use of power weighted shortest path distance functions for clustering high dimensional Euclidean data, under the assumption that the data is drawn from a collection of disjoint low dimensional manifolds. We argue, theoretically and experimentally, that this leads to higher clustering accuracy. We also present a fast algorithm for computing these distances. △ Less

Submitted 3 September, 2019; v1 submitted 30 May, 2019; originally announced May 2019.

Comments: 24 pages. Final version. To appear in Foundations of Data Science

MSC Class: 05C85; 05C80 ACM Class: I.5.3

On the neighborliness of dual flow polytopes of quivers

Authors: Patricio Gallardo, Daniel Mckenzie

Abstract: In this note we investigate under which conditions the dual of the flow polytope (henceforth referred to as the `dual flow polytope') of a quiver is k-neighborly, for generic weights near the canonical weight. We provide a lower bound on k, depending only on the edge connectivity of the underlying graph, for such weights. In the case where the canonical weight is in fact generic, we explicitly det… ▽ More In this note we investigate under which conditions the dual of the flow polytope (henceforth referred to as the `dual flow polytope') of a quiver is k-neighborly, for generic weights near the canonical weight. We provide a lower bound on k, depending only on the edge connectivity of the underlying graph, for such weights. In the case where the canonical weight is in fact generic, we explicitly determine a vertex presentation of the dual flow polytope. Finally, we specialize our results to the case of complete, bipartite quivers where we show that the canonical weight is indeed generic, and are able to provide an improved bound on k. Hence, we are able to produce many new examples of high-dimensional, k-neighborly polytopes. △ Less

Submitted 5 November, 2018; originally announced November 2018.

Comments: 15 pages, one figure

MSC Class: 14M25; 52B20

Compressive Sensing for cut improvement and local clustering

Authors: Ming-Jun Lai, Daniel Mckenzie

Abstract: We show how one can phrase the cut improvement problem for graphs as a sparse recovery problem, whence one can use algorithms originally developed for use in compressive sensing (such as SubspacePursuit or CoSaMP) to solve it. We show that this approach to cut improvement is fast, both in theory and practice and moreover enjoys statistical guarantees of success when applied to graphs drawn from pr… ▽ More We show how one can phrase the cut improvement problem for graphs as a sparse recovery problem, whence one can use algorithms originally developed for use in compressive sensing (such as SubspacePursuit or CoSaMP) to solve it. We show that this approach to cut improvement is fast, both in theory and practice and moreover enjoys statistical guarantees of success when applied to graphs drawn from probabilistic models such as the Stochastic Block Model. Using this new cut improvement approach, which we call ClusterPursuit, as an algorithmic primitive we then propose new methods for local clustering and semi-supervised clustering, which enjoy similar guarantees of success and speed. Finally, we verify the promise of our approach with extensive numerical benchmarking. △ Less

Submitted 25 February, 2020; v1 submitted 17 August, 2018; originally announced August 2018.

Comments: 25 pages. Generalizes and improves upon the earlier versions arxiv: 1808.05780 and arXiv:1708.09477. To appear in SIMODS

MSC Class: 68Q25; 68R10; 68U05; 94A12

Thermodynamics of a Quantum Ising system coupled to a spin bath: Zero Temperature Results

Authors: R. D. McKenzie, P. C. E Stamp

Abstract: We study the effect of coupling a spin bath environment to a system which, at low energies, can be modeled as a quantum Ising system. A field theoretic formalism incorporating both thermal and quantum fluctuations is developed to derive results for the thermodynamic properties and response functions, both for a toy model and for the $LiHoF_4$ system, in which spin-8 electronic spins couple to a sp… ▽ More We study the effect of coupling a spin bath environment to a system which, at low energies, can be modeled as a quantum Ising system. A field theoretic formalism incorporating both thermal and quantum fluctuations is developed to derive results for the thermodynamic properties and response functions, both for a toy model and for the $LiHoF_4$ system, in which spin-8 electronic spins couple to a spin-$7/2$ nuclear spin bath: the phase transition then occurs in a system of electronuclear degrees of freedom, coupled by long-range dipolar interactions. The quantum Ising phase transition still exists, and one hybridized mode of the Ising and bath spins always goes soft at the transition. △ Less

Submitted 3 June, 2018; v1 submitted 15 November, 2017; originally announced November 2017.

Comments: 18 pages, 9 figures

Journal ref: Phys. Rev. B 97, 214430 (2018)

A Compressive Sensing Approach to Community Detection with Applications

Authors: Ming-Jun Lai, Daniel Mckenzie

Abstract: The community detection problem for graphs asks one to partition the n vertices V of a graph G into k communities, or clusters, such that there are many intracluster edges and few intercluster edges. Of course this is equivalent to finding a permutation matrix P such that, if A denotes the adjacency matrix of G, then PAP^T is approximately block diagonal. As there are k^n possible partitions of n… ▽ More The community detection problem for graphs asks one to partition the n vertices V of a graph G into k communities, or clusters, such that there are many intracluster edges and few intercluster edges. Of course this is equivalent to finding a permutation matrix P such that, if A denotes the adjacency matrix of G, then PAP^T is approximately block diagonal. As there are k^n possible partitions of n vertices into k subsets, directly determining the optimal clustering is clearly infeasible. Instead one seeks to solve a more tractable approximation to the clustering problem. In this paper we reformulate the community detection problem via sparse solution of a linear system associated with the Laplacian of a graph G and then develop a two-stage approach based on a thresholding technique and a compressive sensing algorithm to find a sparse solution which corresponds to the community containing a vertex of interest in G. Crucially, our approach results in an algorithm which is able to find a single cluster of size n_0 in O(nlog(n)n_0) operations and all k clusters in fewer than O(n^2ln(n)) operations. This is a marked improvement over the classic spectral clustering algorithm, which is unable to find a single cluster at a time and takes approximately O(n^3) operations to find all k clusters. Moreover, we are able to provide robust guarantees of success for the case where G is drawn at random from the Stochastic Block Model, a popular model for graphs with clusters. Extensive numerical results are also provided, showing the efficacy of our algorithm on both synthetic and real-world data sets. △ Less

Submitted 20 August, 2018; v1 submitted 30 August, 2017; originally announced August 2017.

Comments: 39 pages, 10 figures Version 2, disabled 'showkeys' package. Note that there is an error in the proof of Lemma 5.1. A correct version of this lemma, as well as a greatly improved version of the central algorithm of this paper, is available at: arXiv:1808.05780

An exploration of heating mechanisms in a supra-arcade plasma sheet formed after a coronal mass ejection

Authors: Katharine K. Reeves, Michael S. Freed, David E. McKenzie, Sabrina L. Savage

Abstract: We perform a detailed analysis of the thermal structure of the region above the post-eruption arcade for a flare that occurred on 2011 October 22. During this event, a sheet of hot plasma is visible above the flare loops in the 131 Å bandpass of the Atmospheric Imaging Assembly (AIA) on the {\it Solar Dynamics Observatory (SDO)}. Supra-arcade downflows (SADs) are observed traveling sunward through… ▽ More We perform a detailed analysis of the thermal structure of the region above the post-eruption arcade for a flare that occurred on 2011 October 22. During this event, a sheet of hot plasma is visible above the flare loops in the 131 Å bandpass of the Atmospheric Imaging Assembly (AIA) on the {\it Solar Dynamics Observatory (SDO)}. Supra-arcade downflows (SADs) are observed traveling sunward through the post-eruption plasma sheet. We calculate differential emission measures using the AIA data and derive an emission measure weighted average temperature in the supra-arcade region. In areas where many SADs occur, the temperature of the supra-arcade plasma tends to increase, while in areas where no SADs are observed, the temperature tends to decrease. We calculate the plane-of-sky velocities in the supra-arcade plasma and use them to calculate the potential heating due to adiabatic compression and viscous heating. Ten of the 13 SADs studied have noticeable signatures in both the adiabatic and the viscous terms. The adiabatic heating due to compression of plasma in front of the SADs is on the order of 0.1 - 0.2 MK/s, which is similar in magnitude to the estimated conductive cooling rate. This result supports the notion that SADs contribute locally to the heating of plasma in the supra-arcade region. We also find that in the region without SADs, the plasma cools at a rate slower than the estimated conductive cooling, indicating additional heating mechanisms may act globally to keep the plasma temperature high. △ Less

Submitted 12 January, 2017; originally announced January 2017.

Comments: 21 pages, 9 figures, accepted in ApJ

Analysis of Flows Inside Quiescent Prominences as Captured by Hinode/Solar Optical Telescope

Authors: M. S. Freed, D. E. McKenzie, D. W. Longcope, M. Wilburn

Abstract: Developing an understanding of how magnetic fields can become entangled in a prominence is important for predicting a possible eruption. This work investigates the kinetic energy and vorticity associated with plasma motion residing inside quiescent prominences (QPs). These plasma flow characteristics can be utilized to improve our understanding of how the prominence maintains a stable magnetic fie… ▽ More Developing an understanding of how magnetic fields can become entangled in a prominence is important for predicting a possible eruption. This work investigates the kinetic energy and vorticity associated with plasma motion residing inside quiescent prominences (QPs). These plasma flow characteristics can be utilized to improve our understanding of how the prominence maintains a stable magnetic field configuration. Three different contrast-enhanced solar prominence observations from Hinode/Solar Optical Telescope were used to construct velocity maps -- in the plane of the sky -- via a Fourier local correlation tracking program. The resulting velocities were then used to perform the first ever analysis of the two-dimensional kinetic energy and enstrophy spectra of a prominence. Enstrophy is introduced here as a means of quantifying the vorticity that has been observed in many QPs. The kinetic energy power spectral density (PSD) produced indices ranging from -1.00 to -1.60. There was a consistent anisotropy in the kinetic energy spectrum of all three prominences examined. Examination of the intensity PSD reveals that a different scaling relationship exists between the observed prominence structure and velocity maps. All of the prominences exhibited an inertial range of at least $0.8 \leq k\leq 2.0\; \textrm{rads} \: \textrm{Mm}^{-1}$. Quasi-periodic oscillations were also detected in the centroid of the velocity distributions for one prominence. Additionally, a lower limit was placed on the kinetic energy density ($εいぷしろん\, \sim 0.22-7.04\: \mathrm{km}^{2}\textrm{s}^{-2}$) and enstrophy density ($ωおめが\, \sim 1.43-13.69\: \times 10^{-16} \, \textrm{s}^{-2}$) associated with each prominence. △ Less

Submitted 11 February, 2016; originally announced February 2016.

Comments: 55 pages, 13 figures, 5 tables, Accepted and Published in ApJ on Feb. 10, 2016, The Astrophysical Journal 2016

Local Spectroscopic Characterization of Spin and Layer Polarization in WSe$_2$

Authors: Matthew Yankowitz, Devin McKenzie, Brian J. LeRoy

Abstract: We report scanning tunneling microscopy (STM) and spectroscopy (STS) measurements of monolayer and bilayer WSe$_2$. We measure a band gap of 2.21 $\pm$ 0.08 eV in monolayer WSe$_2$, which is much larger than the energy of the photoluminescence peak indicating a large excitonic binding energy. We additionally observe significant electronic scattering arising from atomic-scale defects. Using Fourier… ▽ More We report scanning tunneling microscopy (STM) and spectroscopy (STS) measurements of monolayer and bilayer WSe$_2$. We measure a band gap of 2.21 $\pm$ 0.08 eV in monolayer WSe$_2$, which is much larger than the energy of the photoluminescence peak indicating a large excitonic binding energy. We additionally observe significant electronic scattering arising from atomic-scale defects. Using Fourier transform STS (FT-STS), we map the energy versus momentum dispersion relations for monolayer and bilayer WSe$_2$. Further, by tracking allowed and forbidden scattering channels as a function of energy we infer the spin texture of both the conduction and valence bands. We observe a large spin-splitting of the valence band due to strong spin-orbit coupling, and additionally observe spin-valley-layer coupling in the conduction band of bilayer WSe$_2$. △ Less

Submitted 29 July, 2015; v1 submitted 1 May, 2015; originally announced May 2015.

Comments: 5 pages, 4 figures

Intrinsic disorder in graphene on transition metal dichalcogenide heterostructures

Authors: Matthew Yankowitz, Stefano Larentis, Kyounghwan Kim, Jiamin Xue, Devin McKenzie, Shengqiang Huang, Marina Paggen, Mazhar N. Ali, Robert J. Cava, Emanuel Tutuc, Brian J. LeRoy

Abstract: The electronic properties of two-dimensional materials such as graphene are extremely sensitive to their environment, especially the underlying substrate. Planar van der Waals bonded substrates such as hexagonal boron nitride (hBN) have been shown to greatly improve the electrical performance of graphene devices by reducing topographic variations and charge fluctuations compared to amorphous insul… ▽ More The electronic properties of two-dimensional materials such as graphene are extremely sensitive to their environment, especially the underlying substrate. Planar van der Waals bonded substrates such as hexagonal boron nitride (hBN) have been shown to greatly improve the electrical performance of graphene devices by reducing topographic variations and charge fluctuations compared to amorphous insulating substrates}. Semiconducting transition metal dichalchogenides (TMDs) are another family of van der Waals bonded materials that have recently received interest as alternative substrates to hBN for graphene as well as for components in novel graphene-based device heterostructures. Additionally, their semiconducting nature permits dynamic gate voltage control over the interaction strength with graphene. Through local scanning probe measurements we find that crystalline defects intrinsic to TMDs induce scattering in graphene which results in significant degradation of the heterostructure quality, particularly compared to similar graphene on hBN devices. △ Less

Submitted 24 November, 2014; originally announced November 2014.

Comments: 18 pages, 7 figures

Journal ref: Nano Lett. 15, 1925-1929 (2015)

Forward Modeling Transient Brightenings and Microflares around an Active Region Observed with Hi-C

Authors: Adam R. Kobelski, David E. McKenzie

Abstract: Small scale flare-like brightenings around active regions are among the smallest and most fundamental of energetic transient events in the corona, providing a testbed for models of heating and active region dynamics. In a previous study, we modeled a large collection of these microflares observed with Hinode/XRT using EBTEL and found that they required multiple heating events, but could not distin… ▽ More Small scale flare-like brightenings around active regions are among the smallest and most fundamental of energetic transient events in the corona, providing a testbed for models of heating and active region dynamics. In a previous study, we modeled a large collection of these microflares observed with Hinode/XRT using EBTEL and found that they required multiple heating events, but could not distinguish between multiple heating events on a single strand, or multiple strands each experiencing a single heating event. We present here a similar study, but with EUV data of Active Region 11520 from the High Resolution Coronal Imager (Hi-C) sounding rocket. Hi- C provides an order of magnitude improvement to the spatial resolution of XRT, and a cooler temperature sensitivity, which combine to provide significant improvements to our ability to detect and model microflare activity around active regions. We have found that at the spatial resolution of Hi-C (approximately 0.3 arcseconds), the events occur much more frequently than expected (57 events detected, only 1 or 2 expected), and are most likely made from strands of order 100 km wide, each of which is impulsively heated with multiple heating events. These findings tend to support bursty reconnection as the cause of the energy release responsible for the brightenings. △ Less

Submitted 22 August, 2014; originally announced August 2014.

Comments: 13 pages, 6 figures, Accepted in ApJ