Showing 1–50 of 174 results for author: Werner, W

System stabilization with policy optimization on unstable latent manifolds

Authors: Steffen W. R. Werner, Benjamin Peherstorfer

Abstract: Stability is a basic requirement when studying the behavior of dynamical systems. However, stabilizing dynamical systems via reinforcement learning is challenging because only little data can be collected over short time horizons before instabilities are triggered and data become meaningless. This work introduces a reinforcement learning approach that is formulated over latent manifolds of unstabl… ▽ More Stability is a basic requirement when studying the behavior of dynamical systems. However, stabilizing dynamical systems via reinforcement learning is challenging because only little data can be collected over short time horizons before instabilities are triggered and data become meaningless. This work introduces a reinforcement learning approach that is formulated over latent manifolds of unstable dynamics so that stabilizing policies can be trained from few data samples. The unstable manifolds are minimal in the sense that they contain the lowest dimensional dynamics that are necessary for learning policies that guarantee stabilization. This is in stark contrast to generic latent manifolds that aim to approximate all -- stable and unstable -- system dynamics and thus are higher dimensional and often require higher amounts of data. Experiments demonstrate that the proposed approach stabilizes even complex physical systems from few data samples for which other methods that operate either directly in the system state space or on generic latent manifolds fail. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 29 pages, 10 figures, 1 table

MSC Class: 37N35; 68T07; 90C30; 93C57; 93D15

MIRI MRS Observations of Beta Pictoris II. The Spectroscopic Case for a Recent Giant Collision

Authors: Christine H. Chen, Cicero X. Lu, Kadin Worthen, David R. Law, B. A. Sargent, Amaya Moro-Martin, G. C. Sloan, Carey M. Lisse, Dan M. Watson, Julien H. Girard, Yiwei Chai, Dean C. Hines, Jens Kammerer, Alexis Li, Marshall Perrin, Laurent Pueyo, Isabel Rebollido, Karl R. Stapelfeldt, Christopher Stark, Michael W. Werner

Abstract: Modeling observations of the archetypal debris disk around $β$ Pic, obtained in 2023 January with the MIRI MRS on board JWST, reveals significant differences compared with that obtained with the IRS on board Spitzer. The bright 5 - 15 $μ$m continuum excess modeled using a $\sim$600 K black body has disappeared. The previously prominent 18 and 23 $μ$m crystalline forsterite emission features, arisi… ▽ More Modeling observations of the archetypal debris disk around $β$ Pic, obtained in 2023 January with the MIRI MRS on board JWST, reveals significant differences compared with that obtained with the IRS on board Spitzer. The bright 5 - 15 $μ$m continuum excess modeled using a $\sim$600 K black body has disappeared. The previously prominent 18 and 23 $μ$m crystalline forsterite emission features, arising from cold dust ($\sim$100 K) in the Rayleigh limit, have disappeared and been replaced by very weak features arising from the hotter 500 K dust population. Finally, the shape of the 10 $μ$m silicate feature has changed, consistent with a shift in the temperature of the warm dust population from $\sim$300 K to $\sim$500 K and an increase in the crystalline fraction of the warm, silicate dust. Stellar radiation pressure may have blown both the hot and the cold crystalline dust particles observed in the Spitzer spectra out of the planetary system during the intervening 20 years between the Spitzer and JWST observations. These results indicate that the $β$ Pic system has a dynamic circumstellar environment, and that periods of enhanced collisions can create large clouds of dust that sweep through the planetary system. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 15 pages, 8 figures, ApJ in press

Deep polytopic autoencoders for low-dimensional linear parameter-varying approximations and nonlinear feedback design

Authors: Jan Heiland, Yongho Kim, Steffen W. R. Werner

Abstract: Polytopic autoencoders provide low-dimensional parametrizations of states in a polytope. For nonlinear PDEs, this is readily applied to low-dimensional linear parameter-varying (LPV) approximations as they have been exploited for efficient nonlinear controller design via series expansions of the solution to the state-dependent Riccati equation. In this work, we develop a polytopic autoencoder for… ▽ More Polytopic autoencoders provide low-dimensional parametrizations of states in a polytope. For nonlinear PDEs, this is readily applied to low-dimensional linear parameter-varying (LPV) approximations as they have been exploited for efficient nonlinear controller design via series expansions of the solution to the state-dependent Riccati equation. In this work, we develop a polytopic autoencoder for control applications and show how it outperforms standard linear approaches in view of LPV approximations of nonlinear systems and how the particular architecture enables higher order series expansions at little extra computational effort. We illustrate the properties and potentials of this approach to computational nonlinear controller design for large-scale systems with a thorough numerical study. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 9 pages, 6 figures, 2 tables

Interpolatory model order reduction of large-scale dynamical systems with root mean squared error measures

Authors: Sean Reiter, Steffen W. R. Werner

Abstract: The root mean squared error is an important measure used in a variety of applications such as structural dynamics and acoustics to model averaged deviations from standard behavior. For large-scale systems, simulations of this quantity quickly become computationally prohibitive. Classical model order reduction techniques attempt to resolve this issue via the construction of surrogate models that em… ▽ More The root mean squared error is an important measure used in a variety of applications such as structural dynamics and acoustics to model averaged deviations from standard behavior. For large-scale systems, simulations of this quantity quickly become computationally prohibitive. Classical model order reduction techniques attempt to resolve this issue via the construction of surrogate models that emulate the root mean squared error measure using an intermediate linear system. However, this approach requires a potentially large number of linear outputs, which can be disadvantageous in the design of reduced-order models. In this work, we consider directly the root mean squared error as the quantity of interest using the concept of quadratic-output models and propose several new model reduction techniques for the construction of appropriate surrogates. We test the proposed methods on a model for the vibrational response of a plate with tuned vibration absorbers. △ Less

Submitted 11 June, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: 8 pages, 2 figures

Parity questions in critical planar Brownian loop-soups (or "where did the free planar bosons go?")

Authors: Matthis Lehmkuehler, Wei Qian, Wendelin Werner

Abstract: The critical two-dimensional Brownian loop-soup is an infinite collection of non-interacting Brownian loops in a planar domain that possesses some combinatorial features related to the notion of indistinguishability of bosons. The properly renormalized occupation time field of this collection of loops is known to be distributed like the properly defined square of a Gaussian free field. In the pres… ▽ More The critical two-dimensional Brownian loop-soup is an infinite collection of non-interacting Brownian loops in a planar domain that possesses some combinatorial features related to the notion of indistinguishability of bosons. The properly renormalized occupation time field of this collection of loops is known to be distributed like the properly defined square of a Gaussian free field. In the present paper, we investigate aspects of the question about how much information these fields provide about the loop-soup. Among other things, we show that the exact set of points that are actually visited by some loops in the loop-soup is not determined by these fields. We further prove that given the fields, a dense family of special points will each have a conditional probability 1/2 of being part of the loop-soup. We also exhibit another instance where the possible decompositions (given the field) into individual loops and excursions can be grouped into two clearly different groups, each having a conditional probability 1/2 of occurring. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Comments: 29 pages

Using $LDL^{T}$ factorizations in Newton's method for solving general large-scale algebraic Riccati equations

Authors: Jens Saak, Steffen W. R. Werner

Abstract: Continuous-time algebraic Riccati equations can be found in many disciplines in different forms. In the case of small-scale dense coefficient matrices, stabilizing solutions can be computed to all possible formulations of the Riccati equation. This is not the case when it comes to large-scale sparse coefficient matrices. In this paper, we provide a reformulation of the Newton-Kleinman iteration sc… ▽ More Continuous-time algebraic Riccati equations can be found in many disciplines in different forms. In the case of small-scale dense coefficient matrices, stabilizing solutions can be computed to all possible formulations of the Riccati equation. This is not the case when it comes to large-scale sparse coefficient matrices. In this paper, we provide a reformulation of the Newton-Kleinman iteration scheme for continuous-time algebraic Riccati equations using indefinite symmetric low-rank factorizations. This allows the application of the method to the case of general large-scale sparse coefficient matrices. We provide convergence results for several prominent realizations of the equation and show in numerical examples the effectiveness of the approach. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 26 pages, 1 figure, 9 tables

MSC Class: 15A24; 49M15; 65F45; 65H10; 93A15

Energy Dissipation of Fast Electrons in Polymethylmetacrylate (PMMA): Towards a Universal Curve for Electron Beam Attenuation in Solids for Energies between ~0 eV and 100 keV

Authors: Wolfgang S. M. Werner, Florian Simperl, Felix Bloedorn, Julian Brunner, Johannes Kero, Alessandra Bellissimo, Olga Ridzel

Abstract: Wolfgang S.M. Werner, Florian Simperl, Felix Bloedorn, Julian Brunner, Johannes Kero, Alessandra Bellissimo and Olga Ridzel Spectroscopy of correlated electron pairs was employed to investigate the energy dissipation process as well as the transport and the emission of low energy electrons on a polymethylmetracylate (PMMA) surface, providing secondary electron (SE) spectra causally related to th… ▽ More Wolfgang S.M. Werner, Florian Simperl, Felix Bloedorn, Julian Brunner, Johannes Kero, Alessandra Bellissimo and Olga Ridzel Spectroscopy of correlated electron pairs was employed to investigate the energy dissipation process as well as the transport and the emission of low energy electrons on a polymethylmetracylate (PMMA) surface, providing secondary electron (SE) spectra causally related to the energy loss of the primary electron. Two groups of electrons are identified in the cascade of slow electrons, corresponding to different stages in the energy dissipation process. For both groups, the characteristic lengths for attenuation due to collective excitations and momentum relaxation are quantified and are found to be distinctly different: l1=(12.0+/-2) Angstroem and l2=(61.5+/-11) Angstroem. The results strongly contradict the commonly employed model of exponential attenuation with the electron inelastic mean free path (IMFP) as characteristic length, but essentially agree with a theory used for decades in astrophysics and neutron transport, albeit with characteristic lengths expressed in units of Angstroems rather than lightyears. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 6 pages 5 figures

Mass Derivation of planets K2-21b and K2-21c from Transit Timing Variations

Authors: Maryame El Moutamid, Kevin B. Stevenson, Billy Quarles, Nikole K. Lewis, Erik Petigura Daniel Fabrycky, Jacob L. Bean, Diana Dragomir, Kristin S. Sotzenvand Michael W. Werner

Abstract: While various indirect methods are used to detect exoplanets, one of the most effective and accurate methods is the transit method, which measures the brightness of a given star for periodic dips when an exoplanet is passing in front of the parent star. For systems with multiple transiting planets, the gravitational perturbations between planets affect their transit times. The difference in transi… ▽ More While various indirect methods are used to detect exoplanets, one of the most effective and accurate methods is the transit method, which measures the brightness of a given star for periodic dips when an exoplanet is passing in front of the parent star. For systems with multiple transiting planets, the gravitational perturbations between planets affect their transit times. The difference in transit times allows a measurement of the planet masses and orbital eccentricities. These parameters help speculating on the formation, evolution and stability of the system. Using Transit Timing Variations (TTVs), we measure the masses and eccentricities of two planets orbiting K2-21, a relatively bright K7 dwarf star. These two planets exhibit measurable TTVs, have orbital periods of about 9.32 days and 15.50 days, respectively, and a period ratio of about 1.66, which is relatively near to the 5:3 mean motion resonance. We report that the inner and outer planets in the K2-21 system have properties consistent with the presence of a hydrogen and helium dominated atmospheres, as we estimate their masses to be 1.59^{+0.52}_{-0.44} M_E and 3.88^{+1.22}_{-1.07} M_E and densities of 0.22^{+0.05}_{-0.04} rho_E and 0.34^{+0.08}_{-0.06} rho_E, respectively (M_E and rho_E are the mass and density of Earth, respectively). Our results show that the inner planet is less dense than the outer planet; one more counter-intuitive exoplanetary system such as Kepler-105, LTT 1445, TOI-175 and Kepler-279 systems. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Adaptive choice of near-optimal expansion points for interpolation-based structure-preserving model reduction

Authors: Quirin Aumann, Steffen W. R. Werner

Abstract: Interpolation-based methods are well-established and effective approaches for the efficient generation of accurate reduced-order surrogate models. Common challenges for such methods are the automatic selection of good or even optimal interpolation points and the appropriate size of the reduced-order model. An approach that addresses the first problem for linear, unstructured systems is the Iterati… ▽ More Interpolation-based methods are well-established and effective approaches for the efficient generation of accurate reduced-order surrogate models. Common challenges for such methods are the automatic selection of good or even optimal interpolation points and the appropriate size of the reduced-order model. An approach that addresses the first problem for linear, unstructured systems is the Iterative Rational Krylov Algorithm (IRKA), which computes optimal interpolation points through iterative updates by solving linear eigenvalue problems. However, in the case of preserving internal system structures, optimal interpolation points are unknown, and heuristics based on nonlinear eigenvalue problems result in numbers of potential interpolation points that typically exceed the reasonable size of reduced-order systems. In our work, we propose a projection-based iterative interpolation method inspired by IRKA for generally structured systems to adaptively compute near-optimal interpolation points as well as an appropriate size for the reduced-order system. Additionally, the iterative updates of the interpolation points can be chosen such that the reduced-order model provides an accurate approximation in specified frequency ranges of interest. For such applications, our new approach outperforms the established methods in terms of accuracy and computational effort. We show this in numerical examples with different structures. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 24 pages, 5 figures, 4 tables

MSC Class: 30E05; 41A30; 65D05; 93A15; 93C80

Journal ref: Adv. Comput. Math., 50(4):79, 2024

Structured interpolation for multivariate transfer functions of quadratic-bilinear systems

Authors: Peter Benner, Serkan Gugercin, Steffen W. R. Werner

Abstract: High-dimensional/high-fidelity nonlinear dynamical systems appear naturally when the goal is to accurately model real-world phenomena. Many physical properties are thereby encoded in the internal differential structure of these resulting large-scale nonlinear systems. The high-dimensionality of the dynamics causes computational bottlenecks, especially when these large-scale systems need to be simu… ▽ More High-dimensional/high-fidelity nonlinear dynamical systems appear naturally when the goal is to accurately model real-world phenomena. Many physical properties are thereby encoded in the internal differential structure of these resulting large-scale nonlinear systems. The high-dimensionality of the dynamics causes computational bottlenecks, especially when these large-scale systems need to be simulated for a variety of situations such as different forcing terms. This motivates model reduction where the goal is to replace the full-order dynamics with accurate reduced-order surrogates. Interpolation-based model reduction has been proven to be an effective tool for the construction of cheap-to-evaluate surrogate models that preserve the internal structure in the case of weak nonlinearities. In this paper, we consider the construction of multivariate interpolants in frequency domain for structured quadratic-bilinear systems. We propose definitions for structured variants of the symmetric subsystem and generalized transfer functions of quadratic-bilinear systems and provide conditions for structure-preserving interpolation by projection. The theoretical results are illustrated using two numerical examples including the simulation of molecular dynamics in crystal structures. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: 29 pages, 7 figures

MSC Class: 30E05; 34K17; 65D05; 93C10; 93A15

Journal ref: Adv. Comput. Math., 50(2):18, 2024

Structured barycentric forms for interpolation-based data-driven reduced modeling of second-order systems

Authors: Ion Victor Gosea, Serkan Gugercin, Steffen W. R. Werner

Abstract: An essential tool in data-driven modeling of dynamical systems from frequency response measurements is the barycentric form of the underlying rational transfer function. In this work, we propose structured barycentric forms for modeling dynamical systems with second-order time derivatives using their frequency domain input-output data. By imposing a set of interpolation conditions, the systems' tr… ▽ More An essential tool in data-driven modeling of dynamical systems from frequency response measurements is the barycentric form of the underlying rational transfer function. In this work, we propose structured barycentric forms for modeling dynamical systems with second-order time derivatives using their frequency domain input-output data. By imposing a set of interpolation conditions, the systems' transfer functions are rewritten in different barycentric forms using different parametrizations. Loewner-like algorithms are developed for the explicit computation of second-order systems from data based on the developed barycentric forms. Numerical experiments show the performance of these new structured data driven modeling methods compared to other interpolation-based data-driven modeling techniques from the literature. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: 27 pages, 5 figures

MSC Class: 41A20; 65D05; 93B15; 93C05; 93C80

Journal ref: Adv. Comput. Math., 50(2):26, 2024

Low-complexity linear parameter-varying approximations of incompressible Navier-Stokes equations for truncated state-dependent Riccati feedback

Authors: Jan Heiland, Steffen W. R. Werner

Abstract: Nonlinear feedback design via state-dependent Riccati equations is well established but unfeasible for large-scale systems because of computational costs. If the system can be embedded in the class of linear parameter-varying (LPV) systems with the parameter dependency being affine-linear, then the nonlinear feedback law has a series expansion with constant and precomputable coefficients. In this… ▽ More Nonlinear feedback design via state-dependent Riccati equations is well established but unfeasible for large-scale systems because of computational costs. If the system can be embedded in the class of linear parameter-varying (LPV) systems with the parameter dependency being affine-linear, then the nonlinear feedback law has a series expansion with constant and precomputable coefficients. In this work, we propose a general method to approximating nonlinear systems such that the series expansion is possible and efficient even for high-dimensional systems. We lay out the stabilization of incompressible Navier-Stokes equations as application, discuss the numerical solution of the involved matrix-valued equations, and confirm the performance of the approach in a numerical example. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 8 pages, 3 figures

Journal ref: IEEE Control Syst. Lett., 7:3012-3017, 2023

Atomic oxygen abundance toward Sagittarius B2

Authors: Dariusz C. Lis, Paul F. Goldsmith, Rolf Güsten, Peter Schilke, Helmut Wiesemeyer, Youngmin Seo, Michael W. Werner

Abstract: A substantial fraction of oxygen in diffuse clouds is unaccounted for by observations and is postulated to be in an unknown refractory form, referred to as unidentified depleted oxygen (UDO), which, depending on the local gas density, may contribute up to 50% of the total oxygen content. Previous Infrared Space Observatory (ISO) observations suggest that a significant fraction of oxygen in even de… ▽ More A substantial fraction of oxygen in diffuse clouds is unaccounted for by observations and is postulated to be in an unknown refractory form, referred to as unidentified depleted oxygen (UDO), which, depending on the local gas density, may contribute up to 50% of the total oxygen content. Previous Infrared Space Observatory (ISO) observations suggest that a significant fraction of oxygen in even denser, translucent clouds may be in atomic form. We have analyzed velocity-resolved archival SOFIA observations of the 63 $μ$m fine-structure [O I] transition toward the high-mass star-forming region Sgr B2(M) in the Central Molecular Zone. The foreground spiral-arm clouds as well as the extended Sgr B2 envelope between the Sun and the background dust continuum source produce multiple [O i] absorption components, spectrally separated in velocity space. The gas-phase atomic oxygen column density in foreground clouds toward Sgr B2 is well correlated with the total hydrogen column density, with an average atomic oxygen abundance of $(2.51 \pm 0.69) \times 10^{-4}$ with respect to hydrogen nuclei. This value is in good agreement with the earlier ISO measurements on the same line of sight, and is about 35% lower than the total interstellar medium oxygen abundance in the low-density warm gas, as measured in the UV. We find no evidence that a significant fraction of the oxygen on the line of sight toward Sagittarius B2 is in the form of UDO. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: 7 pages, 7 figures, and 2 appendices. Astronomy & Astrophysics Letters, in press

Journal ref: A&A 669, L15 (2023)

Extensions of special 3-fields

Authors: Steven Duplij, Wend Werner

Abstract: We investigate finite field extensions of the unital 3-field, consisting of the unit element alone, and find considerable differences to classical field theory. Furthermore, the structure of their automorphism groups is clarified and the respective subfields are determined. In an attempt to better understand the structure of 3-fields that show up here we look at ways in which new unital 3-fields c… ▽ More We investigate finite field extensions of the unital 3-field, consisting of the unit element alone, and find considerable differences to classical field theory. Furthermore, the structure of their automorphism groups is clarified and the respective subfields are determined. In an attempt to better understand the structure of 3-fields that show up here we look at ways in which new unital 3-fields can be obtained from known ones in terms of product structures, one of them the Cartesian product which has no analogue for binary fields. △ Less

Submitted 15 October, 2022; originally announced December 2022.

Comments: 22 pages, amslatex

MSC Class: 17A40; 20N15

Is the hot, dense sub-Neptune TOI-824b an exposed Neptune mantle? Spitzer detection of the hot day side and reanalysis of the interior composition

Authors: Pierre-Alexis Roy, Björn Benneke, Caroline Piaulet, Ian J. M. Crossfield, Laura Kreidberg, Diana Dragomir, Drake Deming, Michael W. Werner, Vivien Parmentier, Jessie L. Christiansen, Courtney D. Dressing, Stephen R. Kane, Farisa Y. Morales

Abstract: The Kepler and TESS missions revealed a remarkable abundance of sub-Neptune exoplanets. Despite this abundance, our understanding of the nature and compositional diversity of sub-Neptunes remains limited, to a large part because atmospheric studies via transmission spectroscopy almost exclusively aimed for low-density sub-Neptunes and even those were often affected by high-altitude clouds. The rec… ▽ More The Kepler and TESS missions revealed a remarkable abundance of sub-Neptune exoplanets. Despite this abundance, our understanding of the nature and compositional diversity of sub-Neptunes remains limited, to a large part because atmospheric studies via transmission spectroscopy almost exclusively aimed for low-density sub-Neptunes and even those were often affected by high-altitude clouds. The recent TESS discovery of the hot, dense TOI-824b ($2.93\,R_\oplus$ and $18.47\,M_\oplus$) opens a new window into sub-Neptune science by enabling the study of a dense sub-Neptune via secondary eclipses. Here, we present the detection of TOI-824b's hot day side via Spitzer secondary eclipse observations in the $3.6$ and $4.5\,\mathrm{μm}$ channels, combined with a reanalysis of its interior composition. The measured eclipse depths (142$^{+57}_{-52}$ and 245$^{+75}_{-77}$ ppm) and brightness temperatures (1463$^{+183}_{-196}$ and 1484$^{+180}_{-202}$ K) indicate a poor heat redistribution ($f>$ 0.49) and a low Bond albedo (A$_{B}<$ 0.26). We conclude that TOI-824b could be an "exposed Neptune mantle": a planet with a Neptune-like water-rich interior that never accreted a hydrogen envelope or that subsequently lost it. The hot day-side temperature is then naturally explained by a high-metallicity envelope re-emitting the bulk of the incoming radiation from the day side. TOI-824b's density is also consistent with a massive rocky core that accreted up to 1% of hydrogen, but the observed eclipse depths favor our high-metallicity GCM simulation to a solar-metallicity GCM simulation with a likelihood ratio of 7:1. The new insights into TOI-824b's nature suggest that the sub-Neptune population may be more diverse than previously thought, with some of the dense hot sub-Neptunes potentially not hosting a hydrogen-rich envelope as generally assumed for sub-Neptunes. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: Accepted for publication in The Astrophysical Journal, 14 pages, 5 figures

A sub-Neptune transiting the young field star HD 18599 at 40 pc

Authors: Jerome P. de Leon, John H. Livingston, James S. Jenkins, Jose I. Vines, Robert A. Wittenmyer, Jake T. Clark, Joshua I. M. Winn, Brett Addison, Sarah Ballard, Daniel Bayliss, Charles Beichman, Björn Benneke, David Anthony Berardo, Brendan P. Bowler, Tim Brown, Edward M. Bryant, Jessie Christiansen, David Ciardi, Karen A. Collins, Kevin I. Collins, Ian Crossfield, Drake Deming, Diana Dragomir, Courtney D. Dressing, Akihiko Fukui , et al. (45 additional authors not shown)

Abstract: Transiting exoplanets orbiting young nearby stars are ideal laboratories for testing theories of planet formation and evolution. However, to date only a handful of stars with age <1 Gyr have been found to host transiting exoplanets. Here we present the discovery and validation of a sub-Neptune around HD 18599, a young (300 Myr), nearby (d=40 pc) K star. We validate the transiting planet candidate… ▽ More Transiting exoplanets orbiting young nearby stars are ideal laboratories for testing theories of planet formation and evolution. However, to date only a handful of stars with age <1 Gyr have been found to host transiting exoplanets. Here we present the discovery and validation of a sub-Neptune around HD 18599, a young (300 Myr), nearby (d=40 pc) K star. We validate the transiting planet candidate as a bona fide planet using data from the TESS, Spitzer, and Gaia missions, ground-based photometry from IRSF, LCO, PEST, and NGTS, speckle imaging from Gemini, and spectroscopy from CHIRON, NRES, FEROS, and Minerva-Australis. The planet has an orbital period of 4.13 d, and a radius of 2.7Rearth. The RV data yields a 3-sigma mass upper limit of 30.5Mearth which is explained by either a massive companion or the large observed jitter typical for a young star. The brightness of the host star (V~9 mag) makes it conducive to detailed characterization via Doppler mass measurement which will provide a rare view into the interior structure of young planets. △ Less

Submitted 14 October, 2022; originally announced October 2022.

Comments: submitted to MNRAS

Context-aware controller inference for stabilizing dynamical systems from scarce data

Authors: Steffen W. R. Werner, Benjamin Peherstorfer

Abstract: This work introduces a data-driven control approach for stabilizing high-dimensional dynamical systems from scarce data. The proposed context-aware controller inference approach is based on the observation that controllers need to act locally only on the unstable dynamics to stabilize systems. This means it is sufficient to learn the unstable dynamics alone, which are typically confined to much lo… ▽ More This work introduces a data-driven control approach for stabilizing high-dimensional dynamical systems from scarce data. The proposed context-aware controller inference approach is based on the observation that controllers need to act locally only on the unstable dynamics to stabilize systems. This means it is sufficient to learn the unstable dynamics alone, which are typically confined to much lower dimensional spaces than the high-dimensional state spaces of all system dynamics and thus few data samples are sufficient to identify them. Numerical experiments demonstrate that context-aware controller inference learns stabilizing controllers from orders of magnitude fewer data samples than traditional data-driven control techniques and variants of reinforcement learning. The experiments further show that the low data requirements of context-aware controller inference are especially beneficial in data-scarce engineering problems with complex physics, for which learning complete system dynamics is often intractable in terms of data and training costs. △ Less

Submitted 18 January, 2023; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: 27 pages, 10 figures

Journal ref: Proc. R. Soc. A: Math. Phys. Eng. Sci., 479(2270):20220506, 2023

A unifying framework for tangential interpolation of structured bilinear control systems

Authors: Peter Benner, Serkan Gugercin, Steffen W. R. Werner

Abstract: In this paper, we consider the structure-preserving model order reduction problem for multi-input/multi-output bilinear control systems by tangential interpolation. We propose a new type of tangential interpolation problem for structured bilinear systems, for which we develop a new structure-preserving interpolation framework. This new framework extends and generalizes different formulations of ta… ▽ More In this paper, we consider the structure-preserving model order reduction problem for multi-input/multi-output bilinear control systems by tangential interpolation. We propose a new type of tangential interpolation problem for structured bilinear systems, for which we develop a new structure-preserving interpolation framework. This new framework extends and generalizes different formulations of tangential interpolation for bilinear systems from the literature and also provides a unifying framework. We then derive explicit conditions on the projection spaces to enforce tangential interpolation in different settings, including conditions for tangential Hermite interpolation. The analysis is illustrated by means of three numerical examples. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: 33 pages, 9 figures

MSC Class: 30E05; 34K17; 65D05; 93A15; 93C10

Journal ref: Numer. Math., 155(3-4): 445-483, 2023

Multi-fidelity robust controller design with gradient sampling

Authors: Steffen W. R. Werner, Michael L. Overton, Benjamin Peherstorfer

Abstract: Robust controllers that stabilize dynamical systems even under disturbances and noise are often formulated as solutions of nonsmooth, nonconvex optimization problems. While methods such as gradient sampling can handle the nonconvexity and nonsmoothness, the costs of evaluating the objective function may be substantial, making robust control challenging for dynamical systems with high-dimensional s… ▽ More Robust controllers that stabilize dynamical systems even under disturbances and noise are often formulated as solutions of nonsmooth, nonconvex optimization problems. While methods such as gradient sampling can handle the nonconvexity and nonsmoothness, the costs of evaluating the objective function may be substantial, making robust control challenging for dynamical systems with high-dimensional state spaces. In this work, we introduce multi-fidelity variants of gradient sampling that leverage low-cost, low-fidelity models with low-dimensional state spaces for speeding up the optimization process while nonetheless providing convergence guarantees for a high-fidelity model of the system of interest, which is primarily accessed in the last phase of the optimization process. Our first multi-fidelity method initiates gradient sampling on higher fidelity models with starting points obtained from cheaper, lower fidelity models. Our second multi-fidelity method relies on ensembles of gradients that are computed from low- and high-fidelity models. Numerical experiments with controlling the cooling of a steel rail profile and laminar flow in a cylinder wake demonstrate that our new multi-fidelity gradient sampling methods achieve up to two orders of magnitude speedup compared to the single-fidelity gradient sampling method that relies on the high-fidelity model alone. △ Less

Submitted 5 December, 2022; v1 submitted 30 May, 2022; originally announced May 2022.

Comments: 28 pages, 4 figures

MSC Class: 37N35; 37N40; 65K10; 90C30; 90C59

Journal ref: SIAM J. Sci. Comput., 45(2):A933-A957, 2023

Trends in Silicates in the $β$ Pictoris Disk

Authors: Cicero X. Lu, Christine H. Chen, B. A. Sargent, Dan M. Watson, Carey M. Lisse, Joel D. Green, Michael L. Sitko, Tushar Mittal, V. Lebouteiller, G. C. Sloan, Isabel Rebollido, Dean C. Hines, Julien H. Girard, Michael W. Werner, Karl R. Stapelfeldt, Winston Wu, Kadin Worthen

Abstract: While beta Pic is known to host silicates in ring-like structures, whether the properties of these silicate dust vary with stellocentric distance remains an open question. We re-analyze the beta Pictoris debris disk spectrum from the Spitzer Infrared Spectrograph (IRS) and a new IRTF/SpeX spectrum to investigate trends in Fe/Mg ratio, shape, and crystallinity in grains as a function of wavelength,… ▽ More While beta Pic is known to host silicates in ring-like structures, whether the properties of these silicate dust vary with stellocentric distance remains an open question. We re-analyze the beta Pictoris debris disk spectrum from the Spitzer Infrared Spectrograph (IRS) and a new IRTF/SpeX spectrum to investigate trends in Fe/Mg ratio, shape, and crystallinity in grains as a function of wavelength, a proxy for stellocentric distance. By analyzing a re-calibrated and re-extracted spectrum, we identify a new 18 micron forsterite emission feature and recover a 23 micron forsterite emission feature with a substantially larger line-to-continuum ratio than previously reported. We find that these prominent spectral features are primarily produced by small submicron-sized grains, which are continuously generated and replenished from planetesimal collisions in the disk and can elucidate their parent bodies' composition. We discover three trends about these small grains: as stellocentric distance increases, (1) small silicate grains become more crystalline (less amorphous), (2) they become more irregular in shape, and (3) for crystalline silicate grains, the Fe/Mg ratio decreases. Applying these trends to beta Pic's planetary architecture, we find that the dust population exterior to the orbits of beta Pic b and c differs substantially in crystallinity and shape. We also find a tentative 3-5 micron dust excess due to spatially unresolved hot dust emission close to the star. From our findings, we infer that the surfaces of large planetesimals are more Fe-rich and collisionally-processed closer to the star but more Fe-poor and primordial farther from the star. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 19 pages, 12 figures, Accepted for Publication in ApJ

On the sample complexity of stabilizing linear dynamical systems from data

Authors: Steffen W. R. Werner, Benjamin Peherstorfer

Abstract: Learning controllers from data for stabilizing dynamical systems typically follows a two step process of first identifying a model and then constructing a controller based on the identified model. However, learning models means identifying generic descriptions of the dynamics of systems, which can require large amounts of data and extracting information that are unnecessary for the specific task o… ▽ More Learning controllers from data for stabilizing dynamical systems typically follows a two step process of first identifying a model and then constructing a controller based on the identified model. However, learning models means identifying generic descriptions of the dynamics of systems, which can require large amounts of data and extracting information that are unnecessary for the specific task of stabilization. The contribution of this work is to show that if a linear dynamical system has dimension (McMillan degree) $n$, then there always exist $n$ states from which a stabilizing feedback controller can be constructed, independent of the dimension of the representation of the observed states and the number of inputs. By building on previous work, this finding implies that any linear dynamical system can be stabilized from fewer observed states than the minimal number of states required for learning a model of the dynamics. The theoretical findings are demonstrated with numerical experiments that show the stabilization of the flow behind a cylinder from less data than necessary for learning a model. △ Less

Submitted 22 July, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

Comments: 29 pages, 4 figures

MSC Class: 65F55; 65P99; 93B52; 93C57; 93D15

Journal ref: Found. Comput. Math., 24(3):955-987, 2024

On-orbit Performance of the Spitzer Space Telescope: Science Meets Engineering

Authors: Michael W. Werner, Patrick J. Lowrance, Tom Roellig, Varoujan Gorjian, Joseph Hunt, C. Matt Bradford, Jessica Krick

Abstract: The Spitzer Space Telescope operated for over 16 years in an Earth-trailing solar orbit, returning not only a wealth of scientific data but, as a by-product, spacecraft and instrument engineering data which will be of interest to future mission planners. These data will be particularly useful because Spitzer operated in an environment essentially identical to that at the L2 LaGrange point where ma… ▽ More The Spitzer Space Telescope operated for over 16 years in an Earth-trailing solar orbit, returning not only a wealth of scientific data but, as a by-product, spacecraft and instrument engineering data which will be of interest to future mission planners. These data will be particularly useful because Spitzer operated in an environment essentially identical to that at the L2 LaGrange point where many future astrophysics missions will operate. In particular, the radiative cooling demonstrated by Spitzer has been adopted by other infrared space missions, from JWST to SPHEREx. This paper aims to facilitate the utility of the Spitzer engineering data by collecting the more unique and potentially useful portions into a single, readily-accessible publication. We avoid discussion of less unique systems, such as the telecom, flight software, and electronics systems and do not address the innovations in mission and science operations which the Spitzer team initiated. These and other items of potential interest are addressed in references supplied in an appendix to this paper. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: 37 pages, 12 figures, Published in SPIE's Journal of Astronomical Telescopes, Instruments, and Systems (JATIS)

Journal ref: Telesc. Instrum. Syst. 8(1), 014002 (2022)

Structured model order reduction for vibro-acoustic problems using interpolation and balancing methods

Authors: Quirin Aumann, Steffen W. R. Werner

Abstract: Vibration and dissipation in vibro-acoustic systems can be assessed using frequency response analysis. Evaluating a frequency sweep on a full-order model can be very costly, so model order reduction methods are employed to compute cheap-to-evaluate surrogates. This work compares structure-preserving model reduction methods based on rational interpolation and balanced truncation with a specific foc… ▽ More Vibration and dissipation in vibro-acoustic systems can be assessed using frequency response analysis. Evaluating a frequency sweep on a full-order model can be very costly, so model order reduction methods are employed to compute cheap-to-evaluate surrogates. This work compares structure-preserving model reduction methods based on rational interpolation and balanced truncation with a specific focus on their applicability to vibro-acoustic systems. Such models typically exhibit a second-order structure and their material properties as well as their excitation may be depending on the driving frequency. We show and compare the effectiveness of all considered methods by applying them to numerical models of vibro-acoustic systems depicting structural vibration, sound transmission, acoustic scattering, and poroelastic problems. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: 31 pages, 18 figures

Journal ref: J. Sound Vib., 543:117363, 2023

A low-rank solution method for Riccati equations with indefinite quadratic terms

Authors: Peter Benner, Jan Heiland, Steffen W. R. Werner

Abstract: Algebraic Riccati equations with indefinite quadratic terms play an important role in applications related to robust controller design. While there are many established approaches to solve these in case of small-scale dense coefficients, there is no approach available to compute solutions in the large-scale sparse setting. In this paper, we develop an iterative method to compute low-rank approxima… ▽ More Algebraic Riccati equations with indefinite quadratic terms play an important role in applications related to robust controller design. While there are many established approaches to solve these in case of small-scale dense coefficients, there is no approach available to compute solutions in the large-scale sparse setting. In this paper, we develop an iterative method to compute low-rank approximations of stabilizing solutions of large-scale sparse continuous-time algebraic Riccati equations with indefinite quadratic terms. We test the developed approach for dense examples in comparison to other established matrix equation solvers, and investigate the applicability and performance in large-scale sparse examples. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 19 pages, 2 figures, 5 tables

Journal ref: Numer. Algorithms, 92(2):1083-1103, 2023

Structured vector fitting framework for mechanical systems

Authors: Steffen W. R. Werner, Ion Victor Gosea, Serkan Gugercin

Abstract: In this paper, we develop a structure-preserving formulation of the data-driven vector fitting algorithm for the case of modally damped mechanical systems. Using the structured pole-residue form of the transfer function of modally damped second-order systems, we propose two possible structured extensions of the barycentric formula of system transfer functions. Integrating these new forms within th… ▽ More In this paper, we develop a structure-preserving formulation of the data-driven vector fitting algorithm for the case of modally damped mechanical systems. Using the structured pole-residue form of the transfer function of modally damped second-order systems, we propose two possible structured extensions of the barycentric formula of system transfer functions. Integrating these new forms within the classical vector fitting algorithm leads to the formulation of two new algorithms that allow the computation of modally damped mechanical systems from data in a least squares fashion. Thus, the learned model is guaranteed to have the desired structure. We test the proposed algorithms on two benchmark models. △ Less

Submitted 18 October, 2021; originally announced October 2021.

Comments: 8 pages, 2 figures

Journal ref: IFAC-Pap., 55(20):163-168, 2022

TOI-674b: an oasis in the desert of exo-Neptunes transiting a nearby M dwarf

Authors: F. Murgas, N. Astudillo-Defru, X. Bonfils, Ian Crossfield, J. M. Almenara, John Livingston, Keivan G. Stassun, Judith Korth, Jaume Orell-Miquel, G. Morello, Jason D. Eastman, Jack J. Lissauer, Stephen R. Kane, Farisa Y. Morales, Michael W. Werner, Varoujan Gorjian, Björn Benneke, Diana Dragomir, Elisabeth C. Matthews, Steve B. Howell, David Ciardi, Erica Gonzales, Rachel Matson, Charles Beichman, Joshua Schlieder , et al. (37 additional authors not shown)

Abstract: We use TESS, Spitzer, ground-based light curves and HARPS spectrograph radial velocity measurements to establish the physical properties of the transiting exoplanet candidate TOI-674b. We perform a joint fit of the light curves and radial velocity time series to measure the mass, radius, and orbital parameters of the candidate. We confirm and characterize TOI-674b, a low-density super-Neptune tran… ▽ More We use TESS, Spitzer, ground-based light curves and HARPS spectrograph radial velocity measurements to establish the physical properties of the transiting exoplanet candidate TOI-674b. We perform a joint fit of the light curves and radial velocity time series to measure the mass, radius, and orbital parameters of the candidate. We confirm and characterize TOI-674b, a low-density super-Neptune transiting a nearby M dwarf. The host star (TIC 158588995, $V = 14.2$ mag, $J = 10.3$ mag) is characterized by its M2V spectral type with $\mathrm{M}_\star=0.420\pm 0.010$ M$_\odot$, $\mathrm{R}_\star = 0.420\pm 0.013$ R$_\odot$, and $\mathrm{T}_{\mathrm{eff}} = 3514\pm 57$ K, and is located at a distance $d=46.16 \pm 0.03$ pc. Combining the available transit light curves plus radial velocity measurements and jointly fitting a circular orbit model, we find an orbital period of $1.977143 \pm 3\times 10^{-6}$ days, a planetary radius of $5.25 \pm 0.17$ $\mathrm{R}_\oplus$, and a mass of $23.6 \pm 3.3$ $\mathrm{M}_\oplus$ implying a mean density of $ρ_\mathrm{p} = 0.91 \pm 0.15$ [g cm$^{-3}$]. A non-circular orbit model fit delivers similar planetary mass and radius values within the uncertainties. Given the measured planetary radius and mass, TOI-674b is one of the largest and most massive super-Neptune class planets discovered around an M type star to date. It is also a resident of the so-called Neptunian desert and a promising candidate for atmospheric characterisation using the James Webb Space Telescope. △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: 22 pages, 16 figures. Accepted for publication in A&A

Journal ref: A&A 653, A60 (2021)

Robust output-feedback stabilization for incompressible flows using low-dimensional $\mathcal{H}_{\infty}$-controllers

Authors: Peter Benner, Jan Heiland, Steffen W. R. Werner

Abstract: Output-based controllers are known to be fragile with respect to model uncertainties. The standard $\mathcal{H}_{\infty}$-control theory provides a general approach to robust controller design based on the solution of the $\mathcal{H}_{\infty}$-Riccati equations. In view of stabilizing incompressible flows in simulations, two major challenges have to be addressed: the high-dimensional nature of th… ▽ More Output-based controllers are known to be fragile with respect to model uncertainties. The standard $\mathcal{H}_{\infty}$-control theory provides a general approach to robust controller design based on the solution of the $\mathcal{H}_{\infty}$-Riccati equations. In view of stabilizing incompressible flows in simulations, two major challenges have to be addressed: the high-dimensional nature of the spatially discretized model and the differential-algebraic structure that comes with the incompressibility constraint. This work demonstrates the synthesis of low-dimensional robust controllers with guaranteed robustness margins for the stabilization of incompressible flow problems. The performance and the robustness of the reduced-order controller with respect to linearization and model reduction errors are investigated and illustrated in numerical examples. △ Less

Submitted 6 October, 2021; v1 submitted 2 March, 2021; originally announced March 2021.

Comments: 20 pages, 4 figures, 3 tables

Journal ref: Comput. Optim. Appl., 82(1):225-249, 2022

TESS Hunt for Young and Maturing Exoplanets (THYME) IV: Three small planets orbiting a 120 Myr-old star in the Pisces--Eridanus stream

Authors: Elisabeth R. Newton, Andrew W. Mann, Adam L. Kraus, John H. Livingston, Andrew Vanderburg, Jason L. Curtis, Pa Chia Thao, Keith Hawkins, Mackenna L. Wood, Aaron C. Rizzuto, Abderahmane Soubkiou, Benjamin M. Tofflemire, George Zhou, Ian J. M. Crossfield, Logan A. Pearce, Karen A. Collins, Dennis M. Conti, Thiam-Guan Tan, Steven Villeneuva, Alton Spencer, Diana Dragomir, Samuel N. Quinn, Eric L. N. Jensen, Kevin I. Collins, Chris Stockdale , et al. (28 additional authors not shown)

Abstract: Young exoplanets can offer insight into the evolution of planetary atmospheres, compositions, and architectures. We present the discovery of the young planetary system TOI 451 (TIC 257605131, Gaia DR2 4844691297067063424). TOI 451 is a member of the 120-Myr-old Pisces--Eridanus stream (Psc--Eri). We confirm membership in the stream with its kinematics, its lithium abundance, and the rotation and U… ▽ More Young exoplanets can offer insight into the evolution of planetary atmospheres, compositions, and architectures. We present the discovery of the young planetary system TOI 451 (TIC 257605131, Gaia DR2 4844691297067063424). TOI 451 is a member of the 120-Myr-old Pisces--Eridanus stream (Psc--Eri). We confirm membership in the stream with its kinematics, its lithium abundance, and the rotation and UV excesses of both TOI 451 and its wide binary companion, TOI 451 B (itself likely an M dwarf binary). We identified three candidate planets transiting in the TESS data and followed up the signals with photometry from Spitzer and ground-based telescopes. The system comprises three validated planets at periods of 1.9, 9.2 and 16 days, with radii of 1.9, 3.1, and 4.1 Earth radii, respectively. The host star is near-solar mass with V=11.0 and H=9.3 and displays an infrared excess indicative of a debris disk. The planets offer excellent prospects for transmission spectroscopy with HST and JWST, providing the opportunity to study planetary atmospheres that may still be in the process of evolving. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: 20 pages, appendix on UV excess

Journal ref: The Astronomical Journal, 2021, Volume 161, Issue 2

Structure-Preserving Model Reduction for Dissipative Mechanical Systems

Authors: Rebekka S. Beddig, Peter Benner, Ines Dorschky, Timo Reis, Paul Schwerdtner, Matthias Voigt, Steffen W. R. Werner

Abstract: Suppressing vibrations in mechanical systems, usually described by second-order dynamical models, is a challenging task in mechanical engineering in terms of computational resources even nowadays. One remedy is structure-preserving model order reduction to construct easy-to-evaluate surrogates for the original dynamical system having the same structure. In our work, we present an overview of recen… ▽ More Suppressing vibrations in mechanical systems, usually described by second-order dynamical models, is a challenging task in mechanical engineering in terms of computational resources even nowadays. One remedy is structure-preserving model order reduction to construct easy-to-evaluate surrogates for the original dynamical system having the same structure. In our work, we present an overview of recently developed structure-preserving model reduction methods for second-order systems. These methods are based on modal and balanced truncation in different variants, as well as on rational interpolation. Numerical examples are used to illustrate the effectiveness of all described methods. △ Less

Submitted 25 January, 2022; v1 submitted 13 October, 2020; originally announced October 2020.

Comments: 21 pages, 8 figures

Journal ref: Lect. Notes Appl. Comput. Mech., 102:209-230, 2024

Structure-Preserving Interpolation for Model Reduction of Parametric Bilinear Systems

Authors: Peter Benner, Serkan Gugercin, Steffen W. R. Werner

Abstract: In this paper, we present an interpolation framework for structure-preserving model order reduction of parametric bilinear dynamical systems. We introduce a general setting, covering a broad variety of different structures for parametric bilinear systems, and then provide conditions on projection spaces for the interpolation of structured subsystem transfer functions such that the system structure… ▽ More In this paper, we present an interpolation framework for structure-preserving model order reduction of parametric bilinear dynamical systems. We introduce a general setting, covering a broad variety of different structures for parametric bilinear systems, and then provide conditions on projection spaces for the interpolation of structured subsystem transfer functions such that the system structure and parameter dependencies are preserved in the reduced-order model. Two benchmark examples with different parameter dependencies are used to demonstrate the theoretical analysis. △ Less

Submitted 22 July, 2020; originally announced July 2020.

Comments: 10 pages, 2 figures

Journal ref: Automatica J. IFAC, 132:109799, 2021

Non-simple conformal loop ensembles on Liouville quantum gravity and the law of CLE percolation interfaces

Authors: Jason Miller, Scott Sheffield, Wendelin Werner

Abstract: We study the structure of the Liouville quantum gravity (LQG) surfaces that are cut out as one explores a conformal loop-ensemble CLE$_{κ'}$ for $κ'$ in $(4,8)$ that is drawn on an independent $γ$-LQG surface for $γ^2=16/κ'$. The results are similar in flavor to the ones from our paper dealing with CLE$_κ$ for $κ$ in $(8/3,4)$, where the loops of the CLE are disjoint and simple. In particular, we… ▽ More We study the structure of the Liouville quantum gravity (LQG) surfaces that are cut out as one explores a conformal loop-ensemble CLE$_{κ'}$ for $κ'$ in $(4,8)$ that is drawn on an independent $γ$-LQG surface for $γ^2=16/κ'$. The results are similar in flavor to the ones from our paper dealing with CLE$_κ$ for $κ$ in $(8/3,4)$, where the loops of the CLE are disjoint and simple. In particular, we encode the combined structure of the LQG surface and the CLE$_{κ'}$ in terms of stable growth-fragmentation trees or their variants, which also appear in the asymptotic study of peeling processes on decorated planar maps. This has consequences for questions that do a priori not involve LQG surfaces: Our previous paper "CLE percolations" described the law of interfaces obtained when coloring the loops of a CLE$_{κ'}$ independently into two colors with respective probabilities $p$ and $1-p$. This description was complete up to one missing parameter $ρ$. The results of the present paper about CLE on LQG allow us to determine its value in terms of $p$ and $κ'$. It shows in particular that CLE$_{κ'}$ and CLE$_{16/κ'}$ are related via a continuum analog of the Edwards-Sokal coupling between FK$_q$ percolation and the $q$-state Potts model (which makes sense even for non-integer $q$ between $1$ and $4$) if and only if $q=4\cos^2(4π/κ')$. This provides further evidence for the long-standing belief that CLE$_{κ'}$ and CLE$_{16/κ'}$ represent the scaling limits of FK$_q$ percolation and the $q$-Potts model when $q$ and $κ'$ are related in this way. Another consequence of the formula for $ρ(p,κ')$ is the value of half-plane arm exponents for such divide-and-color models (a.k.a. fuzzy Potts models) that turn out to take a somewhat different form than the usual critical exponents for two-dimensional models. △ Less

Submitted 28 May, 2021; v1 submitted 25 June, 2020; originally announced June 2020.

Comments: Dedicated to the memory of Harry Kesten. To appear in Probab. Theory Rel. Fields, 30 pages, 18 figures

Secondary Electron Emission by Plasmon Induced Symmetry Breaking in Highly Oriented Pyrolitic Graphite (HOPG)

Authors: Wolfgang S. M. Werner, Vytautas Astašauskas, Philipp Ziegler, Alessandra Bellissimo, Giovanni Stefani, Lukas Linhart, Florian Libisch

Abstract: Two-particle spectroscopy with correlated electron pairs is used to establish the causal link between the secondary electron spectrum, the $(π+σ)-$plasmon peak and the unoccupied band structure of highly oriented pyrolitic graphite. The plasmon spectrum is resolved with respect to the involved interband transitions and clearly exhibits final state effects, in particular due to the energy gap betwe… ▽ More Two-particle spectroscopy with correlated electron pairs is used to establish the causal link between the secondary electron spectrum, the $(π+σ)-$plasmon peak and the unoccupied band structure of highly oriented pyrolitic graphite. The plasmon spectrum is resolved with respect to the involved interband transitions and clearly exhibits final state effects, in particular due to the energy gap between the interlayer resonances along the $Γ$A-direction. The corresponding final state effects can also be identified in the secondary electron spectrum. Interpretation of the results is performed on the basis of density functional theory and tight binding calculations. Excitation of the plasmon perturbs the symmetry of the system and leads to hybridisation of the interlayer resonances with atom-like $σ^*$ bands along the $ΓA$-direction. These hybrid states have a high density of states as well as sufficient mobility along the graphite $c$-axis leading to the sharp $\sim$3\ eV resonance in the spectrum of emitted secondary electrons reported throughout the literature. △ Less

Submitted 9 June, 2020; originally announced June 2020.

Comments: 9 pages, 4 figures

Journal ref: Phys. Rev. Lett. 125, 196603 (2020)

Structure-Preserving Interpolation of Bilinear Control Systems

Authors: Peter Benner, Serkan Gugercin, Steffen W. R. Werner

Abstract: In this paper, we extend the structure-preserving interpolatory model reduction framework, originally developed for linear systems, to structured bilinear control systems. Specifically, we give explicit construction formulae for the model reduction bases to satisfy different types of interpolation conditions. First, we establish the analysis for transfer function interpolation for single-input sin… ▽ More In this paper, we extend the structure-preserving interpolatory model reduction framework, originally developed for linear systems, to structured bilinear control systems. Specifically, we give explicit construction formulae for the model reduction bases to satisfy different types of interpolation conditions. First, we establish the analysis for transfer function interpolation for single-input single-output structured bilinear systems. Then, we extend these results to the case of multi-input multi-output structured bilinear systems by matrix interpolation. The effectiveness of our structure-preserving approach is illustrated by means of various numerical examples. △ Less

Submitted 2 May, 2020; originally announced May 2020.

Comments: 26 pages, 9 figures

Journal ref: Adv. Comput. Math., 47(3):43, 2021

Lecture notes on the Gaussian Free Field

Authors: Wendelin Werner, Ellen Powell

Abstract: The Gaussian Free Field (GFF) in the continuum appears to be the natural generalisation of Brownian motion, when one replaces time by a multidimensional continuous parameter. The goal of these lecture notes is to describe some aspects of the continuum GFF and of its discrete counterpart defined on lattices, with the aim of providing a gentle self-contained introduction to some recent developments… ▽ More The Gaussian Free Field (GFF) in the continuum appears to be the natural generalisation of Brownian motion, when one replaces time by a multidimensional continuous parameter. The goal of these lecture notes is to describe some aspects of the continuum GFF and of its discrete counterpart defined on lattices, with the aim of providing a gentle self-contained introduction to some recent developments on this topic, such as the relation between the continuum GFF, Brownian loop-soups and the Conformal Loop Ensembles CLE(4). This is an updated and expanded version of the notes written by the first author for graduate courses at ETH Zürich. The exercises that are interspersed in the first half of these notes mostly originate from the exercise sheets prepared by the second author for this course in 2018. △ Less

Submitted 22 April, 2021; v1 submitted 9 April, 2020; originally announced April 2020.

Comments: Final version, 145 pages, to appear in the "Cours Spécialisés" series published by the Société Mathématique de France

MORLAB -- The Model Order Reduction LABoratory

Authors: Peter Benner, Steffen W. R. Werner

Abstract: For an easy use of model order reduction techniques in applications, software solutions are needed. In this paper, we describe the MORLAB, Model Order Reduction LABoratory, toolbox as an efficient implementation of model reduction techniques for dense, medium-scale linear time-invariant systems. Giving an introduction to the underlying programming principles of the toolbox, we show the basic idea… ▽ More For an easy use of model order reduction techniques in applications, software solutions are needed. In this paper, we describe the MORLAB, Model Order Reduction LABoratory, toolbox as an efficient implementation of model reduction techniques for dense, medium-scale linear time-invariant systems. Giving an introduction to the underlying programming principles of the toolbox, we show the basic idea of spectral splitting and present an overview about implemented model reduction techniques. Two numerical examples are used to illustrate different use cases of the MORLAB toolbox. △ Less

Submitted 28 February, 2020; originally announced February 2020.

Comments: 17 pages, 6 figures, 5 tables

Journal ref: International Series of Numerical Mathematics, 171:393-415, 2021

On clusters of Brownian loops in d dimensions

Authors: Wendelin Werner

Abstract: We discuss random geometric structures obtained by percolation of Brownian loops, in relation to the Gaussian Free Field, and how their existence and properties depend on the dimension of the ambient space. We formulate a number of conjectures for the cases d=3,4,5 and prove some results when d > 6. We discuss random geometric structures obtained by percolation of Brownian loops, in relation to the Gaussian Free Field, and how their existence and properties depend on the dimension of the ambient space. We formulate a number of conjectures for the cases d=3,4,5 and prove some results when d > 6. △ Less

Submitted 26 February, 2020; originally announced February 2020.

Simple Conformal Loop Ensembles on Liouville Quantum Gravity

Authors: Jason Miller, Scott Sheffield, Wendelin Werner

Abstract: We show that when one draws a simple conformal loop ensemble (CLE$_κ$ for $κ\in (8/3,4)$) on an independent $\sqrtκ$-Liouville quantum gravity (LQG) surface and explores the CLE in a natural Markovian way, the quantum surfaces (e.g., corresponding to the interior of the CLE loops) that are cut out form a Poisson point process of quantum disks. This construction allows us to make direct links betwe… ▽ More We show that when one draws a simple conformal loop ensemble (CLE$_κ$ for $κ\in (8/3,4)$) on an independent $\sqrtκ$-Liouville quantum gravity (LQG) surface and explores the CLE in a natural Markovian way, the quantum surfaces (e.g., corresponding to the interior of the CLE loops) that are cut out form a Poisson point process of quantum disks. This construction allows us to make direct links between CLE on LQG, asymmetric $(4/κ)$-stable processes, and labeled branching trees. The ratio between positive and negative jump intensities of these processes turns out to be $-\cos (4 π/ κ)$, which can be interpreted as a "density" of CLE loops in the CLE on LQG setting. Positive jumps correspond to the discovery of a CLE loop (where the LQG length of the loop is given by the jump size) and negative jumps correspond to the moments where the discovery process splits the remaining to be discovered domain into two pieces. Some consequences are the following: (i) It provides a construction of a CLE on LQG as a patchwork/welding of quantum disks. (ii) It allows to construct the "natural quantum measure" that lives in a CLE carpet. (iii) It enables us to derive some new properties and formulas for SLE processes and CLE themselves (without LQG) such as the exact distribution of the trunk of the general asymmetric SLE$_κ(κ-6)$ processes. The present work deals directly with structures in the continuum and makes no reference to discrete models, but our calculations match those for scaling limits of O(N) models on planar maps with large faces and CLE on LQG. Indeed, our Lévy-tree descriptions are exactly the ones that appear in the study of the large-scale limit of peeling of discrete decorated planar maps such as in recent work of Bertoin, Budd, Curien and Kortchemski. The case of non-simple CLEs on LQG is the topic of another paper. △ Less

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

Comments: 46 pages, final version to appear in Ann. Probab

Frequency- and Time-Limited Balanced Truncation for Large-Scale Second-Order Systems

Authors: Peter Benner, Steffen W. R. Werner

Abstract: Considering the use of dynamical systems in practical applications, often only limited regions in the time or frequency domain are of interest. Therefor, it usually pays off to compute local approximations of the used dynamical systems in the frequency and time domain. In this paper, we consider a structure-preserving extension of the frequency- and time-limited balanced truncation methods to seco… ▽ More Considering the use of dynamical systems in practical applications, often only limited regions in the time or frequency domain are of interest. Therefor, it usually pays off to compute local approximations of the used dynamical systems in the frequency and time domain. In this paper, we consider a structure-preserving extension of the frequency- and time-limited balanced truncation methods to second-order dynamical systems. We give a full overview about the first-order limited balanced truncation methods and extend those methods to second-order systems by using the different second-order balanced truncation formulas from the literature. Also, we present numerical methods for solving the arising large-scale sparse matrix equations and give numerical modifications to deal with the problematic case of second-order systems. The results are then illustrated on three numerical examples. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 30 pages, 10 figures, 9 tables

Journal ref: Linear Algebra Appl., 623:68-103, 2021

Spitzer transit follow-up of planet candidates from the K2 mission

Authors: John H. Livingston, Ian J. M. Crossfield, Michael W. Werner, Varoujan Gorjian, Erik A. Petigura, David R. Ciardi, Courtney D. Dressing, Benjamin J. Fulton, Teruyuki Hirano, Joshua E. Schlieder, Evan Sinukoff, Molly Kosiarek, Rachel Akeson, Charles A. Beichman, Björn Benneke, Jessie L. Christiansen, Bradley M. S. Hansen, Andrew W. Howard, Howard Isaacson, Heather A. Knutson, Jessica Krick, Arturo O. Martinez, Bun'ei Sato, Motohide Tamura

Abstract: We present precision 4.5 $μ$m Spitzer transit photometry of eight planet candidates discovered by the K2 mission: K2-52 b, K2-53 b, EPIC 205084841.01, K2-289 b, K2-174 b, K2-87 b, K2-90 b, and K2-124 b. The sample includes four sub-Neptunes and two sub-Saturns, with radii between 2.6 and 18 $R_\oplus$, and equilibrium temperatures between 440 and 2000 K. In this paper we identify several targets o… ▽ More We present precision 4.5 $μ$m Spitzer transit photometry of eight planet candidates discovered by the K2 mission: K2-52 b, K2-53 b, EPIC 205084841.01, K2-289 b, K2-174 b, K2-87 b, K2-90 b, and K2-124 b. The sample includes four sub-Neptunes and two sub-Saturns, with radii between 2.6 and 18 $R_\oplus$, and equilibrium temperatures between 440 and 2000 K. In this paper we identify several targets of potential interest for future characterization studies, demonstrate the utility of transit follow-up observations for planet validation and ephemeris refinement, and present new imaging and spectroscopy data. Our simultaneous analysis of the K2 and Spitzer light curves yields improved estimates of the planet radii, and multi-wavelength information which help validate their planetary nature, including the previously un-validated candidate EPIC 205686202.01 (K2-289 b). Our Spitzer observations yield an order of magnitude increase in ephemeris precision, thus paving the way for efficient future study of these interesting systems by reducing the typical transit timing uncertainty in mid-2021 from several hours to a dozen or so minutes. K2-53 b, K2-289 b, K2-174 b, K2-87 b, and K2-90 b are promising radial velocity (RV) targets given the performance of spectrographs available today or in development, and the M3V star K2-124 hosts a temperate sub-Neptune that is potentially a good target for both RV and atmospheric characterization studies. △ Less

Submitted 17 January, 2019; originally announced January 2019.

Comments: 21 pages, 8 figures, 9 tables, accepted for publication in AJ

Science Impacts of the SPHEREx All-Sky Optical to Near-Infrared Spectral Survey II: Report of a Community Workshop on the Scientific Synergies Between the SPHEREx Survey and Other Astronomy Observatories

Authors: Olivier Doré, Michael W. Werner, Matthew L. N. Ashby, Lindsey E. Bleem, Jamie Bock, Jennifer Burt, Peter Capak, Tzu-Ching Chang, Jonás Chaves-Montero, Christine H. Chen, Francesca Civano, I. Ilsedore Cleeves, Asantha Cooray, Brendan Crill, Ian J. M. Crossfield, Michael Cushing, Sylvain de la Torre, Tiziana DiMatteo, Niv Dvory, Cora Dvorkin, Catherine Espaillat, Simone Ferraro, Douglas Finkbeiner, Jenny Greene, Jackie Hewitt , et al. (38 additional authors not shown)

Abstract: SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$μ$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$σ$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx… ▽ More SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$μ$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$σ$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop. △ Less

Submitted 24 May, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: 50 pages, 24 figures, more details at http://spherex.caltech.edu

The law of a point process of Brownian excursions in a domain is determined by the law of its trace

Authors: Wei Qian, Wendelin Werner

Abstract: We show the result that is stated in the title of the paper, which has consequences about decomposition of Brownian loop-soup clusters in two dimensions. We show the result that is stated in the title of the paper, which has consequences about decomposition of Brownian loop-soup clusters in two dimensions. △ Less

Submitted 14 December, 2018; v1 submitted 7 February, 2018; originally announced February 2018.

Comments: 24 pages, to appear in Electron.J.Probab

Journal ref: Electronic J. Probab., Volume 23, paper no. 128, 23 pp., 2018

Progress on Experiments towards LWFA-driven Transverse Gradient Undulator-Based FELs

Authors: Axel Bernhard, Veronica Afonso Rodriguez, Stephan Kuschel, Maria Leier, Peter Peiffer, Alexander Saevert, Matthew Schwab, Walter Werner, Christina Widmann, Andreas Will, Anke-Susanne Mueller, Malte Kaluza

Abstract: Free Electron Lasers (FEL) are commonly regarded as the potential key application of laser wakefield accelerators (LWFA). It has been found that electron bunches exiting from state-of-the-art LWFAs exhibit a normalized 6-dimensional beam brightness comparable to those in conventional linear accelerators. Effectively exploiting this beneficial beam property for LWFA-based FELs is challenging due to… ▽ More Free Electron Lasers (FEL) are commonly regarded as the potential key application of laser wakefield accelerators (LWFA). It has been found that electron bunches exiting from state-of-the-art LWFAs exhibit a normalized 6-dimensional beam brightness comparable to those in conventional linear accelerators. Effectively exploiting this beneficial beam property for LWFA-based FELs is challenging due to the extreme initial conditions particularly in terms of beam divergence and energy spread. Several different approaches for capturing, reshaping and matching LWFA beams to suited undulators, such as bunch decompression or transverse-gradient undulator schemes, are currently being explored. In this article the transverse gradient undulator concept will be discussed with a focus on recent experimental achievements. △ Less

Submitted 13 December, 2017; originally announced December 2017.

Comments: NIMA_PROCEEDINGS-D-17-00178

Coupling the Gaussian free fields with free and with zero boundary conditions via common level lines

Authors: Wei Qian, Wendelin Werner

Abstract: We point out a new simple way to couple the Gaussian Free Field (GFF) with free boundary conditions in a two-dimensional domain with the GFF with zero boundary conditions in the same domain: Starting from the latter, one just has to sample at random all the signs of the height gaps on its boundary-touching zero-level lines (these signs are alternating for the zero-boundary GFF) in order to obtain… ▽ More We point out a new simple way to couple the Gaussian Free Field (GFF) with free boundary conditions in a two-dimensional domain with the GFF with zero boundary conditions in the same domain: Starting from the latter, one just has to sample at random all the signs of the height gaps on its boundary-touching zero-level lines (these signs are alternating for the zero-boundary GFF) in order to obtain a free boundary GFF. Constructions and couplings of the free boundary GFF and its level lines via soups of reflected Brownian loops and their clusters are also discussed. Such considerations show for instance that in a domain with an axis of symmetry, if one looks at the overlay of a single usual Conformal Loop Ensemble CLE(3) with its own symmetric image, one obtains the CLE(4)-type collection of level lines of a GFF with mixed zero/free boundary conditions in the half-domain. △ Less

Submitted 3 March, 2018; v1 submitted 13 March, 2017; originally announced March 2017.

Comments: 26 pages. To appear in Comm. Math. Phys

Journal ref: Communications in Mathematical Physics 361 (2018) 53-80

Connection probabilities for conformal loop ensembles

Authors: Jason Miller, Wendelin Werner

Abstract: The goal of the present paper is to explain, based on properties of the conformal loop ensembles CLE$_κ$ (both with simple and non-simple loops, i.e., for the whole range $κ\in (8/3, 8)$) how to derive the connection probabilities in conformal rectangles for a conditioned version of CLE$_κ$ which can be interpreted as a CLE$_κ$ with wired/free/wired/free boundary conditions on four boundary arcs (… ▽ More The goal of the present paper is to explain, based on properties of the conformal loop ensembles CLE$_κ$ (both with simple and non-simple loops, i.e., for the whole range $κ\in (8/3, 8)$) how to derive the connection probabilities in conformal rectangles for a conditioned version of CLE$_κ$ which can be interpreted as a CLE$_κ$ with wired/free/wired/free boundary conditions on four boundary arcs (the wired parts being viewed as portions of to-be-completed loops). In particular, in the case of a conformal square, we prove that the probability that the two wired sides hook up so that they create one single loop is equal to $1/(1 - 2 \cos (4 π/ κ))$. Comparing this with the corresponding connection probabilities for discrete O(N) models for instance indicates that if a dilute O(N) model (respectively a critical FK(q)-percolation model on the square lattice) has a non-trivial conformally invariant scaling limit, then necessarily this scaling limit is CLE$_κ$ where $κ$ is the value in $(8/3, 4]$ such that $-2 \cos (4 π/ κ)$ is equal to $N$ (resp. the value in $[4,8)$ such that $-2 \cos (4π/ κ)$ is equal to $\sqrt {q}$). Our arguments and computations build on the one hand on Dubédat's SLE commutation relations (as developed and used by Dubédat, Zhan or Bauer-Bernard-Kytölä) and on the other hand, on the construction and properties of the conformal loop ensembles and their relation to Brownian loop-soups, restriction measures, and the Gaussian free field, as recently derived in works with Sheffield and with Qian. △ Less

Submitted 28 May, 2018; v1 submitted 9 February, 2017; originally announced February 2017.

Comments: 37 pages, to appear in Comm. in Math. Phys

Journal ref: Communications in Mathematical Physics 362 (2018) 415-453

Stellar & Planetary Parameters for K2's Late Type Dwarf Systems from C1 to C5

Authors: Arturo O. Martinez, Ian J. M. Crossfield, Joshua E. Schlieder, Courtney D. Dressing, Christian Obermeier, John Livingston, Simona Ciceri, Sarah Peacock, Charles A. Beichman, Sébastien Lépine, Kimberly M. Aller, Quadry A. Chance, Erik A. Petigura, Andrew W. Howard, Michael W. Werner

Abstract: The NASA K2 mission uses photometry to find planets transiting stars of various types. M dwarfs are of high interest since they host more short period planets than any other type of main sequence stars and transiting planets around M dwarfs have deeper transits compared to other main sequence stars. In this paper, we present stellar parameters from K and M dwarfs hosting transiting planet candidat… ▽ More The NASA K2 mission uses photometry to find planets transiting stars of various types. M dwarfs are of high interest since they host more short period planets than any other type of main sequence stars and transiting planets around M dwarfs have deeper transits compared to other main sequence stars. In this paper, we present stellar parameters from K and M dwarfs hosting transiting planet candidates discovered by our team. Using the SOFI spectrograph on the European Southern Observatory's New Technology Telescope, we obtained R ~ 1000 J-, H-, and K-band (0.95 - 2.52 microns) spectra of 34 late-type K2 planet and candidate planet host systems and 12 bright K4-M5 dwarfs with interferometrically measured radii and effective temperatures. Out of our 34 late-type K2 targets, we identify 27 of these stars as M dwarfs. We measure equivalent widths of spectral features, derive calibration relations using stars with interferometric measurements, and estimate stellar radii, effective temperatures, masses, and luminosities for the K2 planet hosts. Our calibrations provide radii and temperatures with median uncertainties of 0.059 R_sun (16.09%) and 160 K (4.33%), respectively. We then reassess the radii and equilibrium temperatures of known and candidate planets based on our spectroscopically derived stellar parameters. Since a planet's radius and equilibrium temperature depends on the parameters of its host star, our study provides more precise planetary parameters for planets and candidates orbiting late-type stars observed with K2. We find a median planet radius and an equilibrium temperature of approximately 3R_earth and 500 K, with several systems (K2-18b and K2-72e) receiving near-Earth like levels of incident irradiation. △ Less

Submitted 7 February, 2017; v1 submitted 3 January, 2017; originally announced January 2017.

Comments: 12 pages, 10 figures, 5 tables, accepted to ApJ. Also see related paper by Dressing et al. titled 'Characterizing K2 Candidate Planetary Systems Orbiting Low-Mass Stars I: Classifying Low-Mass Stars Observed During Campaigns 1-7'. Version 2 updates: Updated table 5, stellar parameters, spectra, and planet parameters posted to ExoFOP-K2 (https://exofop.ipac.caltech.edu/k2/), and updated references

Hankel-Norm Approximation of Large-Scale Descriptor Systems

Authors: Peter Benner, Steffen W. R. Werner

Abstract: The Hankel-norm approximation is a model reduction method which provides the best approximation in the Hankel semi-norm. In this paper the computation of the optimal Hankel-norm approximation is generalized to the case of linear time-invariant continuous-time descriptor systems. An efficient algorithm is developed by refining the generalized balanced truncation square root method. For a wide pract… ▽ More The Hankel-norm approximation is a model reduction method which provides the best approximation in the Hankel semi-norm. In this paper the computation of the optimal Hankel-norm approximation is generalized to the case of linear time-invariant continuous-time descriptor systems. An efficient algorithm is developed by refining the generalized balanced truncation square root method. For a wide practical usage, adaptations of the introduced algorithm towards stable computations and sparse systems are made as well as an approach for a projection-free algorithm. To show the approximation behavior of the introduced method, numerical examples are presented. △ Less

Submitted 31 July, 2018; v1 submitted 19 December, 2016; originally announced December 2016.

Comments: 25 pages, 5 figures

Journal ref: Adv. Comput. Math., 46(3):40, 2020

Inductive limits of finite dimensional hermitian symmetric spaces and K-theory

Authors: Dennis Bohle, Wend Werner

Abstract: K-Theory for hermitian symmetric spaces of non-compact type, as developed recently by the authors, allows to put Cartan's classification into a homological perspective. We apply this method to the case of inductive limits of finite dimensional hermitian symmetric spaces. This might be seen as an indication of how much more powerful the homological theory is in comparison to the more classical appr… ▽ More K-Theory for hermitian symmetric spaces of non-compact type, as developed recently by the authors, allows to put Cartan's classification into a homological perspective. We apply this method to the case of inductive limits of finite dimensional hermitian symmetric spaces. This might be seen as an indication of how much more powerful the homological theory is in comparison to the more classical approach. When seen from high above, we follow the path laid out by a similar result in the theory of C*-algebras. Important is a clear picture of the behavior of morphisms between bounded symmetric domains of finite dimensions, which is more complex than in the C*-case, as well as an accessible K-theory. We furthermore have to slightly modify the invariant from our previous work. Roughly, we use traces left by co-root lattices on K-groups, instead of co-roots themselves, which had been used previously. △ Less

Submitted 22 September, 2016; originally announced September 2016.

MSC Class: 17C10; 17C37; 19E99; 32M15; 46G20; 46L08; 46L80; 46M15; 46M40; 46T05; 53B35; 53C35; 58B20; 58B25

Non-simple SLE curves are not determined by their range

Authors: Jason Miller, Scott Sheffield, Wendelin Werner

Abstract: We show that when observing the range of a chordal SLE$_κ$ curve for $κ\in (4,8)$, it is not possible to recover the order in which the points have been visited. We also derive related results about conformal loop ensembles (CLE): (i) The loops in a CLE$_κ$ for $κ\in (4,8)$ are not determined by the CLE$_κ$ gasket. (ii) The continuum percolation interfaces defined in the fractal carpets of conform… ▽ More We show that when observing the range of a chordal SLE$_κ$ curve for $κ\in (4,8)$, it is not possible to recover the order in which the points have been visited. We also derive related results about conformal loop ensembles (CLE): (i) The loops in a CLE$_κ$ for $κ\in (4,8)$ are not determined by the CLE$_κ$ gasket. (ii) The continuum percolation interfaces defined in the fractal carpets of conformal loop ensembles CLE$_κ$ for $κ\in (8/3, 4)$ (we defined these percolation interfaces in previous work, and showed there that they are SLE$_{16/κ}$ curves) are not determined by the CLE$_κ$ carpet that they are defined in. △ Less

Submitted 19 November, 2018; v1 submitted 15 September, 2016; originally announced September 2016.

Comments: Final version. 44 pages, including an appendix on bichordal resampling. To appear in J. Europ. Math. Soc

Journal ref: J. Europ. Math. Soc. 22, 669-716 , 2020

The random pseudo-metric on a graph defined via the zero-set of the Gaussian free field on its metric graph

Authors: Titus Lupu, Wendelin Werner

Abstract: We further investigate properties of the Gaussian free field (GFF) on the metric graph associated to a discrete weighted graph (where the edges of the latter are replaced by continuous line-segments of appropriate length) that has been introduced by the first author. On such a metric graph, the GFF is a random continuous function that generalises one-dimensional Brownian bridges so that one-dimens… ▽ More We further investigate properties of the Gaussian free field (GFF) on the metric graph associated to a discrete weighted graph (where the edges of the latter are replaced by continuous line-segments of appropriate length) that has been introduced by the first author. On such a metric graph, the GFF is a random continuous function that generalises one-dimensional Brownian bridges so that one-dimensional techniques can be used. In the present paper, we define and study the pseudo-metric defined on the metric graph (and therefore also on the discrete graph itself), where the length of a path on the metric graph is defined to be the local time at level zero accumulated by the Gaussian free field along this path. We first derive a pathwise transformation that relates the GFF on the metric graph with the reflected GFF on the metric graph via the pseudo-distance defined by the latter. This is a generalisation of Paul Lévy's result relating the local time at zero of Brownian motion to the supremum of another Brownian motion. We also compute explicitly the distribution of certain functionals of this pseudo-metric and of the GFF. In particular, we point out that when the boundary consists of just two points, the law of the pseudo-distance between them depends solely on the resistance of the network between them. We then discuss questions related to the scaling limit of this pseudo-metric in the two-dimensional case, which should be the conformally invariant way to measure distances between CLE(4) loops introduced and studied by the second author with Wu, and by Sheffield, Watson and Wu. Our explicit laws on metric graphs also lead to new conjectures for related functionals of the continuum GFF on fairly general Riemann surfaces. △ Less

Submitted 5 July, 2017; v1 submitted 21 July, 2016; originally announced July 2016.

Comments: 32 pages, to appear in Probab. Th. and rel. Fields

Journal ref: Probability Theory and Related Fields 171, 775-818, 2018

Science Impacts of the SPHEREx All-Sky Optical to Near-Infrared Spectral Survey: Report of a Community Workshop Examining Extragalactic, Galactic, Stellar and Planetary Science

Authors: Olivier Doré, Michael W. Werner, Matt Ashby, Pancha Banerjee, Nick Battaglia, James Bauer, Robert A. Benjamin, Lindsey E. Bleem, Jamie Bock, Adwin Boogert, Philip Bull, Peter Capak, Tzu-Ching Chang, Jean Chiar, Seth H. Cohen, Asantha Cooray, Brendan Crill, Michael Cushing, Roland de Putter, Simon P. Driver, Tim Eifler, Chang Feng, Simone Ferraro, Douglas Finkbeiner, B. Scott Gaudi , et al. (43 additional authors not shown)

Abstract: SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $μ$m [with R$\sim$41.4] and 4.18 and 5.00 $μ$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circums… ▽ More SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $μ$m [with R$\sim$41.4] and 4.18 and 5.00 $μ$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST. △ Less

Submitted 22 June, 2016; originally announced June 2016.

Comments: Report of the First SPHEREx Community Workshop, http://spherex.caltech.edu/Workshop.html , 84 pages, 28 figures