Random sequential multilayer deposition of different‐sized k mers on a one dimensional infinite substrate

Kinetics of random sequential, irreversible multilayer deposition of macromolecules of two different sizes on a one dimensional infinite lattice is analyzed at the mean field level. A formal solution for the corresponding rate equation is obtained. The Jamming limits and the distribution of gaps of exact sizes are discussed. In the absence of screening, the jamming limits are shown to be the same for all the layers. A detailed analysis for the components differing by one monomer unit is presented. The small and large time behaviors and the dependence of the individual jamming limits of the k mers and (k−1) mers on k and the rate parameters are analyzed.

Quantifying the influence of humic acid adsorption on colloidal microsphere deposition onto iron-oxide-coated sand.

This article describes an approach for quantifying microsphere deposition onto iron-oxide-coated sand under the influence of adsorbed Suwannee River Humic Acid (SRHA). The experimental technique involved a triple pulse injection of model latex microspheres (microspheres) in pulses of (1) microspheres, followed by (2) SRHA, and then (3) microspheres, into a column filled with iron-coated quartz sand as a water-saturated porous medium. A random sequential adsorption model (RSA) simulated the gradual rise in the first (microsphere) breakthrough curve (BTC). Using the same model calibration parameters a dramatic increase in concentration at the start of the second particle BTC, generated after SRHA injection, could be simulated by matching microsphere concentrations to extrapolated RSA output. RSA results and microsphere/SRHA recoveries showed that 1 mg of SRHA could block 5.90 plus or minus 0.14 x 10^9 microsphere deposition sites. This figure was consistent between experiments injecting different SRHA masses, despite contrasting microsphere deposition/release regimes generating the second microsphere BTC.

Quantifying mircoorganism/organic matter interactions in saturated porous media.

Proteins and humic acids are common constituents of waste water. Latex colloids (colloids) acted as surrogates for microorganisms in multiple pulse dynamic column experiments (MPEs) that permitted colloid mobility to be quantified before and after the injection of either BSA (a protein), or Suwannee River humic acid (SRHA).
At low OM coverage colloid breakthrough curves demonstrated both BSA and SRHA reduced colloid deposition rates, but did not affect colloid irreversible deposition mechanisms. By contrast, high levels of SRHA surface coverage not only further reduced the matrix’s ability to attenuate colloids, but also resulted in reversible adsorption of a significant fraction of colloids deposited. Modelling of colloid responses using random sequential adsorption modelling suggested that 1 microgram of SRHA had the same effect as the deposition of 5.90±0.14 x109 colloids; the model suggested that adsorption of the same mass of BSA was equivalent to the deposition of between 7.1x108 and 2.3x109 colloids.
Colloid responses in MPEs where BSA coverage of colloid deposition sites approached saturation demonstrated the sand matrix remained capable of adsorbing colloids. However, in contrast to responses observed in MPEs at low surface coverage, continued colloid injection showed that the sand’s attenuation capacity increased with time, i.e. colloid concentrations declined as more were deposited (filter ripening).
Importance: Study results highlight the contrasting responses that may arise due to the interactions between colloids and OM in porous media. Results not only underscore that colloids can interact differently with various forms of deposited OM, but also that a single type of OM may generate dramatically different responses depending on the degree of surface coverage. The MPE method provides a means of quantifying the influence of OM on microorganism mobility in porous media such as filter beds, which may be used for either drinking water treatment or waste water treatment. In the wider environment study findings have potential to allow more confident predictions of the mobility of sewage derived pathogens discharging to groundwater.

Exact correlation functions in particle-reaction models with immobile particles

Exact results on particle densities as well as correlators in two models of immobile particles, containing either a single species or else two distinct species, are derived. The models evolve following a descent dynamics through pair annihilation where each particle interacts once at most throughout its entire history. The resulting large number of stationary states leads to a non-vanishing configurational entropy. Our results are established for arbitrary initial conditions and are derived via a generating function method. The single-species model is the dual of the 1D zero-temperature kinetic Ising model with Kimball-Deker-Haake dynamics. In this way, both in finite and semi-infinite chains and also the Bethe lattice can be analysed. The relationship with the random sequential adsorption of dimers and weakly tapped granular materials is discussed.

Modellierung der Oberflächenreaktionen von Vesikeln und Proteinen mit Hilfe dynamischer Monte-Carlo-Simulationen

Die vorliegende Arbeit behandelt die Adsorption von Phospholipidvesikeln, Proteinen und Latex Beads an funktionalisierten festkörpergestützten Oberflächen. Der Schwerpunkt der Abhandlung liegt in der Modellierung der Kinetik mit Hilfe eines Simulationsprogramms unter Berücksichtigung des Massentransports in Staupunktgeometrie, der Adsorption und nachfolgender konformativer Änderungen adhärierter Partikel. Das Simulationsprogramm basiert auf einem RSA (random sequential adsorption) Algorithmus und erlaubt es aufgrund zusätzlich implementierter Optionen, z.B. zur Behandlung elektrostatischer Abstoßung, von Spreitprozessen oder von Desorptionsereignissen, den Adsorptionsprozess unter realistischen physikalischen Bedingungen nachzuempfinden und spezifische biologische Fragestellungen zu beantworten. Aus Anpassungen von Simulationen an experimentelle Daten ergeben sich dynamische Parameter, wie z.B. die Transport-, Adsorptions-, Spreit- und die Desorptionsrate. Die experimentellen Daten wurden mit Hilfe der Quarzmikrowaagetechnik (QCM), der Impedanzspektroskopie (IS) und der Rasterkraftmikroskopie (AFM) erhoben. Zusätzlich zur Kinetik gibt die graphische Ausgabe des Simulationsprogramms Aufschluss über die Oberflächenverteilung der adsorbierten und gespreiteten Partikel. Bei der Untersuchung von Systemen, in denen die Adsorption reversibel ist und infolge dessen Desorptionsprozesse eine wichtige Rolle spielen, wurde mit Hilfe des Simulationsprogramms ein völlig neuartiges Sensorkonzept in Aussicht gestellt. Es wurde gezeigt, dass aus der Analyse des Leistungsspektrums der Fluktuationen in der Belegung die Kinetik von Adsorption und Desorption ermittelt werden kann. Dies ist aus sensorischer Sicht interessant, da eine Vielzahl störender Einflüsse, wie z.B: der Massentransport und die elektronische Drift in diesem Fall nicht berücksichtigt werden müssen.

Untersuchung der lateralen Organisation und Phasenumwandlung von festkörperunterstützten Lipidmembranen

Im ersten Teil der Arbeit wurde das Bindungsverhalten von Annexin A1 und Annexin A2t an festkörperunterstützte Lipidmembranen aus POPC und POPS untersucht. Für beide Proteine konnte mit Hilfe der Fluoreszenzmikroskopie gezeigt werden, dass irreversible Bindung nur in Anwesenheit von POPS auftritt. Durch rasterkraftmikroskopische Aufnahmen konnte die laterale Organisation der Annexine auf der Lipidmembran dargestellt werden. Beide Proteine lagern sich in Form lateraler Aggregate (zweidimensionale Domänen) auf der Oberfläche an, außerdem ist der Belegungsgrad und die Größe der Domänen von der Membranzusammensetzung und der Calciumkonzentration abhängig. Mit zunehmendem POPS-Gehalt und Calciumkonzentration steigt der Belegungsgrad an und der mittlere Domänenradius wird kleiner. Diese Ergebnisse konnten in Verbindung mit detaillierten Bindungsstudien des Annexins A1 mit der Quarzmikrowaage verwendet werden, um ein Bindungsmodell auf Basis einer heterogenen Oberfläche zu entwickeln. Auf einer POPC-reichen Matrix findet reversible Adsorption statt und auf POPS-reichen Domänen irreversible Adsorption. Durch die Anpassung von dynamischen Monte Carlo-Simulationen basierend auf einer zweidimensionalen zufälligen sequentiellen Adsorption konnten Erkenntnisse über die Membranstruktur und die kinetischen Ratenkonstanten in Abhängigkeit von der Calciumkonzentration und der Inkubationszeit des Proteins gewonnen werden. Die irreversible Bindung ist in allen Calciumkonzentrationsbereichen schneller als die reversible. Außerdem zeigt die irreversible Adsorption eine deutlich stärkere Abhängigkeit von der Calciumkonzentration. Ein kleinerer Belegungsgrad bei niedrigen Ca2+-Gehalten ist hauptsächlich durch die Abnahme der verfügbaren Bindungsplätze auf der Oberfläche zu erklären. Die gute Übereinstimmung der aus den Monte Carlo-Simulationen erhaltenen Domänenstrukturen mit den rasterkraftmikroskopischen Aufnahmen und die Tatsache, dass sich die simulierten Resonanzfrequenzverläufe problemlos an die experimentellen Kurven aus den QCM-Messungen anpassen ließen, zeigt die gute Anwendbarkeit des entwickelten Simulationsprogramms auf die Adsorption von Annexin A1. Die Extraktion der kinetischen Parameter aus dem zweidimensionalen RSA-Modell ist mit Sicherheit einem einfachen Langmuir-Ansatz überlegen. Bei einem Langmuir-Modell erfolgt eine integrale Erfassung einer einzelnen makroskopischen Geschwindigkeitskonstante, während durch das RSA-Modell eine differenzierte Betrachtung des reversiblen und irreversiblen Bindungsprozesses möglich ist. Zusätzlich lassen sich mikroskopische Informationen über die Oberflächenbeschaffenheit gewinnen. Im zweiten Teil der Arbeit wurde das thermotrope Phasenverhalten von festkörperunterstützten Phospholipidbilayern untersucht. Dazu wurden mikrostrukturierte, frei stehende Membranstreifen präpariert und mit Hilfe der bildgebenden Ellipsometrie untersucht. Dadurch konnten die temperaturabhängigen Verläufe der Schichtdicke und der lateralen Membranausdehnung parallel beobachtet werden. Die ermittelten Phasenübergangstemperaturen von DMPC, diC15PC und DPPC lagen 2 - 3 °C oberhalb der Literaturwerte für vesikuläre Systeme. Außerdem wurde eine deutliche Verringerung der Kooperativität der Phasenumwandlung gefunden, was auf einen großen Einfluss des Substrats bei den festkörperunterstützten Lipidmembranen schließen lässt. Zusätzlich wurde ein nicht systematischer Zusammenhang der Ergebnisse von der Oberflächenpräparation gefunden, der es unabdingbar macht, bei Untersuchungen von festkörperunterstützten Substraten einen internen Standard einzuführen. Bei der Analyse des thermotropen Phasenübergangsverhaltens von DMPC/Cholesterol - Gemischen wurde daher die individuelle Adressierbarkeit der strukturierten Lipidmembranen ausgenutzt und ein Lipidstreifen aus reinem DMPC als Standard verwendet. Auf diese Weise konnte gezeigt werden, dass das für Phospholipide typische Phasenübergangsverhalten ab 30 mol% Cholesterol in der Membran nicht mehr vorhanden ist. Dies ist auf die Bildung einer nur durch höhere Sterole induzierten fluiden Phase mit hoch geordneten Acylketten zurückzuführen. Abschließend konnte durch die Zugabe von Ethanol zu einer mikrostrukturierten DMPC-Membran die Bildung eines interdigitierten Bilayers nachgewiesen werden. Die bildgebende Ellipsometrie ist eine sehr gute Methode zur Untersuchung festkörperunterstützter Lipidmembranen, da sie über ein sehr gutes vertikales und ein ausreichendes laterales Auflösungsvermögen besitzt. Sie ist darin zwar einem Rasterkraftmikroskop noch unterlegen, besitzt dafür aber eine einfachere Handhabung beim Umgang mit Flüssigkeiten und in der Temperierung, eine schnellere Bildgebung und ist als optische Methode nicht-invasiv.

Three-dimensional distribution of cortical synapses: a replicated point pattern-based analysis

The biggest problem when analyzing the brain is that its synaptic connections are extremely complex. Generally, the billions of neurons making up the brain exchange information through two types of highly specialized structures: chemical synapses (the vast majority) and so-called gap junctions (a substrate of one class of electrical synapse). Here we are interested in exploring the three-dimensional spatial distribution of chemical synapses in the cerebral cortex. Recent research has showed that the three-dimensional spatial distribution of synapses in layer III of the neocortex can be modeled by a random sequential adsorption (RSA) point process, i.e., synapses are distributed in space almost randomly, with the only constraint that they cannot overlap. In this study we hypothesize that RSA processes can also explain the distribution of synapses in all cortical layers. We also investigate whether there are differences in both the synaptic density and spatial distribution of synapses between layers. Using combined focused ion beam milling and scanning electron microscopy (FIB/SEM), we obtained three-dimensional samples from the six layers of the rat somatosensory cortex and identified and reconstructed the synaptic junctions. A total volume of tissue of approximately 4500μm3 and around 4000 synapses from three different animals were analyzed. Different samples, layers and/or animals were aggregated and compared using RSA replicated spatial point processes. The results showed no significant differences in the synaptic distribution across the different rats used in the study. We found that RSA processes described the spatial distribution of synapses in all samples of each layer. We also found that the synaptic distribution in layers II to VI conforms to a common underlying RSA process with different densities per layer. Interestingly, the results showed that synapses in layer I had a slightly different spatial distribution from the other layers.

Red blood cell adhesion on a solid/liquid interface

Red blood cells (RBCs), previously fixed with glutaraldehyde, adhere to glass slides coated with fibrinogen. The RBC deposition process on the horizontal glass surface is investigated by analyzing the relative surface covered by the RBCs, as well as the variance of this surface coverage, as a function of the concentration of particles. This study is performed by optical microscopy and image analysis. A model, derived from the classical random sequential adsorption model, has been developed to account for the experimental results. This model highlights the strong influence of the hydrodynamic interactions during the deposition process.

DRIFTS studies on CO and NO adsorption and NO plus CO reaction over Pd2+-substituted CeO2 and Ce0.75Sn0.25O2 catalysts

Sequential adsorption of CO and NO as well as equimolar NO + CO reaction with variation of temperature over Pd2+ ion-substituted CeO2 and Ce0.75Sn0.25O2 supports has been studied by DRIFTS technique. The results are compared with 2 at.% Pd/Al2O3 containing Pd-0. Both linear and bridging Pd-0-CO bands are observed over 2 at.% Pd/Al2O3. But, band positions are shifted to higher frequencies in Ce0.98Pd0.02O2-delta and Ce0.73Sn0.25Pd0.02O2-delta catalysts that could be associated with Pd delta+-CO species. In contrast, a Pd2+-CO band at 2160 cm(-1) is observed upon CO adsorption over Ce0.98Pd0.02O2-delta and Ce0.73Sn0.25Pd0.02O2-delta catalysts pre-adsorbed with NO and a Pd+-CO band at 2120 cm(-1) is slowly developed on Ce(0.73)Srl(0.25)Pd(0.02)O(2-delta) over time. An intense linear Pd-0-NO band at 1750 cm(-1) found upon NO exposure to CO pre-adsorbed 2 at.% Pd/Al2O3 indicates molecular adsorption of NO. On the other hand, a weak Pd2+-NO band at 1850 cm(-1) is noticed after NO exposure to Ce0.98Pd0.02O2-delta catalyst pre-adsorbed with CO indicating dissociative adsorption of NO which is crucial for NO reduction. Pd-0-NO band is initially formed over CO pre-adsorbed Ce0.73Sn0.25Pd0.02O2-delta which is red-shifted over time along with formation of Pd2+-NO band. Several intense bands related to nitrates and nitrites are observed after exposure of NO to fresh as well as CO pre-adsorbed Ce0.98Pd0.02O2-delta and Ce0.73Sn0.25Pd0.02O2-delta catalysts. Ramping the temperature in a DRIFTS cell upon NO and CO adsorption shows the formation of N2O and NCO surface species, and N2O-formation temperature is comparable with the reaction done in a reactor.

Multifractal characterisation and analysis of complex networks

Complex networks have been studied extensively due to their relevance to many real-world systems such as the world-wide web, the internet, biological and social systems. During the past two decades, studies of such networks in different fields have produced many significant results concerning their structures, topological properties, and dynamics. Three well-known properties of complex networks are scale-free degree distribution, small-world effect and self-similarity. The search for additional meaningful properties and the relationships among these properties is an active area of current research. This thesis investigates a newer aspect of complex networks, namely their multifractality, which is an extension of the concept of selfsimilarity. The first part of the thesis aims to confirm that the study of properties of complex networks can be expanded to a wider field including more complex weighted networks. Those real networks that have been shown to possess the self-similarity property in the existing literature are all unweighted networks. We use the proteinprotein interaction (PPI) networks as a key example to show that their weighted networks inherit the self-similarity from the original unweighted networks. Firstly, we confirm that the random sequential box-covering algorithm is an effective tool to compute the fractal dimension of complex networks. This is demonstrated on the Homo sapiens and E. coli PPI networks as well as their skeletons. Our results verify that the fractal dimension of the skeleton is smaller than that of the original network due to the shortest distance between nodes is larger in the skeleton, hence for a fixed box-size more boxes will be needed to cover the skeleton. Then we adopt the iterative scoring method to generate weighted PPI networks of five species, namely Homo sapiens, E. coli, yeast, C. elegans and Arabidopsis Thaliana. By using the random sequential box-covering algorithm, we calculate the fractal dimensions for both the original unweighted PPI networks and the generated weighted networks. The results show that self-similarity is still present in generated weighted PPI networks. This implication will be useful for our treatment of the networks in the third part of the thesis. The second part of the thesis aims to explore the multifractal behavior of different complex networks. Fractals such as the Cantor set, the Koch curve and the Sierspinski gasket are homogeneous since these fractals consist of a geometrical figure which repeats on an ever-reduced scale. Fractal analysis is a useful method for their study. However, real-world fractals are not homogeneous; there is rarely an identical motif repeated on all scales. Their singularity may vary on different subsets; implying that these objects are multifractal. Multifractal analysis is a useful way to systematically characterize the spatial heterogeneity of both theoretical and experimental fractal patterns. However, the tools for multifractal analysis of objects in Euclidean space are not suitable for complex networks. In this thesis, we propose a new box covering algorithm for multifractal analysis of complex networks. This algorithm is demonstrated in the computation of the generalized fractal dimensions of some theoretical networks, namely scale-free networks, small-world networks, random networks, and a kind of real networks, namely PPI networks of different species. Our main finding is the existence of multifractality in scale-free networks and PPI networks, while the multifractal behaviour is not confirmed for small-world networks and random networks. As another application, we generate gene interactions networks for patients and healthy people using the correlation coefficients between microarrays of different genes. Our results confirm the existence of multifractality in gene interactions networks. This multifractal analysis then provides a potentially useful tool for gene clustering and identification. The third part of the thesis aims to investigate the topological properties of networks constructed from time series. Characterizing complicated dynamics from time series is a fundamental problem of continuing interest in a wide variety of fields. Recent works indicate that complex network theory can be a powerful tool to analyse time series. Many existing methods for transforming time series into complex networks share a common feature: they define the connectivity of a complex network by the mutual proximity of different parts (e.g., individual states, state vectors, or cycles) of a single trajectory. In this thesis, we propose a new method to construct networks of time series: we define nodes by vectors of a certain length in the time series, and weight of edges between any two nodes by the Euclidean distance between the corresponding two vectors. We apply this method to build networks for fractional Brownian motions, whose long-range dependence is characterised by their Hurst exponent. We verify the validity of this method by showing that time series with stronger correlation, hence larger Hurst exponent, tend to have smaller fractal dimension, hence smoother sample paths. We then construct networks via the technique of horizontal visibility graph (HVG), which has been widely used recently. We confirm a known linear relationship between the Hurst exponent of fractional Brownian motion and the fractal dimension of the corresponding HVG network. In the first application, we apply our newly developed box-covering algorithm to calculate the generalized fractal dimensions of the HVG networks of fractional Brownian motions as well as those for binomial cascades and five bacterial genomes. The results confirm the monoscaling of fractional Brownian motion and the multifractality of the rest. As an additional application, we discuss the resilience of networks constructed from time series via two different approaches: visibility graph and horizontal visibility graph. Our finding is that the degree distribution of VG networks of fractional Brownian motions is scale-free (i.e., having a power law) meaning that one needs to destroy a large percentage of nodes before the network collapses into isolated parts; while for HVG networks of fractional Brownian motions, the degree distribution has exponential tails, implying that HVG networks would not survive the same kind of attack.

Dual Drug Delivery Microcapsules via Layer-by-Layer Self-Assembly

The integration of hydrophobic and hydrophilic drugs in the polymer microcapsule offers the possibility of developing a new drug delivery system that combines the best features of these two distinct classes of material. Recently, we have reported the encapsulation of an uncharged water-insoluble drug in the polymer membrane. The hydrophobic drug is deposited using a layer-by-layer (LbL) technique, which is based on the sequential adsorption of oppositely charged polyelectrolytes onto a charged substrate. In this paper, we report the encapsulation of two different drugs, which are invariably different in structure and in their solubility in water. We have characterized these dual drug vehicular capsules by confocal laser scanning microscopy, atomic force microscopy, visible microscopy, and transmission electron microscopy. The growth of a thin film on a flat substrate by LbL was monitored by UV−vis spectra. The desorption kinetics of two drugs from the thin film was modeled by a second-order rate model.

LbL-coated microcapsules as systems for encapsulation of optical brightening agent

We report the encapsulation of optical brightening agent (OBA) into hollow microcapsules prepared by the controlled Layer- by-Layer (LbL) self-assembly process, achieved by the sequential adsorption of oppositely charged polyelectrolytes using negatively charged silica template. Loading takes place by spontaneous deposition method which was proved by confocal laser scanning microscopy (CLSM) using rhodamine 6G (Rd6G) as a fluorescent probe. The loading of the OBA into the microcapsules was found to be dependent on the feeding concentration, pH of the medium, and loading temperature. The encapsulation efficiency of OBA decreased on increasing feeding concentration. Maximum loading was observed at pH 4 and amount of OBA loaded decreased with increase in pH. The loaded OBA was released in a sustained manner for 8 h. No degradation of the OBA was observed during the process of encapsulation and release. Polyelectrolyte capsules potentially offer an innovative way of encapsulating large amounts of active materials for a variety of applications. (c) 2012 Wiley Periodicals, Inc. J. Appl. Polym. Sci. 127: 1609-1614, 2013

Estudo das propriedades críticas do processo de contato por par para diferentes atualizações

We study the critical behavior of the one-dimensional pair contact process (PCP), using the Monte Carlo method for several lattice sizes and three different updating: random, sequential and parallel. We also added a small modification to the model, called Monte Carlo com Ressucitamento" (MCR), which consists of resuscitating one particle when the order parameter goes to zero. This was done because it is difficult to accurately determine the critical point of the model, since the order parameter(particle pair density) rapidly goes to zero using the traditional approach. With the MCR, the order parameter becomes null in a softer way, allowing us to use finite-size scaling to determine the critical point and the critical exponents β, ν and z. Our results are consistent with the ones already found in literature for this model, showing that not only the process of resuscitating one particle does not change the critical behavior of the system, it also makes it easier to determine the critical point and critical exponents of the model. This extension to the Monte Carlo method has already been used in other contact process models, leading us to believe its usefulness to study several others non-equilibrium models

Purification and properties of buffalo (Bubalus bubalis) erythrocyte hexokinase

Buffalo erythrocytes contain one isozyme of hexokinase that apparently lacks microheterogeneity as shown by chromatographic properties. A single protein band was detected by means of Western blotting using an antibody raised in rabbits against homogeneous rat brain hexokinase I. The native protein has a molecular weight of 200,000 +/- 2880 by gel filtration. Partial purification of erythrocyte hexokinase by a combination of several procedures, including affinity chromatography, which was previously applied successfully to the purifica tion of other mammalian type I hexokinases, produced a partially purified enzyme that showed several contami nants after SDS-polyacrylamide gel electrophoresis. The affinity of buffalo erythrocyte hexokinase for glucose (K-m = 0.012 +/- 0.001 mM) is lower than most other mammal hexokinases type I. It phosphorylates other sugars, with considerably higher K-m values. This isozyme is able to use MgATP but does not use MgGTP, MgCTP or MgUTP. We used inhibition patterns, obtained with products to elucidate enzyme sequential mechanisms. Our results are clearly in agreement with a random sequential mechanism and in disagreement with an ordered sequential mechanism with either glucose or ATP as the obligatory first substrates. The ADP inhibition was of mixed type with both ATP and glucose as substrates. (C) 1997 Elsevier B.V.

Three-dimensional spatial distribution of synapses in the neocortex: A dual-beam electron microscopy study

In the cerebral cortex, most synapses are found in the neuropil, but relatively little is known about their 3-dimensional organization. Using an automated dual-beam electron microscope that combines focused ion beam milling and scanning electron microscopy, we have been able to obtain 10 three-dimensional samples with an average volume of 180 µm(3) from the neuropil of layer III of the young rat somatosensory cortex (hindlimb representation). We have used specific software tools to fully reconstruct 1695 synaptic junctions present in these samples and to accurately quantify the number of synapses per unit volume. These tools also allowed us to determine synapse position and to analyze their spatial distribution using spatial statistical methods. Our results indicate that the distribution of synaptic junctions in the neuropil is nearly random, only constrained by the fact that synapses cannot overlap in space. A theoretical model based on random sequential absorption, which closely reproduces the actual distribution of synapses, is also presented.