
We study a class of importance sampling methods for stochastic differential
equations (SDEs). A smallnoise analysis is performed, and the results suggest
that a simple symmetrization procedure can significantly improve the
performance of our importance sampling schemes when the noise is not too large.
We demonstrate that this is indeed the case for a number of linear and
nonlinear examples. Potential applications, e.g., data assimilation, are
discussed.

We compare two approaches to the predictive modeling of dynamical systems
from partial observations at discrete times. The first is continuous in time,
where one uses data to infer a model in the form of stochastic differential
equations, which are then discretized for numerical solution. The second is
discrete in time, where one directly infers a discretetime model in the form
of a nonlinear autoregression moving average model. The comparison is performed
in a special case where the observations are known to have been obtained from a
hypoelliptic stochastic differential equation. We show that the discretetime
approach has better predictive skills, especially when the data are relatively
sparse in time. We discuss open questions as well as the broader significance
of the results.

Implicit samplers are algorithms for producing independent, weighted samples
from multivariate probability distributions. These are often applied in
Bayesian data assimilation algorithms. We use Laplace asymptotic expansions to
analyze two implicit samplers in the small noise regime. Our analysis suggests
a symmetrization of the algo rithms that leads to improved (implicit) sampling
schemes at a rel atively small additional cost. Computational experiments
confirm the theory and show that symmetrization is effective for small noise
sampling problems.

This paper considers the problem of designing timedependent, realtime
control policies for controllable nonlinear diffusion processes, with the goal
of obtaining maximallyinformative observations about parameters of interest.
More precisely, we maximize the expected Fisher information for the parameter
obtained over the duration of the experiment, conditional on observations made
up to that time. We propose to accomplish this with a twostep strategy: when
the full state vector of the diffusion process is observable continuously, we
formulate this as an optimal control problem and apply numerical techniques
from stochastic optimal control to solve it. When observations are incomplete,
infrequent, or noisy, we propose using standard filtering techniques to first
estimate the state of the system, then apply the optimal control policy using
the posterior expectation of the state. We assess the effectiveness of these
methods in 3 situations: a paradigmatic bistable model from statistical
physics, a model of action potential generation in neurons, and a model of a
simple ecological system.

Biological information processing is often carried out by complex networks of
interconnected dynamical units. A basic question about such networks is that of
reliability: if the same signal is presented many times with the network in
different initial states, will the system entrain to the signal in a repeatable
way? Reliability is of particular interest in neuroscience, where large,
complex networks of excitatory and inhibitory cells are ubiquitous. These
networks are known to autonomously produce strongly chaotic dynamics  an
obvious threat to reliability. Here, we show that such chaos persists in the
presence of weak and strong stimuli, but that even in the presence of chaos,
intermittent periods of highly reliable spiking often coexist with unreliable
activity. We elucidate the local dynamical mechanisms involved in this
intermittent reliability, and investigate the relationship between this
phenomenon and certain timedependent attractors arising from the dynamics. A
conclusion is that chaotic dynamics do not have to be an obstacle to precise
spike responses, a fact with implications for signal coding in large networks.

Perturbation theory is an important tool in the analysis of oscillators and
their response to external stimuli. It is predicated on the assumption that the
perturbations in question are "sufficiently weak", an assumption that is not
always valid when perturbative methods are applied. In this paper, we identify
a number of concrete dynamical scenarios in which a standard perturbative
technique, based on the infinitesimal phase response curve (PRC), is shown to
give different predictions than the full model. Shearinduced chaos, i.e.,
chaotic behavior that results from the amplification of small perturbations by
underlying shear, is missed entirely by the PRC. We show also that the presence
of "sticky" phasespace structures tend to cause perturbative techniques to
overestimate the frequencies and regularity of the oscillations. The phenomena
we describe can all be observed in a simple 2D neuron model, which we choose
for illustration as the PRC is widely used in mathematical neuroscience.

We review some recent results surrounding a general mechanism for producing
chaotic behavior in periodicallykicked oscillators. The key geometric ideas
are illustrated via a simple linear shear model.

We report the results of a numerical study of nonequilibrium steady states
for a class of Hamiltonian models. In these models of coupled matterenergy
transport, particles exchange energy through collisions with pinneddown
rotating disks. In [Commun. Math. Phys. 262 (2006)], Eckmann and Young studied
1D chains and showed that certain simple formulas give excellent approximations
of energy and particle density profiles.
Keeping the basic mode of interaction in [EckmannYoung], we extend their
prediction scheme to a number of new settings: 2D systems on different
lattices, driven by a variety of boundary (heat bath) conditions including the
use of thermostats. Particleconserving models of the same type are shown to
behave similarly. The second half of this paper examines memory and finitesize
effects, which appear to impact only minimally the profiles of the models
tested in [EckmannYoung]. We demonstrate that these effects can be significant
or insignificant depending on the local geometry. Dynamical mechanisms are
proposed, and in the case of directional bias in particle trajectories due to
memory, correction schemes are derived and shown to give accurate predictions.

We show that Markov couplings can be used to improve the accuracy of Markov
chain Monte Carlo calculations in some situations where the steadystate
probability distribution is not explicitly known. The technique generalizes the
notion of control variates from classical Monte Carlo integration. We
illustrate it using two models of nonequilibrium transport.

We study the reliability of large networks of coupled neural oscillators in
response to fluctuating stimuli. Reliability means that a stimulus elicits
essentially identical responses upon repeated presentations. We view the
problem on two scales: neuronal reliability, which concerns the repeatability
of spike times of individual neurons embedded within a network, and
pooledresponse reliability, which addresses the repeatability of the total
synaptic output from the network. We find that individual embedded neurons can
be reliable or unreliable depending on network conditions, whereas pooled
responses of sufficiently large networks are mostly reliable. We study also the
effects of noise, and find that some types affect reliability more seriously
than others.

This paper concerns the reliability of a pair of coupled oscillators in
response to fluctuating inputs. Reliability means that an input elicits
essentially identical responses upon repeated presentations regardless of the
network's initial condition. Our main result is that both reliable and
unreliable behaviors occur in this network for broad ranges of coupling
strengths, even though individual oscillators are always reliable when
uncoupled. A new finding is that at low input amplitudes, the system is highly
susceptible to unreliable responses when the feedforward and feedback couplings
are roughly comparable. A geometric explanation based on shearinduced chaos at
the onset of phaselocking is proposed.

We study the reliability of phase oscillator networks in response to
fluctuating inputs. Reliability means that an input elicits essentially
identical responses upon repeated presentations, regardless of the network's
initial condition. In this paper, we extend previous results on twocell
networks to larger systems. The first issue that arises is chaos in the absence
of inputs, which we demonstrate and interpret in terms of reliability. We give
a mathematical analysis of networks that can be decomposed into modules
connected by an acyclic graph. For this class of networks, we show how to
localize the source of unreliability, and address questions concerning
downstream propagation of unreliability once it is produced.

Guided by a geometric understanding developed in earlier works of Wang and
Young, we carry out some numerical studies of shearinduced chaos. The settings
considered include periodic kicking of limit cycles, random kicks at Poisson
times, and continuoustime driving by white noise. The forcing of a
quasiperiodic model describing two coupled oscillators is also investigated.
In all cases, positive Lyapunov exponents are found in suitable parameter
ranges when the forcing is suitably directed.

We present the results of a detailed study of energy correlations at steady
state for a 1D model of coupled energy and matter transport. Our aim is to
discover  via theoretical arguments, conjectures, and numerical simulations
 how spatial covariances scale with system size, their relations to local
thermodynamic quantities, and the randomizing effects of heat baths. Among our
findings are that shortrange covariances respond quadratically to local
temperature gradients, and longrange covariances decay linearly with
macroscopic distance. These findings are consistent with exact results for the
simple exclusion and KMP models.

This letter concerns the reliability of coupled oscillator networks in
response to fluctuating inputs. Reliability means that (following a transient)
an input elicits identical responses upon repeated presentations, regardless of
the system's initial condition. Here, we analyze this property for two coupled
oscillators, demonstrating that oscillator networks exhibit both reliable and
unreliable dynamics for broad ranges of coupling strengths. We further argue
that unreliable dynamics are characterized by strange attractors with random
SRB measures, implying that though unreliable, the responses lie on
lowdimensional sets. Finally, we show that 1:1 phase locking in the zeroinput
system corresponds to high susceptibility for unreliable responses. A geometric
explanation is proposed.

The HodgkinHuxley model describes action potential generation in certain
types of neurons and is a standard model for conductancebased, excitable
cells. Following the early work of Winfree and Best, this paper explores the
response of a spontaneously spiking HodgkinHuxley neuron model to a periodic
pulsatile drive. The response as a function of drive period and amplitude is
systematically characterized. A wide range of qualitatively distinct responses
are found, including entrainment to the input pulse train and persistent chaos.
These observations are consistent with a theory of kicked oscillators developed
by Qiudong Wang and LaiSang Young. In addition to general features predicted
by WangYoung theory, it is found that most combinations of drive period and
amplitude lead to entrainment instead of chaos. This preference for entrainment
over chaos is explained by the structure of the HodgkinHuxley phase resetting
curve.

This paper presents a systematic numerical study of the effects of noise on
the invariant probability densities of dynamical systems with varying degrees
of hyperbolicity. It is found that the rate of convergence of invariant
densities in the smallnoise limit is frequently governed by power laws. In
addition, a simple heuristic is proposed and found to correctly predict the
power law exponent in exponentially mixing systems. In systems which are not
exponentially mixing, the heuristic provides only an upper bound on the power
law exponent. As this numerical study requires the computation of invariant
densities across more than 2 decades of noise amplitudes, it also provides an
opportunity to discuss and compare standard numerical methods for computing
invariant probability densities.

We give a new upper bound on the Selberg zeta function for a convex
cocompact Schottky group acting on $ {\mathbb H}^{n+1}$: in strips parallel to
the imaginary axis the zeta function is bounded by $ \exp (C s^\delta) $
where $ \delta $ is the dimension of the limit set of the group. This bound is
more precise than the optimal global bound $ \exp (C s^{n+1}) $, and it gives
new bounds on the number of resonances (scattering poles) of $ \Gamma
\backslash {\mathbb H}^{n+1} $. The proof of this result is based on the
application of holomorphic $ L^2$techniques to the study of the determinants
of the Ruelle transfer operators and on the quasiselfsimilarity of limit
sets. We also study this problem numerically and provide evidence that the
bound may be optimal. Our motivation comes from molecular dynamics and we
consider $ \Gamma \backslash {\mathbb H}^{n+1} $ as the simplest model of
quantum chaotic scattering. The proof of this result is based on the
application of holomorphic $L^2$techniques to the study of the determinants of
the Ruelle transfer operators and on the quasiselfsimilarity of limit sets.

This paper presents numerical evidence that for quantum systems with chaotic
classical dynamics, the number of scattering resonances near an energy $E$
scales like $\hbar^{\frac{D(K_E)+1}{2}}$ as $\hbar\to{0}$. Here, $K_E$ denotes
the subset of the classical energy surface $\{H=E\}$ which stays bounded for
all time under the flow generated by the Hamiltonian $H$ and $D(K_E)$ denotes
its fractal dimension. Since the number of bound states in a quantum system
with $n$ degrees of freedom scales like $\hbar^{n}$, this suggests that the
quantity $\frac{D(K_E)+1}{2}$ represents the effective number of degrees of
freedom in scattering problems.