Showing 1–50 of 51 results for author: Haber, A

Beyond modified Urca: the nucleon width approximation for flavor-changing processes in dense matter

Authors: Mark G. Alford, Alexander Haber, Ziyuan Zhang

Abstract: Flavor-changing charged current ("Urca") processes are of central importance in the astrophysics of neutron stars. Standard calculations approximate the Urca rate as the sum of two contributions, direct Urca and modified Urca. Attempts to make modified Urca calculations more accurate have been impeded by an unphysical divergence at the direct Urca threshold density. In this paper we describe a sys… ▽ More Flavor-changing charged current ("Urca") processes are of central importance in the astrophysics of neutron stars. Standard calculations approximate the Urca rate as the sum of two contributions, direct Urca and modified Urca. Attempts to make modified Urca calculations more accurate have been impeded by an unphysical divergence at the direct Urca threshold density. In this paper we describe a systematically improvable approach where, in the simplest approximation, instead of modified Urca we include an imaginary part of the nucleon mass (nucleon width). The total Urca rate is then obtained via a straightforward generalization of the direct Urca calculation, yielding results that agree with both direct and modified Urca at the densities where those approximations are valid. At low densities, we observe an enhancement of the rate by more than an order of magnitude, with important ramifications for neutron star cooling and other transport properties. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 6 pages, 6 figures

Finite-temperature expansion of the dense-matter equation of state

Authors: Debora Mroczek, Nanxi Yao, Katherine Zine, Jacquelyn Noronha-Hostler, Veronica Dexheimer, Alexander Haber, Elias R. Most

Abstract: In this work we provide a new, well-controlled expansion of the equation of state of dense matter from zero to finite temperatures ($T$), while covering a wide range of charge fractions ($Y_Q$), from pure neutron to isospin symmetric nuclear matter. Our expansion can be used to describe neutron star mergers and core-collapse supernova explosions using as a starting point neutron star observations,… ▽ More In this work we provide a new, well-controlled expansion of the equation of state of dense matter from zero to finite temperatures ($T$), while covering a wide range of charge fractions ($Y_Q$), from pure neutron to isospin symmetric nuclear matter. Our expansion can be used to describe neutron star mergers and core-collapse supernova explosions using as a starting point neutron star observations, while maintaining agreement with laboratory data, in a model independent way. We suggest new thermodynamic quantities of interest that can be calculated from theoretical models or directly inferred by experimental data that can help constrain the finite $T$ equation of state. With our new method, we can quantify the uncertainty in our finite $T$ and $Y_Q$ expansions in a well-controlled manner without making assumptions about the underlying degrees of freedom. We can reproduce results from a microscopic equation of state up to $T=100$ MeV for baryon chemical potential $μみゅー_B\gtrsim 1100$ MeV ($\sim1-2 \ n_{\rm sat}$) within $5\%$ error, with even better results for larger $μみゅー_B$ and/or lower $T$. We investigate the sources of numerical and theoretical uncertainty and discuss future directions of study. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 33 pages, 12 figures

Verbesserung des Record Linkage für die Gesundheitsforschung in Deutschland

Authors: Timm Intemann, Knut Kaulke, Dennis-Kenji Kipker, Vanessa Lettieri, Christoph Stallmann, Carsten O. Schmidt, Lars Geidel, Martin Bialke, Christopher Hampf, Dana Stahl, Martin Lablans, Florens Rohde, Martin Franke, Klaus Kraywinkel, Joachim Kieschke, Sebastian Bartholomäus, Anatol-Fiete Näher, Galina Tremper, Mohamed Lambarki, Stefanie March, Fabian Prasser, Anna Christine Haber, Johannes Drepper, Irene Schlünder, Toralf Kirsten , et al. (5 additional authors not shown)

Abstract: Record linkage means linking data from multiple sources. This approach enables the answering of scientific questions that cannot be addressed using single data sources due to limited variables. The potential of linked data for health research is enormous, as it can enhance prevention, treatment, and population health policies. Due the sensitivity of health data, there are strict legal requirements… ▽ More Record linkage means linking data from multiple sources. This approach enables the answering of scientific questions that cannot be addressed using single data sources due to limited variables. The potential of linked data for health research is enormous, as it can enhance prevention, treatment, and population health policies. Due the sensitivity of health data, there are strict legal requirements to prevent potential misuse. However, these requirements also limit the use of health data for research, thereby hindering innovations in prevention and care. Also, comprehensive Record linkage in Germany is often challenging due to lacking unique personal identifiers or interoperable solutions. Rather, the need to protect data is often weighed against the importance of research aiming at healthcare enhancements: for instance, data protection officers may demand the informed consent of individual study participants for data linkage, even when this is not mandatory. Furthermore, legal frameworks may be interpreted differently on varying occasions. Given both, technical and legal challenges, record linkage for health research in Germany falls behind the standards of other European countries. To ensure successful record linkage, case-specific solutions must be developed, tested, and modified as necessary before implementation. This paper discusses limitations and possibilities of various data linkage approaches tailored to different use cases in compliance with the European General Data Protection Regulation. It further describes requirements for achieving a more research-friendly approach to linking health data records in Germany. Additionally, it provides recommendations to legislators. The objective of this work is to improve record linkage for health research in Germany. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: in German language

State-Robust Observability Measures for Sensor Selection in Nonlinear Dynamic Systems

Authors: Mohamad H. Kazma, Sebastian A. Nugroho, Aleksandar Haber, Ahmad F. Taha

Abstract: This paper explores the problem of selecting sensor nodes for a general class of nonlinear dynamical networks. In particular, we study the problem by utilizing altered definitions of observability and open-loop lifted observers. The approach is performed by discretizing the system's dynamics using the implicit Runge-Kutta method and by introducing a state-averaged observability measure. The observ… ▽ More This paper explores the problem of selecting sensor nodes for a general class of nonlinear dynamical networks. In particular, we study the problem by utilizing altered definitions of observability and open-loop lifted observers. The approach is performed by discretizing the system's dynamics using the implicit Runge-Kutta method and by introducing a state-averaged observability measure. The observability measure is computed for a number of perturbed initial states in the vicinity of the system's true initial state. The sensor node selection problem is revealed to retain the submodular and modular properties of the original problem. This allows the problem to be solved efficiently using a greedy algorithm with a guaranteed performance bound while showing an augmented robustness to unknown or uncertain initial conditions. The validity of this approach is numerically demonstrated on a $H_{2}/O_{2}$ combustion reaction network. △ Less

Submitted 13 July, 2023; originally announced July 2023.

Comments: To Appear in the 62$^{\text{nd}}$ IEEE Conference on Decision and Control (CDC'2023), Singapore, Decemeber 2023

Isospin Equilibration in Neutron Star Mergers

Authors: Mark G. Alford, Alexander Haber, Ziyuan Zhang

Abstract: We analyze the isospin equilibration properties of neutrinoless nuclear ($npe$) matter in the temperature and density range that is relevant to neutron star mergers. Our analysis incorporates neutrino-transparency corrections to the isospin (``beta'') equilibrium condition which become noticeable at $T\gtrsim 1\,$MeV. We find that the isospin relaxation rate rises rapidly as temperature rises, and… ▽ More We analyze the isospin equilibration properties of neutrinoless nuclear ($npe$) matter in the temperature and density range that is relevant to neutron star mergers. Our analysis incorporates neutrino-transparency corrections to the isospin (``beta'') equilibrium condition which become noticeable at $T\gtrsim 1\,$MeV. We find that the isospin relaxation rate rises rapidly as temperature rises, and at $T\approx 5\,$MeV it is comparable to the timescale of the density oscillations that occur immediately after the merger. This produces a resonant peak in the bulk viscosity at $T\approx 5\,$MeV, which causes density oscillations to be damped on the timescale of the merger. Our calculations suggest that isospin relaxation dynamics may also be relevant when neutrinos are treated more accurately via neutrino transport schemes. △ Less

Submitted 12 May, 2024; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 12 pages, 5 figures v2: minor changes in wording, no changes to conclusions or results of v1, version published in PRC

Journal ref: Phys.Rev.C 109 (2024) 5, 055803

Tabulated Equations of State From Models Informed by Chiral Effective Field Theory

Authors: Mark G. Alford, Liam Brodie, Alexander Haber, Ingo Tews

Abstract: We construct four equation of state (EoS) tables, tabulated over a range of temperatures, densities, and charge fractions, relevant for neutron star applications such as simulations of neutron star mergers. The EoS are computed from a relativistic mean-field theory constrained by the pure neutron matter EoS from chiral effective field theory, inferred properties of isospin-symmetric nuclear matter… ▽ More We construct four equation of state (EoS) tables, tabulated over a range of temperatures, densities, and charge fractions, relevant for neutron star applications such as simulations of neutron star mergers. The EoS are computed from a relativistic mean-field theory constrained by the pure neutron matter EoS from chiral effective field theory, inferred properties of isospin-symmetric nuclear matter, and astrophysical observations of neutron star structure. To model nuclear matter at low densities, we attach an EoS that models inhomogeneous nuclear matter at arbitrary temperatures and charge fractions. The four EoS tables we develop are available from the CompOSE EoS repository compose.obspm.fr/eos/297 and gitlab.com/ahaber/qmc-rmf-tables. △ Less

Submitted 17 November, 2023; v1 submitted 16 April, 2023; originally announced April 2023.

Comments: 8 pages, 3 figures v2: version published in Phys. Scr. 98, added finite T chiral EFT comparison

Report number: LA-UR-23-23837

Journal ref: Phys. Scr. 98 125302 (2023)

Theoretical and Experimental Constraints for the Equation of State of Dense and Hot Matter

Authors: Rajesh Kumar, Veronica Dexheimer, Johannes Jahan, Jorge Noronha, Jacquelyn Noronha-Hostler, Claudia Ratti, Nico Yunes, Angel Rodrigo Nava Acuna, Mark Alford, Mahmudul Hasan Anik, Debarati Chatterjee, Katerina Chatziioannou, Hsin-Yu Chen, Alexander Clevinger, Carlos Conde, Nikolas Cruz-Camacho, Travis Dore, Christian Drischler, Hannah Elfner, Reed Essick, David Friedenberg, Suprovo Ghosh, Joaquin Grefa, Roland Haas, Alexander Haber , et al. (35 additional authors not shown)

Abstract: This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutro… ▽ More This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutron stars and their mergers. The validity of different constraints, concerning specific conditions and ranges of applicability, is also provided. △ Less

Submitted 12 June, 2024; v1 submitted 29 March, 2023; originally announced March 2023.

Journal ref: Living Reviews in Relativity (2024) 27:3

Nuclear and Hybrid Equations of State in Light of the Low-Mass Compact Star in HESS J1731-347

Authors: Liam Brodie, Alexander Haber

Abstract: We sample over $500$ relativistic mean-field theories constrained by chiral effective field theory and properties of isospin-symmetric nuclear matter and test them against known stellar structure constraints. This includes a recent mass and radius measurement of a compact object in supernova remnant HESS J1731-347, with an unusually low mass of $M=0.77^{+0.20}_{-0.17}\,M_{\odot}$ and a compact rad… ▽ More We sample over $500$ relativistic mean-field theories constrained by chiral effective field theory and properties of isospin-symmetric nuclear matter and test them against known stellar structure constraints. This includes a recent mass and radius measurement of a compact object in supernova remnant HESS J1731-347, with an unusually low mass of $M=0.77^{+0.20}_{-0.17}\,M_{\odot}$ and a compact radius of $R=10.4^{+0.86}_{-0.78}$ km. We show that none of the sampled nuclear models meet all constraints at the $68\,\%$ credibility level, but that hybrid equations of state with a quark matter inner core and nuclear outer core easily can. This indicates a tension between astrophysical constraints and low-energy nuclear theory. △ Less

Submitted 11 August, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: 7 pages, 3 figures, version to appear in PRC, minor changes

Journal ref: Phys. Rev. C 108, 025806 (2023)

Data-driven Estimation, Tracking, and System Identification of Deterministic and Stochastic Optical Spot Dynamics

Authors: Aleksandar Haber, Michael Krainak

Abstract: Stabilization, disturbance rejection, and control of optical beams and optical spots are ubiquitous problems that are crucial for the development of optical systems for ground and space telescopes, free-space optical communication terminals, precise beam steering systems, and other types of optical systems. High-performance disturbance rejection and control of optical spots require the development… ▽ More Stabilization, disturbance rejection, and control of optical beams and optical spots are ubiquitous problems that are crucial for the development of optical systems for ground and space telescopes, free-space optical communication terminals, precise beam steering systems, and other types of optical systems. High-performance disturbance rejection and control of optical spots require the development of disturbance estimation and data-driven Kalman filter methods. Motivated by this, we propose a unified and experimentally verified data-driven framework for optical-spot disturbance modeling and tuning of covariance matrices of Kalman filters. Our approach is based on covariance estimation, nonlinear optimization, and subspace identification methods. Also, we use spectral factorization methods to emulate optical-spot disturbances with a desired power spectral density in an optical laboratory environment. We test the effectiveness of the proposed approaches on an experimental setup consisting of a piezo tip-tilt mirror, piezo linear actuator, and a CMOS camera. △ Less

Submitted 29 January, 2023; originally announced January 2023.

Comments: 20 pages, 8 figures

Long Range Plan: Dense matter theory for heavy-ion collisions and neutron stars

Authors: Alessandro Lovato, Travis Dore, Robert D. Pisarski, Bjoern Schenke, Katerina Chatziioannou, Jocelyn S. Read, Philippe Landry, Pawel Danielewicz, Dean Lee, Scott Pratt, Fabian Rennecke, Hannah Elfner, Veronica Dexheimer, Rajesh Kumar, Michael Strickland, Johannes Jahan, Claudia Ratti, Volodymyr Vovchenko, Mikhail Stephanov, Dekrayat Almaalol, Gordon Baym, Mauricio Hippert, Jacquelyn Noronha-Hostler, Jorge Noronha, Enrico Speranza , et al. (39 additional authors not shown)

Abstract: Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theo… ▽ More Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theory of dense baryonic matter that connects low- and high-energy nuclear physics, astrophysics, gravitational waves physics, and data science △ Less

Submitted 7 November, 2022; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 70 pages, 3 figures, White Paper for the Long Range Plan for Nuclear Science

Report number: LA-UR-22-31648

The computational and learning benefits of Daleian neural networks

Authors: Adam Haber, Elad Schneidman

Abstract: Dale's principle implies that biological neural networks are composed of neurons that are either excitatory or inhibitory. While the number of possible architectures of such Daleian networks is exponentially smaller than non-Daleian ones, the computational and functional implications of using Daleian networks by the brain are mostly unknown. Here, we use models of recurrent spiking neural networks… ▽ More Dale's principle implies that biological neural networks are composed of neurons that are either excitatory or inhibitory. While the number of possible architectures of such Daleian networks is exponentially smaller than non-Daleian ones, the computational and functional implications of using Daleian networks by the brain are mostly unknown. Here, we use models of recurrent spiking neural networks and rate-based networks to show, surprisingly, that despite the structural limitations on Daleian networks, they can approximate the computation performed by non-Daleian networks to a very high degree of accuracy. Moreover, we find that Daleian networks are more functionally robust to synaptic noise. We then show that unlike non-Daleian networks, Daleian ones can learn efficiently by tuning single neuron features, nearly as well as learning by tuning individual synaptic weights - suggesting a simpler and more biologically plausible learning mechanism. We thus suggest that in addition to architectural simplicity, Dale's principle confers computational and learning benefits for biological networks, and offers new directions for constructing and training biologically-inspired artificial neural networks △ Less

Submitted 12 October, 2022; originally announced October 2022.

Subspace identification of low-dimensional Structural-Thermal-Optical-Performance (STOP) models of reflective optics

Authors: Aleksandar Haber, John E. Draganov, Michael Krainak

Abstract: In this paper, we investigate the feasibility of using subspace system identification techniques for estimating transient Structural-Thermal-Optical Performance (STOP) models of reflective optics. As a test case, we use a Newtonian telescope structure. This work is motivated by the need for the development of model-based data-driven techniques for prediction, estimation, and control of thermal eff… ▽ More In this paper, we investigate the feasibility of using subspace system identification techniques for estimating transient Structural-Thermal-Optical Performance (STOP) models of reflective optics. As a test case, we use a Newtonian telescope structure. This work is motivated by the need for the development of model-based data-driven techniques for prediction, estimation, and control of thermal effects and thermally-induced wavefront aberrations in optical systems, such as ground and space telescopes, optical instruments operating in harsh environments, optical lithography machines, and optical components of high-power laser systems. We estimate and validate a state-space model of a transient STOP dynamics. First, we model the system in COMSOL Multiphysics. Then, we use LiveLink for MATLAB software module to export the wavefront aberrations data from COMSOL to MATLAB. This data is used to test the subspace identification method that is implemented in Python. One of the main challenges in modeling and estimation of STOP models is that they are inherently large-dimensional. The large-scale nature of STOP models originates from the coupling of optical, thermal, and structural phenomena and physical processes. Our results show that large-dimensional STOP dynamics of the considered optical system can be accurately estimated by low-dimensional state-space models. Due to their low-dimensional nature and state-space forms, these models can effectively be used for the prediction, estimation, and control of thermally-induced wavefront aberrations. The developed MATLAB, COMSOL, and Python codes are available online. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: Accepted for publication and presentation at SPIE Optics and Photonics: Optical Engineering + Applications (OP22O), San Diego (2022)

Emergence of microphysical viscosity in binary neutron star post-merger dynamics

Authors: Elias R. Most, Alexander Haber, Steven P. Harris, Ziyuan Zhang, Mark G. Alford, Jorge Noronha

Abstract: In nuclear matter in neutron stars the flavor content (e.g., proton fraction) is subject to weak interactions, establishing flavor ($βべーた$-)equilibrium. During the merger of two neutron stars there can be deviations from this equilibrium. By incorporating Urca processes into general-relativistic hydrodynamics simulations, we study the resulting out-of-equilibrium dynamics during the collision. We pro… ▽ More In nuclear matter in neutron stars the flavor content (e.g., proton fraction) is subject to weak interactions, establishing flavor ($βべーた$-)equilibrium. During the merger of two neutron stars there can be deviations from this equilibrium. By incorporating Urca processes into general-relativistic hydrodynamics simulations, we study the resulting out-of-equilibrium dynamics during the collision. We provide the first direct evidence that microphysical transport effects at late times reach a hydrodynamic regime with a nonzero bulk viscosity, making neutron star collisions intrinsically viscous. Finally, we identify signatures of this process in the post-merger gravitational wave emission. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Report number: INT-PUB-22-017

Journal ref: ApJL 967 L14 (2024)

Relativistic mean-field theories for neutron-star physics based on chiral effective field theory

Authors: Mark G. Alford, Liam Brodie, Alexander Haber, Ingo Tews

Abstract: We describe and implement a procedure for determining the couplings of a Relativistic Mean-Field Theory (RMFT) that is optimized for application to neutron star phenomenology. In the standard RMFT approach, the couplings are constrained by comparing the theory's predictions for symmetric matter at saturation density with measured nuclear properties. The theory is then applied to neutron stars whic… ▽ More We describe and implement a procedure for determining the couplings of a Relativistic Mean-Field Theory (RMFT) that is optimized for application to neutron star phenomenology. In the standard RMFT approach, the couplings are constrained by comparing the theory's predictions for symmetric matter at saturation density with measured nuclear properties. The theory is then applied to neutron stars which consist of neutron-rich matter at densities ranging up to several times saturation density, which allows for additional astrophysical constraints. In our approach, rather than using the RMFT to extrapolate from symmetric to neutron-rich matter and from finite-sized nuclei to uniform matter, we fit the RMFT to properties of uniform pure neutron matter obtained from chiral effective field theory. Chiral effective field theory incorporates the experimental data for nuclei in the framework of a controlled expansion for nuclear forces valid at nuclear densities and enables us to account for theoretical uncertainties when fitting the RMFT. We construct four simple RMFTs that span the uncertainties provided by chiral effective field theory for neutron matter, and are consistent with current astrophysical constraints on the equation of state. Our RMFTs can be used to model the properties of neutron-rich matter across the vast range of densities and temperatures encountered in neutron stars and their mergers. △ Less

Submitted 21 November, 2022; v1 submitted 20 May, 2022; originally announced May 2022.

Comments: V2: Minor corrections, version published in PRC. 12 pages and 5 figures

Report number: LA-UR-22-23972

Journal ref: Phys.Rev.C 106 (2022) 5, 055804

A novel method for adaptive control of deformable mirrors

Authors: Aleksandar Haber

Abstract: For sufficiently wide ranges of applied control signals (control voltages), MEMS and piezoelectric Deformable Mirrors (DMs), exhibit nonlinear behavior. The nonlinear behavior manifests itself in nonlinear actuator couplings, nonlinear actuator deformation characteristics, and in the case of piezoelectric DMs, hysteresis. Furthermore, in a number of situations, DM behavior can change over time, an… ▽ More For sufficiently wide ranges of applied control signals (control voltages), MEMS and piezoelectric Deformable Mirrors (DMs), exhibit nonlinear behavior. The nonlinear behavior manifests itself in nonlinear actuator couplings, nonlinear actuator deformation characteristics, and in the case of piezoelectric DMs, hysteresis. Furthermore, in a number of situations, DM behavior can change over time, and this requires a procedure for updating the DM models on the basis of the observed data. If not properly modeled and if not taken into account when designing control algorithms, nonlinearities, and time-varying DM behavior, can significantly degrade the achievable closed-loop performance of Adaptive Optics (AO) systems. Widely used approaches for DM control are based on pre-estimated linear time-invariant DM models in the form of influence matrices. Often, these models are not being updated during system operation. Consequently, when the nonlinear DM behavior is excited by control signals with wide operating ranges, or when the DM behavior changes over time, the state-of-the-art DM control approaches relying upon linear control methods, might not be able to produce a satisfactory closed-loop performance of an AO system. Motivated by these key facts, we present a novel method for data-driven DM control. Our approach combines a simple open-loop control method with a recursive least squares method for dynamically updating the DM model. The DM model is constantly being updated on the basis of the dynamically changing DM operating points. That is, the proposed method updates both the control actions and the DM model during the system operation. We experimentally verify this approach on a Boston Micromachines MEMS DM with 140 actuators. Preliminary experimental results reported in this manuscript demonstrate good potential for using the developed method for DM control. △ Less

Submitted 9 March, 2022; originally announced March 2022.

Comments: SPIE Proceedings, SPIE BiOS, 2022, San Francisco, California, United States, Proceedings Volume 11969, Adaptive Optics and Wavefront Control for Biological Systems VIII; 1196902

Dual-Update Data-Driven Control of Deformable Mirrors Using Walsh Basis Functions

Authors: Aleksandar Haber, Thomas Bifano

Abstract: In this paper, we develop a novel data-driven method for Deformable Mirror (DM) control. The developed method updates both the DM model and DM control actions that produce desired mirror surface shapes. The novel method explicitly takes into account actuator constraints and couples a feedback control algorithm with an algorithm for recursive estimation of DM influence function models. In addition… ▽ More In this paper, we develop a novel data-driven method for Deformable Mirror (DM) control. The developed method updates both the DM model and DM control actions that produce desired mirror surface shapes. The novel method explicitly takes into account actuator constraints and couples a feedback control algorithm with an algorithm for recursive estimation of DM influence function models. In addition to this, we explore the possibility of using Walsh basis functions for DM control. By expressing the desired and observed mirror surface shapes as sums of Walsh pattern matrices, we formulate the control problem in the 2D Walsh basis domain. We thoroughly experimentally verify the developed approach on a 140-actuator MEMS DM, developed by Boston Micromachines. Our results show that the novel method produces the root-mean-square surface error in the $14-40$ nanometer range. These results can additionally be improved by tuning the control and estimation parameters. The developed approach is also applicable to other DM types, such as for example, segmented DMs. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 25 pages, 11 figures

Beta equilibrium under neutron star merger conditions

Authors: Mark G. Alford, Alexander Haber, Steven P. Harris, Ziyuan Zhang

Abstract: We calculate the nonzero-temperature correction to the beta equilibrium condition in nuclear matter under neutron star merger conditions, in the temperature range $1\,$MeV$ < T \lesssim 5\,$MeV. We improve on previous work by using a consistent description of nuclear matter based on the IUF and SFHo relativistic mean field models. This includes using relativistic dispersion relations for the nucle… ▽ More We calculate the nonzero-temperature correction to the beta equilibrium condition in nuclear matter under neutron star merger conditions, in the temperature range $1\,$MeV$ < T \lesssim 5\,$MeV. We improve on previous work by using a consistent description of nuclear matter based on the IUF and SFHo relativistic mean field models. This includes using relativistic dispersion relations for the nucleons, which we show is essential in these models. We find that the nonzero-temperature correction can be of order $10$ to $20\,$MeV, and plays an important role in the correct calculation of Urca rates, which can be wrong by factors of $10$ or more if it is neglected. △ Less

Submitted 9 October, 2023; v1 submitted 6 August, 2021; originally announced August 2021.

Comments: v2: version published in Universe, minor changes in nomenclature, 20 pages, 14 figures v3: corrected typo in Eq.(C19)

Report number: INT-PUB-21-017

Journal ref: Universe 2021, 7(11), 399

Strangeness-changing Rates and Hyperonic Bulk Viscosity in Neutron Star Mergers

Authors: Mark G. Alford, Alexander Haber

Abstract: In this paper we present a computation of the rates of strangeness-changing processes and the resultant bulk viscosity in matter at the densities and temperatures typical of neutron star mergers. To deal with the high temperature in this environment we go beyond the Fermi surface approximation in our rate calculations and numerically evaluate the full phase space integral. We include processes whe… ▽ More In this paper we present a computation of the rates of strangeness-changing processes and the resultant bulk viscosity in matter at the densities and temperatures typical of neutron star mergers. To deal with the high temperature in this environment we go beyond the Fermi surface approximation in our rate calculations and numerically evaluate the full phase space integral. We include processes where quarks move between baryons via meson exchange: these have generally been omitted in previous analyses but provide the dominant contribution to the rates of strangeness-changing processes and the bulk viscosity. The calculation of these rates is an essential step towards any calculation of dissipation mechanisms in hyperonic matter in mergers. As one application, we calculate the dissipation times for density oscillations at the frequencies seen in merger simulations. We find that hyperon bulk viscosity for temperatures in the MeV regime can probably be neglected in this context, but becomes highly relevant for keV-range temperatures. △ Less

Submitted 6 December, 2020; v1 submitted 10 September, 2020; originally announced September 2020.

Comments: 17 pages, 11 figures v2: corrected computation of susceptibilities, resonant temperature for maximum bulk viscosity shifted to much lower temperatures

Journal ref: Phys. Rev. C 103, 045810 (2021)

COVID-19 Policy Impact Evaluation: A guide to common design issues

Authors: Noah A Haber, Emma Clarke-Deelder, Joshua A Salomon, Avi Feller, Elizabeth A Stuart

Abstract: Policy responses to COVID-19, particularly those related to non-pharmaceutical interventions, are unprecedented in scale and scope. Epidemiologists are more involved in policy decisions and evidence generation than ever before. However, policy impact evaluations always require a complex combination of circumstance, study design, data, statistics, and analysis. Beyond the issues that are faced for… ▽ More Policy responses to COVID-19, particularly those related to non-pharmaceutical interventions, are unprecedented in scale and scope. Epidemiologists are more involved in policy decisions and evidence generation than ever before. However, policy impact evaluations always require a complex combination of circumstance, study design, data, statistics, and analysis. Beyond the issues that are faced for any policy, evaluation of COVID-19 policies is complicated by additional challenges related to infectious disease dynamics and lags, lack of direct observation of key outcomes, and a multiplicity of interventions occurring on an accelerated time scale. The methods needed for policy-level impact evaluation are not often used or taught in epidemiology, and differ in important ways that may not be obvious. The volume and speed, and methodological complications of policy evaluations can make it difficult for decision-makers and researchers to synthesize and evaluate strength of evidence in COVID-19 health policy papers. In this paper, we (1) introduce the basic suite of policy impact evaluation designs for observational data, including cross-sectional analyses, pre/post, interrupted time-series, and difference-in-differences analysis, (2) demonstrate key ways in which the requirements and assumptions underlying these designs are often violated in the context of COVID-19, and (3) provide decision-makers and reviewers a conceptual and graphical guide to identifying these key violations. The overall goal of this paper is to help epidemiologists, policy-makers, journal editors, journalists, researchers, and other research consumers understand and weigh the strengths and limitations of evidence that is essential to decision-making. △ Less

Submitted 16 April, 2021; v1 submitted 3 September, 2020; originally announced September 2020.

Control Node Selection Algorithm for Nonlinear Dynamic Networks

Authors: Aleksandar Haber, Sebastian A. Nugroho, Patricio Torres, Ahmad F. Taha

Abstract: The coupled problems of selecting control nodes and designing control actions for nonlinear network dynamics are fundamental scientific problems with applications in many diverse fields. These problems are thoroughly studied for linear dynamics; however, in spite of a number of open research questions, methods for nonlinear network dynamics are less developed. As observed by various studies, the p… ▽ More The coupled problems of selecting control nodes and designing control actions for nonlinear network dynamics are fundamental scientific problems with applications in many diverse fields. These problems are thoroughly studied for linear dynamics; however, in spite of a number of open research questions, methods for nonlinear network dynamics are less developed. As observed by various studies, the prevailing graph-based controllability approaches for selecting control nodes might result in significantly suboptimal control performance for nonlinear dynamics. Herein we present a new, intuitive, and simple method for simultaneous control node selection and control sequence design for complex networks with nonlinear dynamics. The method is developed by incorporating the control node selection problem into an open-loop predictive control cost function and by solving the resulting mixed-integer optimization problem using a mesh adaptive direct search method. The developed framework is numerically robust and can deal with stiff networks, networks with non-smooth dynamics, as well as with control and actuator constraints. Good numerical performance of the method is demonstrated by testing it on prototypical Duffing oscillator and associative memory networks. The developed codes that can easily be adapted to models of other complex systems are available online. △ Less

Submitted 23 August, 2020; v1 submitted 17 August, 2020; originally announced August 2020.

Comments: accepted for publication in IEEE Control Systems Letters, 6 pages, 4 figures

Joint Sensor Node Selection and State Estimation for Nonlinear Networks and Systems

Authors: Aleksandar Haber

Abstract: State estimation and sensor selection problems for nonlinear networks and systems are ubiquitous problems that are important for the control, monitoring, analysis, and prediction of a large number of engineered and physical systems. Sensor selection problems are extensively studied for linear networks. However, less attention has been dedicated to networks with nonlinear dynamics. Furthermore, wid… ▽ More State estimation and sensor selection problems for nonlinear networks and systems are ubiquitous problems that are important for the control, monitoring, analysis, and prediction of a large number of engineered and physical systems. Sensor selection problems are extensively studied for linear networks. However, less attention has been dedicated to networks with nonlinear dynamics. Furthermore, widely used sensor selection methods relying on structural (graph-based) observability approaches might produce far from optimal results when applied to nonlinear network dynamics. In addition, state estimation and sensor selection problems are often treated separately, and this might decrease the overall estimation performance. To address these challenges, we develop a novel methodology for selecting sensor nodes for networks with nonlinear dynamics. Our main idea is to incorporate the sensor selection problem into an initial state estimation problem. The resulting mixed-integer nonlinear optimization problem is approximately solved using three methods. The good numerical performance of our approach is demonstrated by testing the algorithms on prototypical Duffing oscillator, associative memory, and chemical reaction networks. The developed codes are available online. △ Less

Submitted 21 March, 2021; v1 submitted 7 June, 2020; originally announced June 2020.

Comments: Paper accepted to IEEE Transactions on Network Science and Engineering, 19 pages, 7 figures

DAG With Omitted Objects Displayed (DAGWOOD): A framework for revealing causal assumptions in DAGs

Authors: Noah A Haber, Mollie E Wood, Sarah Wieten, Alexander Breskin

Abstract: Directed acyclic graphs (DAGs) are frequently used in epidemiology as a method to encode causal inference assumptions. We propose the DAGWOOD framework to bring many of those encoded assumptions to the forefront. DAGWOOD combines a root DAG (the DAG in the proposed analysis) and a set of branch DAGs (alternative hidden assumptions to the root DAG). All branch DAGs share a common ruleset, and mus… ▽ More Directed acyclic graphs (DAGs) are frequently used in epidemiology as a method to encode causal inference assumptions. We propose the DAGWOOD framework to bring many of those encoded assumptions to the forefront. DAGWOOD combines a root DAG (the DAG in the proposed analysis) and a set of branch DAGs (alternative hidden assumptions to the root DAG). All branch DAGs share a common ruleset, and must 1) change the root DAG, 2) be a valid DAG, and either 3a) change the minimally sufficient adjustment set or 3b) change the number of frontdoor paths. Branch DAGs comprise a list of assumptions which must be justified as negligible. We define two types of branch DAGs: exclusion branch DAGs add a single- or bidirectional pathway between two nodes in the root DAG (e.g. direct pathways and colliders), while misdirection branch DAGs represent alternative pathways that could be drawn between objects (e.g., creating a collider by reversing the direction of causation for a controlled confounder). The DAGWOOD framework 1) organizes causal model assumptions, 2) reinforces best DAG practices, 3) provides a framework for evaluation of causal models, and 4) can be used for generating causal models. △ Less

Submitted 23 November, 2021; v1 submitted 8 April, 2020; originally announced April 2020.

Optimal shattering of complex networks

Authors: Nicole Balashov, Reuven Cohen, Avieli Haber, Michael Krivelevich, Simi Haber

Abstract: We consider optimal attacks or immunization schemes on different models of random graphs. We derive bounds for the minimum number of nodes needed to be removed from a network such that all remaining components are fragments of negligible size. We obtain bounds for different regimes of random regular graphs, Erdős-Rényi random graphs, and scale free networks, some of which are tight. We show that t… ▽ More We consider optimal attacks or immunization schemes on different models of random graphs. We derive bounds for the minimum number of nodes needed to be removed from a network such that all remaining components are fragments of negligible size. We obtain bounds for different regimes of random regular graphs, Erdős-Rényi random graphs, and scale free networks, some of which are tight. We show that the performance of attacks by degree is bounded away from optimality. Finally we present a polynomial time attack algorithm and prove its optimal performance in certain cases. △ Less

Submitted 5 December, 2019; originally announced December 2019.

Journal ref: Applied Network Science 4, Article number: 99 (2019)

Steady-State Control and Machine Learning of Large-Scale Deformable Mirror Models

Authors: Aleksandar Haber

Abstract: We use Machine Learning (ML) and system identification validation approaches to estimate neural network models of large-scale Deformable Mirrors (DMs) used in Adaptive Optics (AO) systems. To obtain the training, validation, and test data sets, we simulate a realistic large-scale Finite Element (FE) model of a faceplate DM. The estimated models reproduce the input-output behavior of Vector AutoReg… ▽ More We use Machine Learning (ML) and system identification validation approaches to estimate neural network models of large-scale Deformable Mirrors (DMs) used in Adaptive Optics (AO) systems. To obtain the training, validation, and test data sets, we simulate a realistic large-scale Finite Element (FE) model of a faceplate DM. The estimated models reproduce the input-output behavior of Vector AutoRegressive with eXogenous (VARX) input models and can be used for the design of high-performance AO systems. We address the model order selection and overfitting problems. We also provide an FE based approach for computing steady-state control signals that produce the desired wavefront shape. This approach can be used to predict the steady-state DM correction performance for different actuator spacings and configurations. The presented methods are tested on models with thousands of state variables and hundreds of actuators. The numerical simulations are performed on low-cost high-performance graphic processing units and implemented using the TensorFlow machine learning framework. The used codes are available online. The approaches presented in this paper are useful for the design and optimization of high-performance DMs and AO systems. △ Less

Submitted 18 November, 2019; originally announced November 2019.

Comments: 5 pages, 6 figures

Subspace Identification of Temperature Dynamics

Authors: Aleksandar Haber

Abstract: Data-driven modeling and control of temperature dynamics in mechatronics systems and industrial processes are challenging control engineering problems. This is mainly because the temperature dynamics is inherently infinite-dimensional, nonlinear, spatially distributed, and coupled with other physical processes. Furthermore, the dominant time constants are usually long, implying that in practice du… ▽ More Data-driven modeling and control of temperature dynamics in mechatronics systems and industrial processes are challenging control engineering problems. This is mainly because the temperature dynamics is inherently infinite-dimensional, nonlinear, spatially distributed, and coupled with other physical processes. Furthermore, the dominant time constants are usually long, implying that in practice due to various economic and time constraints, we can only collect a relatively small number of data samples that can be used for data-driven modeling. Finally, since sensing and actuation of temperature dynamics are often spatially discrete, special attention needs to be given to sensor (actuator) placement and identifiability problems. Motivated by these challenges, in this manuscript, we consider the problem of data-driven modeling and validation of temperature dynamics. We have developed an experimental setup consisting of a long aluminum bar whose temperature dynamics is influenced by spatially distributed heat actuators and whose temperature is sensed by spatially distributed thermocouples. We address the noise reduction problem and perform step response and nonlinearity analyses. We combine predictor based subspace identification methods with time series analysis methods to identify a multiple-input multiple-output system model. We provide detailed treatments of model structure selection, validation, and residual analysis problems under different modeling and prediction scenarios. Our extensive experimental results show that the temperature dynamics of the experimental setup can be relatively accurately estimated by low-order models. △ Less

Submitted 6 August, 2019; originally announced August 2019.

Comments: 10 pages, 12 figures

Multicomponent Superfluids and Superconductors in Dense Nuclear and Quark Matter

Authors: Alexander Haber

Abstract: Matter at intermediate baryon densities and low temperatures is notoriously hard to tackle theoretically. Whereas lattice methods cannot cover more than rather small densities, perturbative methods are only applicable at much higher densities. The regime of intermediate chemical potential at low temperatures in the QCD-phase diagram is therefore out of reach of first-principle methods, thus we hav… ▽ More Matter at intermediate baryon densities and low temperatures is notoriously hard to tackle theoretically. Whereas lattice methods cannot cover more than rather small densities, perturbative methods are only applicable at much higher densities. The regime of intermediate chemical potential at low temperatures in the QCD-phase diagram is therefore out of reach of first-principle methods, thus we have to rely on stellar objects to investigate dense nuclear and quark matter at low temperatures. Compact stars can serve as an unique laboratory for this regime. Due to their fast rotation and high magnetic field, phenomena like hydrodynamic instabilities and the formation of vortices/flux tubes become of phenomenological interest. In my thesis, I investigate these multicomponent systems in a consistent multi-fluid treatment. By starting from a fieldtheoretical, bosonic model, the phase structure of a two-fluid system, e.g. consisting of superfluid neutrons and superconducting protons, is explored. Consequently, hydrodynamic instabilities, which might serve as trigger for pulsar glitches, are calculated for these two-fluid systems. By incorporating a gauge field and taking into account the charge of one scalar field, the influence of a superfluid on the magnetic phase structure of a superconductor is studied. In the last part, color superconducting quark matter, which can effectively be described as a multicomponent (color-)supercondcutor, is investigated. △ Less

Submitted 29 November, 2018; originally announced November 2018.

Comments: Phd thesis based on: PRD 93 no.2 (2016) 025011 (arXiv:1510.01982), EPJ Web Conf. 137 (2017) 09003 (arXiv:1612.01865), PRD 95 no.11 (2017) 116016 (arXiv:1704.01575) and J.Phys.G45 no.6 (2018) 065001 (arxiv:1712.08587)

Multi-winding flux tubes in CFL quark matter

Authors: Alexander Haber, Andreas Schmitt

Abstract: Color-flavor locked quark matter can be described as a three-component superconductor and thus shows unconventional behavior in the transition regime from type-I to type-II superconductivity. We discuss this behavior by studying magnetic line defects in a Ginzburg-Landau approach, taking into account all possible values of the three winding numbers. After a brief discussion of the defects that inc… ▽ More Color-flavor locked quark matter can be described as a three-component superconductor and thus shows unconventional behavior in the transition regime from type-I to type-II superconductivity. We discuss this behavior by studying magnetic line defects in a Ginzburg-Landau approach, taking into account all possible values of the three winding numbers. After a brief discussion of the defects that include baryon circulation we focus on pure magnetic flux tubes. We show that at strong coupling, relevant for neutron stars, type-II behavior is conceivable and the most preferred configuration has minimal total winding. Only at weak coupling we find a regime where multi-winding flux tubes are preferred, although this regime most likely requires an unrealistically large superconducting gap. △ Less

Submitted 29 November, 2018; originally announced November 2018.

Comments: 12 pages, 3 figures, proceedings contribution to "XIIIth Quark Confinement and the Hadron Spectrum", 31 July - 6 August 2018, Maynooth University, Ireland

Sparsity Preserving Optimal Control of Discretized PDE Systems

Authors: Aleksandar Haber, Michel Verhaegen

Abstract: We focus on the problem of optimal control of large-scale systems whose models are obtained by discretization of partial differential equations using the Finite Element (FE) or Finite Difference (FD) methods. The motivation for studying this pressing problem originates from the fact that the classical numerical tools used to solve low-dimensional optimal control problems are computationally infeas… ▽ More We focus on the problem of optimal control of large-scale systems whose models are obtained by discretization of partial differential equations using the Finite Element (FE) or Finite Difference (FD) methods. The motivation for studying this pressing problem originates from the fact that the classical numerical tools used to solve low-dimensional optimal control problems are computationally infeasible for large-scale systems. Furthermore, although the matrices of large-scale FE or FD models are usually sparse banded or highly structured, the optimal control solution computed using the classical methods is dense and unstructured. Consequently, it is not suitable for efficient centralized and distributed real-time implementations. We show that the a priori (sparsity) patterns of the exact solutions of the generalized Lyapunov equations for FE and FD models are banded matrices. The a priori pattern predicts the dominant non-zero entries of the exact solution. We furthermore show that for well-conditioned problems, the a priori patterns are not only banded but also sparse matrices. On the basis of these results, we develop two computationally efficient methods for computing sparse approximate solutions of generalized Lyapunov equations. Using these two methods and the inexact Newton method, we show that the solution of the generalized Riccati equation can be approximated by a banded matrix. This enables us to develop a novel computationally efficient optimal control approach that is able to preserve the sparsity of the control law. We perform extensive numerical experiments that demonstrate the effectiveness of our approach. △ Less

Submitted 20 February, 2018; v1 submitted 16 January, 2018; originally announced January 2018.

Comments: 33 pages, 8 figures, accepted to Computer Methods in Applied Mechanics and Engineering (to be published)

Journal ref: Computer Methods in Applied Mechanics and Engineering, 335 (2018), 610-630

New color-magnetic defects in dense quark matter

Authors: Alexander Haber, Andreas Schmitt

Abstract: Color-flavor locked (CFL) quark matter expels color-magnetic fields due to the Meissner effect. One of these fields carries an admixture of the ordinary abelian magnetic field and therefore flux tubes may form if CFL matter is exposed to a magnetic field, possibly in the interior of neutron stars or in quark stars. We employ a Ginzburg-Landau approach for three massless quark flavors, which takes… ▽ More Color-flavor locked (CFL) quark matter expels color-magnetic fields due to the Meissner effect. One of these fields carries an admixture of the ordinary abelian magnetic field and therefore flux tubes may form if CFL matter is exposed to a magnetic field, possibly in the interior of neutron stars or in quark stars. We employ a Ginzburg-Landau approach for three massless quark flavors, which takes into account the multi-component nature of color superconductivity. Based on the weak-coupling expressions for the Ginzburg-Landau parameters, we identify the regime where CFL is a type-II color superconductor and compute the radial profiles of different color-magnetic flux tubes. Among the configurations without baryon circulation we find a new solution that is energetically preferred over the flux tubes previously discussed in the literature in the parameter regime relevant for compact stars. Within the same setup, we also find a new defect in the 2SC phase, namely magnetic domain walls, which emerge naturally from the previously studied flux tubes if a more general ansatz for the order parameter is used. Color-magnetic defects in the interior of compact stars allow for sustained deformations of the star, potentially strong enough to produce detectable gravitational waves. △ Less

Submitted 6 August, 2020; v1 submitted 22 December, 2017; originally announced December 2017.

Comments: 29 pages, 5 figures, v2: small modifications in the text, version accepted for publication in Journal of Physics G, v3: sign error in Sec. V in front of the charge $\tilde{q}_{83}$ corrected (all results unaffected since error was exactly compensated by negative winding number $n_3$, which is positive in the new version)

State observation and sensor selection for nonlinear networks

Authors: Aleksandar Haber, Ferenc Molnar, Adilson E. Motter

Abstract: A large variety of dynamical systems, such as chemical and biomolecular systems, can be seen as networks of nonlinear entities. Prediction, control, and identification of such nonlinear networks require knowledge of the state of the system. However, network states are usually unknown, and only a fraction of the state variables are directly measurable. The observability problem concerns reconstruct… ▽ More A large variety of dynamical systems, such as chemical and biomolecular systems, can be seen as networks of nonlinear entities. Prediction, control, and identification of such nonlinear networks require knowledge of the state of the system. However, network states are usually unknown, and only a fraction of the state variables are directly measurable. The observability problem concerns reconstructing the network state from this limited information. Here, we propose a general optimization-based approach for observing the states of nonlinear networks and for optimally selecting the observed variables. Our results reveal several fundamental limitations in network observability, such as the trade-off between the fraction of observed variables and the observation length on one side, and the estimation error on the other side. We also show that owing to the crucial role played by the dynamics, purely graph- theoretic observability approaches cannot provide conclusions about one's practical ability to estimate the states. We demonstrate the effectiveness of our methods by finding the key components in biological and combustion reaction networks from which we determine the full system state. Our results can lead to the design of novel sensing principles that can greatly advance prediction and control of the dynamics of such networks. △ Less

Submitted 16 October, 2017; v1 submitted 16 June, 2017; originally announced June 2017.

Comments: Matches publication version to appear in IEEE Transactions on Control of Network Systems. 28 pages and 13 figures

Journal ref: IEEE Trans. Control Netw. Syst. 5(2), 694 (2018)

Critical magnetic fields in a superconductor coupled to a superfluid

Authors: Alexander Haber, Andreas Schmitt

Abstract: We study a superconductor that is coupled to a superfluid via density and derivative couplings. Starting from a Lagrangian for two complex scalar fields, we derive a temperature-dependent Ginzburg-Landau potential, which is then used to compute the phase diagram at nonzero temperature and external magnetic field. This includes the calculation of the critical magnetic fields for the transition to a… ▽ More We study a superconductor that is coupled to a superfluid via density and derivative couplings. Starting from a Lagrangian for two complex scalar fields, we derive a temperature-dependent Ginzburg-Landau potential, which is then used to compute the phase diagram at nonzero temperature and external magnetic field. This includes the calculation of the critical magnetic fields for the transition to an array of magnetic flux tubes, based on an approximation for the interaction between the flux tubes. We find that the transition region between type-I and type-II superconductivity changes qualitatively due to the presence of the superfluid: the phase transitions at the upper and lower critical fields in the type-II regime become first order, opening the possibility of clustered flux tube phases. These flux tube clusters may be realized in the core of neutron stars, where superconducting protons are expected to be coupled to superfluid neutrons. △ Less

Submitted 19 June, 2017; v1 submitted 5 April, 2017; originally announced April 2017.

Comments: 29 pages, 7 figures; v2: references and clarifications added, version to appear in Phys. Rev. D

Journal ref: Phys. Rev. D 95, 116016 (2017)

Mixing of charged and neutral Bose condensates at nonzero temperature and magnetic field

Authors: Alexander Haber, Andreas Schmitt

Abstract: It is expected that in the interior of compact stars a proton superconductor coexists with and couples to a neutron superfluid. Starting from a field-theoretical model for two complex scalar fields - one of which is electrically charged - we derive a Ginzburg-Landau potential which includes entrainment between the two fluids and temperature effects from thermal excitations of the two scalar fields… ▽ More It is expected that in the interior of compact stars a proton superconductor coexists with and couples to a neutron superfluid. Starting from a field-theoretical model for two complex scalar fields - one of which is electrically charged - we derive a Ginzburg-Landau potential which includes entrainment between the two fluids and temperature effects from thermal excitations of the two scalar fields and the gauge field. The Ginzburg-Landau description is then used for an analysis of the phase structure in the presence of an external magnetic field. In particular, we study the effect of the superfluid on the flux tube phase by computing the various critical magnetic fields and deriving an approximation for the flux tube interaction. As a result, we point out differences to the naive expectations from an isolated superconductor, for instance the existence of a first-order flux tube onset, resulting in a more complicated phase structure in the region between type-I and type-II superconductivity. △ Less

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

Comments: 11 pages, contribution to the proceedings of "Quark Confinement and the Hadron Spectrum XII", Aug 28 - Sep 4, 2016, in Thessaloniki, Greece; v2: small corrections in Eqs. (9) and (13), references added

Instabilities in relativistic two-component (super)fluids

Authors: Alexander Haber, Andreas Schmitt, Stephan Stetina

Abstract: We study two-fluid systems with nonzero fluid velocities and compute their sound modes, which indicate various instabilities. For the case of two zero-temperature superfluids we employ a microscopic field-theoretical model of two coupled bosonic fields, including an entrainment coupling and a non-entrainment coupling. We analyse the onset of the various instabilities systematically and point out t… ▽ More We study two-fluid systems with nonzero fluid velocities and compute their sound modes, which indicate various instabilities. For the case of two zero-temperature superfluids we employ a microscopic field-theoretical model of two coupled bosonic fields, including an entrainment coupling and a non-entrainment coupling. We analyse the onset of the various instabilities systematically and point out that the dynamical two-stream instability can only occur beyond Landau's critical velocity, i.e., in an already energetically unstable regime. A qualitative difference is found for the case of two normal fluids, where certain transverse modes suffer a two-stream instability in an energetically stable regime if there is entrainment between the fluids. Since we work in a fully relativistic setup, our results are very general and of potential relevance for (super)fluids in neutron stars and, in the non-relativistic limit of our results, in the laboratory. △ Less

Submitted 25 January, 2016; v1 submitted 7 October, 2015; originally announced October 2015.

Comments: 25 pages, 7 figures, v2: minor modifications; references added; parts of Sec III.B moved into new appendix B; version published in Phys. Rev. D

Journal ref: Phys. Rev. D 93, 025011 (2016)

Integration of Heterogeneous Modeling Languages via Extensible and Composable Language Components

Authors: Arne Haber, Markus Look, Antonio Navarro Perez, Bernhard Rumpe, Steven Völkel, Andreas Wortmann

Abstract: Effective model-driven engineering of complex systems requires to appropriately describe different specific system aspects. To this end, efficient integration of different heterogeneous modeling languages is essential. Modeling language integaration is onerous and requires in-depth conceptual and technical knowledge and ef- fort. Traditional modeling lanugage integration approches require language… ▽ More Effective model-driven engineering of complex systems requires to appropriately describe different specific system aspects. To this end, efficient integration of different heterogeneous modeling languages is essential. Modeling language integaration is onerous and requires in-depth conceptual and technical knowledge and ef- fort. Traditional modeling lanugage integration approches require language engineers to compose monolithic language aggregates for a specific task or project. Adapting these aggregates cannot be to different contexts requires vast effort and makes these hardly reusable. This contribution presents a method for the engineering of grammar-based language components that can be independently developed, are syntactically composable, and ultimately reusable. To this end, it introduces the concepts of language aggregation, language embed- ding, and language inheritance, as well as their realization in the language workbench MontiCore. The result is a generalizable, systematic, and efficient syntax-oriented composition of languages that allows the agile employment of modeling languages efficiently tailored for individual software projects. △ Less

Submitted 15 September, 2015; originally announced September 2015.

Comments: 12 pages, 11 figures. Proceedings of the 3rd International Conference on Model-Driven Engineering and Software Development. Angers, Loire Valley, France, pp. 19-31, 2015

Journal ref: Proceedings of the 3rd International Conference on Model-Driven Engineering and Software Development. Angers, Loire Valley, France, pp. 19-31, 2015

Baryon onset in a magnetic field

Authors: Alexander Haber, Florian Preis, Andreas Schmitt

Abstract: The critical baryon chemical potential for the onset of nuclear matter is a function of the vacuum mass and the binding energy. Both quantities are affected by an external magnetic field. We show within two relativistic mean-field models - including magnetic catalysis, but omitting the anomalous magnetic moment - that a magnetic field increases both the vacuum mass and the binding energy. For suff… ▽ More The critical baryon chemical potential for the onset of nuclear matter is a function of the vacuum mass and the binding energy. Both quantities are affected by an external magnetic field. We show within two relativistic mean-field models - including magnetic catalysis, but omitting the anomalous magnetic moment - that a magnetic field increases both the vacuum mass and the binding energy. For sufficiently large magnetic fields, the effect on the vacuum mass dominates and as a result the critical baryon chemical potential is increased. △ Less

Submitted 19 December, 2014; originally announced December 2014.

Comments: 5 pages, 3 figures, to appear in the proceedings of "Quark Confinement and the Hadron Spectrum XI", St Petersburg, September 8-12, 2014

Invariance property of wave scattering through disordered media

Authors: Romain Pierrat, Philipp Ambichl, Sylvain Gigan, Alexander Haber, Rémi Carminati, Stefan Rotter

Abstract: A fundamental insight in the theory of diffusive random walks is that the mean length of trajectories traversing a finite open system is independent of the details of the diffusion process. Instead, the mean trajectory length depends only on the system's boundary geometry and is thus unaffected by the value of the mean free path. Here we show that this result is rooted on a much deeper level than… ▽ More A fundamental insight in the theory of diffusive random walks is that the mean length of trajectories traversing a finite open system is independent of the details of the diffusion process. Instead, the mean trajectory length depends only on the system's boundary geometry and is thus unaffected by the value of the mean free path. Here we show that this result is rooted on a much deeper level than that of a random walk, which allows us to extend the reach of this universal invariance property beyond the diffusion approximation. Specifically, we demonstrate that an equivalent invariance relation also holds for the scattering of waves in resonant structures as well as in ballistic, chaotic or in Anderson localized systems. Our work unifies a number of specific observations made in quite diverse fields of science ranging from the movement of ants to nuclear scattering theory. Potential experimental realizations using light fields in disordered media are discussed. △ Less

Submitted 18 December, 2014; v1 submitted 25 September, 2014; originally announced September 2014.

Comments: 10 pages, 4 figures, appendix included; open access journal publication at PNAS

Journal ref: PNAS 111, 17765 (2014)

MontiArc - Architectural Modeling of Interactive Distributed and Cyber-Physical Systems

Authors: Arne Haber, Jan Oliver Ringert, Bernhard Rumpe

Abstract: This report presents MontiArc, a modeling language for the description of Component & Connector architectures. A component is a unit executing computations and/or storing data. Information flow between components is modeled via unidirectional connectors connecting typed, directed ports of the interfaces of components. Language features of the ADL MontiArc include hierarchical decomposition of comp… ▽ More This report presents MontiArc, a modeling language for the description of Component & Connector architectures. A component is a unit executing computations and/or storing data. Information flow between components is modeled via unidirectional connectors connecting typed, directed ports of the interfaces of components. Language features of the ADL MontiArc include hierarchical decomposition of components, subtyping by structural inheritance, component type definitions and reference declarations for reuse, generic component types and configurable components, syntactic sugar for connectors, and controlled implicit creation of connections and subcomponent declarations. This technical report gives an overview of the MontiArc language and is a reference for the MontiArc grammar intended to enable reuse and extension of MontiArc and MontiArc related tools. MontiArc is implemented using the DSL framework MontiCore. Available tools include an editor with syntax highlighting and code completion as well as a simulation framework with a Java code generator. △ Less

Submitted 22 September, 2014; originally announced September 2014.

Comments: 102 pages, RWTH Aachen University, Technical Report. AIB-2012-03. February 2012

Report number: AIB-2012-03

Towards Architectural Programming of Embedded Systems

Authors: Arne Haber, Jan Oliver Ringert, Bernhard Rumpe

Abstract: Integrating architectural elements with a modern programming language is essential to ensure a smooth combination of architectural design and programming. In this position statement, we motivate a combination of architectural description for distributed, asynchronously communicating systems and Java as an example for such an integration. The result is an ordinary programming language, that exhibit… ▽ More Integrating architectural elements with a modern programming language is essential to ensure a smooth combination of architectural design and programming. In this position statement, we motivate a combination of architectural description for distributed, asynchronously communicating systems and Java as an example for such an integration. The result is an ordinary programming language, that exhibits architecture, data structure and behavior within one view. Mappings or tracing between different views is unnecessary. A prototypical implementation of a compiler demonstrates the possibilities and challenges of architectural programming. △ Less

Submitted 8 September, 2014; originally announced September 2014.

Comments: 10 pages, 6 figures. Tagungsband des Dagstuhl-Workshop MBEES: Modellbasierte Entwicklung eingebetteter Systeme VI. Informatik-Bericht 2010-01, Institut für Software Systems Engineering, TU Braunschweig, 2010

Delta Modeling for Software Architectures

Authors: Arne Haber, Holger Rendel, Bernhard Rumpe, Ina Schaefer

Abstract: Architectural modeling is an integral part of modern software development. In particular, diverse systems benefit from precise architectural models since similar components can often be reused between different system variants. However, during all phases of diverse system development, system variability has to be considered and modeled by appropriate means. Delta modeling is a language-independent… ▽ More Architectural modeling is an integral part of modern software development. In particular, diverse systems benefit from precise architectural models since similar components can often be reused between different system variants. However, during all phases of diverse system development, system variability has to be considered and modeled by appropriate means. Delta modeling is a language-independent approach for modeling system variability. A set of diverse systems is represented by a core system and a set of deltas specifying modifications to the core system. In this paper, we give a first sketch of how to apply delta modeling in MontiArc, an existing architecture description language, in order to obtain an integrated modeling language for architectural variability. The developed language, MontiArc, allows the modular modeling of variable software architectures and supports proactive as well as extractive product line development. △ Less

Submitted 8 September, 2014; originally announced September 2014.

Comments: 10 pages, 6 figures. > Tagungsband des Dagstuhl-Workshop MBEES: Modellbasierte Entwicklung eingebetteter Systeme VII, fortiss GmbH München, February 2011

Hierarchical Variability Modeling for Software Architectures

Authors: Arne Haber, Holger Renel, Bernhard Rumpe, Ina Schaefer, Frank van der Linden

Abstract: Hierarchically decomposed component-based system development reduces design complexity by supporting distribution of work and component reuse. For product line development, the variability of the components to be deployed in different products has to be represented by appropriate means. In this paper, we propose hierarchical variability modeling which allows specifying component variability integr… ▽ More Hierarchically decomposed component-based system development reduces design complexity by supporting distribution of work and component reuse. For product line development, the variability of the components to be deployed in different products has to be represented by appropriate means. In this paper, we propose hierarchical variability modeling which allows specifying component variability integrated with the component hierarchy and locally to the components. Components can contain variation points determining where components may vary. Associated variants define how this variability can be realized in different component configurations. We present a meta model for hierarchical variability modeling to formalize the conceptual ideas. In order to obtain an implementation of the proposed approach together with tool support, we extend the existing architectural description language MontiArc with hierarchical variability modeling. We illustrate the presented approach using an example from the automotive systems domain. △ Less

Submitted 8 September, 2014; originally announced September 2014.

Comments: 10 pages, 9 figures. Proceedings of International Software Product Lines Conference (SPLC 2011), IEEE Computer Society, August 2011

Delta-oriented Architectural Variability Using MontiCore

Authors: Arne Haber, Thomas Kutz, Holger Rendel, Bernhard Rumpe, Ina Schaefer

Abstract: Modeling of software architectures is a fundamental part of software development processes. Reuse of software components and early analysis of software topologies allow the reduction of development costs and increases software quality. Integrating variability modeling concepts into architecture description languages (ADLs) is essential for the development of diverse software systems with high dema… ▽ More Modeling of software architectures is a fundamental part of software development processes. Reuse of software components and early analysis of software topologies allow the reduction of development costs and increases software quality. Integrating variability modeling concepts into architecture description languages (ADLs) is essential for the development of diverse software systems with high demands on software quality. In this paper, we present the integration of delta modeling into the existing ADL MontiArc. Delta modeling is a language-independent variability modeling approach supporting proactive, reactive and extractive product line development. We show how ?-MontiArc, a language for explicit modeling of architectural variability based on delta modeling, is implemented as domain-specific language (DSL) using the DSL development framework MontiCore. We also demonstrate how MontiCore's language reuse mechanisms provide efficient means to derive an implementation of ?-MontiArc tool implementation. We evaluate ?-Monti-Arc by comparing it with annotative variability modeling. △ Less

Submitted 8 September, 2014; originally announced September 2014.

Comments: 10 pages, 9 figures. ECSA '11 5th European Conference on Software Architecture: Companion Volume, ACM New York, NY, USA, Article No. 6, September 2011

Towards a Family-based Analysis of Applicability Conditions in Architectural Delta Models

Authors: Arne Haber, Thomas Kutz, Holger Rendel, Bernhard Rumpe, Ina Schaefer

Abstract: Modeling variability in software architectures is a fundamental part of software product line development. ?-MontiArc allows describing architectural variability in a modular way by a designated core architecture and a set of architectural delta models modifying the core architecture to realize other architecture variants. Delta models have to satisfy a set of applicability conditions for the defi… ▽ More Modeling variability in software architectures is a fundamental part of software product line development. ?-MontiArc allows describing architectural variability in a modular way by a designated core architecture and a set of architectural delta models modifying the core architecture to realize other architecture variants. Delta models have to satisfy a set of applicability conditions for the definedness of the architectural variants. The applicability conditions can in principle be checked by generating all possible architecture variants, which requires considering the same intermediate architectures repeatedly. In order to reuse previously computed architecture variants, we propose a family-based analysis of the applicability conditions using the concept of inverse deltas. △ Less

Submitted 8 September, 2014; originally announced September 2014.

Comments: 10 pages, 7 figures. Variability for You Proceedings of VARY InternationalWorkshop affiliated with ACM/IEEE 14th International Conference on Model Driven Engineering Languages and Systems (MODELS'11), IT University Technical Report Series TR-2011-144

Report number: TR-2011-144

Evolving Delta-oriented Software Product Line Architectures

Authors: Arne Haber, Holger Renel, Bernhard Rumpe, Ina Schaefer

Abstract: Diversity is prevalent in modern software systems. Several system variants exist at the same time in order to adapt to changing user requirements. Additionally, software systems evolve over time in order to adjust to unanticipated changes in their application environment. In modern software development, software architecture modeling is an important means to deal with system complexity by architec… ▽ More Diversity is prevalent in modern software systems. Several system variants exist at the same time in order to adapt to changing user requirements. Additionally, software systems evolve over time in order to adjust to unanticipated changes in their application environment. In modern software development, software architecture modeling is an important means to deal with system complexity by architectural decomposition. This leads to the need of architectural description languages that can represent spatial and temporal variability. In this paper, we present delta modeling of software architectures as a uniform modeling formalism for architectural variability in space and in time. In order to avoid degeneration of the product line model under system evolution, we present refactoring techniques to maintain and improve the quality of the variability model. Using a running example from the automotive domain, we evaluate our approach by carrying out a case study that compares delta modeling with annotative variability modeling. △ Less

Submitted 8 September, 2014; originally announced September 2014.

Comments: 26 pages, 6 figures

Journal ref: Large-Scale Complex IT Systems. Development, Operation and Management, 17th Monterey Workshop 2012, Oxford, UK, March 19-21, 2012. LNCS Vol. 7539, pp. 183-208

Magnetic catalysis in nuclear matter

Authors: Alexander Haber, Florian Preis, Andreas Schmitt

Abstract: A strong magnetic field enhances the chiral condensate at low temperatures. This so-called magnetic catalysis thus seeks to increase the vacuum mass of nucleons. We employ two relativistic field-theoretical models for nuclear matter, the Walecka model and an extended linear sigma model, to discuss the resulting effect on the transition between vacuum and nuclear matter at zero temperature. In both… ▽ More A strong magnetic field enhances the chiral condensate at low temperatures. This so-called magnetic catalysis thus seeks to increase the vacuum mass of nucleons. We employ two relativistic field-theoretical models for nuclear matter, the Walecka model and an extended linear sigma model, to discuss the resulting effect on the transition between vacuum and nuclear matter at zero temperature. In both models we find that the creation of nuclear matter in a sufficiently strong magnetic field becomes energetically more costly due to the heaviness of magnetized nucleons, even though it is also found that nuclear matter is more strongly bound in a magnetic field. Our results are potentially important for dense nuclear matter in compact stars, especially since previous studies in the astrophysical context have always ignored the contribution of the magnetized Dirac sea and thus the effect of magnetic catalysis. △ Less

Submitted 10 December, 2014; v1 submitted 1 September, 2014; originally announced September 2014.

Comments: 19 pages, 5 figures, v3: version to appear in Phys. Rev. D

Journal ref: Phys. Rev. D 90, 125036 (2014)

Engineering Delta Modeling Languages

Authors: Arne Haber, Katrin Hölldobler, Carsten Kolassa, Markus Look, Klaus Müller, Bernhard Rumpe, Ina Schaefer

Abstract: Delta modeling is a modular, yet flexible approach to capture spatial and temporal variability by explicitly representing the differences between system variants or versions. The conceptual idea of delta modeling is language-independent. But, in order to apply delta modeling for a concrete language, so far, a delta language had to be manually developed on top of the base language leading to a larg… ▽ More Delta modeling is a modular, yet flexible approach to capture spatial and temporal variability by explicitly representing the differences between system variants or versions. The conceptual idea of delta modeling is language-independent. But, in order to apply delta modeling for a concrete language, so far, a delta language had to be manually developed on top of the base language leading to a large variety of heterogeneous language concepts. In this paper, we present a process that allows deriving a delta language from the grammar of a given base language. Our approach relies on an automatically generated language extension that can be manually adapted to meet domain-specific needs. We illustrate our approach using delta modeling on a textual variant of statecharts. △ Less

Submitted 25 August, 2014; originally announced August 2014.

Comments: 10 pages, 8 figures. Proceedings of the 17th International Software Product Line Conference, Tokyo, September 2013, pp.22-31, ACM, 2013

ACM Class: D.2; D.2.2; D.2.3

Journal ref: Proceedings of the 17th International Software Product Line Conference, Tokyo, September 2013, pp.22-31, ACM, 2013

First-Class Variability Modeling in Matlab/Simulink

Authors: Arne Haber, Carsten Kolassa, Peter Manhart, Pedram Mir Seyed Nazari, Bernhard Rumpe, Ina Schaefer

Abstract: Modern cars exist in an vast number of variants. Thus, variability has to be dealt with in all phases of the development process, in particular during model-based development of software-intensive functionality using Matlab/Simulink. Currently, variability is often encoded within a functional model leading to so called 150%-models which easily become very complex and do not scale for larger produc… ▽ More Modern cars exist in an vast number of variants. Thus, variability has to be dealt with in all phases of the development process, in particular during model-based development of software-intensive functionality using Matlab/Simulink. Currently, variability is often encoded within a functional model leading to so called 150%-models which easily become very complex and do not scale for larger product lines. To counter these problems, we propose a modular variability modeling approach for Matlab/Simulink based on the concept of delta modeling [8, 9, 24]. A functional variant is described by a delta encapsulating a set of modifications. A sequence of deltas can be applied to a core product to derive the desired variant. We present a prototypical implementation, which is integrated into Matlab/Simulink and offers graphical editing of delta models. △ Less

Submitted 25 August, 2014; originally announced August 2014.

Comments: 8 pages, 8 figures. Workshop on Variability Modelling of Software-intensive Systems. VaMoS 2013 in Pisa, Italy

ACM Class: D.2.6; D.2.2

Journal ref: Proceedings of the Seventh International Workshop on Variability Modelling of Software-intensive Systems, 23.-25.1.2013, pp. 11-18, ACM, New York, NY, USA. 2013

Sparse solution of the Lyapunov equation for large-scale interconnected systems

Authors: Aleksandar Haber, Michel Verhaegen

Abstract: We consider the problem of computing an approximate banded solution of the continuous-time Lyapunov equation $\underline{A}\underline{X}+\underline{X}\underline{A}^{T}=\underline{P}$, where the coefficient matrices $\underline{A}$ and $\underline{P}$ are large, symmetric banded matrices. The (sparsity) pattern of $\underline{A}$ describes the interconnection structure of a large-scale interconnect… ▽ More We consider the problem of computing an approximate banded solution of the continuous-time Lyapunov equation $\underline{A}\underline{X}+\underline{X}\underline{A}^{T}=\underline{P}$, where the coefficient matrices $\underline{A}$ and $\underline{P}$ are large, symmetric banded matrices. The (sparsity) pattern of $\underline{A}$ describes the interconnection structure of a large-scale interconnected system. Recently, it has been shown that the entries of the solution $\underline{X}$ are spatially localized or decaying away from a banded pattern. We show that the decay of the entries of $\underline{X}$ is faster if the condition number of $\underline{A}$ is smaller. By exploiting the decay of entries of $\underline{X}$, we develop two computationally efficient methods for approximating $\underline{X}$ by a banded matrix. For a well-conditioned and sparse banded $\underline{A}$, the computational and memory complexities of the methods scale linearly with the state dimension. We perform extensive numerical experiments that confirm this, and that demonstrate the effectiveness of the developed methods. The methods proposed in this paper can be generalized to (sparsity) patterns of $\underline{A}$ and $\underline{P}$ that are more general than banded matrices. The results of this paper open the possibility for developing computationally efficient methods for approximating the solution of the large-scale Riccati equation by a sparse matrix. △ Less

Submitted 16 May, 2016; v1 submitted 18 August, 2014; originally announced August 2014.

Comments: Accepted in Automatica, Final version, 16 pages, 10 figures

Sparse approximate inverses of Gramians and impulse response matrices of large-scale interconnected systems

Authors: Aleksandar Haber, Michel Verhaegen

Abstract: In this paper we show that inverses of well-conditioned, finite-time Gramians and impulse response matrices of large-scale interconnected systems described by sparse state-space models, can be approximated by sparse matrices. The approximation methodology established in this paper opens the door to the development of novel methods for distributed estimation, identification and control of large-sca… ▽ More In this paper we show that inverses of well-conditioned, finite-time Gramians and impulse response matrices of large-scale interconnected systems described by sparse state-space models, can be approximated by sparse matrices. The approximation methodology established in this paper opens the door to the development of novel methods for distributed estimation, identification and control of large-scale interconnected systems. The novel estimators (controllers) compute local estimates (control actions) simply as linear combinations of inputs and outputs (states) of local subsystems. The size of these local data sets essentially depends on the condition number of the finite-time observability (controllability) Gramian. Furthermore, the developed theory shows that the sparsity patterns of the system matrices of the distributed estimators (controllers) are primarily determined by the sparsity patterns of state-space matrices of large-scale systems. The computational and memory complexity of the approximation algorithms are $O(N)$, where $N$ is the number of local subsystems of the interconnected system. Consequently, the proposed approximation methodology is computationally feasible for interconnected systems with an extremely large number of local subsystems. △ Less

Submitted 11 May, 2014; originally announced May 2014.

Comments: 11 pages, 4 figures. Submitted to IEEE Transactions on Automatic Control

Subspace identification of large-scale interconnected systems

Authors: Aleksandar Haber, Michel Verhaegen

Abstract: We propose a decentralized subspace algorithm for identification of large-scale, interconnected systems that are described by sparse (multi) banded state-space matrices. First, we prove that the state of a local subsystem can be approximated by a linear combination of inputs and outputs of the local subsystems that are in its neighborhood. Furthermore, we prove that for interconnected systems with… ▽ More We propose a decentralized subspace algorithm for identification of large-scale, interconnected systems that are described by sparse (multi) banded state-space matrices. First, we prove that the state of a local subsystem can be approximated by a linear combination of inputs and outputs of the local subsystems that are in its neighborhood. Furthermore, we prove that for interconnected systems with well-conditioned, finite-time observability Gramians (or observability matrices), the size of this neighborhood is relatively small. On the basis of these results, we develop a subspace identification algorithm that identifies a state-space model of a local subsystem from the local input-output data. Consequently, the developed algorithm is computationally feasible for interconnected systems with a large number of local subsystems. Numerical results confirm the effectiveness of the new identification algorithm. △ Less

Submitted 14 February, 2014; v1 submitted 19 September, 2013; originally announced September 2013.

Comments: revised version, conditionally accepted for publication in IEEE Transaction on Automatic Control

Two-Dimensional Helioseismic Power, Phase, and Coherence Spectra of {\it Solar Dynamics Observatory} Photospheric and Chromospheric Observables

Authors: Rachel Howe, Kiran Jain, Richard S. Bogart, Deborah A. Haber, Charles S. Baldner

Abstract: While the {\it Helioseismic and Magnetic Imager} (HMI) onboard the {\it Solar Dynamics Observatory} (SDO) provides Doppler velocity [$V$], continuum intensity [$I_C$], and line-depth [$Ld$] observations, each of which is sensitive to the five-minute acoustic spectrum, the {\it Atmospheric Imaging Array} (AIA) also observes at wavelengths -- specifically the 1600 and 1700 Angstrom bands -- that are… ▽ More While the {\it Helioseismic and Magnetic Imager} (HMI) onboard the {\it Solar Dynamics Observatory} (SDO) provides Doppler velocity [$V$], continuum intensity [$I_C$], and line-depth [$Ld$] observations, each of which is sensitive to the five-minute acoustic spectrum, the {\it Atmospheric Imaging Array} (AIA) also observes at wavelengths -- specifically the 1600 and 1700 Angstrom bands -- that are partly formed in the upper photosphere and have good sensitivity to acoustic modes. In this article we consider the characteristics of the spatio--temporal Fourier spectra in AIA and HMI observables for a 15-degree region around NOAA Active Region 11072. We map the spatio--temporal-power distribution for the different observables and the HMI Line Core [$I_L$], or Continuum minus Line Depth, and the phase and coherence functions for selected observable pairs, as a function of position and frequency. Five-minute oscillation power in all observables is suppressed in the sunspot and also in plage areas. Above the acoustic cut-off frequency, the behaviour is more complicated: power in HMI $I_C$ is still suppressed in the presence of surface magnetic fields, while power in HMI $I_L$ and the AIA bands is suppressed in areas of surface field but enhanced in an extended area around the active region, and power in HMI $V$ is enhanced in a narrow zone around strong-field concentrations and suppressed in a wider surrounding area. The relative phase of the observables, and their cross-coherence functions, are also altered around the active region. These effects may help us to understand the interaction of waves and magnetic fields in the different layers of the photosphere, and will need to be taken into account in multi-wavelength local helioseismic analysis of active regions. △ Less

Submitted 8 August, 2012; originally announced August 2012.

Comments: 18 pages, 15 figures, to be published in Solar Physics

Journal ref: 2012 Solar Phys. 281, 533