Showing 1–50 of 145 results for author: Hof, M

A flexible model for Record Linkage

Authors: Kayané Robach, Stéphanie van der Pas, Mark van de Wiel, Michel H. Hof

Abstract: Combining data from various sources empowers researchers to explore innovative questions, for example those raised by conducting healthcare monitoring studies. However, the lack of a unique identifier often poses challenges. Record linkage procedures determine whether pairs of observations collected on different occasions belong to the same individual using partially identifying variables (e.g. bi… ▽ More Combining data from various sources empowers researchers to explore innovative questions, for example those raised by conducting healthcare monitoring studies. However, the lack of a unique identifier often poses challenges. Record linkage procedures determine whether pairs of observations collected on different occasions belong to the same individual using partially identifying variables (e.g. birth year, postal code). Existing methodologies typically involve a compromise between computational efficiency and accuracy. Traditional approaches simplify this task by condensing information, yet they neglect dependencies among linkage decisions and disregard the one-to-one relationship required to establish coherent links. Modern approaches offer a comprehensive representation of the data generation process, at the expense of computational overhead and reduced flexibility. We propose a flexible method, that adapts to varying data complexities, addressing registration errors and accommodating changes of the identifying information over time. Our approach balances accuracy and scalability, estimating the linkage using a Stochastic Expectation Maximisation algorithm on a latent variable model. We illustrate the ability of our methodology to connect observations using large real data applications and demonstrate the robustness of our model to the linking variables quality in a simulation study. The proposed algorithm FlexRL is implemented and available in an open source R package. △ Less

Submitted 9 July, 2024; originally announced July 2024.

The 23.01 release of Cloudy

Authors: Chamani M. Gunasekera, Peter A. M. van Hoof, Marios Chatzikos, Gary J. Ferland

Abstract: We announce the C23.01 update of Cloudy. This corrects a simple coding error, present since $\sim$ 1990, in one routine that required a conversion from the line-center to the mean normalization of the Ly$αあるふぁ$ optical depth. This affects the destruction of H I Ly$αあるふぁ$ by background opacities. Its largest effect is upon the Ly$αあるふぁ$ intensity in high-ionization dusty clouds, where the predicted intensity i… ▽ More We announce the C23.01 update of Cloudy. This corrects a simple coding error, present since $\sim$ 1990, in one routine that required a conversion from the line-center to the mean normalization of the Ly$αあるふぁ$ optical depth. This affects the destruction of H I Ly$αあるふぁ$ by background opacities. Its largest effect is upon the Ly$αあるふぁ$ intensity in high-ionization dusty clouds, where the predicted intensity is now up to three times stronger. Other properties that depend on Ly$αあるふぁ$ destruction, such as grain infrared emission, change in response. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 4 pages, 1 figure

MELCHIORS: The Mercator Library of High Resolution Stellar Spectroscopy

Authors: P. Royer, T. Merle, K. Dsilva, S. Sekaran, H. Van Winckel, Y. Frémat, M. Van der Swaelmen, S. Gebruers, A. Tkachenko, M. Laverick, M. Dirickx, G. Raskin, H. Hensberge, M. Abdul-Masih, B. Acke, M. L. Alonso, S. Bandhu Mahato, P. G. Beck, N. Behara, S. Bloemen, B. Buysschaert, N. Cox, J. Debosscher, P. De Cat, P. Degroote , et al. (49 additional authors not shown)

Abstract: Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundr… ▽ More Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundred objects from the southern hemisphere (UVES-POP) . We aim to extend the sample, offering a finer coverage of effective temperatures and surface gravity with a uniform collection of spectra obtained in the northern hemisphere. Between 2010 and 2020, we acquired several thousand echelle spectra of bright stars with the Mercator-HERMES spectrograph located in the Roque de Los Muchachos Observatory in La Palma, whose pipeline offers high-quality data reduction products. We have also developed methods to correct for the instrumental response in order to approach the true shape of the spectral continuum. Additionally, we have devised a normalisation process to provide a homogeneous normalisation of the full spectral range for most of the objects. We present a new spectral library consisting of 3256 spectra covering 2043 stars. It combines high signal-to-noise and high spectral resolution over the entire range of effective temperatures and luminosity classes. The spectra are presented in four versions: raw, corrected from the instrumental response, with and without correction from the atmospheric molecular absorption, and normalised (including the telluric correction). △ Less

Submitted 5 November, 2023; originally announced November 2023.

Comments: 17 pages, 18 figures Preview and access to the library: https://www.royer.se/melchiors.html

Journal ref: A&A 681, A107 (2024)

A framework for interpretation and testing of sparse canonical correlations

Authors: Nuria Senar, Mark van de Wiel, Aeilko Zwinderman, Michel Hof

Abstract: In clinical and biomedical research, multiple high-dimensional datasets are nowadays routinely collected from omics and imaging devices. Multivariate methods, such as Canonical Correlation Analysis (CCA), integrate two (or more) datasets to discover and understand underlying biological mechanisms. For an explorative method like CCA, interpretation is key. We present a sparse CCA method based on so… ▽ More In clinical and biomedical research, multiple high-dimensional datasets are nowadays routinely collected from omics and imaging devices. Multivariate methods, such as Canonical Correlation Analysis (CCA), integrate two (or more) datasets to discover and understand underlying biological mechanisms. For an explorative method like CCA, interpretation is key. We present a sparse CCA method based on soft-thresholding that produces near-orthogonal components, allows for browsing over various sparsity levels, and permutation-based hypothesis testing. Our soft-thresholding approach avoids tuning of a penalty parameter. Such tuning is computationally burdensome and may render unintelligible results. In addition, unlike alternative approaches, our method is less dependent on the initialisation. We examined the performance of our approach with simulations and illustrated its use on real cancer genomics data from drug sensitivity screens. Moreover, we compared its performance to Penalised Matrix Analysis (PMA), which is a popular alternative of sparse CCA with a focus on yielding interpretable results. Compared to PMA, our method offers improved interpretability of the results, while not compromising, or even improving, signal discovery. he software and simulation framework are available at https://github.com/nuria-sv/toscca. △ Less

Submitted 4 October, 2023; v1 submitted 3 October, 2023; originally announced October 2023.

JWST observations of the Ring Nebula (NGC 6720): I. Imaging of the rings, globules, and arcs

Authors: R. Wesson, Mikako Matsuura, Albert A. Zijlstra, Kevin Volk, Patrick J. Kavanagh, Guillermo García-Segura, I. McDonald, Raghvendra Sahai, M. J. Barlow, Nick L. J. Cox, Jeronimo Bernard-Salas, Isabel Aleman, Jan Cami, Nicholas Clark, Harriet L. Dinerstein, K. Justtanont, Kyle F. Kaplan, A. Manchado, Els Peeters, Griet C. Van de Steene, Peter A. M. van Hoof

Abstract: We present JWST images of the well-known planetary nebula NGC 6720 (the Ring Nebula), covering wavelengths from 1.6$μみゅー$m to 25 $μみゅー$m. The bright shell is strongly fragmented with some 20 000 dense globules, bright in H$_2$, with a characteristic diameter of 0.2 arcsec and density $n_{\rm H} \sim 10^5$-$10^6$ cm$^{-3}$. The shell contains a thin ring of polycyclic aromatic hydrocarbon (PAH) emission.… ▽ More We present JWST images of the well-known planetary nebula NGC 6720 (the Ring Nebula), covering wavelengths from 1.6$μみゅー$m to 25 $μみゅー$m. The bright shell is strongly fragmented with some 20 000 dense globules, bright in H$_2$, with a characteristic diameter of 0.2 arcsec and density $n_{\rm H} \sim 10^5$-$10^6$ cm$^{-3}$. The shell contains a thin ring of polycyclic aromatic hydrocarbon (PAH) emission. H$_2$ is found throughout the shell and in the halo. H$_2$ in the halo may be located on the swept-up walls of a biconal polar flow. The central cavity is shown to be filled with high ionization gas and shows two linear structures. The central star is located 2 arcsec from the emission centroid of the cavity and shell. Linear features (`spikes') extend outward from the ring, pointing away from the central star. Hydrodynamical simulations are shown which reproduce the clumping and possibly the spikes. Around ten low-contrast, regularly spaced concentric arc-like features are present; they suggest orbital modulation by a low-mass companion with a period of about 280 yr. A previously known much wider companion is located at a projected separation of about 15 000 au; we show that it is an M2-M4 dwarf. The system is therefore a triple star. These features, including the multiplicity, are similar to those seen in the Southern Ring Nebula (NGC 3132) and may be a common aspect of such nebulae. △ Less

Submitted 21 August, 2023; v1 submitted 17 August, 2023; originally announced August 2023.

Comments: 25 pages, 23 figures. Submitted to Monthly Notices of the Royal Astronomical Society. Corrected typos in metadata

The heart of Sakurai's Object revealed by ALMA

Authors: Daniel Tafoya, Peter A. M. van Hoof, Jesus A. Toala, Griet Van de Steene, Suzanna Randall, Ramlal Unnikrishnan, Stefan Kimeswenger, Marcin Hajduk, Daniela Barria, Albert Zijlstra

Abstract: We present high angular-resolution observations of Sakurai's object using the Atacama Large Millimeter Array, shedding new light on its morpho-kinematical structure. The millimetre continuum emission, observed at an angular resolution of 20 milliarcsec (corresponding to 70 AUえーゆー), reveals a bright compact central component whose spectral index indicates that it composed of amorphous carbon dust. Base… ▽ More We present high angular-resolution observations of Sakurai's object using the Atacama Large Millimeter Array, shedding new light on its morpho-kinematical structure. The millimetre continuum emission, observed at an angular resolution of 20 milliarcsec (corresponding to 70 AUえーゆー), reveals a bright compact central component whose spectral index indicates that it composed of amorphous carbon dust. Based on these findings, we conclude that this emission traces the previously suggested dust disc observed in mid-infrared observations. Therefore, our observations provide the first direct imaging of such a disc. The H$^{12}$CN($J$=4$\rightarrow$3) line emission, observed at an angular resolution of 300 milliarcsec (corresponding to 1000 AUえーゆー), displays bipolar structure with a north-south velocity gradient. From the position-velocity diagram of this emission we identify the presence of an expanding disc and a bipolar molecular outflow. The inclination of the disc is determined to be $i$=72$^\circ$. The derived values for the de-projected expansion velocity and the radius of the disc are $v_{\rm exp}$=53 km s$^{-1}$ and $R$=277 AUえーゆー, respectively. On the other hand, the de-projected expansion velocity of the bipolar outflow detected in the H$^{12}$CN($J$=4$\rightarrow$3) emission of approximately 1000 km s$^{-1}$. We propose that the molecular outflow has an hourglass morphology with an opening angle of around 60$^{\circ}$. Our observations unambiguously show that an equatorial disc and bipolar outflows formed in Sakurai's object in less than 30 years after the born-again event occurred, providing important constraints for future modelling efforts of this phenomenon. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: 9 pages, 5 figures. Accepted for publication as a Letter in Astronomy and Astrophysics

The 2023 release of Cloudy

Authors: Marios Chatzikos, Stefano Bianchi, Francesco Camilloni, Priyanka Chakraborty, Chamani M. Gunasekera, Francisco Guzmán, Jonathan S. Milby, Arnab Sarkar, Gargi Shaw, Peter A. M. van Hoof, Gary J. Ferland

Abstract: We describe the 2023 release of the spectral synthesis code Cloudy. Since the previous major release, migrations of our online services motivated us to adopt git as our version control system. This change alone led us to adopt an annual release scheme, accompanied by a short release paper, the present being the inaugural. Significant changes to our atomic and molecular data have improved the accur… ▽ More We describe the 2023 release of the spectral synthesis code Cloudy. Since the previous major release, migrations of our online services motivated us to adopt git as our version control system. This change alone led us to adopt an annual release scheme, accompanied by a short release paper, the present being the inaugural. Significant changes to our atomic and molecular data have improved the accuracy of Cloudy predictions: we have upgraded our instance of the Chianti database from version 7 to 10; our H- and He-like collisional rates to improved theoretical values; our molecular data to the most recent LAMDA database, and several chemical reaction rates to their most recent UDfA and KiDA values. Finally, we describe our progress on upgrading Cloudy's capabilities to meet the requirements of the X-ray microcalorimeters aboard the upcoming XRISM and Athena missions, and outline future development that will make Cloudy of use to the X-ray community. △ Less

Submitted 11 August, 2023; originally announced August 2023.

Comments: 18 pages, 10 figures, accepted for publication in RMxAA

Indexing and Partitioning the Spatial Linear Model for Large Data Sets

Authors: Jay M. Ver Hoef, Michael Dumelle, Matt Higham, Erin E. Peterson, Daniel J. Isaak

Abstract: We consider four main goals when fitting spatial linear models: 1) estimating covariance parameters, 2) estimating fixed effects, 3) kriging (making point predictions), and 4) block-kriging (predicting the average value over a region). Each of these goals can present different challenges when analyzing large spatial data sets. Current research uses a variety of methods, including spatial basis fun… ▽ More We consider four main goals when fitting spatial linear models: 1) estimating covariance parameters, 2) estimating fixed effects, 3) kriging (making point predictions), and 4) block-kriging (predicting the average value over a region). Each of these goals can present different challenges when analyzing large spatial data sets. Current research uses a variety of methods, including spatial basis functions (reduced rank), covariance tapering, etc, to achieve these goals. However, spatial indexing, which is very similar to composite likelihood, offers some advantages. We develop a simple framework for all four goals listed above by using indexing to create a block covariance structure and nearest-neighbor predictions while maintaining a coherent linear model. We show exact inference for fixed effects under this block covariance construction. Spatial indexing is very fast, and simulations are used to validate methods and compare to another popular method. We study various sample designs for indexing and our simulations showed that indexing leading to spatially compact partitions are best over a range of sample sizes, autocorrelation values, and generating processes. Partitions can be kept small, on the order of 50 samples per partition. We use nearest-neighbors for kriging and block kriging, finding that 50 nearest-neighbors is sufficient. In all cases, confidence intervals for fixed effects, and prediction intervals for (block) kriging, have appropriate coverage. Some advantages of spatial indexing are that it is available for any valid covariance matrix, can take advantage of parallel computing, and easily extends to non-Euclidean topologies, such as stream networks. We use stream networks to show how spatial indexing can achieve all four goals, listed above, for very large data sets, in a matter of minutes, rather than days, for an example data set. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Marginal Inference for Hierarchical Generalized Linear Mixed Models with Patterned Covariance Matrices Using the Laplace Approximation

Authors: Jay M. Ver Hoef, Eryn Blagg, Michael Dumelle, Philip M. Dixon, Dale L. Zimmerman, Paul Conn

Abstract: Using a hierarchical construction, we develop methods for a wide and flexible class of models by taking a fully parametric approach to generalized linear mixed models with complex covariance dependence. The Laplace approximation is used to marginally estimate covariance parameters while integrating out all fixed and latent random effects. The Laplace approximation relies on Newton-Raphson updates,… ▽ More Using a hierarchical construction, we develop methods for a wide and flexible class of models by taking a fully parametric approach to generalized linear mixed models with complex covariance dependence. The Laplace approximation is used to marginally estimate covariance parameters while integrating out all fixed and latent random effects. The Laplace approximation relies on Newton-Raphson updates, which also leads to predictions for the latent random effects. We develop methodology for complete marginal inference, from estimating covariance parameters and fixed effects to making predictions for unobserved data, for any patterned covariance matrix in the hierarchical generalized linear mixed models framework. The marginal likelihood is developed for six distributions that are often used for binary, count, and positive continuous data, and our framework is easily extended to other distributions. The methods are illustrated with simulations from stochastic processes with known parameters, and their efficacy in terms of bias and interval coverage is shown through simulation experiments. Examples with binary and proportional data on election results, count data for marine mammals, and positive-continuous data on heavy metal concentration in the environment are used to illustrate all six distributions with a variety of patterned covariance structures that include spatial models (e.g., geostatistical and areal models), time series models (e.g., first-order autoregressive models), and mixtures with typical random intercepts based on grouping. △ Less

Submitted 4 May, 2023; originally announced May 2023.

On a Canonical Distributed Controller in the Behavioral Framework

Authors: Tom R. V. Steentjes, Mircea Lazar, Paul M. J. Van den Hof

Abstract: Control in a classical transfer function or state-space setting typically views a controller as a signal processor: sensor outputs are mapped to actuator inputs. In behavioral system theory, control is simply viewed as interconnection; the interconnection of a plant with a controller. In this paper we consider the problem of control of interconnected systems in a behavioral setting. The behavioral… ▽ More Control in a classical transfer function or state-space setting typically views a controller as a signal processor: sensor outputs are mapped to actuator inputs. In behavioral system theory, control is simply viewed as interconnection; the interconnection of a plant with a controller. In this paper we consider the problem of control of interconnected systems in a behavioral setting. The behavioral setting is especially fit for modelling interconnected systems, because it allows for the interconnection of subsystems without imposing inputs and outputs. We introduce a so-called canonical distributed controller that implements a given interconnected behavior that is desired, provided that necessary and sufficient conditions hold true. The controller design can be performed in a decentralized manner, in the sense that a local controller only depends on the local system behavior. Regularity of interconnections is an important property in behavioral control that yields feedback interconnections. We provide conditions under which the interconnection of this distributed controller with the plant is regular. Furthermore, we show that the interconnections of subsystems of the canonical distributed controller are regular if and only if the interconnections of the plant and desired behavior are regular. △ Less

Submitted 10 April, 2023; v1 submitted 12 October, 2022; originally announced October 2022.

Recombination of hot ionized nebulae: The old planetary nebula around V4334 Sgr (Sakurai's star)

Authors: Martin Reichel, Stefan Kimeswenger, Peter A. M. van Hoof, Albert A. Zijlstra, Daniela Barria, Marcin Hajduk, Griet C. Van de Steene, Daniel Tafoya

Abstract: After becoming ionized, low-density astrophysical plasmas will begin a process of slow recombination. Models for this still have significant uncertainties. The recombination cannot normally be observed in isolation, because the ionization follows the evolutionary time scale of the ionizing source. Laboratory experiments are unable to reach the appropriate conditions because of the required very lo… ▽ More After becoming ionized, low-density astrophysical plasmas will begin a process of slow recombination. Models for this still have significant uncertainties. The recombination cannot normally be observed in isolation, because the ionization follows the evolutionary time scale of the ionizing source. Laboratory experiments are unable to reach the appropriate conditions because of the required very long time scales. The extended nebula around the very late helium flash (VLTP) star V4334 Sgr provides a unique laboratory for this kind of study. The sudden loss of the ionizing UV radiation after the VLTP event has allowed the nebula to recombine free from other influences. More than 290 long slit spectra taken with FORS1/2 at the ESO VLT between 2007 and 2022 are used to follow the time evolution of lines of H, He, N, S, O, Ar. Hydrogen and helium lines, representing most of the ionized mass, do not show significant changes. A small increase is seen in [N II] (+2.8 %/yr; significance 2.7 sigma), while we see a decrease in [O III] (-1.96 %/yr; 2.0 sigma). The [S II] lines show a change of +3.0 %/yr; 1.6 sigma). The lines of [S III] and of Ar III] show no significant change. For [S III], the measurement differs from the predicted decrease by 4.5 sigma. A possible explanation is that the fraction of [S IV] and higher is larger than expected. Such an effect could provide a potential solution for the sulfur anomaly in planetary nebulae. △ Less

Submitted 8 September, 2022; originally announced September 2022.

Comments: 13 pages, 8 figures, 3 tables; Accepted for Publication in ApJ

Learning linear modules in a dynamic network with missing node observations

Authors: Karthik R. Ramaswamy, Giulio Bottegal, Paul M. J. Van den Hof

Abstract: In order to identify a system (module) embedded in a dynamic network, one has to formulate a multiple-input estimation problem that necessitates certain nodes to be measured and included as predictor inputs. However, some of these nodes may not be measurable in many practical cases due to sensor selection and placement issues. This may result in biased estimates of the target module. Furthermore,… ▽ More In order to identify a system (module) embedded in a dynamic network, one has to formulate a multiple-input estimation problem that necessitates certain nodes to be measured and included as predictor inputs. However, some of these nodes may not be measurable in many practical cases due to sensor selection and placement issues. This may result in biased estimates of the target module. Furthermore, the identification problem associated with the multiple-input structure may require determining a large number of parameters that are not of particular interest to the experimenter, with increased computational complexity in large-sized networks. In this paper, we tackle these problems by using a data augmentation strategy that allows us to reconstruct the missing node measurements and increase the accuracy of the estimated target module. To this end, we develop a system identification method using regularized kernel-based methods coupled with approximate inference methods. Keeping a parametric model for the module of interest, we model the other modules as Gaussian Processes (GP) with a kernel given by the so-called stable spline kernel. An Empirical Bayes (EB) approach is used to estimate the parameters of the target module. The related optimization problem is solved using an Expectation-Maximization (EM) method, where we employ a Markov-chain Monte Carlo (MCMC) technique to reconstruct the unknown missing node information and the network dynamics. Numerical simulations on dynamic network examples illustrate the potentials of the developed method. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 17 pages, 7 Figures. Draft paper submitted to Automatica

Informativity conditions for data-driven control based on input-state data and polyhedral cross-covariance noise bounds

Authors: Tom R. V. Steentjes, Mircea Lazar, Paul M. J. Van den Hof

Abstract: Modeling and control of dynamical systems rely on measured data, which contains information about the system. Finite data measurements typically lead to a set of system models that are unfalsified, i.e., that explain the data. The problem of data-informativity for stabilization or control with quadratic performance is concerned with the existence of a controller that stabilizes all unfalsified sys… ▽ More Modeling and control of dynamical systems rely on measured data, which contains information about the system. Finite data measurements typically lead to a set of system models that are unfalsified, i.e., that explain the data. The problem of data-informativity for stabilization or control with quadratic performance is concerned with the existence of a controller that stabilizes all unfalsified systems or achieves a desired quadratic performance. Recent results in the literature provide informativity conditions for control based on input-state data and ellipsoidal noise bounds, such as energy or magnitude bounds. In this paper, we consider informativity of input-state data for control where noise bounds are defined through the cross-covariance of the noise with respect to an instrumental variable; bounds that were introduced originally as a noise characterization in parameter bounding identification. The considered cross-covariance bounds are defined by a finite number of hyperplanes, which induce a (possibly unbounded) polyhedral set of unfalsified systems. We provide informativity conditions for input-state data with polyhedral cross-covariance bounds for stabilization and $\mathcal{H}_2$/$\mathcal{H}_\infty$ control through vertex/half-space representations of the polyhedral set of unfalsified systems. △ Less

Submitted 18 February, 2022; originally announced February 2022.

SSNbayes: An R package for Bayesian spatio-temporal modelling on stream networks

Authors: Edgar Santos-Fernandez, Jay M. Ver Hoef, James M. McGree, Daniel J. Isaak, Kerrie Mengersen, Erin E. Peterson

Abstract: Spatio-temporal models are widely used in many research areas from ecology to epidemiology. However, most covariance functions describe spatial relationships based on Euclidean distance only. In this paper, we introduce the R package SSNbayes for fitting Bayesian spatio-temporal models and making predictions on branching stream networks. SSNbayes provides a linear regression framework with multipl… ▽ More Spatio-temporal models are widely used in many research areas from ecology to epidemiology. However, most covariance functions describe spatial relationships based on Euclidean distance only. In this paper, we introduce the R package SSNbayes for fitting Bayesian spatio-temporal models and making predictions on branching stream networks. SSNbayes provides a linear regression framework with multiple options for incorporating spatial and temporal autocorrelation. Spatial dependence is captured using stream distance and flow connectivity while temporal autocorrelation is modelled using vector autoregression approaches. SSNbayes provides the functionality to make predictions across the whole network, compute exceedance probabilities and other probabilistic estimates such as the proportion of suitable habitat. We illustrate the functionality of the package using a stream temperature dataset collected in Idaho, USA. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Excitation allocation for generic identifiability of linear dynamic networks with fixed modules

Authors: H. J. Dreef, S. Shi, X. Cheng, M. C. F. Donkers, P. M. J. Van den Hof

Abstract: Identifiability of linear dynamic networks requires the presence of a sufficient number of external excitation signals. The problem of allocating a minimal number of external signals for guaranteeing generic network identifiability has been recently addressed in the literature. Here we will extend that work by explicitly incorporating the situation that some network modules are known, and thus are… ▽ More Identifiability of linear dynamic networks requires the presence of a sufficient number of external excitation signals. The problem of allocating a minimal number of external signals for guaranteeing generic network identifiability has been recently addressed in the literature. Here we will extend that work by explicitly incorporating the situation that some network modules are known, and thus are fixed in the parametrized model set. The graphical approach introduced earlier is extended to this situation, showing that the presence of fixed modules reduces the required number of external signals. An algorithm is presented that allocates the external signals in a systematic fashion. △ Less

Submitted 13 May, 2022; v1 submitted 22 January, 2022; originally announced January 2022.

Journal ref: IEEE Control Systems Letters, Vol. 6, pp. 2587-2592, 2022

Identifiability in Dynamic Acyclic Networks with Partial Excitations and Measurements

Authors: Xiaodong Cheng, Shengling Shi, Ioannis Lestas, Paul M. J. Van den Hof

Abstract: This paper deals with dynamic networks in which the causality relations between the vertex signals are represented by linear time-invariant transfer functions (modules). Considering an acyclic network where only a subset of its vertices are measured and a subset of the vertices are excited, we explore conditions under which all the modules are identifiable on the basis of measurement data. Two suf… ▽ More This paper deals with dynamic networks in which the causality relations between the vertex signals are represented by linear time-invariant transfer functions (modules). Considering an acyclic network where only a subset of its vertices are measured and a subset of the vertices are excited, we explore conditions under which all the modules are identifiable on the basis of measurement data. Two sufficient conditions are presented in the paper, where the first condition concerns an identifiability analysis that needs to be performed for each vertex, while the second condition, based on the concept of tree/anti-tree covering, results from a graphical synthesis approach to allocate actuators and sensors in an acyclic network for achieving generic identifiability. △ Less

Submitted 5 April, 2024; v1 submitted 19 January, 2022; originally announced January 2022.

First images of the molecular gas around a born-again star revealed by ALMA

Authors: Daniel Tafoya, Jesús A. Toalá, Ramlal Unnikrishnan, Wouter H. T. Vlemmings, Martín A. Guerrero, Stefan Kimeswenger, Peter A. M. van Hoof, Luis A. Zapata, Sandra P. Treviño-Morales, Janis B. Rodríguez-González

Abstract: Born-again stars allow probing stellar evolution in human timescales and provide the most promising path for the formation of hydrogen-deficient post-asymptotic giant branch objects, but their cold and molecular components remain poorly explored. Here we present ALMA observations of V605 Aql that unveil for the first time the spatio-kinematic distribution of the molecular material associated to a… ▽ More Born-again stars allow probing stellar evolution in human timescales and provide the most promising path for the formation of hydrogen-deficient post-asymptotic giant branch objects, but their cold and molecular components remain poorly explored. Here we present ALMA observations of V605 Aql that unveil for the first time the spatio-kinematic distribution of the molecular material associated to a born-again star. Both the continuum and molecular line emission exhibit a clumpy ring-like structure with a total extent of $\approx$1$^{\prime\prime}$ in diameter. The bulk of the molecular emission is interpreted as being produced in a radially-expanding disk-like structure with an expansion velocity v$_{\rm exp}$$\sim$90 km s$^{-1}$ and an inclination $i$$\approx$60$^{\circ}$ with respect to the line-of-sight. The observations also reveal a compact high-velocity component, v$_{\rm exp}$$\sim$280 km s$^{-1}$, that is aligned perpendicularly to the expanding disk. This component is interpreted as a bipolar outflow with a kinematical age $τたう$$\lesssim$20 yr, which could either be material that is currently being ejected from V605 Aql, or it is being dragged from the inner parts of the disk by a stellar wind. The dust mass of the disk is in the range $M_{\rm dust}$$\sim$0.2-8$\times$10$^{-3}$ M$_{\odot}$, depending on the dust absorption coefficient. The mass of the CO is $M_{\rm CO}$$\approx$1.1$\times10^{-5}$ $M_{\odot}$, which is more than three orders of magnitude larger than the mass of the other detected molecules. We estimate a $^{12}$C/$^{13}$C ratio of 5.6$\pm$0.6, which is consistent with the single stellar evolution scenario in which the star experienced a very late thermal pulse instead of a nova-like event as previously suggested. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 9 pages, 5 Figures; accepted to ApJL

On data-driven control: informativity of noisy input-output data with cross-covariance bounds

Authors: Tom R. V. Steentjes, Mircea Lazar, Paul M. J. Van den Hof

Abstract: In this paper we develop new data informativity based controller synthesis methods that extend existing frameworks in two relevant directions: a more general noise characterization in terms of cross-covariance bounds and informativity conditions for control based on input-output data. Previous works have derived necessary and sufficient informativity conditions for noisy input-state data with quad… ▽ More In this paper we develop new data informativity based controller synthesis methods that extend existing frameworks in two relevant directions: a more general noise characterization in terms of cross-covariance bounds and informativity conditions for control based on input-output data. Previous works have derived necessary and sufficient informativity conditions for noisy input-state data with quadratic noise bounds via an S-procedure. Although these bounds do not capture cross-covariance bounds in general, we show that the S-procedure is still applicable for obtaining non-conservative conditions on the data. Informativity-conditions for stability, $\mathcal{H}_\infty$ and $\mathcal{H}_2$ control are developed, which are sufficient for input-output data and also necessary for input-state data. Simulation experiments illustrate that cross-covariance bounds can be less conservative for informativity, compared to norm bounds typically employed in the literature. △ Less

Submitted 23 March, 2022; v1 submitted 28 November, 2021; originally announced November 2021.

Comments: In the previous version of this arXiv article and the corresponding L-CSS article, DOI: 10.1109/LCSYS.2021.3139526, there was a typographical error in equation (21) and (25). The equations contain $P_-$, which should be $Y_-$. This typo is corrected in this version

Tensor Product Scheme for Computing Bound States of the Quantum Mechanical Three-Body Problem

Authors: Jonas Thies, Moritz Travis Hof, Matthias Zimmermann, Maxim Efremov

Abstract: We develop a computationally and numerically efficient method to calculate binding energies and corresponding wave functions of quantum mechanical three-body problems in low dimensions. Our approach exploits the tensor structure of the multidimensional stationary Schrödinger equation, being expressed as a discretized linear eigenvalue problem. In one spatial dimension, we solve the three-body prob… ▽ More We develop a computationally and numerically efficient method to calculate binding energies and corresponding wave functions of quantum mechanical three-body problems in low dimensions. Our approach exploits the tensor structure of the multidimensional stationary Schrödinger equation, being expressed as a discretized linear eigenvalue problem. In one spatial dimension, we solve the three-body problem with the help of iterative methods. Here the application of the Hamiltonian operator is represented by dense matrix-matrix products. In combination with a newly-designed preconditioner for the Jacobi-Davidson QR, our highly accurate tensor method offers a significantly faster computation of three-body energies and bound states than other existing approaches. For the two-dimensional case, we additionally make use of a hybrid distributed/shared memory parallel implementation to calculate the corresponding three-body energies. Our novel method is of high relevance for the analysis of few-body systems and their universal behavior, which is only governed by the particle masses, overall symmetries, and the spatial dimensionality. Our results have straightforward applications for ultracold atomic gases that are widespread and nowadays utilized in quantum sensors. △ Less

Submitted 2 March, 2022; v1 submitted 3 November, 2021; originally announced November 2021.

H, He-like recombination spectra IV; H, He-like recombination spectra IV: clarification and refinement of methodology for $l$-changing collisions

Authors: N. R. Badnell, F. Guzmán, S. Brodie, R. J. R. Williams, P. A. M. van Hoof, M. Chatzikos, G. J. Ferland

Abstract: Precise spectral diagnostic modelling of H~{\sc i} and He~{\sc ii} recombination spectra can constrain theoretical models which describe many astrophysical environments. Simple analytic expressions are of interest for collisional $l$-changing rate coefficients that are used by large-scale population modelling codes. We review, clarify and improve-upon the modified Pengelly \& Seaton formulae of Gu… ▽ More Precise spectral diagnostic modelling of H~{\sc i} and He~{\sc ii} recombination spectra can constrain theoretical models which describe many astrophysical environments. Simple analytic expressions are of interest for collisional $l$-changing rate coefficients that are used by large-scale population modelling codes. We review, clarify and improve-upon the modified Pengelly \& Seaton formulae of Guzmán \etal We show that the recent poor results for it shown by Vrinceanu \etal are due to their misinterpretation of its usage. We also detail efficient numerical algorithms which should enable the full quantum mechanical expression for such rate coefficients to be used much more routinely by modelling codes. We illustrate with some collisional-radiative population modelling for hydrogen. △ Less

Submitted 20 August, 2021; originally announced August 2021.

Comments: accepted in MNRAS

A scalable multi-step least squares method for network identification with unknown disturbance topology

Authors: Stefanie J. M. Fonken, Karthik R. Ramaswamy, Paul M. J. Van den Hof

Abstract: Identification methods for dynamic networks typically require prior knowledge of the network and disturbance topology, and often rely on solving poorly scalable non-convex optimization problems. While methods for estimating network topology are available in the literature, less attention has been paid to estimating the disturbance topology, i.e., the (spatial) noise correlation structure and the n… ▽ More Identification methods for dynamic networks typically require prior knowledge of the network and disturbance topology, and often rely on solving poorly scalable non-convex optimization problems. While methods for estimating network topology are available in the literature, less attention has been paid to estimating the disturbance topology, i.e., the (spatial) noise correlation structure and the noise rank in a filtered white noise representation of the disturbance signal. In this work we present an identification method for dynamic networks, in which an estimation of the disturbance topology precedes the identification of the full dynamic network with known network topology. To this end we extend the multi-step Sequential Linear Regression and Weighted Null Space Fitting methods to deal with reduced rank noise, and use these methods to estimate the disturbance topology and the network dynamics in the full measurement situation. As a result, we provide a multi-step least squares algorithm with parallel computation capabilities and that rely only on explicit analytical solutions, thereby avoiding the usual non-convex optimizations involved. Consequently we consistently estimate dynamic networks of Box Jenkins model structure, while keeping the computational burden low. We provide a consistency proof that includes path-based data informativity conditions for allocation of excitation signals in the experimental design. Numerical simulations performed on a dynamic network with reduced rank noise clearly illustrate the potential of this method. △ Less

Submitted 25 May, 2022; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: 17 pages, 4 figures, accepted and published in Automatica Volume 141, July 2022

Journal ref: Volume 141, July 2022, 110295

Planet Hunters TESS III: two transiting planets around the bright G dwarf HD 152843

Authors: Nora L. Eisner, Belinda A. Nicholson, Oscar Barragán, Suzanne Aigrain, Chris Lintott, Laurel Kaye, Baptiste Klein, Grant Miller, Jake Taylor, Norbert Zicher, Lars A. Buchhave, Douglas A. Caldwell, Jonti Horner, Joe Llama, Annelies Mortier, Vinesh M. Rajpaul, Keivan Stassun, Avi Sporer, Andrew Tkachenko, Jon M. Jenkins, David W. Latham, George R. Ricker, Sara Seager, Joshua N. Winn, Safaa Alhassan , et al. (15 additional authors not shown)

Abstract: We report on the discovery and validation of a two-planet system around a bright (V = 8.85 mag) early G dwarf (1.43 $R_{\odot}$, 1.15 $M_{\odot}$, TOI 2319) using data from NASA's Transiting Exoplanet Survey Satellite (TESS). Three transit events from two planets were detected by citizen scientists in the month-long TESS light curve (sector 25), as part of the Planet Hunters TESS project. Modellin… ▽ More We report on the discovery and validation of a two-planet system around a bright (V = 8.85 mag) early G dwarf (1.43 $R_{\odot}$, 1.15 $M_{\odot}$, TOI 2319) using data from NASA's Transiting Exoplanet Survey Satellite (TESS). Three transit events from two planets were detected by citizen scientists in the month-long TESS light curve (sector 25), as part of the Planet Hunters TESS project. Modelling of the transits yields an orbital period of \Pb\ and radius of $3.41 _{ - 0.12 } ^ { + 0.14 }$ $R_{\oplus}$ for the inner planet, and a period in the range 19.26-35 days and a radius of $5.83 _{ - 0.14 } ^ { + 0.14 }$ $R_{\oplus}$ for the outer planet, which was only seen to transit once. Each signal was independently statistically validated, taking into consideration the TESS light curve as well as the ground-based spectroscopic follow-up observations. Radial velocities from HARPS-N and EXPRES yield a tentative detection of planet b, whose mass we estimate to be $11.56 _{ - 6.14 } ^ { + 6.58 }$ $M_{\oplus}$, and allow us to place an upper limit of $27.5$ $M_{\oplus}$ (99 per cent confidence) on the mass of planet c. Due to the brightness of the host star and the strong likelihood of an extended H/He atmosphere on both planets, this system offers excellent prospects for atmospheric characterisation and comparative planetology. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: Accepted for publication in MNRAS (15 pages, 10 figures, 3 tables)

Identification of diffusively coupled linear networks through structured polynomial models

Authors: E. M. M., Kivits, Paul M. J. Van den Hof

Abstract: Physical dynamic networks most commonly consist of interconnections of physical components that can be described by diffusive couplings. These diffusive couplings imply that the cause-effect relationships in the interconnections are symmetric and therefore physical dynamic networks can be represented by undirected graphs. This paper shows how prediction error identification methods developed for l… ▽ More Physical dynamic networks most commonly consist of interconnections of physical components that can be described by diffusive couplings. These diffusive couplings imply that the cause-effect relationships in the interconnections are symmetric and therefore physical dynamic networks can be represented by undirected graphs. This paper shows how prediction error identification methods developed for linear time-invariant systems in polynomial form can be configured to consistently identify the parameters and the interconnection structure of diffusively coupled networks. Further, a multi-step least squares convex optimization algorithm is developed to solve the nonconvex optimization problem that results from the identification method. △ Less

Submitted 1 July, 2022; v1 submitted 3 June, 2021; originally announced June 2021.

A frequency domain approach for local module identification in dynamic networks

Authors: Karthik R. Ramaswamy, Péter Zoltán Csurcsia, Johan Schoukens, Paul M. J. Van den Hof

Abstract: In classical approaches of dynamic network identification, in order to identify a system (module) embedded in a dynamic network, one has to formulate a Multi-input-Single-output (MISO) identification problem that requires identification of a parametric model for all the modules constituting the MISO setup including (possibly) the noise model, and determine their model order. This requirement leads… ▽ More In classical approaches of dynamic network identification, in order to identify a system (module) embedded in a dynamic network, one has to formulate a Multi-input-Single-output (MISO) identification problem that requires identification of a parametric model for all the modules constituting the MISO setup including (possibly) the noise model, and determine their model order. This requirement leads to model order selection steps for modules that are of no interest to the experimenter which increases the computational complexity for large-sized networks. Also, identification using a parametric noise model (like BJ method) can suffer from local minima, however neglecting the noise model has its impact on the variance of the estimates. In this paper, we provide a two-step identification approach to avoid these problems. The first step involves performing a non-parametric indirect approach for a MISO identification problem to get the non-parametric frequency response function (FRF) estimates and its variance as a function of frequency. In the second step, the estimated non-parametric FRF of the target module is smoothed using a parametric frequency domain estimator with the estimated variance from the previous step as the non-parametric noise model. The developed approach is practical with weak assumptions on noise, uses the available toolbox, requires a parametric model only for the target module of interest, and uses a non-parametric noise model to reduce the variance of the estimates. Numerical simulations illustrate the potentials of the introduced method in comparison with the classical identification methods. △ Less

Submitted 23 May, 2021; originally announced May 2021.

Comments: 8 pages, 5 figures, Submitted for publication in Automatica, 22 May 2021

A Necessary Condition for Network Identifiability with Partial Excitation and Measurement

Authors: Xiaodong Cheng, Shengling Shi, Ioannis Lestas, Paul M. J. Van den Hof

Abstract: This paper considers dynamic networks where vertices and edges represent manifest signals and causal dependencies among the signals, respectively. We address the problem of how to determine if the dynamics of a network can be identified when only partial vertices are measured and excited. A necessary condition for network identifiability is presented, where the analysis is performed based on ident… ▽ More This paper considers dynamic networks where vertices and edges represent manifest signals and causal dependencies among the signals, respectively. We address the problem of how to determine if the dynamics of a network can be identified when only partial vertices are measured and excited. A necessary condition for network identifiability is presented, where the analysis is performed based on identifying the dependency of a set of rational functions from excited vertices to measured ones. This condition is further characterised by using an edge-removal procedure on the associated bipartite graph. Moreover, on the basis of necessity analysis, we provide a necessary and sufficient condition for identifiability in circular networks. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Guaranteed $\mathcal{H}_\infty$ performance analysis and controller synthesis for interconnected linear systems from noisy input-state data

Authors: Tom R. V. Steentjes, Mircea Lazar, Paul M. J. Van den Hof

Abstract: The increase in available data and complexity of dynamical systems has sparked the research on data-based system performance analysis and controller design. Recent approaches can guarantee performance and robust controller synthesis based on noisy input-state data of a single dynamical system. In this paper, we extend a recent data-based approach for guaranteed performance analysis to distributed… ▽ More The increase in available data and complexity of dynamical systems has sparked the research on data-based system performance analysis and controller design. Recent approaches can guarantee performance and robust controller synthesis based on noisy input-state data of a single dynamical system. In this paper, we extend a recent data-based approach for guaranteed performance analysis to distributed analysis of interconnected linear systems. We present a new set of sufficient LMI conditions based on noisy input-state data that guarantees $\mathcal{H}_\infty$ performance and have a structure that lends itself well to distributed controller synthesis from data. Sufficient LMI conditions based on noisy data are provided for the existence of a dynamic distributed controller that achieves $\mathcal{H}_\infty$ performance. The presented approach enables scalable analysis and control of large-scale interconnected systems from noisy input-state data sets. △ Less

Submitted 7 April, 2021; v1 submitted 26 March, 2021; originally announced March 2021.

Bayesian spatio-temporal models for stream networks

Authors: Edgar Santos-Fernandez, Jay M. Ver Hoef, Erin E. Peterson, James McGree, Daniel Isaak, Kerrie Mengersen

Abstract: Spatio-temporal models are widely used in many research areas including ecology. The recent proliferation of the use of in-situ sensors in streams and rivers supports space-time water quality modelling and monitoring in near real-time. A new family of spatio-temporal models is introduced. These models incorporate spatial dependence using stream distance while temporal autocorrelation is captured u… ▽ More Spatio-temporal models are widely used in many research areas including ecology. The recent proliferation of the use of in-situ sensors in streams and rivers supports space-time water quality modelling and monitoring in near real-time. A new family of spatio-temporal models is introduced. These models incorporate spatial dependence using stream distance while temporal autocorrelation is captured using vector autoregression approaches. Several variations of these novel models are proposed using a Bayesian framework. The results show that our proposed models perform well using spatio-temporal data collected from real stream networks, particularly in terms of out-of-sample RMSPE. This is illustrated considering a case study of water temperature data in the northwestern United States. △ Less

Submitted 14 February, 2022; v1 submitted 5 March, 2021; originally announced March 2021.

Comments: 30 pages, 10 figs

Optical Near-Field Electron Microscopy

Authors: Raphaël Marchand, Radek Šachl, Martin Kalbáč, Martin Hof, Rudolf Tromp, Mariana Amaro, Sense J. van der Molen, Thomas Juffmann

Abstract: Imaging dynamical processes at interfaces and on the nanoscale is of great importance throughout science and technology. While light-optical imaging techniques often cannot provide the necessary spatial resolution, electron-optical techniques damage the specimen and cause dose-induced artefacts. Here, Optical Near-field Electron Microscopy (ONEM) is proposed, an imaging technique that combines non… ▽ More Imaging dynamical processes at interfaces and on the nanoscale is of great importance throughout science and technology. While light-optical imaging techniques often cannot provide the necessary spatial resolution, electron-optical techniques damage the specimen and cause dose-induced artefacts. Here, Optical Near-field Electron Microscopy (ONEM) is proposed, an imaging technique that combines non-invasive probing with light, with a high spatial resolution read-out via electron optics. Close to the specimen, the optical near-fields are converted into a spatially varying electron flux using a planar photocathode. The electron flux is imaged using low energy electron microscopy, enabling label-free nanometric resolution without the need to scan a probe across the sample. The specimen is never exposed to damaging electrons. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Journal ref: Phys. Rev. Applied 16, 014008 (2021)

Single module identifiability in linear dynamic networks with partial excitation and measurement

Authors: Shengling Shi, Xiaodong Cheng, Paul M. J. Van den Hof

Abstract: Identifiability of a single module in a network of transfer functions is determined by whether a particular transfer function in the network can be uniquely distinguished within a network model set, on the basis of data. Whereas previous research has focused on the situations that all network signals are either excited or measured, we develop generalized analysis results for the situation of parti… ▽ More Identifiability of a single module in a network of transfer functions is determined by whether a particular transfer function in the network can be uniquely distinguished within a network model set, on the basis of data. Whereas previous research has focused on the situations that all network signals are either excited or measured, we develop generalized analysis results for the situation of partial measurement and partial excitation. As identifiability conditions typically require a sufficient number of external excitation signals, this work introduces a novel network model structure such that excitation from unmeasured noise signals is included, which leads to less conservative identifiability conditions than relying on measured excitation signals only. More importantly, graphical conditions are developed to verify global and generic identifiability of a single module based on the topology of the dynamic network. Depending on whether the input or the output of the module can be measured, we present four identifiability conditions which cover all possible situations in single module identification. These conditions further lead to synthesis approaches for allocating excitation signals and selecting measured signals, to warrant single module identifiability. In addition, if the identifiability conditions are satisfied, indirect identification methods are developed to provide a consistent estimate of the module. All the obtained results are also extended to identifiability of multiple modules in the network. △ Less

Submitted 20 December, 2021; v1 submitted 21 December, 2020; originally announced December 2020.

Improved Fe II emission line models for AGN using new atomic datasets

Authors: A. Sarkar, G. J. Ferland, M. Chatzikos, F. Guzmán, P. A. M. van Hoof, R. T. Smyth, C. A. Ramsbottom, F. P. Keenan, C. P. Ballance

Abstract: Understanding the Fe II emission from Active Galactic Nuclei (AGN) has been a grand challenge for many decades. The rewards from understanding the AGN spectra would be immense, involving both quasar classification schemes such as "Eigenvector 1" and tracing the chemical evolution of the cosmos. Recently, three large Fe II atomic datasets with radiative and electron collisional rates have become av… ▽ More Understanding the Fe II emission from Active Galactic Nuclei (AGN) has been a grand challenge for many decades. The rewards from understanding the AGN spectra would be immense, involving both quasar classification schemes such as "Eigenvector 1" and tracing the chemical evolution of the cosmos. Recently, three large Fe II atomic datasets with radiative and electron collisional rates have become available. We have incorporated these into the spectral synthesis code Cloudy and examine predictions using a new generation of AGN Spectral Energy Distribution (SED), which indicates that the UV emission can be quite different depending on the dataset utilized. The Smyth et al dataset better reproduces the observed Fe II template of the I ZW 1 Seyfert galaxy in the UV and optical regions, and we adopt these data. We consider both thermal and microturbulent clouds and show that a microturbulence of $\approx$ 100 km/s reproduces the observed shape and strength of the so-called Fe II "UV bump". Comparing our predictions with the observed Fe II template, we derive a typical cloud density of $10^{11}$ cm$^{-3}$ and photon flux of $10^{20}$ cm$^{-2}$ s$^{-1}$, and show that these largely reproduce the observed Fe II emission in the UV and optical. We calculate the $I$(Fe II)/$I$(Mg II) emission-line intensity ratio using our best-fitting model and obtain log($I$(Fe II)/$I$(Mg II)) $\sim$ 0.7, suggesting many AGNs have a roughly solar Fe/Mg abundance ratio. Finally, we vary the Eddington ratio and SED shape as a step in understanding the Eigenvector 1 correlation. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Data-driven distributed control: Virtual reference feedback tuning in dynamic networks

Authors: Tom R. V. Steentjes, Mircea Lazar, Paul M. J. Van den Hof

Abstract: In this paper, the problem of synthesizing a distributed controller from data is considered, with the objective to optimize a model-reference control criterion. We establish an explicit ideal distributed controller that solves the model-reference control problem for a structured reference model. On the basis of input-output data collected from the interconnected system, a virtual experiment setup… ▽ More In this paper, the problem of synthesizing a distributed controller from data is considered, with the objective to optimize a model-reference control criterion. We establish an explicit ideal distributed controller that solves the model-reference control problem for a structured reference model. On the basis of input-output data collected from the interconnected system, a virtual experiment setup is constructed which leads to a network identification problem. We formulate a prediction-error identification criterion that has the same global optimum as the model-reference criterion, when the controller class contains the ideal distributed controller. The developed distributed controller synthesis method is illustrated on an academic example network of nine subsystems and the influence of the controller interconnection structure on the achieved closed-loop performance is analyzed. △ Less

Submitted 27 October, 2020; originally announced October 2020.

A Bayesian method for inference of effective connectivity in brain networks for detecting the Mozart effect

Authors: Rik J. C. van Esch, Shengling Shi, Antoine Bernas, Svitlana Zinger, Albert P. Aldenkamp, Paul M. J. Van den Hof

Abstract: Several studies claim that listening to Mozart music affects cognition and can be used to treat neurological conditions like epilepsy. Research into this Mozart effect has not addressed how dynamic interactions between brain networks, i.e. effective connectivity, are affected. The Granger-causality analysis is often used to infer effective connectivity. First, we investigate if a new method, Bayes… ▽ More Several studies claim that listening to Mozart music affects cognition and can be used to treat neurological conditions like epilepsy. Research into this Mozart effect has not addressed how dynamic interactions between brain networks, i.e. effective connectivity, are affected. The Granger-causality analysis is often used to infer effective connectivity. First, we investigate if a new method, Bayesian topology identification, can be used as an alternative. Both methods are evaluated on simulation data, where the Bayesian method outperforms the Granger-causality analysis in the inference of connectivity graphs of dynamic networks, especially for short data lengths. In the second part, the Bayesian method is extended to enable the inference of changes in effective connectivity between groups of subjects. Next, we apply both methods to fMRI scans of 16 healthy subjects, who were scanned before and after exposure to Mozart's sonata K448 at least 2 hours a day for 7 days. Here, we investigate if the effective connectivity of the subjects significantly changed after listening to Mozart music. The Bayesian method detected changes in effective connectivity between networks related to cognitive processing and control: First, in the connection from the central executive to the superior sensori-motor network. Second, in the connection from the posterior default mode to the fronto-parietal right network. Finally, in the connection from the anterior default mode to the dorsal attention network, but only in a subgroup of subjects with a longer listening duration. Only in this last connection an effect was found by the Granger-causality analysis. △ Less

Submitted 15 October, 2020; originally announced October 2020.

The cooling-down central star of the planetary nebula SwSt\,1: a late thermal pulse in a massive post-AGB star?

Authors: Marcin Hajduk, Helge Todt, Wolf-Rainer Hamann, Karolina Borek, Peter A. M. van Hoof, Albert A. Zijlstra

Abstract: SwSt 1 (PN G001.5-06.7) is a bright and compact planetary nebula containing a late [WC]-type central star. Previous studies suggested that the nebular and stellar lines are slowly changing with time. We studied new and archival optical and ultraviolet spectra of the object. The [OIII] 4959 and 5007 A to $\mathrm{H}βべーた$ line flux ratios decreased between about 1976 and 1997/2015. The stellar spectrum… ▽ More SwSt 1 (PN G001.5-06.7) is a bright and compact planetary nebula containing a late [WC]-type central star. Previous studies suggested that the nebular and stellar lines are slowly changing with time. We studied new and archival optical and ultraviolet spectra of the object. The [OIII] 4959 and 5007 A to $\mathrm{H}βべーた$ line flux ratios decreased between about 1976 and 1997/2015. The stellar spectrum also shows changes between these epochs. We modeled the stellar and nebular spectra observed at different epochs. The analyses indicate a drop of the stellar temperature from about 42 kK to 40.5 kK between 1976 and 1993. We do not detect significant changes between 1993 and 2015. The observations show that the star performed a loop in the H-R diagram. This is possible when a shell source is activated during its post-AGB evolution. We infer that a late thermal pulse (LTP) experienced by a massive post-AGB star can explain the evolution of the central star. Such a star does not expand significantly as the result of the LTP and does not became a born-again red giant. However, the released energy can remove the tiny H envelope of the star. △ Less

Submitted 30 September, 2020; originally announced September 2020.

Comments: 18 pages, 15 figures, published in MNRAS

Journal ref: 2020, Volume 498, pp. 1205-1220

Generic identifiability of subnetworks in a linear dynamic network: the full measurement case

Authors: Shengling Shi, Xiaodong Cheng, Paul M. J. Van den Hof

Abstract: Identifiability conditions for single or multiple modules in a dynamic network specify under which conditions the considered modules can be uniquely recovered from the second-order statistical properties of the measured signals. Conditions for generic identifiability of multiple modules, i.e. a subnetwork, are developed for the situation that all node signals are measured and excitation of the net… ▽ More Identifiability conditions for single or multiple modules in a dynamic network specify under which conditions the considered modules can be uniquely recovered from the second-order statistical properties of the measured signals. Conditions for generic identifiability of multiple modules, i.e. a subnetwork, are developed for the situation that all node signals are measured and excitation of the network is provided by both measured excitation signals and unmeasured disturbance inputs. Additionally, the network model set is allowed to contain non-parametrized modules that are fixed, and e.g. reflect modules of which the dynamics are known to the user. The conditions take the form of path-based conditions on the graph of the network model set. Based on these conditions, synthesis results are formulated for allocating external excitation signals to achieve generic identifiability of particular subnetworks. If there are a sufficient number of measured external excitation signals, the formulated results give rise to a generalized indirect type of identification algorithm that requires only the measurement of a subset of the node signals in the network. △ Less

Submitted 26 October, 2021; v1 submitted 4 August, 2020; originally announced August 2020.

Learning linear modules in a dynamic network using regularized kernel-based methods

Authors: Karthik R. Ramaswamy, Giulio Bottegal, Paul M. J. Van den Hof

Abstract: In order to identify one system (module) in an interconnected dynamic network, one typically has to solve a Multi-Input-Single-Output (MISO) identification problem that requires identification of all modules in the MISO setup. For application of a parametric identification method this would require estimating a large number of parameters, as well as an appropriate model order selection step for a… ▽ More In order to identify one system (module) in an interconnected dynamic network, one typically has to solve a Multi-Input-Single-Output (MISO) identification problem that requires identification of all modules in the MISO setup. For application of a parametric identification method this would require estimating a large number of parameters, as well as an appropriate model order selection step for a possibly large scale MISO problem, thereby increasing the computational complexity of the identification algorithm to levels that are beyond feasibility. An alternative identification approach is presented employing regularized kernel-based methods. Keeping a parametric model for the module of interest, we model the impulse response of the remaining modules in the MISO structure as zero mean Gaussian processes (GP) with a covariance matrix (kernel) given by the first-order stable spline kernel, accounting for the noise model affecting the output of the target module and also for possible instability of systems in the MISO setup. Using an Empirical Bayes (EB) approach the target module parameters are estimated through an Expectation-Maximization (EM) algorithm with a substantially reduced computational complexity, while avoiding extensive model structure selection. Numerical simulations illustrate the potentials of the introduced method in comparison with the state-of-the-art techniques for local module identification. △ Less

Submitted 26 January, 2021; v1 submitted 13 May, 2020; originally announced May 2020.

Comments: 15 pages, 7 figures, Submitted for publication in Automatica, 12 May 2020. Final version of paper submitted on 06 January 2021 (To appear in Automatica)

A Linear Mixed Model Formulation for Spatio-Temporal Random Processes with Computational Advances for the Separable and Product-Sum Covariances

Authors: Michael Dumelle, Jay M. Ver Hoef, Claudio Fuentes, Alix Gitelman

Abstract: We describe spatio-temporal random processes using linear mixed models. We show how many commonly used models can be viewed as special cases of this general framework and pay close attention to models with separable or product-sum covariances. The proposed linear mixed model formulation facilitates the implementation of a novel algorithm using Stegle eigendecompositions, a recursive application of… ▽ More We describe spatio-temporal random processes using linear mixed models. We show how many commonly used models can be viewed as special cases of this general framework and pay close attention to models with separable or product-sum covariances. The proposed linear mixed model formulation facilitates the implementation of a novel algorithm using Stegle eigendecompositions, a recursive application of the Sherman-Morrison-Woodbury formula, and Helmert-Wolf blocking to efficiently invert separable and product-sum covariance matrices, even when every spatial location is not observed at every time point. We show our algorithm provides noticeable improvements over the standard Cholesky decomposition approach. Via simulations, we assess the performance of the separable and product-sum covariances and identify scenarios where separable covariances are noticeably inferior to product-sum covariances. We also compare likelihood-based and semivariogram-based estimation and discuss benefits and drawbacks of both. We use the proposed approach to analyze daily maximum temperature data in Oregon, USA, during the 2019 summer. We end by offering guidelines for choosing among these covariances and estimation methods based on properties of observed data. △ Less

Submitted 2 May, 2020; originally announced May 2020.

Comments: 43 pages (including an Appendix) and 8 figures

Journal ref: Spatial Staistics, Volume 43, 2021

Current and future development of the photoionization code Cloudy

Authors: P. A. M. van Hoof, G. C. Van de Steene, F. Guzmán, M. Dehghanian, M. Chatzikos, G. J. Ferland

Abstract: The gas that is present in the interstellar medium is usually very far removed from (local) thermodynamic equilibrium, and in some cases may also not be in a steady-state equilibrium with its surroundings. The physics of this material is complex and one needs a sophisticated numerical code to study it. For this purpose the open-source photoionization code Cloudy was created. It models the physical… ▽ More The gas that is present in the interstellar medium is usually very far removed from (local) thermodynamic equilibrium, and in some cases may also not be in a steady-state equilibrium with its surroundings. The physics of this material is complex and one needs a sophisticated numerical code to study it. For this purpose the open-source photoionization code Cloudy was created. It models the physical state of the gas and predicts the spectrum that it emits. Cloudy is continually being developed to improve the treatment of the microphysical processes and the database of fundamental data that it uses. In this paper we will discuss how we are developing the code to improve our high-density predictions by implementing better collisional-radiative models for all ions. We will also briefly discuss the experimental mode in Cloudy to model gas that is not in steady-state equilibrium and present a preliminary model of recombining gas in a planetary nebula that is on the cooling track. We finish with a short discussion of how we are speeding up the code by using parallelization. △ Less

Submitted 13 February, 2020; originally announced February 2020.

Comments: 12 pages, 5 figures. This paper was presented as an invited lecture at the 12th Serbian Conference on Spectral Line Shapes in Astrophysics in Vrdnik, Serbia

Journal ref: Contributions of the Astronomical Observatory Skalnaté Pleso, 2020, vol. 50, pp 32-43

Scalable distributed and decentralized $\mathscr{H}_2$ controller synthesis for interconnected linear discrete-time systems

Authors: Tom R. V. Steentjes, Mircea Lazar, Paul M. J. Van den Hof

Abstract: The current limitation in the synthesis of distributed $\mathscr{H}_2$ controllers for linear interconnected systems is scalability due to non-convex or unstructured synthesis conditions. In this paper we develop convex and structured conditions for the existence of a distributed $\mathscr{H}_2$ controller for discrete-time interconnected systems with an interconnection structure that corresponds… ▽ More The current limitation in the synthesis of distributed $\mathscr{H}_2$ controllers for linear interconnected systems is scalability due to non-convex or unstructured synthesis conditions. In this paper we develop convex and structured conditions for the existence of a distributed $\mathscr{H}_2$ controller for discrete-time interconnected systems with an interconnection structure that corresponds to an arbitrary graph. Neutral interconnections and a storage function with a block-diagonal structure are utilized to attain coupling conditions that are of a considerably lower computational complexity compared to the corresponding centralized $\mathscr{H}_2$ controller synthesis problem. Additionally, the developed conditions are adapted for the corresponding decentralized $\mathscr{H}_2$ controller synthesis problem with fixed supply functions for the interconnections. The effectiveness and scalability of the developed distributed $\mathscr{H}_2$ controller synthesis method is demonstrated for small- to large-scale oscillator networks on a cycle graph. △ Less

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

Comments: Changes in this version include: overview of dissipativity-based results for interconnected systems is removed, result on the existence of a decentralized controller is included and a simulation example for the illustration of scalability replaces the previous simulation example

SSNdesign -- an R package for pseudo-Bayesian optimal and adaptive sampling designs on stream networks

Authors: Alan R. Pearse, James M. McGree, Nicholas A. Som, Catherine Leigh, Jay M. Ver Hoef, Paul Maxwell, Erin E. Peterson

Abstract: Streams and rivers are biodiverse and provide valuable ecosystem services. Maintaining these ecosystems is an important task, so organisations often monitor the status and trends in stream condition and biodiversity using field sampling and, more recently, autonomous in-situ sensors. However, data collection is often costly and so effective and efficient survey designs are crucial to maximise info… ▽ More Streams and rivers are biodiverse and provide valuable ecosystem services. Maintaining these ecosystems is an important task, so organisations often monitor the status and trends in stream condition and biodiversity using field sampling and, more recently, autonomous in-situ sensors. However, data collection is often costly and so effective and efficient survey designs are crucial to maximise information while minimising costs. Geostatistics and optimal and adaptive design theory can be used to optimise the placement of sampling sites in freshwater studies and aquatic monitoring programs. Geostatistical modelling and experimental design on stream networks pose statistical challenges due to the branching structure of the network, flow connectivity and directionality, and differences in flow volume. Thus, unique challenges of geostatistics and experimental design on stream networks necessitates the development of new open-source software for implementing the theory. We present SSNdesign, an R package for solving optimal and adaptive design problems on stream networks that integrates with existing open-source software. We demonstrate the mathematical foundations of our approach, and illustrate the functionality of SSNdesign using two case studies involving real data from Queensland, Australia. In both case studies we demonstrate that the optimal or adaptive designs outperform random and spatially balanced survey designs. The SSNdesign package has the potential to boost the efficiency of freshwater monitoring efforts and provide much-needed information for freshwater conservation and management. △ Less

Submitted 1 December, 2019; originally announced December 2019.

Comments: Main document: 18 pages, 7 figures Supp Info A: 11 pages, 0 figures Supp Info B: 24 pages, 6 figures Supp Info C: 3 pages, 0 figures

Allocation of Excitation Signals for Generic Identifiability of Linear Dynamic Networks

Authors: Xiaodong Cheng, Shengling Shi, Paul M. J. Van den Hof

Abstract: A recent research direction in data-driven modeling is the identification of dynamic networks, in which measured vertex signals are interconnected by dynamic edges represented by causal linear transfer functions. The major question addressed in this paper is where to allocate external excitation signals such that a network model set becomes generically identifiable when measuring all vertex signal… ▽ More A recent research direction in data-driven modeling is the identification of dynamic networks, in which measured vertex signals are interconnected by dynamic edges represented by causal linear transfer functions. The major question addressed in this paper is where to allocate external excitation signals such that a network model set becomes generically identifiable when measuring all vertex signals. To tackle this synthesis problem, a novel graph structure, referred to as \textit{directed pseudotree}, is introduced, and the generic identifiability of a network model set can be featured by a set of disjoint directed pseudotrees that cover all the parameterized edges of an \textit{extended graph}, which includes the correlation structure of the process noises. Thereby, an algorithmic procedure is devised, aiming to decompose the extended graph into a minimal number of disjoint pseudotrees, whose roots then provide the appropriate locations for excitation signals. Furthermore, the proposed approach can be adapted using the notion of \textit{anti-pseudotrees} to solve a dual problem, that is to select a minimal number of measurement signals for generic identifiability of the overall network, under the assumption that all the vertices are excited. △ Less

Submitted 15 February, 2021; v1 submitted 10 October, 2019; originally announced October 2019.

A wind-based unification model for NGC 5548: spectral holidays, non-disk emission, and implications for changing-look quasars

Authors: M. Dehghanian, G. J. Ferland, B. M. Peterson, G. A. Kriss, K. T. Korista, M. Chatzikos, F. Guzman, N. Arav, G. De Rosa, M. R. Goad, M. Mehdipour, P. A. M. Van Hoof

Abstract: The 180-day Space Telescope and Optical Reverberation Mapping campaign on NGC 5548 discovered an anomalous period, the broad-line region (BLR) holiday, in which the emission lines decorrelated from the continuum variations. This is important since the correlation between the continuum-flux variations and the emission-line response is the basic assumption for black hole (BH) mass determinations thr… ▽ More The 180-day Space Telescope and Optical Reverberation Mapping campaign on NGC 5548 discovered an anomalous period, the broad-line region (BLR) holiday, in which the emission lines decorrelated from the continuum variations. This is important since the correlation between the continuum-flux variations and the emission-line response is the basic assumption for black hole (BH) mass determinations through reverberation mapping. During the BLR holiday, the high-ionization intrinsic absorption lines also decorrelated from the continuum as a result of variable covering factor of the line of sight (LOS) obscurer. The emission lines are not confined to the LOS, so this does not explain the BLR holiday. If the LOS obscurer is a disk wind, its streamlines must extend down to the plane of the disk and the base of the wind would lie between the BH and the BLR, forming an equatorial obscurer. This obscurer can be transparent to ionizing radiation, or can be translucent, blocking only parts of the SED, depending on its density. An emission-line holiday is produced if the wind density increases only slightly above its transparent state. Both obscurers are parts of the same wind, so they can have associated behavior in a way that explains both holidays. A very dense wind would block nearly all ionizing radiation, producing a Seyfert 2 and possibly providing a contributor to the changing-look AGN phenomenon. Disk winds are very common and we propose that the equatorial obscurers are too, but mostly in a transparent state. △ Less

Submitted 22 August, 2019; v1 submitted 20 August, 2019; originally announced August 2019.

Comments: Accepted for publication in ApJL

A local direct method for module identification in dynamic networks with correlated noise

Authors: Karthik R. Ramaswamy, Paul M. J. Van den Hof

Abstract: The identification of local modules in dynamic networks with known topology has recently been addressed by formulating conditions for arriving at consistent estimates of the module dynamics, under the assumption of having disturbances that are uncorrelated over the different nodes. The conditions typically reflect the selection of a set of node signals that are taken as predictor inputs in a MISO… ▽ More The identification of local modules in dynamic networks with known topology has recently been addressed by formulating conditions for arriving at consistent estimates of the module dynamics, under the assumption of having disturbances that are uncorrelated over the different nodes. The conditions typically reflect the selection of a set of node signals that are taken as predictor inputs in a MISO identification setup. In this paper an extension is made to arrive at an identification setup for the situation that process noises on the different node signals can be correlated with each other. In this situation the local module may need to be embedded in a MIMO identification setup for arriving at a consistent estimate with maximum likelihood properties. This requires the proper treatment of confounding variables. The result is a set of algorithms that, based on the given network topology and disturbance correlation structure, selects an appropriate set of node signals as predictor inputs and outputs in a MISO or MIMO identification setup. Three algorithms are presented that differ in their approach of selecting measured node signals. Either a maximum or a minimum number of measured node signals can be considered, as well as a preselected set of measured nodes. △ Less

Submitted 31 October, 2020; v1 submitted 2 August, 2019; originally announced August 2019.

Comments: Submitted for publication in IEEE Transactions on Automatic Control, 2 August 2019. arXiv admin note: text overlap with arXiv:1809.07502. Final version of paper scheduled to appear in IEEE Transactions on Automatic Control, Vol. 66, No. 11, November 2021

MSC Class: 93A99

Bayesian topology identification of linear dynamic networks

Authors: Shengling Shi, Giulio Bottegal, Paul M. J. Van den Hof

Abstract: In networks of dynamic systems, one challenge is to identify the interconnection structure on the basis of measured signals. Inspired by a Bayesian approach in [1], in this paper, we explore a Bayesian model selection method for identifying the connectivity of networks of transfer functions, without the need to estimate the dynamics. The algorithm employs a Bayesian measure and a forward-backward… ▽ More In networks of dynamic systems, one challenge is to identify the interconnection structure on the basis of measured signals. Inspired by a Bayesian approach in [1], in this paper, we explore a Bayesian model selection method for identifying the connectivity of networks of transfer functions, without the need to estimate the dynamics. The algorithm employs a Bayesian measure and a forward-backward search algorithm. To obtain the Bayesian measure, the impulse responses of network modules are modeled as Gaussian processes and the hyperparameters are estimated by marginal likelihood maximization using the expectation-maximization algorithm. Numerical results demonstrate the effectiveness of this method. △ Less

Submitted 14 March, 2019; originally announced March 2019.

H-, He-like recombination spectra III: $n$-changing collisions in highly-excited Rydberg states and their impact on the radio, IR and optical recombination lines

Authors: F. Guzmán, M. Chatzikos, P. A. M. van Hoof, Dana S. Balser, M. Dehghanian, N. R. Badnell, G. J. Ferland

Abstract: At intermediate to high densities, electron (de-)excitation collisions are the dominant process for populating or depopulating high Rydberg states. In particular, the accurate knowledge of the energy changing ($n$-changing) collisional rates is determinant for predicting the radio recombination spectra of gaseous nebula. The different datasets present in the literature come either from impact para… ▽ More At intermediate to high densities, electron (de-)excitation collisions are the dominant process for populating or depopulating high Rydberg states. In particular, the accurate knowledge of the energy changing ($n$-changing) collisional rates is determinant for predicting the radio recombination spectra of gaseous nebula. The different datasets present in the literature come either from impact parameter calculations or semi-empirical fits and the rate coefficients agree within a factor of two. We show in this paper that these uncertainties cause errors lower than 5% in the emission of radio recombination lines (RRL) of most ionized plasmas of typical nebulae. However, in special circumstances where the transitions between Rydberg levels are amplified by maser effects, the errors can increase up to 20%. We present simulations of the optical depth and H$nαあるふぁ$ line emission of Active Galactic Nuclei (AGN) Broad Line Regions (BLRs) and the Orion Nebula Blister to showcase our findings. △ Less

Submitted 18 March, 2019; v1 submitted 13 March, 2019; originally announced March 2019.

Comments: 16 pages (accepted by MNRAS, replaced)

The Belgian repository of fundamental atomic data and stellar spectra (BRASS) II. Quality assessment of atomic data for unblended lines in FGK stars

Authors: M. Laverick, A. Lobel, P. Royer, T. Merle, C. Martayan, P. A. M. van Hoof, M. Van der Swaelmen, M. David, H. Hensberge, E. Thienpont

Abstract: Fundamental atomic transition parameters, such as oscillator strengths and wavelengths, play a key role in modelling and understanding the chemical composition of stars in the universe. Despite the significant work under way to produce these parameters for many ions, uncertainties in these parameters remain large and can limit the accuracy of chemical abundance determinations. The Belgian reposito… ▽ More Fundamental atomic transition parameters, such as oscillator strengths and wavelengths, play a key role in modelling and understanding the chemical composition of stars in the universe. Despite the significant work under way to produce these parameters for many ions, uncertainties in these parameters remain large and can limit the accuracy of chemical abundance determinations. The Belgian repository of fundamental atomic data and stellar spectra (BRASS) aims to provide a large systematic and homogeneous quality assessment of the atomic data available for quantitative spectroscopy. BRASS shall compare synthetic spectra against extremely high quality observed spectra, at a resolution of ~85000 and signal-noise ratios of ~1000, for around 20 bright BAFGK spectral type stars, in order to evaluate the atomic data available for over a thousand potentially useful spectral lines. A large-scale homogeneous selection of atomic lines is performed by synthesising theoretical spectra of literature atomic lines, for FGK-type stars including the Sun, resulting in a selection of 1091 theoretically deep and unblended lines, in the wavelength range 4200-6800~Å, which may be suitable for quality assessment. Astrophysical log(gf) values are determined for the 1091 transitions using two commonly employed methods. The agreement of these log(gf) values are used to select well-behaving lines for quality assessment. 845 atomic lines were found to be suitable for quality assessment, of which 408 were found to be robust against systematic differences between analysis methods. Around 54% of the quality-assessed lines were found to have at least one literature log(gf) value in agreement with our derived values, though the remaining values can disagree by as much as 0.5 dex. Only 38% of FeI lines were found to have sufficiently accurate log(gf) values, increasing to ~70-75% for the remaining Fe-group lines. △ Less

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

Comments: 15 pages + 35 appendix pages, 20 figures, 6 tables. Full appendix tables available via CDS

Journal ref: A&A 624, A60 (2019)

Abstractions of linear dynamic networks for input selection in local module identification

Authors: Harm H. M. Weerts, Jonas Linder, Martin Enqvist, Paul M. J. Van den Hof

Abstract: In abstractions of linear dynamic networks, selected node signals are removed from the network, while keeping the remaining node signals invariant. The topology and link dynamics, or modules, of an abstracted network will generally be changed compared to the original network. Abstractions of dynamic networks can be used to select an appropriate set of node signals that are to be measured, on the b… ▽ More In abstractions of linear dynamic networks, selected node signals are removed from the network, while keeping the remaining node signals invariant. The topology and link dynamics, or modules, of an abstracted network will generally be changed compared to the original network. Abstractions of dynamic networks can be used to select an appropriate set of node signals that are to be measured, on the basis of which a particular local module can be estimated. A method is introduced for network abstraction that generalizes previously introduced algorithms, as e.g. immersion and the method of indirect inputs. For this abstraction method it is shown under which conditions on the selected signals a particular module will remain invariant. This leads to sets of conditions on selected measured node variables that allow identification of the target module. △ Less

Submitted 9 April, 2020; v1 submitted 2 January, 2019; originally announced January 2019.

Comments: 17 pages, 7 figures. Paper to appear in Automatica, Vol. 117, July 2020

Space Telescope and Optical Reverberation Mapping Project. X. Understanding the Absorption-Line Holiday in NGC 5548

Authors: M. Dehghanian, G. J. Ferland, G. A. Kriss, B. M. Peterson, S. Mathur, M. Mehdipour, F. Guzman, M. Chatzikos, P. A. M. Van Hoof, R. J. R. Williams, N. Arav, A. J. Barth, M. C. Bentz, S. Bisogni, W. N. Brandt, D. M. Crenshaw, E. Dalla Bonta, G. De Rosa, M. M. Fausnaugh, J. M. Gelbord, M. R. Goad, A. Gupta, Keith Horne, J. Kaastra, C. Knigge , et al. (5 additional authors not shown)

Abstract: The flux variations in the emission lines in active galactic nuclei (AGNs) are driven by variations in the ionizing continuum flux --which are usually reflected in the observable UV-optical continuum. The "Reverberation mapping" technique measures the delay between line and continuum variations to determine the size of the line emitting region, this is the basis for measurements of the central bla… ▽ More The flux variations in the emission lines in active galactic nuclei (AGNs) are driven by variations in the ionizing continuum flux --which are usually reflected in the observable UV-optical continuum. The "Reverberation mapping" technique measures the delay between line and continuum variations to determine the size of the line emitting region, this is the basis for measurements of the central black hole mass in AGNs. The Space Telescope and Optical Reverberation Mapping Project (AGN STORM) on NGC 5548 in 2014 is the most intensive multi-wavelength AGN monitoring campaign ever. For most of the campaign, the emission-line variations followed changes in the continuum with a time lag, as expected. However, the lines varied independently of the observed UV-optical continuum during a 60 -- 70 day "holiday." To understand this remarkable phenomenon, we study the intrinsic absorption lines present in NGC 5548. We identify a novel cycle that reproduces the absorption line variability and thus identify the physics that allows the holiday to occur. In our model, variations in this obscurer's line-of-sight covering factor modify the soft X-ray continuum. This leads to changes in the ionization of helium gas in the broad-line region. Ionizing radiation produced by recombining helium then affects the ionization of other species as observed during the AGN STORM holiday. It is likely that any other model which selectively changes the soft X-ray part of the continuum during the holiday can also explain the anomalous emission line behavior observed. △ Less

Submitted 24 April, 2019; v1 submitted 30 December, 2018; originally announced December 2018.

Comments: This is the version accepted by ApJ. Will be published shortly

Comparing Spatial Regression to Random Forests for Large Environmental Data Sets

Authors: Eric W. Fox, Jay M. Ver Hoef, Anthony R. Olsen

Abstract: Environmental data may be "large" due to number of records, number of covariates, or both. Random forests has a reputation for good predictive performance when using many covariates with nonlinear relationships, whereas spatial regression, when using reduced rank methods, has a reputation for good predictive performance when using many records that are spatially autocorrelated. In this study, we c… ▽ More Environmental data may be "large" due to number of records, number of covariates, or both. Random forests has a reputation for good predictive performance when using many covariates with nonlinear relationships, whereas spatial regression, when using reduced rank methods, has a reputation for good predictive performance when using many records that are spatially autocorrelated. In this study, we compare these two techniques using a data set containing the macroinvertebrate multimetric index (MMI) at 1859 stream sites with over 200 landscape covariates. A primary application is mapping MMI predictions and prediction errors at 1.1 million perennial stream reaches across the conterminous United States. For the spatial regression model, we develop a novel transformation procedure that estimates Box-Cox transformations to linearize covariate relationships and handles possibly zero-inflated covariates. We find that the spatial regression model with transformations, and a subsequent selection of significant covariates, has cross-validation performance slightly better than random forests. We also find that prediction interval coverage is close to nominal for each method, but that spatial regression prediction intervals tend to be narrower and have less variability than quantile regression forest prediction intervals. A simulation study is used to generalize results and clarify advantages of each modeling approach. △ Less

Submitted 26 December, 2018; originally announced December 2018.

Local module identification in dynamic networks with correlated noise: the full input case

Authors: Paul M. J. Van den Hof, Karthik R. Ramaswamy, Arne G. Dankers, Giulio Bottegal

Abstract: The identification of local modules in dynamic networks with known topology has recently been addressed by formulating conditions for arriving at consistent estimates of the module dynamics, typically under the assumption of having disturbances that are uncorrelated over the different nodes. The conditions typically reflect the selection of a set of node signals that are taken as predictor inputs… ▽ More The identification of local modules in dynamic networks with known topology has recently been addressed by formulating conditions for arriving at consistent estimates of the module dynamics, typically under the assumption of having disturbances that are uncorrelated over the different nodes. The conditions typically reflect the selection of a set of node signals that are taken as predictor inputs in a MISO identification setup. In this paper an extension is made to arrive at an identification setup for the situation that process noises on the different node signals can be correlated with each other. In this situation the local module may need to be embedded in a MIMO identification setup for arriving at a consistent estimate with maximum likelihood properties. This requires the proper treatment of confounding variables. The result is an algorithm that, based on the given network topology and disturbance correlation structure, selects an appropriate set of node signals as predictor inputs and outputs in a MISO or MIMO identification setup. As a first step in the analysis, we restrict attention to the (slightly conservative) situation where the selected output node signals are predicted based on all of their in-neighbor node signals in the network. △ Less

Submitted 25 March, 2019; v1 submitted 20 September, 2018; originally announced September 2018.

Comments: Extended version of paper submitted to the 58th IEEE Conf. Decision and Control, Nice, 2019

Suppression of Dielectronic Recombination Due to Finite Density Effects II: Analytical Refinement and Application to Density-dependent Ionization Balances and AGN Broad-line Emission

Authors: D. Nikolić, T. W. Gorczyca, K. T. Korista, M. Chatzikos, G. J. Ferland, F. Guzmán, P. A. M. van Hoof, R. J. R. Williams, N. R. Badnell

Abstract: We present improved fits to our treatment of suppression of dielectronic recombination at intermediate densities. At low densities, most recombined excited states eventually decay to the ground state, and therefore the total dielectronic recombination rate to all levels is preserved. At intermediate densities, on the other hand, collisions can lead to ionization of higher-lying excited states, the… ▽ More We present improved fits to our treatment of suppression of dielectronic recombination at intermediate densities. At low densities, most recombined excited states eventually decay to the ground state, and therefore the total dielectronic recombination rate to all levels is preserved. At intermediate densities, on the other hand, collisions can lead to ionization of higher-lying excited states, thereby suppressing the dielectronic recombination rate. The improved suppression factors presented here, although highly approximate, allow summed recombination rate coefficients to be used to intermediate densities. There have been several technical improvements to our previously presented fits. For H- through B-like ions the activation log densities have been adjusted to better reproduce existing data. For B-, C-, Al-, and Si-like ions secondary autoionization is now included. The treatment of density discontinuity in electron excitations out of ground state H-, He-, and Ne-like ions has been improved. These refined dielectronic recombination suppression factors are used in the most recent version of the plasma simulation code Cloudy. We show how the ionization and emission spectrum change when this physics is included. Although these suppression factors improve the treatment of intermediate densities, they are highly approximate and are not a substitution for a complete collisional-radiative model of the ionization balance. △ Less

Submitted 22 August, 2018; originally announced August 2018.