Research Article | Open Access
Stochastic Signatures of Phase Space Decomposition
We explore the consequences of metrically decomposing a finite phase space, modeled as a d-dimensional lattice, into disjoint subspaces (lattices). Ergodic flows of a test particle undergoing an unbiased random walk are characterized by implementing the theory of finite Markov processes. Insights drawn from number theory are used to design the sublattices, the roles of lattice symmetry and system dimensionality are separately considered, and new lattice invariance relations are derived to corroborate the numerical accuracy of the calculated results. We find that the reaction efficiency in a finite system is strongly dependent not only on whether the system is compartmentalized, but also on whether the overall reaction space of the microreactor is further partitioned into separable reactors. We find that the reaction efficiency in a finite system is strongly dependent not only on whether the system is compartmentalized, but also on whether the overall reaction space of the microreactor is further partitioned into separable reactors. The sensitivity of kinetic processes in nanoassemblies to the dimensionality of compartmentalized reaction spaces is quantified.
To provide a physical motivation for the present study, understanding the factors influencing self-assembly in nanophase materials is a major experimental and theoretical challenge [1, 2]. Modeling even the early stages of self-assembly already introduces fundamental questions, for example, whether the process is reducible to a sequence of elementary steps occurring under equilibrium conditions [3, 4]. The role of different system geometries in influencing the efficiency of reactions in a finite, compartmentalized system  is an ongoing problem, one aspect of which is explored quantitatively here.
Before discussing the technical aspects of our work, we illustrate the main idea with a simple example. Consider a dimensional lattice of sites. One possible representation of such a lattice would be a system (see Figure 1). But, another representation that comprises sites would be the disjoint pair of lattices, a lattice with sites, and an lattice with sites, linked at a common corner site. To monitor ergodic flows on the assigned phase space(s), we place a sink (trap) at one of the corner vertices of the lattice and at the common vertex of the two sublattices. The transit time of a test particle undergoing unbiased, random displacements until, eventually, it is (irreversibly) localized at the trap can be determined by formulating and then solving numerically the stochastic master equation  for the problem. By extracting the first moment of the underlying probability distribution function, specified by the set of site-specific mean walklengths , the overall mean number of steps before localization, a signature of the mean transit time, can be obtained.
One fully expects to get different values of for these two lattice representations, but the question is “how different?” Also, the problem was illustrated above using (three) lattices. How do the results change if one or more of the disjoint spaces are lattices, that is, not approximately “square?” Is the difference in reaction efficiency magnified or suppressed? And, further, do these differences persist if one considers phase spaces of higher dimension? It is to the above questions that the present contribution is directed.
Fundamental results in number theory can be used to motivate the choices of lattices for which the metric decomposability of phase space can be studied. For dimension , one can take advantage of Fermat’s last theorem, namely, that one can find three positive integers , , and that satisfy the equation for for a given choice of . In , the decomposition is not unique; for any given , there may be (and usually is) more than one choice of that satisfies the theorem.
In three dimensions, however, the possibility of finding multiple lattices with such a degeneracy is much more limited, in fact, just the set of Hardy-Ramanujan numbers. The first is ; the second “taxi cab” number is These lattices are, practically speaking, as far as one can go. The next smallest Hardy-Ramanujan number is the last factorization one of several possible. Given the computational demands in implementing the theory of finite Markov processes, these lattices are just too large; hence, our results for dimension presented in Section 3 are (much) more limited than for .
The following section specifies the lattice statistical parameters needed to develop the Markovian theory. Section 3 presents the results of a series of calculations for and dimensional lattices. In Section 4, site-specific walklength data are used to deduce two new invariance relations, valid for finite lattices, which must be satisfied exactly to confirm that the computational results are numerically exact. In the final section, we comment on the relevance of our results to self-assembly in nanosystems.
We consider finite planar lattices, each site of which is of valence (connectivity) , and finite dimensional lattices, each site of which is of valence . At any interior site of the lattice, the test particle can move with equal a priori probability in any one of directions. If the particle initiates its motion at an edge site on a lattice, it can move in one of three directions on the lattice, but if it attempts to step “out of” the lattice, it is reset at the starting position. If the particle happens to be at a corner site, it can move to two adjacent lattice sites, remaining on the lattice, or is “reset” (stalls) at the corner site twice, rather than leaving the lattice. A similar protocol is implemented on lattices, where a particle on an interior site can move in directions, on a face site in directions with one reset, on a edge site in directions with two resets, and on a corner site in directions with three resets.
In all calculations here, the sink is placed at one corner site. Previous studies  have documented that moving the trap away from a centrosymmetric location will increase the walklength. However, the parent lattice or the lattices that result when one deconstructs a given initial lattice may not have a unique centrosymmetric site; for definiteness, we anchor the trap at a corner site. All calculations were carried out for nearest-neighbor displacements only and in the absence of any governing potentials; both restrictions can be lifted in further work.
a1729 = (52 + 112 + 712)/ 3. |
Displayed in Table 1 are representative results for the dimensional case. In every instance, one expects and finds that the results for are systematically larger for the “parent” lattice than for the corresponding two disjoint lattices. The more “symmetric” and “similar” are the pair of disjoint lattices, the more efficient is the underlying ergodic process. Computing the percent difference between results for the composite lattice versus two disjoint (but linked) lattices, one finds that the difference is always greater than 10%, and frequently much greater (in the range from 10%–30%, depending on lattice symmetry).
In Table 2, we find that the difference in results for dimensions is much less severe. For , the percentage is less than 4%. Hence, dimensionality appears to play the critical role in cases studied here. Note that there is a companion, dimensional result in Table 2 which was obtained by taking advantage of the decomposition of into the average of the greatest member in each pair of Brown numbers , and . Once again, the difference between the (here) three disjoint lattices and the parent lattices is ~15%.
4. Invariance Relations
Over 40 years ago, Montroll and Weiss  reported for lattices subject to periodic boundary conditions an exact, analytic result for the average number of steps taken by a random walker from a site which is a first nearest-neighbor to a (deep) trap. Further, they discussed how one would proceed to derive similar analytic expressions starting from sites second, third, th nearest neighbor(s) to the trap, emphasizing that the derived expressions will depend on the structure of the lattice.
In previous work, we pursued the Montroll-Weiss program and obtained additional invariance relations for hexagonal, square planar, and cubic lattices subject to periodic boundary conditions, each having a single, centro-symmetric, deep trap [8–10]. Since satisfaction of lattice invariance relations is an “acid test” on the accuracy of numerical calculations of the site-specific (and overall mean ), we used the methods described in [8–10] to develop invariance relations for finite lattices.
The corresponding first Montroll-Weiss invariance relation for finite lattices holds for a trap positioned either at an interior site or a surface site. The second invariance relation obtained for finite lattices is for lattices with a trap at a corner site and for lattices with a trap at a corner site.
It is important to stress that the right-hand side of each of the above expressions is an integer. Hence, calculating the in integer format, constructing , and , and comparing the results obtained with the (integer) results predicted by the right-hand side of the exact formulae above confirms (or not) the numerical accuracy of our calculations. Table 3 gives the results for , and for all the lattices considered in this study and provides the desired confirmation.
In this contribution we have studied ergodic flows of a random walker undergoing unbiased displacements in a positional phase space represented by a host lattice, and characterized quantitatively the consequences of different metric decompositions of a given, finite parent lattice. As a model for studying the early stages of self-assembly of nanoparticles, the results here complement those reported in [3, 4] where the phase (reaction) space was defined by a host, periodic lattice. Specifically, by breaking the translational symmetry of the reaction space, we find that the reaction efficiency in a finite system is strongly dependent not only on whether the system is compartmentalized, but also on whether the overall reaction space of the microreactor is further decoupled into separable reactors.
Finally, the influence of system dimensionality on the efficiency of kinetic processes carried out in compartmentalized microreactors is quantified. We find that more significant differences are encountered when a reaction space is deconstructed than a reaction space. Phrased in the language of experimental studies on nanoparticle kinetics, our results suggest that reactions carried out on articulated surfaces of a compartmentalized nanoparticle are likely to influence the efficiency of reaction to a greater extent than if the same reaction occurs in the partitioned interior of a nanoassembly.
One of the authors (J. J. Kozak) gratefully appreciates conversations with and the technical assistance of Amelia E. Pawlak, DePaul University. The work of Professor R. A. Garza-López was supported by the Hirsh Research Initiation Grant, the Howard Hughes Medical Institute Research Program, and the Summer Undergraduate Research Program from the Pomona College.
- A. N. Goldstein, Handbook of Nanophase Materials, Dekker, New York, NY, USA, 1977.
- P. Jensen, “Growth of nanostructures by cluster deposition : a review,” Reviews of Modern Physics, vol. 71, no. 5, pp. 1695–1735, 1999.
- J. J. Kozak, C. Nicolis, and G. Nicolis, “Modeling the early stages of self-assembly in nanophase materials,” Journal of Chemical Physics, vol. 126, no. 15, Article ID 154701, 2007.
- J. J. Kozak and G. Nicolis, “Modeling the early stages of self-assembly in nanophase materials. II. Role of symmetry and dimensionality,” Journal of Chemical Physics, vol. 134, no. 6, Article ID 064701, 8 pages, 2011.
- J. J. Kozak, “Chemical reactions and reaction efficiency in compartmentalized systems,” Advances in Chemical Physics, vol. 115, pp. 245–406, 2000.
- R. A. Garza-López, P. Bouchard, G. Nicolis, M. Sleutel, J. Brzezinski, and J. J. Kozak, “Kinetics of docking in postnucleation stages of self-assembly,” Journal of Chemical Physics, vol. 128, no. 11, Article ID 114701, 2008.
- E. W. Montroll and G. H. Weiss, “Random walks on lattices. II,” Journal of Mathematical Physics, vol. 6, no. 2, pp. 167–181, 1965.
- R. A. Garza-López and J. J. Kozak, “Invariance relations for random walks on hexagonal lattices,” Chemical Physics Letters, vol. 371, no. 3-4, pp. 365–370, 2003.
- R. A. Garza-López and J. J. Kozak, “Invariance relations for random walks on square-planar lattices,” Chemical Physics Letters, vol. 406, no. 1–3, pp. 38–43, 2005.
- R. A. Garza-López, A. Linares, A. Yoo, G. Evans, and J. J. Kozak, “Invariance relations for random walks on simple cubic lattices,” Chemical Physics Letters, vol. 421, no. 1–3, pp. 287–294, 2006.
Copyright © 2012 John J. Kozak and Roberto A. Garza-López. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.