Boundary conditions, probability distributions

For the solution of real tasks, depending on the concrete setup of the problem, either the forward or the backward Kolmogorov equation may be used. If the one-dimensional probability density with known initial distribution deserves needs to be determined, then it is natural to use the forward Kolmogorov equation. Contrariwise, if it is necessary to calculate the distribution of the mean first passage time as a function of initial state xo, then one should use the backward Kolmogorov equation. Let us now focus at the time on Eq. (2.6) as much widely used than (2.7) and discuss boundary conditions and methods of solution of this equation. 

So, Eq. (3.14) with boundary conditions is the equation for eigenfunction X (x) of the nth order. For X0(x), Eq. (3.14) will be an equation for stationary probability distribution with zero eigenvalue y0 = 0, and for X (x) the equation will have the following form ... 

The Transition Probability. Suppose we have a Brownian particle located at an initial instant of time at the point xo, which corresponds to initial delta-shaped probability distribution. It is necessary to find the probability Qc,d(t,xo) = Q(t,xo) of transition of the Brownian particle from the point c 0 Q(t,xo) = W(x, t) dx + Jrf+ X W(x, t) dx. The considered transition probability Q(t,xo) is different from the well-known probability to pass an absorbing boundary. Here we suppose that c and d are arbitrary chosen points of an arbitrary potential profile (x), and boundary conditions at these points may be arbitrary W(c, t) > 0, W(d, t) > 0. 

Smoluchowski, who worked on the rate of coagulation of colloidal particles, was a pioneer in the development of the theory of diffusion-controlled reactions. His theory is based on the assumption that the probability of reaction is equal to 1 when A and B are at the distance of closest approach (Rc) ( absorbing boundary condition ), which corresponds to an infinite value of the intrinsic rate constant kR. The rate constant k for the dissociation of the encounter pair can thus be ignored. As a result of this boundary condition, the concentration of B is equal to zero on the surface of a sphere of radius Rc, and consequently, there is a concentration gradient of B. The rate constant for reaction k (t) can be obtained from the flux of B, in the concentration gradient, through the surface of contact with A. This flux depends on the radial distribution function of B, p(r, t), which is a solution of Fick s equation... 

An elementary treatment of a three-level system under pulsed excitation was given in Sec. II.C. Pollack treats the steady-state condition and the turn-off condition as well. These cases are quite interesting and deserve further discussion. Figure 51 shows the three-level system he analyzed. The quantities Ni(f), N2(f), and N3(f)are the electron-distribution functions (populations) and a9 b, and c represent the total transition probabilities per unit time. The boundary condition Nx + N2 + N = JV =constant is assumed. 

The solution of eqn. (44) for a coulomb potential with boundary conditions (45) and (46) for either initial conditions (48) or (49) has only been developed in recent years. Hong and Noolandi [72] showed that the solution of the Debye—Smoluchowski equation is related to the Mathieu equation. Many of the details of their analysis are discussed in the Appendix A, Sect. 4, and the Appendix eqn. (A.21) is the Green s function (fundamental solution), which is the probability that a reactant B is at r given that it was initially at r0. This equation is developed as the Laplace transform. To obtain the density of interest p(r, ), with either condition, the Green s function has to be averaged over the initial distribution, as in eqn. (A.12), and the Laplace transform inverted. Alternatively, the density p(r, ) can be found from the inverse Laplace transform of the linear combination of independent solutions (A.17) which satisfy the boundary and initial conditions. This is shown in Fig. 10. For a Boltzmann initial condition, Hong and Noolandi [72] found... 

Furthermore, the initial and outer boundary conditions are effectively identical [eqns. (3), (4) and (165)] as are also the partially reflecting boundary conditions [eqns. (46) and (165)]. This can be shown by substituting p by exp — p p in the boundary conditions (165). Consequently, the relationship between the survival probability of an ion-pair at a time t0 after they were formed at time t and separation r and the density distribution of an initial (time t0) homogeneous distribution of the majority ion species around the minority ionic species, p(r, f f0), is an identity. 

Exercise. The equation (6.5) (without boundaries) does not have in general an equilibrium distribution of the separable form JX( )7(y). Hence the conditional probability P( y) for given y is not the same as the marginal probability P( ). For the case U(y) = j(x>oy2 show explicitly that this is so. 

We emphasize that the question of stability of a CA under small random perturbations is in itself an important unsolved problem in the theory of fluctuations [92-94] and the difficulties in solving it are similar to those mentioned above. Thus it is unclear at first glance how an analogy between these two unsolved problems could be of any help. However, as already noted above, the new method for statistical analysis of fluctuational trajectories [60,62,95,112] based on the prehistory probability distribution allows direct experimental insight into the almost deterministic dynamics of fluctuations in the limit of small noise intensity. Using this techique, it turns out to be possible to verify experimentally the existence of a unique solution, to identify the boundary condition on a CA, and to find an accurate approximation of the optimal control function. 

To find the boundary conditions on the CA, we analyze the prehistory probability distribution Ph(q, t qf, f/) of the escape trajectories. The corresponding distribution is shown in the Fig. 16. It can be inferred by the inspection of how the ridge of the most probable escape path merges the CA that most of the escape trajectories pass close to the saddle cycle of the period 5 embedded into the CA. 

The stationary solution of Eq. (12) with the boundary condition Eq. (13) gives the asymptotic distribution of protein coordinates which coincides with the length distribution of peptides p(x) (because probability of a strand to be cut is independent of its coordinate). The problem in Eqs. (12), (13) admits the stationary solution ... 

Concerning the boundary conditions of this problem, we can have various situations (i) in the first situation, the probabilities are null but not the probability gradients at z = 0 zero. For example, for a negative speed Vk(0,x), the particle is not in the stochastic space of displacement. However, at z = 0, we have a maximum probability for the output of the particle from the stochastic displacement space. Indeed, the flux of the characteristic probability must be a maximum and, consequently, dPk(0, x)/dz = 0 (ii) we have a similar situation at z = (iii) in other situations we can have uniformly distributed probabilities at the input in the stochastic displacement space then we can write the following expression ... 

At the aquifer scale, the most important contribution of age tracers is probably reduction in the nonuniqueness of numerical models. Lack of uniqueness stems, among other things, from inadequate knowledge of the distribution of hydraulic properties within groundwater systems and from poor constraints on boundary conditions (Konikow and Bredehoeft, 1992 Maloszewski and Zuber, 1993). Commonly, groundwater flow... 

Equation [16] is known as the Liouville equation and is, in fact, a statement of the conservation of the phase space probability density. Indeed, it can be seen that the Liouville equation takes the form of a continuity equation for a flow field on the phase space satisfying the incompressibility condition dIdT F = 0. Thus, given an initial phase space distribution function /(F, 0) and some appropriate boundary conditions on the phase space satisfied by f, Eq. [16] can be used to determine /(F, t) at any time t later. 

Fig. 13. Illustration of the grand-canonical simulation technique for temperature UbT/e = 1.68 and p = pcoex- A cuboidal system geometry 13.8cr x 13.8cr x 27.6cr is used with periodic boundary conditions in all three directions. The solid line corresponds to the negative logarithm of the probability distribution, P p), in the grand canonical ensemble. The two minima correspond to the coexisting phases and the arrows on the p axis mark their densities. The height of the plateau yields an accurate estimate for the interfacial tension, yLV- The dashed line is a parabohc fit in the vicinity of the liquid phase employed to determine the compressibihty. Representative system configurations are sketched schematically. From [62]...

The notation DX(t) means a path integral (summation over all trajectories, X(t)), that starts at Xq and ends at X,. The action depends on the boundary coordinates as well as on the trajectory. At the limit of a small time step the distribution of the trajectories, X(t), is sharply peaked at the exact trajectory. In fact, in the numerical examples discussed below we only consider exact trajectories. The effects of large steps and the corresponding errors on the conditional probability will be the topics of future work. 

In contrast, a system in contact with a thermal bath (constant-temperature, constant-volume ensemble) can be in a state of all energies, from zero to arbitrary large energies however, the state probability is different. The distribution of the probabilities is obtained under the assumption that the system plus the bath constimte a closed system. The imposed temperature varies linearly from start-temp to end-temp. The main techniques used to keep the system at a given temperature are velocity rescaling. Nose, and Nos Hoover-based thermostats. In general, the Nose-Hoover-based thermostat is known to perform better than other temperature control schemes and produces accurate canonical distributions. The Nose-Hoover chain thermostat has been found to perform better than the single thermostat, since the former provides a more flexible and broader frequency domain for the thermostat [29]. The canonical ensemble is the appropriate choice when conformational searches of molecules are carried out in vacuum without periodic boundary conditions. 

We now restrict ourselves to periodic functions V x) of x, so that we may impose periodic boundary conditions. Thus directly proportional to the equilibrium distribution. The probability density P may now be written in terms of the so-called biorthogonal expansion given by Morse and Feschbach [61]... 

For the computer simulation, a two-dimensional matrix of 60 x 60 with a hexagonal packing of elements was randomly filled by 1 and 0 (corresponding to Dad and Had atoms formed at the sites where appropriate reactions took place) with a deuterium content equal to that detected experimentally. A numerical value of the sum of n-n elements (0-6) was used to calculate the probability of H2, HD, and D2 formation for each element. The content of H2, HD, and D2 was calculated in the central part of the matrix of 30 x 30 elements to avoid the edge effects during both simultaneous irreversible recombination of the elements into dimmers and rearrangements of the matrix (alternatively, periodic boundary conditions could be applied). The calculated ratio H2% HD% D2% for the matrix randomly filled by Had and Dad at D2 mol% = 50.17 was found to be 25.02 49.61 25.37 confirming a proper statistical distribution of atoms in the initial matrix. 

This implies that Vv(r, t) can be used to describe an exponentially decaying probability distribution. The time-independent partial wave solutions above a threshold obey scattering boundary conditions. They are thus proportional to the regular solution tpe(E, r) a origin and tend to infinity like Jost solutions f E,r) ... 

---