abstract: We study the passage of active and passive granular particles through a bottleneck under gravitational bias. The grains are pharmaceutical capsules with nearly spherocylindrical shapes on a vibrating table. The vibrating plate is slightly tilted in order to break the in-plane symmetry and to give particles a motivation to move in a preferential direction on the plate. The passage through a narrow gate with openings comparable to the grain length is studied using video imaging, and particle positions and velocities are extracted from the recorded frames. We compare the behaviour of asymmetrically filled, active capsules with symmetrically filled, passive ones.
author: Mahdieh Mohammadi$^{1,2}$, Kirsten Harth$^1$, Dmitry Puzyrev$^1$, Tina Hanselka$^1$,
Torsten Trittel$^1$, and Ralf Stannarius$^1$
$^1$ Institute of Physics, Otto von Guericke University Magdeburg,
39106 Magdeburg, Germany
$^2$ Department of Physics, Institute for Advanced Studies in Basic Sciences (IASBS),
Zanjan, Iran
date: 2024-10-02
title: A simple model system for self-propelled particles passing a bottleneck

# Introduction

The evacuation of a crowd of people from a room, egress of animals from a corral and discharge of granular material from a silo are, for different reasons, very active, interrelated topics of current research. Seeking the similarities and differences in these processes promises a more universal insight into their underlying mechanisms. On the one hand, granular systems can surely not mimic the complex, individual behavior of spontaneously decision-making people or animals, but may nevertheless reproduce fundamental features of egress studies. On the other hand, conducting experiments or field studies with humans, or with animals, is restricted by ethical principles and laws, requiring tremendous efforts and costs. There are usually restrictions to safe situations (no panic) with comparatively few participants, or one has to refer to data of coincidentally recorded events. The outcome is influenced by psychological components which are hard to determine quantitatively, and it usually depends on social, cultural and professional background, as well as the age of the test persons . Computer simulations, succeeding the works of Helbing et al. , often capture only a fraction of the complexity, or may require substantial computation effort. This motivates studies of less complex, inanimate objects that reproduce some important features of these processes. The properties of such mechanical objects are easier to control, and experiments can be readily handled and reproduced in a laboratory. One such example are vibration-driven granular particles, e. g. . The usefulness of such self-propelled rod-like objects to mimic the behavior of biological systems in vitro has been comprehensively discussed recently in a review article by Bär et al. . A good overview of recent studies of pedestrian dynamics can be found in Refs. . A more general review of the dynamics of systems of active particles was provided by Bechinger et al. . For a comparison of living and granular systems see, for example, Refs. .

We will first recollect some similarities between the passage of living objects through narrow outlets and the discharge of grains from a silo. Emptying of containers through their orifices can follow two different dynamic scenarios : In the first one, which occurs at sufficiently large outlet widths, the outflow is practically continuous until the container is practically empty, without considerable flow rate fluctuations. The second one is characterized by transient or permanent blockages of the outlet, usually undesired. Such clogging or strong fluctuations of the flow rate cause severe problems, e.g. in construction, agricultural, or pharmaceutical industries. In crowds of humans or animals, an even more severe consequence is a build-up of high pressure near the outlet , which may lead to casualties and unresolvable blockages of the door.

In both biological and granular systems, not only the prevention of clogs, but also a minimization of the total time needed for emptying the room or container (the egress time) and the avoidance of strong pressure fluctuations are of substantial importance. For people or animals moving through constrictions, the actual *motivation* of the individuals is expected to strongly affect the pressure near the outlet, as qualitatively confirmed by Pastor et al. . Similar observations were made in experiments with animals such as, e. g., sheep, mice or ants , and in numerical simulations  . Obviously, this intrinsic motivation to pass the outlet is hard to quantify or even to control in experiments with living creatures. For sufficiently wide exits, without congestions, the outflow rate (in persons per door width $d$ and time$t$) is commonly assumed to be independent of the exit width  (but see also ).

Whereas in living systems, the motivation to pass an outlet is set by a certain stimulus of the individuals, studies of granular material discharging from hoppers and silos typically rely exclusively on gravity forcing the grains toward the exit. As long as a certain fill height of the container is exceeded, the outflow rate is independent of the fill level and satisfactorily described by the Beverloo equation . This applies in the free flow regime as well as during avalanches between clogs . It predicts an outflow rate that is proportional to $d^{1.5}$ in a 2D system.

The pressure at the bottom of the silo saturates above a certain height related to the container diameter, because forces are transferred to the side walls (Janssen effect) . The actual pressure depends on the geometry of the container near the outlet . The outflow and clogging characteristics are affected by inter-particle friction, particle shape and deformability. A prominent feature of living systems that distinguishes them from grains in a classical hopper experiment is their activity and their intrinsic ability to perform deliberate rearrangements of their positions and orientations in situations perceived unfavorable.

Some features of biological systems can be simulated by using (almost) frictionless, deformable constituents, such as emulsion drops  or viscoelastic hydrogel spheres . Then, permanent clogging occurs only at substantially smaller opening widths compared to hard frictional grains. Interestingly, slow particle rearrangements in the container can evoke intermittent flow with non-permanent congestions , and the pressure near the outlet can vary substantially in time - a similarity to animals and pedestrians .

Here, we report experiments with elongated grains with specially prepared mass distributions, which are positioned on a flat plate mounted on an electro-mechanical shaker. We drive these spherocylinders by vertical vibrations of the plate. The driving frequency and amplitude set the diffusion parameters of the spherocylinders. At the same time, we break the in-plane symmetry of the shaker dish by tilting the setup by a few degrees, defining a preferred direction of particle motion. Particles are initially placed in a restricted area with a single exit. The 'motivation' to pass this bottleneck exit is controlled by the tilt angle of the setup that sets an effective gravitational force towards the exit. Experiments are conducted with symmetric grains that exhibit purely diffusive dynamics, and asymmetrically prepared, self-propelled grains. We compare the egress behavior of the active grains to the equally shaped but purely diffusive ones. The behavior of this purely mechanical model system is compared to features of colloidal systems and the biological systems described above.

Before the experiment is described in detail, it is necessary to add a few remarks to differentiate our experiment from experiments with vibrated hoppers. Those, at first glance, appear to have some resemblance to the present setup: Our sample container, viewed from above, looks quite similar to a quasi two-dimensional hopper. Indeed, the particles essentially remain in a monolayer. The important difference is that gravity does not act primarily in the hopper plane. In our setup, this plane is almost horizontal, nearly perpendicular to the gravitational force. There is still some small effective gravitation in the direction of the orifice when the shaker is tilted, which is needed for the experiment. However, the effective gravity is too small to overcome friction, so that the particles do not slide downward when the shaker is off.

Several studies of vibrated hoppers probed the stability and the destruction of clogs formed at the orifice, denoted as unjamming. Some other studies applied vibrations during silo discharge in order to characterize either the inhibition or the support of discharge rates . Regarding the present study, the consequences of the vibrations are considerable different. In all of these earlier experiments, the grains are driven downward towards the orifice by the gravitational force, and the vibrations in the same direction break arches and prevent clogging. In our experiments, the oscillations are perpendicular to the slope of the container and serve only as the energy supply to generate diffusive or directed particle motion.

Similar to typical biological situations, the particles we use are substantially elongated. The prolate shape of the grains strongly affects their diffusion characteristics . Shape-anisotropy can also affect the packing in denser ensembles (see e.g. Refs. ). In particular, elongated constituents tend to prefer alignment towards the exit . Elongated shapes can suppress motility-induced phase separation . Caging effects cause a decrease of perpendicular diffusion compared to axial diffusion .

This manuscript is organized as follows: We first describe the experimental setup, the particles and the image evaluation methods. Then, we characterize the dynamics of individual active or passive particles on a horizontal plate, in dependence of the excitation parameters. After this, the outflow characteristics of the grains are characterized for specific driving parameters and tilt. Representative orientational properties near the orifice are shown. Finally, we discuss our results in the context of typical exit scenarios from the literature.

# Experimental setup, materials and methods

## Particle preparation

We employed two types of particles with identical shape and size: 'active' and 'passive' capsules. The particles are excited with a vibrating table, and their dynamic behaviour differs qualitatively in that the 'passive' particles show an anisotropic diffusive behaviour with preferential motion along their long axes, but without preferred sense of direction, while the 'active' particles exhibit an additional directed motion along their long axes. Details are discussed in Sec. .

The shells of all particles are commercial pharmaceutical spherocylindrical capsules of 21 mm length and 7 mm diameter (Figure ). The artificial active particles contain a copper-plated lead bullet (air rifle ammunition) of 0.48 g mass which is asymmetrically glued inside the capsule at one end. The unequal mass distribution is the basis for the self propulsion on the shaken plane. Since this bullet has a smaller diameter than the capsule, it breaks the rotational symmetry of the particle so that sideways rolling of the capsule is inhibited. Active capsules have two colors, white and blue, in order to distinguish their heavier head (white) from the lighter tail (blue) in image processing.

Passive capsules were prepared with same size and total mass as the active ones, but there are two fundamental differences to distinguish them from the active ones: First, the passive particles contain a thin aluminum bar of 15 mm length which is glued symmetrically along the major axis of the capsules. Their motion is symmetric with respect to the long axis of the capsule. Second, the passive particles are mono colored (red). They appear symmetric in image processing.

## Particle excitation and observation

The particles are excited by vertical vibrations of a shaker dish with flat bottom and a slightly elevated border in order to prevent them from leaving the dish. Figure a shows the side view of the experimental setup, without the cameras. The shaker is driven by a sinusoidal wave form, whose amplitude and frequencies can be adjusted in a custom-made Labview program. In order to provide an appropriate excitation for the particles, we adjusted the acceleration to values slightly larger than the gravitational acceleration ($g=9.81$ m/s$^2$) to 1.4 $g$, 1.6 $g$ and 1.8 $g$. The corresponding vibration amplitudes were in the range of a few dozen micrometers, the excitation frequencies in the range between 50 and 100 Hz. Below the above given minimum excitation in our experiment, the motion of the particles is very slow. Above 1.8 $g$, the excitation becomes too strong, the particles begin to jump, to stand vertically frequently and topple over in the tilted hopper, and to leave the horizontal plate.

In order to observe the behaviour of the particles near a bottleneck we manufactured trianglular plates with different geometries serving as sample containers (artificial 2D 'hopper'). Figure b shows the top view of the shaker with the hopper compartment filled with particles. This metal sheet provides the confinement in our experiment. Seen from above, this sheet has funnel shape, its three edges are bent upward to serve as barriers for the particles. At one side where two edges approach each other, they leave a gap which serves as the bottleneck the particles have to pass.

To motivate the particles to move through the gap, the whole setup including the shaker and cameras is slightly tilted. We use a car jack to adjust angles of a few degrees with an accuracy of approximately 0.15$^\circ$. The desired slopes in different experiments were between 2.5$^\circ$ and 4$^\circ$. At larger tilt angles, particles tend to jump and show a tendency to rise with one end above neighbors. At angles smaller than 2$^\circ$, there is no noticeable trend of the capsules to move downward within the experimental time of approximate one minute.

This compartment is fixed onto the vibrating dish, with a few millimeters height gap between them. Then, the particles leaving the hopper through the gap drop down onto the dish and do not block the orifice after they passed it. The top edge of the compartment is 190 mm wide, its vertical distance to the orifice is 110 mm. The opening sizes were chosen as 15 mm, 18 mm and 21 mm in individual experiments. Figure shows the top view of the compartment initially containing the particles. The direction of the tilt towards the orifice breaks the symmetry for particle motion in the plane and defines the mean preferential direction of motion of both active and passive particles.

The particles were initially filled into the compartment while the shaker was at rest. The two images in Fig.  show typical initial configurations of the compartment filled with active and passive capsules. Then, the shaker was switched on and the emptying of the silo during excitation was recorded with optical video imaging. In order to characterize the individual particles, we have additionally performed an analysis of single capsules moving on a vibrating horizontal plate with the same excitation parameters but zero tilt. These measurements show the qualitatively different kinematic behaviour of both species.

The collective bottleneck passage experiments were recorded with a Canon EOS customer camera at a frame rate of 60 fps with a resolution of 4 pixels/mm. The setup was illuminated with diffuse light in order to avoid reflexes at the particle surfaces which are undesired in image processing. The diffusion experiments with single capsules were performed with a XiaoYi 4K camera with a resolution of 26 pixel/mm at 100 fps. The same camera and frame rate were used to record close-ups of the orifice to determine the time intervals between the passages of active particles. The typical duration of a discharge experiment was of the order of one minute. The video shown in the supplemental information shows downgraded footage of a typical experiment.

## Image processing and analysis

In order to obtain the trajectories of the capsules, one first needs to detect the particles by means of image processing and analysis. The detection of single particles in a video frame can be performed by separation of colors in the image. For a passive particle, its center of mass and orientation are then easily extracted using MATLAB as segmented region properties. For an active particle, one can find centroids of colored (blue and white) regions and the orientation of the vector between them. Moreover, if particles in the image are lying horizontally on the plate and do not overlap in the camera view, simultaneous detection of multiple capsules is possible. The regions corresponding to each single passive particle are easily separated from each other, so the detection of multiple particles is easy in this case.

In order to reliably find the positions and orientations of multiple active particles, one has to perform an additional step to match image regions corresponding to the blue and white capsule parts of one and the same capsule. The optimal assignment problem can be solved so that all vector lengths between matched blue and white regions' centroids correspond to the parameters of real particles. This avoids a situation where parts of different capsules are matched. Figure a) demonstrates the result of the detection of active particles, zoomed in from one exemplary video frame.

After detecting particles on each video frame, the Crocker-Grier algorithm is applied to find their trajectories. While tracking works well for a single particle as well as multiple non-overlapping particles, it is problematic closer to the orifice, where capsules often jump, assume tilted positions or even partially overlap. This is true especially for higher excitation amplitudes and smaller orifice sizes.

In the situations where particle tracking is problematic, one can adopt another method which allows to extract the times when a particle leaves the silo (passage times). A similar technique was previously used for the extraction of data in experiments with sheep and humans and soft grains . Here, from each video frame a thin image slice from a line across the outlet is sampled. If the frame rate is sufficiently high, one obtains a long space-time image with blue and white regions corresponding to each capsule which passes the orifice by stacking these slices. It suffices to focus on the more easily detectable white capsule parts, while they are oriented back from the orifice for the most particles leaving the silo. Moreover, it is mentioned in Ref.  that the exact part of the object that is chosen as a reference to measure the passage time does not substantially affect the statistics. Figure b shows a short cutout from this image which corresponds to 400 frames (4 seconds). By finding the rightmost (latest) borders of the segmented white regions, marked by red lines in Fig. b), one gets the frame numbers and hence exact times when capsules pass the outlet. From these data we obtain the number of particles remaining in the silo at each instant of time during the discharge, as well as the times between subsequent particles leaving it. Additionally, one can extract the orientations of the capsules.

The identification of passive capsules can be done analogously. As mentioned before, these capsules are completely red. In order to detect red particles, the red channel is extracted with MATLAB and then subtracted from the specified background (same background subtraction for all images). Finally, we analyze the properties of the detected objects such as the centers of mass, positions, orientations and areas for each frame. As the passive particles show practically no overlaps and do not leave the hopper plane, this simple particle detection is sufficient to access their egress and orientation behaviour.

# Results

## Single particle kinematics

For the determination of the single particle motion, trajectories of isolated grains were recorded and their mean square displacements (MSD) and mean square rotational displacements were measured. Figure shows the evolution of the distribution of displacements with progressing delay time for isolated active particles on a horizontal vibrating plate . It is clearly evident that there is a drift along the capsule axis, and a slight broadening of the distribution by additional random, diffusive contributions. In contrast, the distribution of displacements perpendicular to the capsule axis remains symmetric with respect to the two sides. It is a pure diffusion process without drift. The mean value of this distribution remains close to zero.

Figure shows typical mean square displacements for the active particles. The displacement along the capsule axis is nearly two orders of magnitude larger than that in perpendicular direction. Moreover, the MSD grows with the square of $\Delta t$. Here, we have chosen a short period of 2 seconds where the change of the direction of motion by rotational diffusion is still small. The displacements refer to the initial capsule orientation, which can change by a few degrees within the evaluation period $\Delta t$. The diffusion perpendicular to the capsule axis is mainly caused by such fluctuations of the long axis direction, therefore it is also proportional to $\Delta t^2$.

Particle trajectories keep their current directions for a few seconds, over some persistence length of the order of the capsule length (see ). On long time scales ($\Delta t\gg 2$ s), purely diffusive behavior takes over, which is well known for active grains .

For the passive particles, the distributions of both the displacements along and perpendicular to the particle axis are symmetric with respect to positive and negative displacements. The passive particles diffuse much more slowly than the active ones, as seen in Fig. . Along their axis direction, they diffuse more than one order of magnitude faster than perpendicular to it. The MSD curves parallel to the long axis grows linearly with the delay time $\Delta t$. The perpendicular diffusion has a stronger dependence on $\Delta t$, since diffusion parallel to the axis is 'mixed in' due to rotational diffusion (see paragraph on active particles). Typical distributions of displacements in dependence of $\Delta t$ are shown in the Appendix A in Fig. . This appendix also shows more results for active capsules at different excitation strengths. The essential message from the single-particle measurements is that we have qualitatively distinct behaviour of both types of capsules that are quantitatively reproducible. On the basis of these results, we investigate now how the differences in the individual particle motions influence the collective behavior in the bottleneck passage experiment.

We conclude from these measurements that the asymmetrically loaded capsules are useful models for active matter, and a comparison of their individual and collective behavior with that of the symmetrically loaded, passive particles will yield useful insight on the role of activity in several experimental situations like passing bottlenecks, self-organization of collective motion (swarms), sedimentation, and others.

## Collective passage through narrow bottlenecks

First, we checked how the fill level of our container affects the discharge characteristics. Figure shows two examples of the outflow characteristics of active grains. We selected the $d=18$ mm orifice data at 1.6 $g$ excitation strength and varied the tilt angle. It is seen that individual discharges vary statistically around an average value (red curves) for each tilt angle. The discharge curves are on average rather linear as long as enough grains remain in the container. This means that the mean outflow rate is roughly constant until about 2/3 of the initial filling has left the container. The last $\approx 15$ particles discharge slightly slower than the first ones. In particular, the outflow rate drops significantly when less than 10 capsules remain in the container. Some of these remaining grains are stuck head-on at the side walls in positions that they keep for a long time when there are no other capsules colliding with them.

An interesting aspect is the dependence of the discharge rate on the orifice size. There are different predictions for pedestrians or animals passing a gate, where the rate $\dot N$ of individuals exiting is more or less linearly related to the gate width $d$, $\dot N\propto d^\beta$ with $\beta =1$. In contrast, the outflow rate of granular material, in reasonable approximation described by Beverloo's law, goes with the exponent $\beta=3/2$ for 2D systems. This can be derived from a simple geometric arguments: the product of the orifice width and the velocity of particles passing the exit. When one assumes that the granular particles in a silo fall from an arch with approximate height $d$ above the orifice, the exit velocity is $\propto \sqrt{d}$, which yields $\beta=3/2$. In our system of passive particles (Fig.), we find rates of 2.5/s, 2.9/s and 3.4/s for the three orifice widths 15 mm, 18 mm and 21 mm, respectively, from the linear slopes of the averaged graphs when the container is well-filled. This yields an exponent $\beta\approx 1$. This is a quite reasonable result: The width of the orifice determines the number of particles passing per time, and they all pass with roughly the same velocity since they are not in free fall. They are subject to a continuous drift down the slope. The surprising result is found for the active particles. Here, the exponent is $\beta= 2.2\pm 0.2$ for the three mean discharge curves in Fig. . This behavior is qualitatively different from that of passive capsules.

A key feature of passages of living matter through bottlenecks is their delay time between individuals passing the outlet . We extract this characteristics directly from our experiments. Each graph shown below is constructed from several independent discharge runs, the active capsule data, for example, from 50 individual runs each for each parameter combination. The delay times $\tau_i$ between two objects subsequently passing the exit are distributed according to a density $p(\tau)$ with $\int_0^\infty   p(\tau) d\tau=1$. One can define a mean delay time $\bar\tau$ under the condition that the integral $\bar\tau=\int_0^\infty \tau p(\tau) d\tau$ is finite, otherwise one needs to resort to other measures describing the overall discharge characteristics. The probability of delays being longer than a time $\tau$ is given by the complementary cumulative distribution function (CDF) $$G(\tau) = \int_\tau^\infty p(\tau') d\tau'    .
    \label{eq:taudef}$$ with $G(0)=1$ and $G(\infty)=0$. This cumulative distribution of experimental data is much easier to handle than $p(\tau)$. We have no theoretical prediction of the actual form of $G(\tau)$ but it turns out that the empirical ansatz $$G(\tau) = \frac{1}{1+(\tau/\tau_0)^\gamma}
    \label{eq:tau2}$$ with two adjustable parameters $\tau_0$ and $\gamma$ is well suited to describe all experimental data with desired accuracy. The asymptotic tail of this distribution yields a power law dependence $G \propto \tau^{-\gamma}$ which can be directly compared to data from other studies . We note that the form of $p(\tau)$ is then not as simple, $$p(\tau) = \frac{\gamma}{\tau_0}\frac{(\tau/\tau_0)^{\gamma-1}}{[1+(\tau/\tau_0)^\gamma]^2}
    \label{eq:tau3}$$ An advantage of the form of the cumulative function is that the integral $$\int_0^\infty G(\tau)d\tau = \int_0^\infty \tau p(\tau) d\tau = \bar\tau
    \label{eq:tau_mean}$$ exists for $\gamma > 1$ and yields the analytical expression $$\bar\tau  = \frac{\pi}{\gamma \sin \pi/\gamma} \,\tau_0$$ for the mean delay time $\bar\tau$ as a function of the fit parameters $\tau_0$ and $\gamma$. The very good agreement between the experimental $G(\tau)$ values and the empirical fit function Eq. () is demonstrated in Appendix B.

The passage experiments were performed for three different opening sizes, three different excitation strengths and several tilt angles. In order to restrict the data in this paper to the essentials, we present as examples an excitation strength dependence at fixed tilt and orifice (Fig. ), an orifice size dependence with fixed tilt and excitation (Fig. ), and a tilt angle dependence at fixed orifice and excitation parameters (Fig. ).

The influence of different excitation strengths on the distribution of delay times between passages of individual particles is shown in Fig. a for active and passive capsules. It is seen that the influence of a higher or lower excitation strength on the global shapes of the curves is limited to a shift to longer or shorter times, respectively. With larger excitation amplitudes, the curves are shifted to shorter delay times $\tau_0$ but the slopes remain the same within experimental uncertainty. This means in particular that the exponents $\gamma$ of the fit curves, reflected in the slopes in the log-log plot for large $\tau$, are practically independent of the excitation strength. These $\gamma$ are shown for different orifice sizes and excitation strengths in Fig. b. The fit parameters $\tau_0$ are systematically smaller for the passive particles than for the active ones at given $d$ and $\Gamma$ and $\theta$ for all experiments. It is also seen that the exponents for passive and active capsules are substantially different. The average coefficient for active grains is $\gamma_{\rm a} =1.57 \pm 0.1$, whereas the average for passive grains is $\gamma_{\rm p} =2.58 \pm 0.05$. Both coefficients are large enough to guarantee the convergence of the integral in Eq. (). The corresponding mean delay times $\bar\tau$ are $\approx 0.58~\tau_0$ for our active and $\approx 0.36~\tau_0$ for our passive particles.

Figure  shows the tilt angle dependence of the delay distributions $G(\tau)$. The orifice width and excitation strength were kept constant. It is obvious that the slopes at large $\tau$ are only slightly influenced by the variation of the effective gravitation ('motivation'). A change of the tilt angle $\theta$ from $2.5^\circ$ to $4^\circ$ reduces $\gamma$ by only 10% for active and 7% for passive capsules. On the other hand, the larger tilt angle reduces $\tau_0$ by 65% for the active grains and 33% for the passive ones.

# Discussion and summary

The investigation of symmetrically and asymmetrically loaded spherocylindrical capsules on a vibrating plate can be used to explore the characteristics of self-propelled objects with qualitatively different features. We have found that capsules equipped with an asymmetric weight perform a directed motion along their long axis with velocities between 5 mm/s and 10 mm/s. The directed motion is in the direction of the heavier part. Velocities of individual capsules vary within a narrow range. For excitation strengths between 1.4 $g$ and 1.8 $g$ acceleration amplitude (70 Hz frequency), an increase of the excitation strength causes a faster motion. Above 1.8 $g$, the particles start to lift their lighter halves from the vibrating plate, and the propagation velocity drops. At the same time, the capsules undergo fluctuations of their orientations on the scale of a few seconds, so that the persistence length of their individual trajectories is of the order of a few particle lengths. This behavior reflects some of the features of bacteria, worms or ants in the search for food. When we break the in-plane symmetry of the vibrating plate by a slight tilt of $2^\circ$ to about $5^\circ$, the active motion of the particles is superimposed by a drift down the slope. There is a small tendency of individual active capsules to orient with the heavier part uphill so that the active motion counteracts the downward drift. When the active particles are enclosed in the container at large packing fractions, and they have contact to their neighbors, this tendency is less pronounced (see video in ). Thus, we assume that this effect is not the primary reason for the observed unique discharge characteristics of the container.

The symmetrically loaded, passive particles perform no directed motion, they show a purely diffusive behaviour, where the diffusion along the capsule axis is dominating and diffusion perpendicular to it can be practically neglected. They reflect the typical behavior of inanimate matter. The comparison of both types of objects may help to identify some statistical consequences of biological and granular particle activity.

Besides the motion of individual particles, among the interesting problems is the interaction of multiple particles, alignment and local order as a function of the packing density, the sedimentation of active particles, and the passage around obstacles or through bottlenecks. Our study was focused on the latter problem. The first surprising observation is the discharge rate vs. orifice width. Passive capsules behave just as intuitively predicted, with the discharge rate roughly proportional to the orifice width. This means that the capsules leave the orifice at some average velocity that does not depend on the orifice width. Such a behaviour is comparable to many egress situations with living objects. In contrast, the outflow rate of the active capsules follows an orifice width dependence with at least square power. The exponent determined from the experiments is, independent of excitation strength and tilt angle, $\beta=2.2$ within 10% experimental uncertainty. The reason for this is not obvious and may require further experimental studies and, primarily, appropriate numerical simulations.

The second difference between the discharge characteristics of our active and passive particles is the distribution of delay times $\tau$ between subsequent passages of particles. The cumulative functions $G(\tau)$ of passive capsules are characterized by an exponent $\gamma_{\rm p} = 2.58\pm 0.1$ of the long-time tail, whereas the active capsules exhibit a much weaker exponent of $\gamma_{\rm a} = 1.57\pm 0.15$. This means that their probability of long delays is considerably higher than for the passive particles. Compared to animals passing a gate or pedestrians exiting through a door , the exponents $\gamma$ in our experiments are considerably smaller, i. e. the probability of long-term congestions is much larger for our active capsules. One of the reasons is that the activity is not directed towards the outlet like in usual egress experiments with living matter, but the individual particle velocities are determined by the 'erratic' orientations of the individuals. This behavior can at best be compared to random motion of lower organisms or a crowd in disordered panic.

In vibrated hoppers, Janda et al. determined time lapses between subsequent grain passages and found tails with exponents near $\tau^{-2}$ in the $p(\tau)$ graphs (corresponding to asymptotic tails $\propto \tau^{-1}$ in the cumulative distributions). In some of these experiments, the integral to calculate the mean time lapse diverges and one needs an alternative characterization of the dynamic state of the system. An experiment with 'bristlebots', self-propelled elongated mechanical toys , yielded exponents that are also much smaller than those found in our study, close to the critical value of $\gamma_{\rm c} =1$ where $\bar \tau$ diverges.

Finally, suspensions of particles passing bottlenecks show a similar intermittent behavior, but with tails consistently flatter than $\tau^{-2}$ in $p(\tau)$  ($\gamma <1$). With respect to the intermittent character of the flow, the present system thus represents an intermediate dynamics between those of living objects and those of suspensions in bottleneck geometries.

The results obtained for all these experimental parameters are surprisingly robust. Changes of the tilt angle (the 'incentive') and the excitation parameters do essentially affect the time scales of the egress, but leave the scaling exponents $\gamma$ and $\beta$ practically unaffected. This suggests that these dynamical parameters are determined essentially by the particle dynamics, not by the geometrical or dynamical details. The opening size merely determines the outflow rate but not the shapes of the characteristic curves. This makes us confident that many of the above results can be generalized.

# Acknowledgements

We cordially thank C. Wagner for leaving the shaker device to us. The authors acknowledge DLR for funding within project EVA (50WM2048). M. M. acknowledges a research grant from the Iranian Ministry of Science, Research and Technology for funding part of her stay in Magdeburg.

# Appendix A: Distributions of displacements for active and passive particles

Figure shows the distribution of translational and rotational displacements of the passive particles, measured from 11 individual trajectories. It is obvious that the distribution functions broaden continuously, leading to the mean values shown in Fig. above, but the spatial displacements are statistically symmetric with respect to positive and negative steps. Note that the diffusion along the particle axis is much faster than in perpendicular direction, by a factor of approximately five. The rotational distribution also broadens continuously, but it has in general a small bias towards one of the two rotation senses. This small asymmetry is not considered relevant here, each individual capsule has a slight bias in one or the other direction.

Figure shows the distribution of translational displacements of the active particles along their axis, measured from 11 individual trajectories each. The distribution functions broaden continuously with larger delay $\Delta t$, but most prominent is the shift of the funcion in the direction of the front end of the particle. This shift becomes faster with stronger excitation. The displacement perpendicular to the active axis are comparable to those of the passive capsules.

# Appendix B: Fit of the cumulative distribution $G(\tau)$

Figure demonstrates the quality of the fit function, Eq. (), for the determination of characteristic parameters describing the experiment. The experimental data are well reproduced by the empirical fit. There is a slight deviation for very long $\tau$ which arises primarily from a few particles dwelling in the container when it is nearly emptied. One can get an impression of the distribution function $p(\tau)$ using the fit parameters and Eq. (). It agrees with the experimental data satisfactorily, but naturally the experimental results for $p(\tau)$ scatter much more than those of its integral $G(\tau)$.


abstract: # Abstract
  .
  We fabricate artificial molecules composed of heavy atom lead on a van der Waals crystal. Pb atoms templated on a honeycomb charge-order superstructure of IrTe$_2$ form clusters ranging from dimers to heptamers including benzene-shaped ring hexamers. Tunneling spectroscopy and electronic structure calculations reveal the formation of unusual relativistic molecular orbitals within the clusters. The spin-orbit coupling is essential both in forming such Dirac electronic states and stabilizing the artificial molecules by reducing the adatom-substrate interaction. Lead atoms are found to be ideally suited for a maximized relativistic effect. This work initiates the use of novel two dimensional orderings to guide the fabrication of artificial molecules of unprecedented properties.
author: Jae Whan Park$^{1,\dagger}$; Hyo Sung Kim$^{1,2,\dagger}$; Thomas Brumme$^3$; Thomas Heine$^{3,4,5}$; Han Woong Yeom$^{1,2}$
title: Artificial Relativistic Molecules

# Introduction

Artificial atomic clusters or lattices supported on solid surfaces can exhibit tailored electronic, magnetic, and topological properties of fundamental and technological importance. For example, atomic-scale chains and clusters disclosed the quantum confinement of electrons , topological edge modes , superlattice Dirac bands , flat bands , atomic scale spin interactions , and topological defects . The direct atom-by-atom manipulation and the self-assembly of atoms or molecules are two major approaches to realize such atomic scale chains, clusters, and finite lattices. Both methods, however, have their own limitations under given interatomic interactions and fabricating energetically and kinetically unfavorable structures or assemblies has been a huge challenge .

In case of the self assembly of supported clusters, such limitations may be overcome by templates, which provide unusual growth environments to produce otherwise unfavored molecular structures. Step edges and one-dimensional surface superstructures were actively used as such templates on metal and Si surfaces . The growth of metal atoms on these substrates leads to unusually anisotropic clusters. On the other hand, nanosized two dimensional superstructures were suggested as novel templates, such as the moiré structure of hexagonal boron nitride , metal-organic frameworks of DNA , and metal-organic honeycomb networks on metal surfaces . However, while the confinement of metal atoms and molecules on these templates were demonstrated, the formation of artificial molecules was not reported yet.

In this work we employ a novel template-guided atomic self-assembly technique and take advantage of the strong relativistic effects in lead atoms that effectively create long-ranged interatomic bonds. This allows us to produce a series of unprecedented artificial molecules with relativistic, so called Dirac, electronic orbitals. The honeycomb charge-order superstructure of the van der Waals (vdW) crystal IrTe$_2$ cages various Pb clusters ranging from dimers to heptamers including those in interesting benzene-like hexagonal rings. Atomically-resolved spectroscopy and electronic structure calculations reveal clearly the formation of molecular orbitals by direct orbital overlap between neighboring Pb adatoms. The spin-orbit coupling (SOC) of Pb atoms and of the substrate is shown to enhance the interatomic interaction to drive the formation of unprecedentedly relativistic molecules at an unusually large interatomic distance. This case exemplifies the chemical stabilization of nanoscale assemblies uniquely by relativistic effects and opens the potential of novel superstructures on vdW crystals for the fabrication of otherwise unfavored molecules with unprecedented properties.

# Results

## Charge order superstructure hosts Pb clusters

Figure 1a shows Pb adatoms deposited on the honeycomb charge-order superstructure at the surface of a IrTe$_2$ vdW crystal. Below 180 K, a honeycomb superstructure forms with a period of 2.2 nm by the spontaneous ordering of 5$d^{3+}$ and 5$d^{4+}$ valence electrons of Ir (with the buckling of the Te layer accompanied) and competes with various stripe charge orders . Pb adatoms sit at hollow sites surrounded by three Te atoms (Fig. 1b and Supplementary Fig. 2) and, very selectively with yet unclear reasons, inside of each honeycomb, which are above Ir atoms with the 5$d^{4+}$ valency (see Supplementary Fig. 1). Density-functional theory (DFT) calculations including SOC consistently predict the four-fold hollow site adsorption with an adsorption energy of -1.93 eV against the three-fold hollow site of -1.74 eV and Te on-top site of -1.16 eV. Clusters ranging from one to seven Pb atoms are confined within a honeycomb unitcell with various lateral configurations. The Pb-Pb interactomic distance is dictated by the substrate with two choices of $\sqrt{3}$a$_0$ or 2a$_0$ (a$_0$, the IrTe$_2$ lattice constant of 3.5 Å ). The 2a$_0$-spaced dimers (Fig. 1c) are more frequently observed than those of $\sqrt{3}$a$_0$ Pb-Pb distance (Fig. 1d). As trimers, a heterogeneous trimer with a Pb-Pb spacing both of $\sqrt{3}$a$_0$ and 2a$_0$ or a compact one of 2a$_0$ or $\sqrt{3}$a$_0$ are formed (Fig. 1e and Supplementary Fig. 4c). Figures 1f and 1g show a tetramer formed by two 2a$_0$ dimers and one with $\sqrt{3}$a$_0$ and 2a$_0$ dimers, respectively. We also find benzene-like hexagonal ring clusters as shown in Figs. 1h, 1i, and 1j among a few other interesting clusters. These configurations are made only by 2a$_0$ building blocks and are the largest clusters confined within the hexagonal unit of the template.

## Relativistic molecular orbitals in Pb dimers

Electronic energy levels of Pb clusters are revealed by scanning tunneling spectroscopy (STS). For a monomer (Fig. 2c), we observe three main spectral features at 1.40, 1.71, and 1.93 eV for Pb 6$p$ valence electrons. Our calculation shows that a Pb adatom is partly ionized by donating 6$p$ electrons into the substrate to shift the 6$p$ states to unoccupied states (Fig. 2o and Fig. 2w). As shown in Fig. 2w, the relativistic effect of SOC splits the 6$p$ states into 6$p_{1/2}$ and 6$p_{3/2}$ Dirac orbitals and the latter splits further due to the substrate-induced splitting of in-plane ($p_{xy}$) and out-of-plane ($p_z$) orbitals. Consequently, three main peaks correspond to the relativistic $p$ orbitals of $p_{1/2}$ ($j$=1/2, $m_j$ = $\pm$`<!-- -->`{=html}1/2) and $p_{3/2}$ ($j$=3/2, $m_j$ = $\pm$`<!-- -->`{=html}3/2, $\pm$`<!-- -->`{=html}1/2), respectively.

In case of a dimer with a shorter Pb-Pb distance, an energy splitting (0.23 eV) of molecular bonding and antibonding states is observed clearly for the 6$p_{1/2}$ state (Fig. 2j). In a longer dimer, the splitting is much smaller, about 0.10 eV, as the orbital overlap reduces (Fig. 2f). Our DFT calculation predicts (Figs. 2r and 2V) the lowest molecular orbital splittings of 0.05 and 0.18 eV for the longer and shorter dimers, respectively, in reasonable agreement with the experiment. The present DFT calculation seems to underestimate the molecular bonding interaction by about 0.05 eV. These energy scale is consistent with what are expected for a freestanding Pb-Pb dimer (Supplementary Fig. 6). This result clearly indicates the molecular orbital formation within the Pb clusters at an unusually large Pb-Pb distance of 7 Å, which is observed consistently for all cluster configurations identified (Supplementary Fig. 7).

The molecular levels of a Pb dimer are well explained by the simple tight-binding interaction of Dirac atomic orbitals perturbed weakly by the substrate as discussed above. That is, three relativistic $p$ orbitals split simply into the combination of molecular orbitals in gerade ($g$) and ungerade ($u$) symmetries (see Fig. 2w). For example, the two lowest peaks correspond to the combination of the $p_{1/2}$ orbitals in $g$ and $u$ symmetries. The resulting Dirac molecular orbitals are clearly differentiated from their scalar-relativistic counterparts of $\sigma$ and $\pi$ molecular orbitals. The $g$ ($u$) orbital consists of one-third $\sigma$-bond ($\sigma$-antibond) and two-third $\pi$-antibond ($\pi$-bond) character, indicating a substantial mixing between $\sigma$- and $\pi$-bonds due to the full relativity effect of SOC. The Dirac basis and the interaction with substrate atoms also breaks the mirror symmetry perpendicular to the dimer axis. This is reasonably confirmed in the experiment (Figs. 2h and 2i) and in the calculation (Figs. 2t and 2u). Further details of the interaction with the substrate will be discussed below, which become important for the understanding of the relativistic effect on the adsorption energy.

The direct and real space observation of relativistic molecular orbitals is exceptional, since most of previous experiments for very few heavy element molecules relied on spectroscopy techniques and, thus, the relativistic effects were mainly discussed in theoretical aspects . Although a Pb dimer has a much longer bond length (7 Å) than other relativistic molecules and there is a substantial adatom-substrate interaction, the molecular splitting mainly arises from the direct overlap of relativistic $p$ orbitals. This is distinguished from the substrate-mediated long-range interaction of other metal clusters on metal surfaces .

## Mechanism of the interaction between Pb atoms

It is important to point out that the attractive interaction between Pb adatoms is only due to SOC and that the 2a$_0$ distance imposed by the substrate is optimal for the formation of a Pb dimer (Fig. 2x) in our calculations. This is consistent with the observation of the higher population of 2a$_0$ dimers. The SOC energetics for the larger Pb molecules, from trimer to heptamer (Fig. 2y and Supplementary Fig. 10), is also consistent with the experimental observation of larger Pb molecules with exclusively 2a$_0$ nearest-neighbor distance. The 2a$_0$ distance is optimized by two competing energy contributions, the ionic (or dipole) repulsion and the orbital overlap energy gain of adatom-substrate hybridized states. The comparison of the calculated differential charges with and without SOC reveals two sophisticated SOC effects for those contributions; (i) a reduced donation of $p$ electrons to the substrate to reduce the ionic repulsion and (ii) an enhanced orbital overlap of the Pb-Te hybridized states (see Supplementary Fig. 10). The former is due to the lowering of the atomic energy level of Dirac orbitals by SOC. It is also partly due to the charge redistribution of the IrTe$_2$ substrate by SOC, namely the charge transfer from the $d_{z^2}$ states of Ir atom to Te atoms. The substantial interaction of Pb with the substrate Te orbitals provides an extra energy gain for the antibonding level of the $p_{1/2}$ ($m_j$ = $\pm$`<!-- -->`{=html}3/2) orbital of the $\sqrt{3}$a$_0$ dimer (see Supplementary FIG. 9b). That is, the interwined effect of the substrate and the relativity determines the stability of the Pb molecules. In these respects, Pb is ideal since it has very strong relativistic effects (in comparison, for example, to Sn, see Supplementary Fig. 11) and a proper interaction strength with the substrate (in comparison to Tl, see Supplementary Fig. 12). Of course, the formation of Pb molecules would also be affected by the size and shape of the charge-order honeycomb template.

## Benzene-like ring-shaped Dirac molecules

While a consistent interaction mechanism holds for all 15 different Pb molecules identified (Supplementary Fig. 4), particularly interesting ones are those based on a six-fold hexagonal ring (Fig. 3). The benzene-like ring molecule of Pb$_6$ with a 2a$_0$ Pb-Pb distance exhibit the weak molecular splitting similar to a 2a$_0$ dimer (Fig. 3e). The intriguing molecular orbital formation is clearly imaged by the $dI/dV$ maps: (i) the broken AB sublattice symmetry for the lowest energy state (Fig. 3b), (ii) a Kekulé-like distortion for the second state (Fig. 3c), and (iii) a mirror-symmetry-broken two-fold-symmetric bonding feature for the third state (Fig. 3d). Our calculation reproduces well the spatial characteristics of these molecular orbitals, which are indeed due to the relativistic character and the interaction with the substrate (Fig. 3l-n and Supplementary Fig.13). Note, however, the remaining energetic discrepancy between experiment and calculation. The present calculations underestimate the molecular bonding interaction and the calculated spectral features are almost degenerate at around 1.3 eV for the benzene-like molecule. We think that this difference is at least partly due to the limitation of our model in taking into account of the substrate electronic states in its charge-ordered correlated state . Another interesting molecule is the Pb pentamer (Figs. 3f-3j), whose Dirac molecular orbitals are also reasonably well simulated (Figs. 3p-3t). Such type of an individual open ring molecule cannot be synthesized. The Pb$_5$ molecule features the edge states at both truncated ends between the bonding and antibonding levels of three Dirac $p$ orbitals (see Supplementary Fig. 14). It indicates the ring-type Pb molecules are circularly delocalized electron systems. The formation of a filled-benzene-ring molecule is also highly unusal (Fig. 1j), where the interatomic overlap is much enhanced by the central Pb atom with unique bonding configurations (Supplementary Figs.15-17). This molecule can be compared with a rare example of a planar B$_7$ molecule except for the strong SOC. While not accessible by the present experiments, the benzene-type Dirac molecular orbitals fabricated here have unique and interesting spin configurations (see Supplementary Fig. 13), which may be exploited further. It becomes very obvious that the present artificial Pb molecules introduce unprecedented molecular configurations combined with the strong SOC.

## Summary

We fabricate various artificial Pb molecules utilizing a novel honeycomb template of the charge-ordered superstructure in IrTe$_2$. The molecular orbital formation is clearly observed in the Pb clusters from dimers to unusual benzene-ring-type molecules. The spin-orbit coupling is found to play crucial roles in both the condensation of Pb atoms and the molecular orbital formation, bringing them into an unprecedented regime of relativistic Dirac molecules. Beyond the relativistic chemistry, the unusual bonding and spin structure of these relativistic molecular orbitals may be combined with the novel electronic properties of the substrate such as the charge ordering and the emerging superconductivity to lead to a new type of a quantum system. In addition to the versatile technique of atom-by-atom manipulation, where the interatomic distance may be varied, the use of various types of 2D superstructures as templates would definitely extend the potential of molecular or cluster self assemblies. Interesting 2D superstructures can include twisted bilayer graphene and domain-wall or twin-boundary networks of transition metal dichalcogenides . The templated self assembly of magnetic atoms on such substrates can be very interesting due to the proximity to the novel 2D electronic states of substrates.

# Methods

## Sample preparation

Single crystals of IrTe$_2$ were grown by Te flux using pre-sintered IrTe$_2$ polycrystals as reported previously. Samples were cleaved in a vacuum better than 5$\times$`<!-- -->`{=html}10$^{-10}$ torr at room temperature. Pb atoms are grown on cleaved IrTe$_2$ surfaces at room temperature. The Pb atom tends to cover the IrTe$_2$ surface up to  0.1 ML. Over this coverage, this Pb growth mode depends on the deposition temperature. For example, when the Pb grows at room temperature, Pb forms a uniform film. However, at low temperature deposition, depending on the coverage and annealing condition, we can found different heights of the Pb islands.

## STM measurement

All the STM measurements were performed with a commercial ultrahigh vacuum cryogenic STM (Specs, Germany) in the constant-current mode with PtIr tips at 4.3 K. The differential conductance, $dI/dV$, was measured using the lock-in detection with a modulation of 1.17 kHz.

## DFT calculations

We perform the relativistic DFT calculations using the Vienna $ab$ $initio$ simulation package within the generalized gradient approximation of the Perdew-Burke-Ernerhof type and the projector augmented-wave method . The IrTe$_2$ surface is modeled by a (5$\times$`<!-- -->`{=html}5) and (7$\times$`<!-- -->`{=html}7) supercell with a single IrTe$_2$ layer and a vacuum spacing of about 23.6 Å. The calculated value 3.839 Å$$ is used as the lattice constant of IrTe$_2$. We use a plane-wave basis with a cutoff of 211 eV for expansion of electronic wave functions and a 15$\times$`<!-- -->`{=html}15$\times$`<!-- -->`{=html}1 $k$-point mesh for the 1$\times$`<!-- -->`{=html}1 Brillouin-zone integrations. To avoid the unexpected substrate distortion of the IrTe$_2$, only adsorbed Pb atoms are relaxed until the residual force components are within 0.02 eV/Å. Here the Fermi level was shifted down by 0.70 eV because to overcome the inherent discrepancy between experiment and theory due to unknown charge-ordered structure, doping effect of the Pb adsorbates and the IrTe$_2$ sample and/or the inaccuracy of DFT band-energy calculations.

## Data availability

The authors declare that the data supporting the findings of this study are available within the article and its Supplementary Information.

## Acknowledgements

J.W.P., H.S.K., and H.W.Y. were supported by Institute for Basic Science (Grant No. IBS-R014-D1). The IrTe$_2$ crystal was provided by S. W. Cheong in Rutgers University. J. Lee is appreciated for his help in data processing. T.B. and T.H. used the ZIH Dresden computational resources.

## Author contributions

J.W.P. performed the DFT calculations. H.S.K. carried out the STM measurements. These two authors equally contributed in this paper. J.W.P. and H.W.Y. analysed the data and wrote the manuscript with the comments of all other authors. H.W.Y. supervised the research.

## Competing interests:

The authors declare no competing interests.


abstract: Attenuation of a surface acoustic wave is used as a highly sensitive and non-invasive probe of persistent photoconductivity effects in ZnCdSe/ZnSe quantum wells. These effects are observed over long time-scales exceeding several minutes at low temperatures. By varying the optical excitation energy and power and temperature we show that these effects arise from carriers photogenerated by interband excitation which are trapped in random potential fluctuations in the quantum wells related to compositional fluctuations. Effects related to defect levels in the bandgap can be excluded and a transition of the conduction mechanism with temperature from a hopping to a percolation regime is observed. The transition temperature observed for our quantum well material is strongly reduced compared to bulk crystals. This indicates a superior structural quality giving rise to only weak potential fluctuation of $\lesssim 3$ meV.
author: D. A. Fuhrmann; H. J. Krenner; A. Wixforth; A. Curran; K. A. Prior; R. J. Warburton; J. Ebbecke
bibliography: references.bib
date: 2024-10-01
title: Non-invasive probing of persistent conductivity in high quality ZnCdSe / ZnSe quantum wells using surface acoustic waves

## Introduction

Surface acoustic waves (SAW) are a highly versatile tool to probe and manipulate optical and electronic properties of materials. Examples in the field of low-dimensional semiconductor heterostructures include charge conveyance and spectroscopy of high-mobility electron systems and optically active systems as quantum wells (QWs) and quantum dots (QDs) . However, most of these investigations focused on the most advanced III-V semiconductor compounds due to their piezoelectricity which allows for the use of interdigital transducers (IDTs) to excite and detect SAW. Despite the fact that most II-VI semiconductors are also weakly piezoelectric still SAW-based studies are very limited .
One important mechanism which has large impact on the electrical transport properties in semiconductors is the persistent photoconductivity effect (PPC). This persistence of free charge carriers in semiconductors long after the removal of the exciting light source can by far exceed typical carrier lifetimes and leads to a pronounced change of conductivity in the material. PPC has been studied in the AlGaAs and other III-V material systems in great detail for many years . In the case of AlGaAs a thermal barrier prevents the recombination of the DX-centers at low temperatures which gives rise to the PPC effect. Similar PPC effects have also been observed in II-VI semiconductors such as ZnCdSe, ZnMgSSe, ZnMgSe and CdTe . In inhomogeneous semiconductors and heterostructures its origin can be of a different nature. In particular for bulk ZnCdSe spatial separation of charge carriers in random local potential fluctuations (RLPF) can occur. These RLPF result from compositional fluctuations during material growth and have been found to be the dominant mechanism giving rise to PPC . One striking consequence is the transition of the conduction mechanism from a localized to a percolation regime at a critical temperature. Surprisingly, a similar transition temperatures were found for both bulk and QWs grown by molecular beam epitaxy indicating that in these materials RLPFs are introduced inherently and cannot be reduced by introducing heterointerfaces.
In this paper we present a study of PPC in a ZnCdSe/ZnSe multiple quantum well (MQW) structure without directly measuring the conductivity by electrical means. Instead we apply a method based on surface acoustic SAWs, which for low SAW power represents a non-invasive and highly sensitive probe of the conductivity. The underlying mechanism is the attenuation of the SAW by mobile carriers which can be measured with high precision. Moreover, we find that for our QW material the transition temperature between the localized and percolation transport regime is significantly reduced to $T_C\lesssim 40$ K compared to values of $T_C^{bulk} \sim 120$ K reported for bulk material . This finding clearly shows the superior quality of our material compared to previous reports and indicates that inhomogeneities resulting in RLPFs are substantially reduced.

## Experimental setup

Our sample consists of a ZnCdSe/ZnSe heterostructure grown by molecular beam epitaxy (MBE). Five optically active, 4 nm wide Zn$_{0.9}$Cd$_{0.1}$Se QWs, separated by 8 nm ZnSe barriers are embedded between two 63 nm ZnSe cladding layers. After growth, the II-VI semiconductor heterostructure was removed from the substrate and transferred onto the strongly piezoelectric LiNbO$_3$ chip using an epitaxial lift-off method . The transferred film is positioned between two IDTs on the host LiNbO$_3$ substrate. Due to the increased coupling coefficient, $K_{\rm eff}$, of LiNbO$_3$ compared to Zn(Cd)Se this hybrid sample layout efficiently gives access to higher piezoelectric fields in the semiconductor material which e.g. was used to study SAW mediated exciton dissociation . A microscope image of one part of the sample is shown in the inset to Fig. . In the upper part the IDT with an aperture of 0.58 mm allows for the excitation and/or detection of SAW at $f_0= 115$ MHz and $f_3= 3\cdot f_0$. The SAW launched by the IDT propagates through the semiconductor film (length 2.6 mm) to a second IDT on the opposite side. The total length of the SAW delay line is 5.4 mm and the transmitted power was measured using a conventional network analyzer. The sample itself was mounted on the coldfinger of an optical helium-flow cryostat equipped with high-frequency connections to the IDTs. For photoluminescence experiments carriers were photogenerated by a laser diode (Laser-1, $\hbar\omega_{ex}=3.18$ eV, $P_{max}= 30$ mW) focused to $\sim 500~\mu$m diameter spot using a $5\times$ microscope objective. The emission from the sample was collected by the same objective, sent to a 0.5 m grating monochromator where it is spectrally analyzed and detected by a $l{\rm N_2}$ cooled charge-coupled device. For excitation energy dependent experiments Laser-1 was replaced by white light spectrally filtered by a monochromator. For optical manipulation of the PPC a second laser diode (Laser-2, $\hbar\omega_{ex}= 1.84$ eV, $P_{max}= 35$ mW) was used.

## SAW on a ZnCdSe/ZnSe-LiNbO$_3$ hybrid

Surface acoustic waves are attenuated in the presence of free charge carriers close to the sample surface. In a piezoelectric material the periodic mechanical deformation is accompanied by electric fields which induce currents dissipating energy. Whilst for high SAW intensities ($P_{SAW}>$ 20 dBm) a modulation of the conduction and valence band edges leads to a dissociation of optically generated excitons, thus quenching the photoluminescence , in the regime of low powers $P_{SAW}= -20$ dBm, used in all presented measurements, the photoluminescence signal of the QW is not yet altered or quenched. The SAW propagation hence acts only as a weak perturbation of the system and a measurement of the SAW transmission in this power range provides a non-invasive probe of the conductivity. The transmitted intensity between the two IDTs is given by $I =I_0 e^{-s\Gamma}$, with $s$ being the length of the interaction region and $\Gamma$ the absorption coefficient. For a two-dimensional system and small SAW amplitudes, the SAW has only a negligible effect on the charge density $\rho(\textbf{r},t)$, therefore, the for a given sheet conductivity, $\sigma$, and absorption coefficient, $\Gamma$, can be written as $$\Gamma=\dfrac{K^{2}_{\rm eff}\cdot k}{2} \cdot\dfrac{\dfrac{\sigma}{\sigma_{m}}}{1+\left(\dfrac{\sigma}{\sigma_{m}}\right)^{2}}.
    \label{eq:gamma}$$ In this equation $K_{\rm eff}$ denotes the effective coupling constant, a measure of the piezoelectric strength of the material, $k$ the SAW wave vector and $\sigma_m$ the so-called characteristic conductivity. From a closer examination of equation () one can see that the absorption coefficient as a function of conductivity increases linearly for $\sigma<\sigma_{m}$, reaches a maximum at $\sigma=\sigma_m$ and subsequently decreases proportional to $\sigma^{-1}$ with further increasing conductivity. The characteristic conductivity $\sigma_{m}$ is given by $\sigma_{m}=\epsilon_{0}(1+\epsilon_{hyb})v$ and is thus dependent on the speed of sound $v$ in the material along the respective crystal direction and the relative dielectric constant of the hybrid system. $\epsilon_{hyb}$ can be determined for a given layer sequence and depends on the SAW wave vector and, therefore, the different harmonics used in our experiments exhibit different values for the characteristic conductivities $$\begin{aligned}
    \sigma_{m}(\mathrm{115MHz})&\approx &1.04\cdot 10^{-6}\,\Omega ^{-1},\\
    \sigma_{m}(\mathrm{345MHz})&\approx &8.75\cdot 10^{-7}\,\Omega ^{-1}.
    \label{eq:eqn1}
\end{aligned}$$ At these critical conductivities pronounced minima in the transmitted SAW power for both frequencies are expected. The small values for $\sigma_{m}$ indicate the high sensitivity of the SAW to measure small conductivities. In the following the measured parameter will be $s\cdot\Gamma$, the length of the interaction region $s$ multiplied by the absorption coefficient, referred to as the SAW attenuation. While the upper bound of $s$ is determined by the size of the epitaxial lift-off film (2.6 mm) the effective value of $s$ over which changes in the conductivity are detected is mainly limited to the excitation spot size of $\sim 0.5$ mm. Since we keep the spot size constant in our experiments $s\cdot\Gamma$ is a direct measure for the conductivity which are linked by equation .

Figure shows a typical measurement of attenuation at the frequency of the third harmonic ($f_{3}= 345$ MHz) at low temperatures ($T= 35$ K) as a function of time. At the onset of illumination ($t=-10$ s) the attenuation instantly rises. Due to the creation of free charge carriers, the transmission of the SAW is attenuated by $\approx$ 35 dB. At $t=0$ s the laser is turned off, after which the attenuation rises to a maximum, before it very slowly decays on a seconds timescale. This behavior is consistent with equation . In particular, at the maximum of the attenuation the conductivity is equal to the characteristic conductivity $\sigma_m$.

The increase of the attenuation during illumination as a function of temperature is shown in Fig. . For $T>150$ K the transmission of the SAW is not damped and no attenuation is measured. As this is the same temperature above which no photoluminescence could be detected for our sample , it hints towards charge carriers from band-to-band transitions as the origin of the PPC. The attenuation of the third harmonic shows a non-linear dependence on the conductivity for $T<70$ K as described by equation (). At this temperature the conductivity is very close to the characteristic conductivity $\sigma_{m}\approx 8.75\cdot 10^{-7}\Omega ^{-1}$ of the third harmonic. Since the value of the characteristic conductivity of the first harmonic is higher, the attenuation remains directly proportional to the conductivity over the entire range of temperatures $T>10$ K. Unless stated otherwise we restrict our analysis to data obtained at $f_1$ due to the linear dependency of the conductivity on the SAW attenuation. For $T>75$ K the measured values of the third SAW harmonic match those of the first multiplied by three (see inset to Fig. ). This indicates that the slight difference of $\sigma_{m}$(345 MHz) and $\sigma_{m}$(115 MHz) expresses itself in the attenuation for conductivities in the range of $\sigma_{m}$ and shows very nicely that the attenuation is a relative measure of sheet conductivity.

## Persistent photoconductivity in ZnCdSe quantum wells

As shown in Fig. , the SAW propagation is considerably affected during and long after illumination. This is a rather striking discovery, because in an ideal system in which only charge neutral excitons are generated, the dominant loss mechanism would be radiative recombination. This process takes place on much faster timescales which are not affected by the small amplitudes of the SAW used in these experiments. In order to exclude optically active impurities giving rise to levels in the bandgap or surface states, the SAW attenuation was investigated for different photon energies. In this experiment Laser-1 was replaced by white light which was spectrally filtered by our monochromator.

Fig. shows the SAW attenuation of the first harmonic during illumination at $T= 60$ K as a function the photon energy. At this temperature the attenuation is directly proportional to the conductivity, as indicated by the temperature dependence of the first SAW harmonic shown in Fig. and, moreover, the low optical powers used in this experiment ensure that the conductivity is much less $\sigma_m$. Photon energies less or equal than the free exciton recombination energy (a PL spectrum is shown in the figure for comparison), $\hbar\omega\approx$ 2.63 eV, have no effect on the SAW transmission. This excludes activation of defects or impurities with levels within the band gap as the mechanism responsible for the observed effect as it was observed in a previous study on II-VI QWs grown by MBE . For energies larger than the QW emission two regions can be distinguished: (i) For energies ranging between 2.63 eV and $\sim2.78$ eV charge carriers are photogenerated solely in the ZnCdSe-QW. (ii) For energies exceeding the bandgap of the ZnSe barrier material, $E_g(\mathrm{ZnSe})\approx 2.78$ eV , carriers are also generated outside the QW region. The increase of the signal at this absorption edge indicates that some of the carriers excited in the barrier material relax into the QWs levels.

To obtain further insight in the underlying mechanism we performed temperature dependent SAW-attenuation and PL experiments under illumination by Laser-1 which are summarized in Fig. . In the SAW-attenuation experiments the third harmonic was used \[Fig. (a)\] and the laser was switched off at $t = 0$ s. Clearly, the time transients undergo a pronounced change as the temperature is raised. In particular the peak in the SAW-attenuation at which the conductivity becomes comparable to $\sigma_m$ shifts in time and its temporal width changes. At a temperature of $T\sim 30-40$ K both the delay and width of this conductivity peak are maximum. For lower and higher temperatures this peak vanishes and, in addition, a slow and fast decay of the attenuation signal are observed, respectively. Since the observed maximum is linked to a shift of $\sigma$ from values larger to values smaller than $\sigma_m$ its width and contrast reflect the recombination mechanism and rate of the charge carriers in the system. Though $\sigma (t=0$ s) is not the same for all temperatures the temporal width of the attenuation peak is directly proportional to the recombination rate and these measurements suggest a minimum recombination rate at $T\sim 30$ K. Thus a critical temperature $T_C$ is exceeded in the range between for $25<T<35$ K above which the recombination rate increases rapidly with temperature and exhibits an activated temperature dependence. This critical temperature corresponds to a thermal energy of $k_BT_C=2.6$ meV which suggests an activation energy in the same range. Since the excitation energy dependence of the PPC indicates that the underlying mechanism is governed by carriers photogenerated across the bandgap we can use PL spectroscopy to obtain further information over the same temperature range. In Fig. (b) we compare the QW emission at $T=10$, 30 and 60 K which are below, close to and above the temperature at which the transition of the SAW attenuation transients is observed. For clarity the corresponding SAW transients at these temperatures are highlighted in Fig. (a). A close examination of the PL signal at $T=10$ K shows two contributions of similar intensity which are split by $\sim 6$ meV. The features at lower and higher energy arise from emission of bound (bX) and free excitons (fX), respectively. Whilst for $T= 30$ K bX is suppressed compared to fX but still well resolved at $T=60$ K only fX emission is detected. For a detailed analysis of this effect in the temperature range between $T=10$ and 60 K we plot the relative intensities of bX and fX given by $I_{rel,\mathrm{bX/fX}}=\dfrac{I_{\mathrm{bX/fX}}}{I_{\mathrm{bX}}+I_{\mathrm{fX}}}$ as symbols in Fig (c). This data can be described by a thermally activated conversion process from bound (bX) into free excitons (fX): $$\begin{aligned}
I_{rel,\mathrm{bX}}(T)&=& \dfrac{I_\mathrm{bX}}{I_{\mathrm{bX}}+I_{\mathrm{fX}}}\nonumber\\ ~&=&1-\left[\exp\left({\dfrac{-\Delta E}{k_BT}}\right)+\alpha\right]\label{eq:fit1}\\
I_{rel,\mathrm{fX}}(T)&=& \dfrac{I_\mathrm{fX}}{I_{\mathrm{bX}}+I_{\mathrm{fX}}}\nonumber\\~&=&\exp\left({\dfrac{-\Delta E}{k_BT}}\right)+\alpha .
\label{eq:fit2}
\end{aligned}$$ Here, $\Delta E$ is the exciton localization energy and $\alpha$ is a constant which reflects the number of free excitons for $T\rightarrow 0$ K for a fixed generation rate i.e. optical pump power. The corresponding fits of equations and are plotted as solid lines in Fig. (c). From these fits we obtain a localization energy $\Delta E =  4.3$ meV, a value which is comparable to and consistent with the energy splitting between the bX and fX peaks. Moreover, the thermal energy $k_BT_C=2.6$ meV at which the transition is observed in the SAW-attenuation is also in good agreement. Since both SAW-attenuation and PL spectroscopy show a transition in the same temperature range we assume that also for the PPC measured via the SAW-attenuation localization of carriers is the dominant mechanism.

In a semiconductor alloy disorder can lead to compositional variations and the resulting variations of the conduction and valence band edges give rise to RLPF. These RLPF can act as local traps for charge carriers as shown schematically in Fig. . For our ZnCdSe MQW, local chemical fluctuations i.e. of the Cd-content can create potential minima for either electrons or holes, which exhibit an average energy barrier $E_{rec}$ which is a measure for the disorder in the system. Moreover, to transfer carriers between RLPF or release excitons bound to an RLPF the corresponding recombination energy $E_{rec}$ has to be provided for example thermally which gives rise to the observed temperature dependence: For low temperatures charge carriers are localized in the potential minima, thus the resulting conductivity is very low. The separation leads to a reduced wavefunction overlap of electrons and holes inhibiting radiative recombination. This consequently results in an increase of carrier lifetime giving rise to PPC. With increasing thermal energy or temperature, the electron conduction is mediated by hopping transport, since only neighboring minima are accessible. The recombination rate is mainly determined by the distribution of charge carriers in the localization centers, i.e. their wavefunction overlap at the time the illumination is turned off. At a critical temperature $T_C$, a transition of the conduction mechanism occurs : Electrons can now percolate through the network of accessible potential minima, whereas holes remain localized because due to their higher effective mass. The probability for electron trapping is increasing for energetically deeper localization centers where they experience minimal potential energy and the recombination rate reaches a minimum, due to the reduced number of holes in the vicinity. With further increasing temperature, the recombination rate rises until $E_{rec}$ becomes comparable to the thermal energy. From this temperature on PPC breaks down and is no longer observed.
As discussed in the previous paragraph, the critical temperature $T_C$ at which the transition of the transport mechanism occurs depends on the depth of the energy minima. These themselves are directly connected to disorder in the material composition and thus the sample quality. Therefore, a small value of $T_C$ indicates that there is a weak and shallow potential modulation due to RLPF. In previous PPC studies on II-VI semiconductor compounds RLPF have been identified in bulk material to be the dominant mechanism at play. The reported values of $T_C>120$ K were found for bulk crystals and MBE-grown QWs . These studies indicate that the introduction of a heterointerface as in a QWs does not necessarily result in a reduction of RLPF. Comparing the findings of reference to our work we do not observe PPC effects for sub-bandgap excitation which indicates that RLPF are the dominant mechanism. Moreover, from the presented data we can deduce a value of $T_ C<40$ K for our sample. This strong reduction of the critical temperature and the corresponding shallow band edge modulation of $E_{rec}\lesssim 3$ meV underlines the high quality of our films. This value is comparable to the activation energy between bound and free excitons of $\sim 4.3$ meV. The small deviation could arise from different localization energies for individual carriers and bound electron-hole pairs.

## Manipulation and saturation of the PPC

Another way to transfer carriers trapped in RLPF is by optical excitation using light with an energy larger than $E_{rec}$ as shown schematically in Fig. (b). We demonstrated in Fig. that for photon energies below the effective bandgap of the QW no extra carriers are generated and, therefore, no PPC can build up. However, if carriers are already present in the system and trapped in RLPF they can be activated giving rise to a response in the PPC. We want to note that the ability to manipulate the PPC is a clear indication that the underlying mechanism in our II-VI QW-system has a different origin e.g. in AlGaAs where PPC is dominated by DX-centers. To confirm this assumption a second laser (Laser-2) with an energy $\hbar\omega_{ex}= 1.84$ eV *less* than the effective bandgap of the QW was used for excitation. Due to its smaller energy no excitons are generated by this laser via interband absorption. Fig. shows in the main panel the attenuation of the first SAW harmonic over time for a temperature of $T= 45$ K, which ensures a linear dependency on the conductivity. At $t=0$ s the above bandgap excitation $\hbar\omega_{ex}= 3.18$ eV of Laser-1 is turned off (gray shaded area) and the SAW attenuation decays slowly as described before. At $t= 3.5$ s Laser-2 is turned on and a sharp increase of the SAW attenuation is observed. This increase can be understood and explained by an intraband activation of trapped electrons and holes as shown schematically in Fig. (b). After activation these carriers now contribute to the sheet conductivity until they recombine. This gives rise to the observed decay of the conductivity even though Laser-2 remains on. This is furthermore consistent with the fact that no additional electrons and holes are photogenerated and only trapped carriers are activated which are lost after recombination. Clearly, the time constant of the decay is faster with Laser-2 switched on which was also observed using an infrared light emitting diode instead of Laser-2. To further confirm this additional activation effect we performed a steady state experiment shown in the inset of Fig. . Here, Laser-1 is on all the time and at $t = 10$ s also Laser-2 is used to excite the sample. When Laser-2 is switched on, the SAW attenuation $s\Gamma$ drops by $\sim 25\%$ since additional activation and subsequent recombination of carriers trapped in RPFL occurs. After Laser-2 is switched off again the original attenuation level is recovered.

To further support our model based on RLPF we performed SAW-attenuation experiments under continuous excitation using Laser-1 in which the optical pump power was varied and thus, the generation rate of excitons. Assuming that the rates of radiative decay of excitions and recombination of charges trapped in RLPF are independent on the optical pump power, a steady state occupation of RLPF and number of mobile carriers will build up for a given optical power. The measured attenuation as a function of the excitation intensity is shown as symbols in Fig. . With increasing laser intensity the SAW attenuation of the third harmonic and, therefore, the sheet conductivity rises and at about ${I}\approx 0.1\cdot I_{0}$, the sheet conductivity matches the characteristic conductivity. The observed behavior can be fitted in this excitation power range using $$s\Gamma \left(I/I_{0}\right) = \cfrac{K_{\rm eff}^{2}\cdot s\cdot k}{2}\cfrac{\dfrac{\sigma( I/I_{0})}{\sigma_{m}}}{1+\left(
    \dfrac{\sigma( {I}/{I_{0}})}{\sigma_{m}}\right)^{2}}
    \label{eq:fit-f}$$ and assuming for $\sigma$ a single exponential increase given by $\sigma({I}/I_{0})\propto \exp({I}/I_{0})$. Clearly, the result of this fit (solid line in Fig. ) agrees well over a wide range of intensities. The initial increase of the conductivity can be well understood by assuming a Gaussian variation of the Cd-content from the nominal value of 0.1 and, thus, a Gaussian distribution of $E_{rec}$. Furthermore, this distribution can be correlated with the distribution of the electron kinetic energy in the RLPF: At $T=60$ K the conduction is in the percolation state and, therefore, at low laser intensities the small number of generated electrons is mainly trapped by deep localization centers with an increased Cd-content, minimizing their kinetic energy. With increasing excitation intensity and electron density, electrons have to populate the shallower localization centers with reduced Cd-content. This leads to a strong increase of total kinetic energy of the trapped electrons and explains the exponential rise of sheet conductivity. The overall good agreement between the experimental data and the fit over a wide range of excitation power proves both the validity of equation () and the assumption made for $\sigma({I}/{I_{0}})$. In addition, the observed saturation of the signal is a clear indication to a finite number of localization site available which are fully populated for laser intensities $> 0.2\cdot I_0$.
At this point we want to note that the applied SAW-technique is capable to sensitively detect conductivity at very low optical pump powers. In particular we want to note that the minimum excitation power required to detect PL signals from our samples was in the order of $\sim0.05\cdot I_0$ setting the lower bound at which this invasive method can be applied. In contrast, at this particular pump power already a pronounced SAW-attenuation of $s\Gamma = 6$ is detected. This attenuation level is a factor of 20 - 30 larger than minimum attenuation levels of $0.2\pm0.1$ determined from the baseline and noise in Fig. . In addition to this high sensitivity our SAW-attenuation method does not rely on high SAW powers. In particular it is capable to detect in the low SAW power regime where no influence on the PL signal is detected or charge conveyance effects occur. Thus, this method does not perturb the carrier system in study and resembles a non-invasive probe of the conductivity. In addition, we have shown that despite of the non-linear dependence on the conductivity given by the SAW-response on the conductivity can be analyzed and fitted accurately \[Fig. \] .

## Conclusions

In conclusion we studied the attenuation of a SAW in a ZnCdSe/ZnSe-LiNbO$_3$-hybrid. Under optical excitation we observe a variation of the SAW attenuation due to persistent photoconductivity effects in the ternary II-VI alloy. This PPC arises from carriers trapped in and activated from random potential fluctuations due to compositional fluctuations in the material. By performing temperature and excitation energy and power dependent experiments impurity or defect related effects can be excluded. It is shown that the carriers trapped in RLPF can be activated optically and thermally. In the temperature dependent experiments a clear transition from a hopping to a percolation regime is observed which also manifests itself in a suppression of bound exciton emission. We find that the critical temperature below which localization occurs is reduced from $T=120$ K for bulk crystals to $T<40$ K in high qualitiy QWs grown by MBE. The corresponding shallow potential modulation of $E_{rec}=\lesssim 3$ meV underlines the high structural quality achievable in MBE-grown II-VI heterostructures. Furthermore, direct comparison between the applied techniques demonstrates that compared to PL spectroscopy the SAW-based method yields a very versatile, non-invasive tool for sensitive conductivity measurements which can be applied to other material systems than the one studied in this work.
This work was supported in parts by the German federal government as part of the Cluster of Excellence *Nanosystems Initiative Munich* (NIM).