skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Brownian bridges for stochastic chemical processes—An approximation method based on the asymptotic behavior of the backward Fokker–Planck equation
Title: Brownian bridges for stochastic chemical processes—An approximation method based on the asymptotic behavior of the backward Fokker–Planck equation
A Brownian bridge is a continuous random walk conditioned to end in a given region by adding an effective drift to guide paths toward the desired region of phase space. This idea has many applications in chemical science where one wants to control the endpoint of a stochastic process—e.g., polymer physics, chemical reaction pathways, heat/mass transfer, and Brownian dynamics simulations. Despite its broad applicability, the biggest limitation of the Brownian bridge technique is that it is often difficult to determine the effective drift as it comes from a solution of a Backward Fokker–Planck (BFP) equation that is infeasible to compute for complex or high-dimensional systems. This paper introduces a fast approximation method to generate a Brownian bridge process without solving the BFP equation explicitly. Specifically, this paper uses the asymptotic properties of the BFP equation to generate an approximate drift and determine ways to correct (i.e., re-weight) any errors incurred from this approximation. Because such a procedure avoids the solution of the BFP equation, we show that it drastically accelerates the generation of conditioned random walks. We also show that this approach offers reasonable improvement compared to other sampling approaches using simple bias potentials.  more » « less
Award ID(s):
2126230
PAR ID:
10366344
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
156
Issue:
18
ISSN:
0021-9606
Page Range / eLocation ID:
184108
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, IEEE Transactions on Automatic Control)
    We consider particles obeying Langevin dynamics while being at known positions and having known velocities at the two end-points of a given interval. Their motion in phase space can be modeled as an Ornstein–Uhlenbeck process conditioned at the two end-points—a generalization of the Brownian bridge. Using standard ideas from stochastic optimal control we construct a stochastic differential equation (SDE) that generates such a bridge that agrees with the statistics of the conditioned process, as a degenerate diffusion. Higher order linear diffusions are also considered. In general, a time-varying drift is sufficient to modify the prior SDE and meet the end-point conditions. When the drift is obtained by solving a suitable differential Lyapunov equation, the SDE models correctly the statistics of the bridge. These types of models are relevant in controlling and modeling distribution of particles and the interpolation of density functions. 
    more » « less
  2. ; (, The Journal of Chemical Physics)
    We present a new method to sample conditioned trajectories of a system evolving under Langevin dynamics based on Brownian bridges. The trajectories are conditioned to end at a certain point (or in a certain region) in space. The bridge equations can be recast exactly in the form of a non-linear stochastic integro-differential equation. This equation can be very well approximated when the trajectories are closely bundled together in space, i.e., at low temperature, or for transition paths. The approximate equation can be solved iteratively using a fixed point method. We discuss how to choose the initial trajectories and show some examples of the performance of this method on some simple problems. This method allows us to generate conditioned trajectories with a high accuracy. 
    more » « less
  3. ; ; (, Inverse problems)
    This paper is concerned with the mathematical analysis of an inverse random source problem for the time fractional diffusion equation, where the source is driven by a fractional Brownian motion. Given the random source, the direct problem is to study the stochastic time fractional diffusion equation. The inverse problem is to determine the statistical properties of the source from the expectation and variance of the final time data. For the direct problem, we show that it is well-posed and has a unique mild solution under a certain condition. For the inverse problem, the uniqueness is proved and the instability is characterized. The major ingredients of the analysis are based on the properties of the Mittag–Leffler function and the stochastic integrals associated with the fractional Brownian motion. 
    more » « less
  4. ; ; ; (, Probability Theory and Related Fields)
    We give an explicit description of a family of jointly invariant measures of the KPZ equation singled out by asymptotic slope conditions. These are couplings of Brownian motions with drift, and can be extended to a cadlag process indexed by all real drift parameters. We name this process the KPZ horizon (KPZH). As a corollary, we resolve a recent conjecture by showing the existence of a random, countably infinite dense set of drift values at which the Busemann process of the KPZ equation is discontinuous. This signals instability, and shows the failure of the one force–one solution principle and the existence of at least two extremal semi-infinite polymer measures in the exceptional directions. The low-temperature limit of the KPZH is the stationary horizon (SH), the unique jointly invariant measure of the KPZ fixed point under the same slope conditions. The high-temperature limit of the KPZH is a coupling of Brownian motions that differ by linear shifts, which is jointly invariant under the Edwards–Wilkinson fixed point. 
    more » « less
  5. ; ; (, AIChE Journal)
    Abstract Using linear operator techniques, we demonstrate an efficient method for investigating rare events in stochastic processes. Specifically, we examine contained trajectories, which are continuous random walks that only leave a specified region of phase space after a set period of time . We show that such trajectories can be efficiently generated through the use of a Brownian Bridge, derived via the solution to the Backward Fokker–Planck (BFP) equation. Using linear operator techniques, we place the BFP operator in self‐adjoint form and show that in the asymptotic limit , the set of paths contained in a specified region is equivalent to paths on a modified potential energy landscape that is related to the dominant eigenfunction of the self‐adjoint BFP operator. We demonstrate this idea on several example problems, one of which is the Graetz problem, where one is interested in the survival time of a particle diffusing in tube flow. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email