Percolation thresholds of two-dimensional continuum systems of rectangles

Lattice model for percolation on a plane of partially aligned sticks with length dispersity

A lattice-based model for continuum percolation is applied to the case of randomly located, partially aligned sticks with unequal lengths in 2D which are allowed to cross each other. Results are obtained for the critical number of sticks per unit area at the percolation threshold in terms of the distributions over length and orientational angle and are compared with findings from computer simulations. Consistent with findings from computer simulations, our model shows that the percolation threshold is (i) elevated by increasing degrees of alignment for a fixed length distribution, and (ii) lowered by increasing degrees of length dispersity for a fixed orientational distribution. The impact of length dispersity is predicted to be governed entirely by the first and second moments of the stick length distribution, and the threshold is shown to be quite sensitive to particulars of the orientational distribution function.

Continuum percolation of two-dimensional adaptive dynamics systems

The percolation phase transition of a continuum adaptive neuron system with homeostasis is investigated. In order to maintain their average activity at a particular level, each neuron (represented by a disk) varies its connection radius until the sum of overlapping areas with neighboring neurons (representing the overall connection strength in the network) has reached a fixed target area for each neuron. Tuning the two key parameters in the model, i.e., the density defined as the number of neurons (disks) per unit area and the sum of the overlapping area of each disk with its adjacent disks, can drive the system into the critical percolating state. These two parameters are inversely proportional to each other at the critical state, and the critical filling factors are fixed about 0.7157, which is much less than the case of the continuum percolation with uniform disks. It is also confirmed that the critical exponents in this model are the same as the two-dimensional standard lattice percolation. Although the critical state is relatively more sensitive and exhibits long-range spatial correlation, local fluctuations do not propagate in a long-range manner through the system by the adaptive dynamics, which renders the system overall robust against perturbations.

Percolation of straight slots on a square grid

This work is devoted to the emergence of a connected network of slots (cracks) on a square grid. Accordingly, extensive Monte Carlo simulations and finite-size scaling analysis have been conducted to study the site percolation of straight slots with length l measured in the number of elementary cells of the grid with the edge size L. A special focus was made on the dependence of the percolation threshold p_{C}(l,L) on the slot length l varying in the range 1≤l≤L-2 for the square grids with edge size in the range 50≤L≤1000. In this way, we found that p_{C}(l,L) strongly decreases with increase of l, whereas the variations of p_{C}(l=const,L) with the variation of ratio l/L are very small. Consequently, we acquire the functional dependencies of the critical filling factor and percolation strength on the slot length. Furthermore, we established that the slot percolation model interpolates between the site percolation on square lattice (l=1) and the continuous percolation of widthless sticks (l→∞) aligned in two orthogonal directions. In this regard, we note that the critical number of widthless sticks per unit area is larger than in the case of randomly oriented sticks. Our estimates for the critical exponents indicate that the slot percolation belongs to the same universality class as standard Bernoulli percolation.

Numerical modeling of the effects of the shape and aspect ratio of 3D curved fiber on the percolation threshold and electrical conductivity of conductive polymer composites

For designing conductive polymer composites (CPCs), understanding how the fiber curvature affects the percolation behavior of curved conductive fibers is essential for determining the effective electrical conductivity σeff of the CPCs. In this work, CPCs were considered as a polymer matrix filled with the random packing of overlapped curved spherocylinders. The geometries of the curved spherocylinders were defined, and inter-curved spherocylinder contact-detecting and system-spanning fiber cluster searching algorithms were developed. The finite-size-scaling method was used to explore how the aspect ratio α and bending central angle θ of a curved spherocylinder affect the percolation threshold ϕc of an overlapped curved spherocylinder system in 3D space. The findings suggest that ϕc decreases as α increases and increases initially before declining as θ increases. An empirical approximation formula was proposed to quantify the effect of the curved spherocylinder's morphology, characterized by the dimensionless excluded volume Vdex of the curved spherocylinder, on ϕc. The new rigorous bound for ϕc of the soft-curved spherocylinder system was further proposed. A random resistor network model was constructed, and the reliability of this model was validated by comparing the simulations and published data. Finally, a fitting formula was developed to assess the impacts of the normalized reduced density (η - ηc)/ηc and Vdex on the σeff of CPCs. A distinct linear correlation between σeff and (η - ηc)/ηc was constructed, denoted as σeff ∼ [(η - ηc)/ηc]t(α,θ). An empirical approximation model was proposed to establish the relationship between the fiber shape and conductivity exponent t. Our study may provide a theoretical hint for the design of CPCs.

Percolation threshold and electrical conductivity of conductive polymer composites filled with curved fibers in two-dimensional space

Quantifying the influence of fiber curvature on the percolation behavior of flexible conductive fiber and further on the electrical conductivity of conductive polymer composites (CPCs) is crucial for the design of CPCs. This study considers CPCs as a random packing of soft curved discorectangles (CDCRs) in a polymer matrix. The geometry of CDCR is developed, and an inter-CDCR contact detection algorithm is used to generate a random packing structure of CDCRs. The effects of aspect ratio α and bending central angles θ of CDCR on the percolation threshold ϕc of the overlapped CDCR system in a two-dimensional plane are then investigated using the finite-size scaling method. The result reveals that ϕc decreases monotonically as α grows and increases monotonically as θ rises. A shape-independent power law formula, denoted as ϕc = 2.2015 A-0.8172dex is developed to quantify the relationship between the Adex and ϕc. A comparison of our numerical simulations, published data, and predictions verifies the reliability and universality of the fitting model. Subsequently, a resistor network searching algorithm (RNSA) is proposed to construct the random resistor network model (RRNM). A power law model, denoted as is developed to evaluate the effects of the normalized reduced density (η - ηc)/ηc on the effective conductivity σeff of CPC. Comparing our predictions with data from the literature and our simulation verifies the reliability of our RNSA and the fitting model. This paper's methodology and findings may provide a theoretical hint for the CPC's design.

Random nanowire networks: Identification of a current-carrying subset of wires using a modified wall follower algorithm

We mimic random nanowire networks by the homogeneous, isotropic, and random deposition of conductive zero-width sticks onto an insulating substrate. The number density (the number of objects per unit area of the surface) of these sticks is supposed to exceed the percolation threshold, i.e., the system under consideration is a conductor. To identify any current-carrying part (the backbone) of the percolation cluster, we have proposed and implemented a modification of the well-known wall follower algorithm-one type of maze solving algorithm. The advantage of the modified algorithm is its identification of the whole backbone without visiting all the edges. The complexity of the algorithm depends significantly on the structure of the graph and varies from O(sqrt[N_{V}]) to Θ(N_{V}). The algorithm has been applied to backbone identification in networks with different number densities of conducting sticks. We have found that (i) for number densities of sticks above the percolation threshold, the strength of the percolation cluster quickly approaches unity as the number density of the sticks increases; (ii) simultaneously, the percolation cluster becomes identical to its backbone plus simplest dead ends, i.e., edges that are incident to vertices of degree 1. This behavior is consistent with the presented analytical evaluations.

Percolation threshold of curved linear objects

In this study, we investigate the percolation threshold of curved linear objects, describing them as quadratic Bézier curves. Using Monte Carlo simulations, we calculate the critical number densities of the curves with different curviness. We also obtain the excluded area of the curves. When an excluded area is given, we can find the critical number density of the curves with arbitrary curviness. Apparent conductivity exponents are computed for the curves, and these values are found to be analogous to that of sticks in the percolative region for a junction resistance dominant system. These results can be used to analyze the optoelectrical performance of metal nanowire films because the high-aspect-ratio metal nanowires can be easily curved during coating.

Percolation thresholds for discorectangles: Numerical estimation for a range of aspect ratios

Using Monte Carlo simulation, we have studied the percolation of discorectangles. Also known as stadiums or two-dimensional spherocylinders, a discorectangle is a rectangle with semicircles at a pair of opposite sides. Scaling analysis was performed to obtain the percolation thresholds in the thermodynamic limits. We found that (i) for the two marginal aspect ratios ɛ=1 (disc) and ɛ→∞ (stick) the percolation thresholds coincide with known values within the statistical error and (ii) for intermediate values of ɛ the percolation threshold lies between the percolation thresholds for ellipses and rectangles and approaches the latter as the aspect ratio increases.

Relaxation in two-dimensional suspensions of rods as driven by Brownian diffusion

Relaxation in a two-dimensional suspensions containing rods was studied by using dynamic Monte Carlo simulations. An off-lattice model with continuous positional and orientational degrees of freedom was considered. The initial state was produced by using a random sequential adsorption model. During the relaxation, the rods underwent translational and rotational Brownian motions. The simulations were run at different values of number density ρ (the number of rods per unit area) and of the initial orientation order parameter S_{i}. The rods were assumed to have core-shell structures. The evolutions of both the connectivity and the order parameter have revealed different relaxation behavior.

Anisotropy in electrical conductivity of films of aligned intersecting conducting rods

Numerical simulations by means of the Monte Carlo method have been performed to study the electrical properties of a two-dimensional composite filled with rodlike particles. The main goal was to study the effect of the alignment of such rods on the anisotropy of its electrical conductivity. A continuous model was used. In this model, the rods have zero-width (i.e., infinite aspect ratio) and they may intersect each other. To involve both the low conductive host matrix and highly conductive fillers (rods) in the consideration, a discretization algorithm based on the use of a supporting mesh was applied. The discretization is equivalent to the substitution of rods with the polyominoes. Once discretized, the Frank-Lobb algorithm was applied to evaluate the electrical conductivity. Our main findings are (i) the alignment of the rods essentially affects the electrical conductivity and its anisotropy, (ii) the discrete nature of computer simulations is crucial. For slightly disordered system, high electrical anisotropy was observed at small filler content, suggesting a method to enable the production of optically transparent and highly anisotropic conducting films.

Anisotropy in electrical conductivity of two-dimensional films containing aligned nonintersecting rodlike particles: Continuous and lattice models

The electrical conductivity of two-dimensional films filled with rodlike particles (rods) was simulated by the Monte Carlo method. The main attention has been paid to the investigation of the effect of the rod alignment on the electrical properties of the films. Both continuous and lattice approaches were used. Intersections of particles were forbidden. Our main findings are (i) both models demonstrate similar behaviors, (ii) at low concentration of rods, both approaches lead to the same dependencies of the electrical conductivity on the concentration of the rods, (iii) the alignment of the rods essentially affects the electrical conductivity, (iv) at some concentrations of partially aligned rods, the films may be conducting only in one direction, and (v) the films may simultaneously be both highly transparent and electrically anisotropic.

Composites Formed from Thermoresponsive Polymers and Conductive Nanowires for Transient Electronic Systems

The disintegration of transient electronic systems after a preprogrammed time or a particular stimulus (e.g., water, light, or temperature) is fundamentally linked to the properties and behavior of the materials used for their construction. Herein, we demonstrate that polymers exhibiting lower critical solution temperature (LCST) behavior can work as thermoresponsive substrates for circuitry and that these materials can be coupled with conductive nanowires to form a transient electronics platform with unique, irreversible temperature-responsive behavior. The transient systems formed from composites of LCST polymers and conductive nanowires exhibit stable electrical performance in solution (T_solution > LCST) for over 24 h until a cooling stimulus triggers a rapid (within 5 min) and gigantic (3-4 orders of magnitude) transition in electrical conductance due to polymer dissolution. Using a parylene mask, we are able to fabricate thermoresponsive electrical components, such as conductive traces and parallel-plate capacitors, demonstrating the versatility of this material and patterning technique. With this unique stimulus-responsive transient system and polymers with LCSTs above room temperature (e.g., poly(N-isopropylacrylamide), methyl cellulose), we have developed a platform in which a circuit requires a source of heat to remain viable and is destroyed and vanishes once this heat source is lost.

Experimental percolation studies of random networks

This report establishes an experimental method of studying electrically percolating networks at a higher resolution than previously implemented. This method measures the current across a conductive sheet as a function of time as elliptical pores are cut into the sheet. This is done utilizing a Universal Laser System X2-600 100 W CO₂ laser system with a 76 × 46 cm² field and 394 dpc (dots/cm) resolution. This laser can cut a random system of elliptical pores into a conductive sheet with a potential voltage applied across it and measures the current versus time. This allows for experimental verification of a percolation threshold as a function of the ellipse's aspect ratio (minor/major diameter). We show that as an ellipse's aspect ratio approaches zero, the percolation threshold approaches one. The benefit of this method is that it can experimentally measure the effect of removing small pores, as well as pores with complex geometries, such as an asterisk from a conductive sheet.

Continuum percolation of congruent overlapping spherocylinders

Continuum percolation of randomly orientated congruent overlapping spherocylinders (composed of cylinder of height H with semispheres of diameter D at the ends) with aspect ratio α=H/D in [0,∞) is studied. The percolation threshold ϕ_{c}, percolation transition width Δ, and correlation-length critical exponent ν for spherocylinders with α in [0, 200] are determined with a high degree of accuracy via extensive finite-size scaling analysis. A generalized excluded-volume approximation for percolation threshold with an exponent explicitly depending on both aspect ratio and excluded volume for arbitrary α values in [0,∞) is proposed and shown to yield accurate predictions of ϕ_{c} for an extremely wide range of α in [0, 2000] based on available numerical and experimental data. We find ϕ_{c} is a universal monotonic decreasing function of α and is independent of the effective particle size. Our study has implications in percolation theory for nonspherical particles and composite material design.

Exceptional stiffening in composite fiber networks

We study the small strain elastic behavior of composite athermal fiber networks constructed by adding stiffer fibers to a cross-linked base network. We observe that if the base network is in the affine deformation regime, the composite behaves similar to a fiber-reinforced continuum. When the base network is in the nonaffine deformation regime, the stiffness of the composite increases by orders of magnitude upon the addition of a small fraction of stiff fibers. The increase is not gradual, but rather occurs in two steps. Of these, one is associated with the stiffness percolation of the network of added fibers. The other, which occurs at very small fractions of stiff fibers, is due to the percolation of perturbation zones, or "interphases," induced in the base network by the stiff fibers, regions where the energy is stored mostly in the axial deformation mode. Their size controls the stiffening transition and depends on base network parameters and the length of added fibers. It is also shown that the perturbation field introduced in the base network by the presence of a stiff fiber is much longer ranged than in the case when the fiber is tied to a continuum of same modulus with the base network.

Conductivity scaling in supercritical percolation of nanoparticles--not a power law

The power-law behavior widely observed in supercritical percolation systems of conductive nanoparticles may merely be a phenomenological approximation to the true scaling law not yet discovered. In this work, we derive a comprehensive yet simple scaling law and verify its extensive applicability to various experimental and numerical systems. In contrast to the power law which lacks theoretical backing, the new scaling law is explanatory and predictive, and thereby helpful to gain more new insights into percolation systems of conductive nanoparticles.

The structure of adsorbed cyclic chains

In order to determine the structure of polymer films formed of cyclic chains (rings) we developed and studied a simple coarse-grained model. Our main goal was to check how the percolation and jamming thresholds in such a system were related to the thresholds obtained for linear flexible chains system, i.e., how the geometry of objects influenced both thresholds. All atomic details were suppressed and polymers were represented as a sequence of identical beads and the chains were embedded to a square lattice (a strictly 2D model). The system was athermal and the excluded volume was the only potential introduced. A random sequential adsorption algorithm was chosen to determine the properties of a polymer monolayer. It was shown that the percolation threshold of cyclic chains was considerably higher than those of linear flexible chains while the jamming thresholds for both chain architectures are very similar. The shape of adsorbed cyclic chains was found to be more prolate when compared to average single chain.