Nonlinear Kinetic Parameter Identification through Map Inversion

Characterization of a Driven Two-Level Quantum System by Supervised Learning

We investigate the extent to which a two-level quantum system subjected to an external time-dependent drive can be characterized by supervised learning. We apply this approach to the case of bang-bang control and the estimation of the offset and the final distance to a given target state. For any control protocol, the goal is to find the mapping between the offset and the distance. This mapping is interpolated using a neural network. The estimate is global in the sense that no a priori knowledge is required on the relation to be determined. Different neural network algorithms are tested on a series of data sets. We show that the mapping can be reproduced with very high precision in the direct case when the offset is known, while obstacles appear in the indirect case starting from the distance to the target. We point out the limits of the estimation procedure with respect to the properties of the mapping to be interpolated. We discuss the physical relevance of the different results.

A closed-loop identification protocol for nonlinear dynamical systems

A previous work introduced an optimal identification (OI) technique for reliably extracting model parameters of biochemical reaction systems from tailored laboratory experiments. The notion of optimality enters through seeking an external control in the laboratory producing data that leads to minimum uncertainties in the identified parameter distributions. A number of algorithmic and operational improvements are introduced in this paper to OI, aiming to build a more practical and efficient closed-loop identification protocol/procedure (CLIP) for nonlinear dynamical systems. The improvements in CLIP include (a) inversion cost function modification to preferably search for the upper and lower boundaries of the parameter distributions consistent with the observed data, (b) dynamic search range updating of the unknown parameters to better exploit the information from the prior iterative experiments, (c) replacing the control genetic algorithm by the simplex method to enable better balance between operational cost and inversion quality, and (d) utilizing virtual sensitivity optimization techniques to further reduce the laboratory costs. The workings of CLIP utilizing these new algorithms are illustrated in indentifying a simulated tRNA proofreading model, and the results demonstrate enhanced performance of CLIP in terms of algorithmic reliability and efficiency.

Extraction of Parameters and Their Error Distributions from Cyclic Voltammograms Using Bootstrap Resampling Enhanced by Solution Maps: Computational Study

The conventional determination of model parameter errors in least-squares regression of experimental cyclic voltammetric data assumes validity of local approximations (e.g., linearization) in the parameter space and normal distributions of the data and parameter errors. Such assumptions may not always be satisfied in practice. Bootstrap resampling techniques present a more universally applicable approach to error estimation, which until now has not been used in cyclic voltammetric studies, owing to the high costs of the required voltammogram simulations. We demonstrate that the burden of computing voltammograms can be significantly reduced by the use of high-dimensional model representation (HDMR) solution mapping techniques, thereby making it feasible to apply the bootstrap data analysis in cyclic voltammetry. We perform computational experiments with bootstrap resampling, enhanced by HDMR maps, for a typical cyclic voltammetric model (i.e., the Eqrev Cirr Eqrev reaction mechanism at a planar macroelectrode under semi-infinite, pure diffusion transport conditions). The experiments reveal that the bootstrap distributions of the estimated parameters provide a satisfactory quantification of the parameter errors and can also be used for detecting statistical correlations of the parameters.

Efficient chemical kinetic modeling through neural network maps

An approach to modeling nonlinear chemical kinetics using neural networks is introduced. It is found that neural networks based on a simple multivariate polynomial architecture are useful in approximating a wide variety of chemical kinetic systems. The accuracy and efficiency of these ridge polynomial networks (RPNs) are demonstrated by modeling the kinetics of H(2) bromination, formaldehyde oxidation, and H(2)+O(2) combustion. RPN kinetic modeling has a broad range of applications, including kinetic parameter inversion, simulation of reactor dynamics, and atmospheric modeling.

Estimation of Molecular Properties by High-Dimensional Model Representation

Additivity models have been widely employed to approximate unknown molecular properties based on previously measured or calculated data for similar molecules. This paper proposes an improved formulation of additivity, which is based on high-dimensional model representation (HDMR). HDMR is a general function-mapping technique that expresses the output of a multivariate system in terms of a hierarchy of cooperative effects among its input variables. HDMR rests on the general observation that, for many physical systems, only relatively low-order input variable cooperativity is significant. A molecule is expressed as a multivariate system by defining binary-valued input variables corresponding to the presence or absence of a chemical bond, with the molecular property as the output. Conventional additivity decomposes a molecular property into contributions from nonoverlapping subcomponents of fixed size. On the other hand, HDMR decomposes a molecular property into the exact contributions from the full hierarchy of its variable-sized subcomponents and contains additivity as a special case. The complete hierarchical structure of HDMR can in many cases lead to a much more accurate estimate than conventional additivity. Also, when full group additivity is not possible, HDMR gives an expression for a lower-order approximation for the missing group additivity value, greatly expanding the scope of HDMR compared to additivity. The component terms in an HDMR approximation have well-defined physical significance. Moreover, HDMR gives an exact expression for the truncation error in any given HDMR approximation, also with a well-defined physical significance. The HDMR model is tested for the enthalpy of formation of a broad range of organic molecules, and its advantages over additivity are illustrated.

Methods for parameter identification in oscillatory networks and application to cortical and thalamic 600 Hz activity

Directed information transfer in the human brain occurs presumably by oscillations. As of yet, most approaches for the analysis of these oscillations are based on time-frequency or coherence analysis. The present work concerns the modeling of cortical 600 Hz oscillations, localized within the Brodmann Areas 3b and 1 after stimulation of the nervus medianus, by means of coupled differential equations. This approach leads to the so-called parameter identification problem, where based on a given data set, a set of unknown parameters of a system of ordinary differential equations is determined by special optimization procedures. Some suitable algorithms for this task are presented in this paper. Finally an oscillatory network model is optimally fitted to the data taken from ten volunteers.