Solution for the density distribution

It is worth noting that a simple solution of eqn. (10) is p(r,t) = constant 

Only one other general solution exists. Two methods may be used to solve a partial differential equation such as the diffusion equation, or wave equation separation of variables or Laplace transformation (Carslaw and Jaeger [26] Crank [27]). The Laplace transformation route is often easier, especially if the inversion of the Laplace transform can be found in standard tables [28]. The Laplace transform of a function of time, (t), is defined as 

Note that the initial condition appears. Replacing the dependent variable 1 

The arbitrary constants are not functions of distance, but may depend on s. The initial condition has already been used, but the boundary conditions [eqns. (4) and (5)] have not been used yet. Their Laplace transforms are, respectively 

Standard tables of inverse Laplace transforms enable this to be inverted very easily giving 

Substituting eqn. (13) into eqn. (12) in the limit as r°o shows that a = 0. Considering the inner boundary condition (14) where r = i in 

---

Y. However, no general solution to Eq. (32) exists due in large part to the ditfieulty of obtaining explicit expressions for the density distribution of electrostatie eharge around the solute molecule, Y y, which are both general and rigorous. 

Smoluchowski theory [29, 30] and its modifications fonu the basis of most approaches used to interpret bimolecular rate constants obtained from chemical kinetics experiments in tenus of difhision effects [31]. The Smoluchowski model is based on Brownian motion theory underlying the phenomenological difhision equation in the absence of external forces. In the standard picture, one considers a dilute fluid solution of reactants A and B with [A] [B] and asks for the time evolution of [B] in the vicinity of A, i.e. of the density distribution p(r,t) = [B](rl)/[B] 2i ] r(t))l ] Q ([B] is assumed not to change appreciably during the reaction). The initial distribution and the outer and inner boundary conditions are chosen, respectively, as... 

Each spray-dried droplet forms a single particle whose size is determined by the droplet size, the dissolved solids of the feed solution, and the density of the resulting solid particle. For a given formulation and process, both the solid content and density of the powder remain constant within a batch and from batch to batch therefore, the distribution of the primary particle size is determined by the droplet size distribution. A narrowly distributed particle size can be achieved with a well-designed atomizer and controlled process parameters. 

Usually, experiments are performed with steady-state photolysis or radiolysis of the solution and the yield of scavenger products determined optically or by ESR methods. There is no direct interest in the actual time evolution of the density or recombination (survival) probability. Consequently, the creation of ion-pairs may be pictured as occurring at a constant rate, say 1 s 1, from time t0 = 0 to infinity. The steady-state ion-pair density distribution, which arises when dp/dt = 0, is the balance between continuous creation of ion pairs at a rate Is-1, recombination and scavenging. Removing the instantaneous creation of an ion-pair at time t = t0 (i.e. removing the 6(f — f0) in the source term), means that ion-pairs were continuously formed from time t = — 00 to t. At long times, f > — oo the density distribution is independent of t and, of course, t0. Let pss(r cs r0) = /i p(r, t cs t0, 0)d 0 be the steady-state ion-pair density distribution for ion pairs continuously formed at r0, and note d/dt J" f pd 0 = 0. The diffusion equation (169) becomes... 

The Laplace Transform relationship between the density distribution for creation, recombination and scavenging and the time-dependent (or s-dependent) density distribution for creation, recombination and escape is a logical consequence of the nature of eqn. (169). The solution of eqn. (169) was shown by Mozumder [364] to be... 

Solution The data on numbers of particles in each particle range given in Table El.3 can be converted to relative frequencies per unit of particle size as given in Table El. 4. The histogram for the relative frequency per unit of particle size for the data is plotted in Fig. El.2 the histogram yields a total area of bars equal to unity. Superimposed on the histogram is the density function for the normal distribution based on Eqs. (1.24) and (1.30). For this distribution, the values for do and ad are evaluated as 0.342 and 0.181, respectively. Also included in the figure is the density function for the log-normal distribution based on Eq. (1.32a). For this distribution, the values for In doi and od are evaluated as —1.209 and 0.531, respectively. 

In practice AFGS is related to both the surface potential and to the flat-band voltage Vf. On the other hand, Vf is related to the virtual gate, which takes into account the presence of the solution as it has been shown in some details in Eqs. (3) and (4). A derivation of Id that does not take into account one of the two hypotheses for the density-of-states distribution considered in Ki-shida el al. (1983) could be more difficult nevertheless, a dependence on gate voltage will always be present and in any case could always be determined by experimental procedure. Also in this case, the 7d control will be exercised by the virtual gate, which contains, in ultimate analysis, the information on the ion concentration in solution. 

The solution to Poisson s equation for the depletion layer is discussed further in Chapter 9. The hatched region in Fig. 4.15 represents the gap states which change their charge state in depletion and so contribute to p(x). When there is a continuous distribution of gap states, p(x) is a spatially varying quantity. For the simpler case of a shallow donor-like level, the space charge equals the donor density N-a and the solution for the dependence of capacitance on applied bias, is (see Section 9.1.1)... 

Consider two parallel planar dissimilar ion-penetrable membranes 1 and 2 at separation h immersed in a solution containing a symmetrical electrolyte of valence z and bulk concentration n. We take an x-axis as shown in Fig. 13.2 [7-9]. We denote by Ni and Zi, respectively, the density and valence of charged groups in membrane 1 and by N2 and Z2 the corresponding quantities of membrane 2. Without loss of generality we may assume that Zj > 0 and Z2 may be either positive or negative and that Eq. (13.1) holds. The Poisson-Boltzmann equations (13.2)-(13.4) for the potential distribution j/(x) are rewritten in terms of the scaled potential y = zeif/IkT as... 

With the above, a formal set of equations Is given, the elaboration of which requiring a solution for the problem that the recurrent relationships p p - p p, . .. diverge. Relatively simple densities, or distribution functions, are converted into more complex ones. A "closure" is needed to "stop this explosion". A number of such closures have been proposed, all involving an assumption of which the rigour has to be tested. Most of these write three-body interactions in terms of three two-body Interactions, weighted in some way. A well known example is Kirkwood s superposition closure, which reads ... 

Table 44 The magnetochiral birefringence, An [in 10 cm /(T g)], per unit magnetic field and density for various chiral molecules. The theoretical method is either the Hartree-Fock or the B3LYP density-functional method, and the calculations were performed for either the gas phase or a solution. For the last two molecules, not only the most stable isomer was considered but also mixtures of six stable structures with populations according to a Boltzmann distribution. [a] is the specific optical rotation [in deg/(cm dm g)]. All results are from ref. 101... 

Some attempts to employ much less rigorous empirical schemes in the estimation of the electronic distribution for the considered structure and its correlation with spectral parameters such as chemical shifts were suggested by Bangov. The charge densities of each candidate structure may be calculated by a fast empirical scheme based on either full or partial equalization of orbital or atomicelectronegativity and further correlated with the chemical shifts from the spectrum of the query structure. The structure providing the best correlation is considered to be the correct solution for the unknown compound. However, this approach, as subsequently discussed, also yields rather vague results. 

---