Refine
Year of publication
Language
- English (1103) (remove)
Keywords
- optimal control (27)
- stability (14)
- integer programming (11)
- Stochastic programming (9)
- finite elements (9)
- mixed integer programming (9)
- Hamiltonian matrix (8)
- finite element method (8)
- model reduction (8)
- state constraints (8)
Persistence of rogue waves in extended nonlinear Schrödinger equations: Integrable Sasa-Satsuma case
(2012)
We present the lowest order rogue wave solution of the Sasa-Satsuma equation (SSE) which is one of the integrable extensions of the nonlinear Schrödinger equation (NLSE). In contrast to the Peregrine solution of the NLSE, it is significantly more involved and contains polynomials of fourth order rather than second order in the corresponding expressions. The correct limiting case of Peregrine solution appears when the extension parameter of the SSE is reduced to zero.
When simulating isolated resonators, the application of transparent boundary conditions causes the approximated spectrum to be polluted with spurious solutions. Distinguishing these artificial solutions from solutions with a physical meaning is often difficult and requires a priori knowledge of the spectrum or the expected field distribution of resonant states. We present an implementation of the pole condition that distinguishes between incoming and outgoing waves by the location of the poles of their Laplace transform as transparent boundary condition. This implementation depends on one tuning parameter. We will use the sensitivity of the computed solutions to perturbations of this parameter as a means to identify spurious solutions. To obtain global statements, we will combine this technique with a convergence monitor for the boundary condition.
The pole condition approach for deriving transparent boundary conditions is extended to the time-dependent, two-dimensional case. Non-physical modes of the solution are identified by the position of poles of the solution's spatial Laplace transform in the complex plane. By requiring the Laplace transform to be analytic on some
problem dependent complex half-plane, these modes can be
suppressed. The resulting algorithm computes a finite number of coefficients of a series expansion of the Laplace transform, thereby providing an approximation to the exact boundary condition. The resulting error decays super-algebraically with the number of coefficients, so relatively few additional degrees of freedom are
sufficient to reduce the error to the level of the discretization error in the interior of the computational domain. The approach shows good results for the Schroedinger and the drift-diffusion equation
but, in contrast to the one-dimensional case, exhibits instabilities for the wave and Klein-Gordon equation. Numerical examples are shown that demonstrate the good performance in the former and the instabilities in the latter case.
The hypergraph assignment problem (HAP) is the generalization of assignments
from directed graphs to directed hypergraphs. It serves, in particular,
as a universal tool to model several train composition rules in vehicle rotation
planning for long distance passenger railways. We prove that even for problems
with a small hyperarc size and hypergraphs with a special partitioned structure
the HAP is NP-hard and APX-hard. Further, we present an extended integer
linear programming formulation which implies, e. g., all clique inequalities.
We consider a semilinear parabolic equation subject to a nonlinear dynamical boundary condition that is related to the so-calles Wentzell boundary condition. First, we prove the existence and uniqueness of global solutions as well as the existence of a global attractor. Then we derive a suitable Lojasiewicz-Simon-type inequality to show the convergence of global solutions to single steady states as time tends to infinity under the assumption that the nonlinear terms $f$, $g$ are real analytic. Moreover, we provide an estimate for the convergence rate.
We investigate a distributed optimal control problem for a phase field
model of Cahn-Hilliard type. The model describes two-species phase segregation
on an atomic lattice under the presence of diffusion; it has been introduced recently in
[4], on the basis of the theory developed in [15], and consists of a system of two
highly nonlinearly coupled PDEs. For this reason, standard arguments of optimal control theory do not apply
directly, although the control constraints and the cost functional are of standard type.
We show that the problem admits a solution, and we derive the first-order
necessary conditions of optimality.
A nonlocal quasilinear multi-phase system with nonconstant specific heat and heat conductivity
(2012)
In this paper, we prove the existence
and global boundedness from above for a solution to an
integrodifferential model for nonisothermal multi-phase
transitions under nonhomogeneous third type boundary conditions.
The system couples a quasilinear internal energy balance
ruling the evolution of the absolute temperature with a vectorial
integro-differential inclusion governing the vectorial
phase-parameter dynamics. The specific heat and the heat
conductivity k are allowed to depend both on the order parameter
$\chi$ and on the absolute temperature $\teta$ of the system, and
the convex component of the free energy may or may not be
singular. Uniqueness and continuous data dependence are
also proved under additional assumptions.
This paper is devoted to an optimal control problem of Maxwell's equations in the presence of pointwise state constraints. The control is given by a divergence--free three--dimensional vector function representing an applied current density. To cope with the divergence--free constraint on the control, we consider a vector potential ansatz. Due to the lack of regularity of the control--to--state mapping, existence of Lagrange multipliers cannot be guaranteed. We regularize the optimal control problem by penalizing the pointwise state constraints. Optimality conditions for the regularized problem can be derived straightforwardly. It also turns out that the solution of the regularized problem enjoys higher regularity which then allows us to establish its convergence towards the solution of the unregularized problem. The second part of the paper focuses on the numerical analysis of the regularized optimal control problem. Here the state and the control are discretized by N\'ed\'elec's curl--conforming edge elements. Employing the higher regularity property of the optimal control, we establish an a priori error estimate for the discretization error in the $\boldsymbol{H}(\bold{curl})$--norm. The paper ends by numerical results including a numerical verification of our theoretical results.
An optimal control problem arising in the context of 3D electromagnetic induction heating is investigated. The state equation is given by a quasilinear stationary heat equation coupled with a semilinear time-harmonic eddy current equation. The temperature-dependent electrical conductivity and the presence of pointwise inequality state-constraints represent the main challenge of the paper. In the first part of the paper, the existence and regularity of the state are addressed. The second part of the paper deals with the analysis of the corresponding linearized equation. Some sufficient conditions are presented which guarantee the solvability of the linearized system. The final part of the paper is concerned with the optimal control. The aim of the optimization is to find the optimal voltage such that a desired temperature can be achieved optimally. The corresponding first-order necessary optimality condition is presented.
This paper is concerned with a PDE-constrained optimization problem of induction heating, where the state equations consist of 3D time--dependent heat equations coupled with 3D time--harmonic eddy current equations. The control parameters are given by finite real numbers representing applied alternating voltages which enter the eddy current equations via impressed current. The optimization problem is to find optimal voltages so that, under certain constraints on the voltages and the temperature, a desired temperature can be optimally achieved. As there are finitely many control parameters but the state constraint has to be satisfied in an infinite number of points, the problem belongs to a class of semi--infinite programming problems. We present a rigorous analysis of the optimization problem and a numerical strategy based on our theoretical result.
We consider a shape implant design problem that arises in the context of facial surgery. We introduce a reformulation as an optimal control problem, where the control acts as a boundary force. The state is modelled as a minimizer of a polyconvex hyperelastic energy functional. We show existence of optimal solutions and derive - on a formal level - first order optimality conditions. Finally, preliminary numerical results are presented.
Markov state models of molecular kinetics (MSMs), in which the long-time statistical dynamics
of a molecule is approximated by a Markov chain on a discrete partition of configuration space, have
seen widespread use in recent years. This approach has many appealing characteristics compared
to straightforward molecular dynamics simulation and analysis, including the potential to mitigate
the sampling problem by extracting long-time kinetic information from short trajectories and the
ability to straightforwardly calculate expectation values and statistical uncertainties of various
stationary and dynamical molecular observables. In this article, we summarize the current state of
the art in generation and validation of MSMs and give some important new results. We describe
an upper bound for the approximation error made by modeling molecular dynamics with an MSM
and we show that this error can be made arbitrarily small with surprisingly little effort. In contrast
to previous practice, it becomes clear that the best MSM is not obtained by the most metastable
discretization, but the MSM can be much improved if non-metastable states are introduced near
the transition states. Moreover, we show that it is not necessary to resolve all slow processes
by the state space partitioning, but individual dynamical processes of interest can be resolved
separately. We also present an efficient estimator for reversible transition matrices and a robust
test to validate that an MSM reproduces the kinetics of the molecular dynamics data.
Dynamical averages based on functionals of dynamical trajectories, such as time-correlation func-
tions, play an important role in determining kinetic or transport properties of matter. At temperatures
of interest, the expectations of these quantities are often dominated by contributions from rare events,
making the precise calculation of these quantities by molecular dynamics simulation difficult. Here,
we present a reweighting method for combining simulations from multiple temperatures (or from
simulated or parallel tempering simulations) to compute an optimal estimate of the dynamical prop-
erties at the temperature of interest without the need to invoke an approximate kinetic model (such as
the Arrhenius law). Continuous and differentiable estimates of these expectations at any temperature
in the sampled range can also be computed, along with an assessment of the associated statistical
uncertainty. For rare events, aggregating data from multiple temperatures can produce an estimate
of the desired precision at greatly reduced computational cost compared with simulations conducted
at a single temperature. Here, we describe use of the method for the canonical (NVT) ensemble us-
ing four common models of dynamics (canonical distribution of Hamiltonian trajectories, Andersen
thermostatting, Langevin, and overdamped Langevin or Brownian dynamics), but it can be applied to
any thermodynamic ensemble provided the ratio of path probabilities at different temperatures can be
computed. To illustrate the method, we compute a time-correlation function for solvated terminally-
blocked alanine peptide across a range of temperatures using trajectories harvested using a modified
parallel tempering protocol.
Dynamical fingerprints of macromolecules obtained from experiments often seem to indicate two- or
three state kinetics while simulations typically reveal a more complex picture. Markov state models of
molecular conformational dynamics can be used to predict these dynamical fingerprints and to reconcile
experiment with simulation. This is illustrated on two model systems: a one-dimensional energy surface
and a four-state model of a protein folding equilibrium. We show that (i) there might be no process
which corresponds to our notion of folding, (ii) often the experiment will be insensitive to some of the
processes present in the system, (iii) with a suitable combination the observable and initial conditions in
a relaxation experiment one can selectively measure specific processes. Furthermore, our method can be
used to design experiments such that specific processes appear with large amplitudes. We demonstrate
that for a fluorescence quenching experiment of the MR121-G9-W peptide.
Protein-ligand interactions are essential for nearly all biological processes, and yet the bio-
physical mechanism that enables potential binding partners to associate before specific binding
occurs remains poorly understood. Fundamental questions include which factors influence the
formation of protein-ligand encounter complexes, and whether designated association path-
ways exist. In this article we introduce a computational approach to systematically analyze
the complete ensemble of association pathways and to thus investigate these questions. This
approach is employed here to study the binding of a phosphate ion to the Escherichia coli
Phosphate Binding Protein. Various mutants of the protein are considered and their effects
on binding free energy profiles, association rates and association pathway distributions are
quantified. The results reveal the existence of two anion attractors, i.e. regions that initially
attract negatively charged particles and allow them to be efficiently screened for phosphate
which is specifically bound subsequently. Point mutations that affect the charge on these
attractors modulate their attraction strength and speed up association to a factor of 10 of
the diffusion limit and thus change the association pathways of the phosphate ligand. It is
demonstrated that a phosphate that pre-binds to such an attractor neutralizes its attraction
effect to the environment, making the simultaneous association of a second phosphate ion
unlikely. Our study suggests ways how structural properties can be used to tune molecular
association kinetics so as to optimize the efficiency of binding, and highlights the importance
of kinetic properties.
While studies of protein-ligand association have mostly focused on the native complex and its
stability (binding affinity), relatively little attention has been paid on the association process
that precedes the formation of the complex. Here we review approaches to study the kinet-
ics of association and association mechanisms, i.e. the probability distribution of association
pathways. Selected methods are described that allow these properties to be calculated quan-
titatively from simulation models. We summarize some applications of these methods and
finally propose a model mechanism by which proteins may efficiently screen potential ligands
for those that can be natively bound.
In this paper, we present a Gaussian Markov random field (GMRF) model for the transition
matrices (TMs) of Markov chains (MCs) by assuming the existence of a neighborhood relationship
between states, and develop the maximum a posteriori (MAP) estimators under different obser-
vation conditions. Unlike earlier work on TM estimation, our method can make full use of the
similarity between different states to improve the estimated accuracy, and the estimator can be
performed very efficiently by solving a convex programming problem. In addition, we discuss the
parameter choice of the proposed model, and introduce a Monte Carlo cross validation (MCCV)
method. The numerical simulations of a diffusion process are employed to show the effectiveness
of the proposed models and algorithms.
Optimal Identification of Semi-Rigid Domains in Macromolecules from Molecular Dynamics Simulation
(2012)
Biological function relies on the fact that biomolecules can switch between different conformations and aggregation states.
Such transitions involve a rearrangement of parts of the biomolecules involved that act as dynamic domains. The reliable
identification of such domains is thus a key problem in biophysics. In this work we present a method to identify semi-rigid
domains based on dynamical data that can be obtained from molecular dynamics simulations or experiments. To this end
the average inter-atomic distance-deviations are computed. The resulting matrix is then clustered by a constrained
quadratic optimization problem. The reliability and performance of the method are demonstrated for two artificial peptides.
Furthermore we correlate the mechanical properties with biological malfunction in three variants of amyloidogenic
transthyretin protein, where the method reveals that a pathological mutation destabilizes the natural dimer structure of the
protein. Finally the method is used to identify functional domains of the GroEL-GroES chaperone, thus illustrating the
efficiency of the method for large biomolecular machines.
Discrete-state Markov (or master equation) models provide a useful simplified representation for
characterizing the long-time statistical evolution of biomolecules in a manner that allows direct
comparison with experiments as well as the elucidation of mechanistic pathways for an inherently
stochastic process. A vital part of meaningful comparison with experiment is the characterization of
the statistical uncertainty in the predicted experimental measurement, which may take the form of
an equilibrium measurement of some spectroscopic signal, the time-evolution of this signal following
a perturbation, or the observation of some statistic (such as the correlation function) of the equilib-
rium dynamics of a single molecule. Without meaningful error bars (which arise due to the finite
quantity of data used to construct the model), there is no way to determine whether the deviations
between model and experiment are statistically meaningful. Previous work has demonstrated that
a Bayesian method that enforces microscopic reversibility can be used to characterize the correlated
uncertainties in state-to-state transition probabilities (and functions thereof) for a model inferred from
molecular simulation data. Here, we extend this approach to include the uncertainty in observables
that are functions of molecular conformation (such as surrogate spectroscopic signals) characteriz-
ing each state, permitting the full statistical uncertainty in computed spectroscopic experiments to be
assessed. We test the approach in a simple model system to demonstrate that the computed uncer-
tainties provide a useful indictor of statistical variation, and then apply it to the computation of the
fluorescence autocorrelation function measured for a dye-labeled peptide previously studied by both
experiment and simulation.
Markov (state) models (MSMs) have attracted a lot of interest recently as they (1) can probe
long-term molecular kinetics based on short-time simulations, (2) offer a way to analyze great
amounts of simulation data with relatively little subjectivity of the analyst, (3) provide insight into
microscopic quantities such as the ensemble of transition pathways, and (4) allow simulation data
to be reconciled with measurement data in a rigorous and explicit way. Here we sketch our current
perspective of Markov models and explain in short their theoretical basis and assumptions. We
describe transition path theory which allows the entire ensemble of protein folding pathways to be
investigated and that combines naturally with Markov models. Experimental observations can be
naturally linked to Markov models with the dynamical fingerprint theory, by which experimentally
observable timescales can be equipped with an understanding of the structural rearrangement
processes that take place at these timescales. The concepts of this paper are illustrated by a
simple kinetic model of protein folding.