Spontaneous Synchronization of Two Bistable Pyridine-Furan Nanosprings Connected by an Oligomeric Bridge

1. Introduction

Spontaneous synchronization is a phenomenon in which a number of coupled dynamical systems self-organize to behave consistently over time. Synchronization phenomena are observed in various physical, chemical, and biological systems over a wide range of spatial and temporal scales (see, for example, Neiman et al. , VanWiggeren and Roy , Jung et al. , Gandhimathi et al. , Kenfack and Singh or Muralidharan and Ali [1,2,3,4,5,6]). The concept of synchronization is especially important in the design of nanoscale devices since there is often a need to amplify the weak output signals from individual molecular-size functional units while maintaining the high sensitivity of the device to weak stimuli.

A canonical example of a spontaneously synchronized system is a set of mechanical or electric oscillators that become coupled with each other via interacting forces.[7,8,9,10,11]. Another canonical example is the spontaneous synchronization of spontaneous vibrations (SV) of coupled bistable systems, i.e. synchronization of the noise-activated random transitions between the states of individual bistable systems. [12]. Unlike an oscillatory system that has eigenfrequency, SV of a bistable system is characterized by a wide frequency distribution[13,14,15]. Nevertheless, SV of coupled bistable systems can also exhibit spontaneous synchronization [12,16] despite the fact that the transitions between the states of the bistable systems remain random.

In our recent papers Avetisov et al. [17], Markina et al. [18], Avetisov et al. [19]) we have shown that nano-sized oligomeric structures stabilized by short-range low-energy interactions, such as weak hydrogen bonds, hydrophilic-hydrophobic interactions, and $\pi$ stacking, could exhibit bistability and SV caused by fluctuations coming from the ambient heat bath. Using molecular dynamic simulations, we demonstrated that short pyridine-furan oligomers a few nanometers in size (named oligo-PF springs), which helical shape was stabilized in water by the $\pi -\pi$ interactions of aromatic groups, exhibited bistability, SV, and the stochastic resonance (SR) effects (see Avetisov et al. [19]). All these effects are characteristic of the classical example of the springs with nonlinear elasticity known as Duffing oscillators. As described in (see Avetisov et al. [19]), these effects appeared in oligo-PF springs due to the contribution of $\pi -\pi$ interactions into the spring stretching energy.

Bistable oligo-PF springs, being two-state dynamic systems, are attractive for applications because the thermal fluctations present at room temperature are sufficient for them to work. However, experimenting with such a small spring and its implementation as an operational unit is challenging since it might require detection extremly weak input signals and sophisticated single-molecule manipulations. It seems desirable to amplify a weak response of an individual bistable nanospring but retain those inviting characteristics that are due to the nanosize of the springs, such as bistability and spontaneous vibrations activated by normal thermal noise from the surrounding thermal bath. This is a non-trivial task because a simple lengthening of an PF-nanospring up to several turns makes it’s vibrational dynamics multimode due to extra degrees of freedom of movements of the turn relative each other. As a result, the spring bistability itself becomes ambiguous. In addition, an increase in the size of the nanospring, even maintaning bistability, would lead to a rise in the activation barrier for SV, and the finest sensitivity of the nanospring to weak perturbations could be lost. Therefore, if one wants to scale the finest characteristics of bistability of individual nanosprings, then coupling and synchronization of the nanosprings seems like a suitable solution. The question is what kind of coupling can lead to the synchronization of SV of oligo-PF springs.

In this paper, we present the simulations of long-term molecular dynamics of two bistable oligo-PF springs in water, which are coupled to each other by an oligo-pyrrole bridge. We show that such coupling leads to significant synchronization of the SV of the nanosprings. We also show that the synchronization effect is enhanced in the SR regime of the PF-spring vibrations. The paper is organized as follows. In Section 1, we introduce the theory of spontaneous synchronization of bistable systems as a footing for computer simulation studies of the dynamics of coupled pyridine-furan nanosprings and provide details of the computer model and simulations. In Section 2, we present the computer simulation data and discuss the synchronization effects in the regimes of SV and SR. The paper ends with the discussions and conclusions sections, Section 3 and Section 4.

2. Materials and Methods

2.1. Theoretical modelling: Langevin approximation

Here we introduce the theoretical underpinning of the bistable system’s synchronization to align the computer simulation studies presented below with the concept of synchronization. First, let us consider the overdamped regime of two coupled identical Duffing oscillators firstly discussed by Neiman in [12]. The two-dimensional dynamics of the system in dimensionless units is described by Langevin equations:

$$\left\{\begin{array}{c}\frac{dx}{dt}=-\frac{dU(x,t)}{dx}+ϵ{\xi }_1\left(t\right)+G_1(x,y)\hfill \\ \frac{dy}{dt}=-\frac{dU(y,t)}{dy}+ϵ{\xi }_2\left(t\right)+G_2(x,y),\hfill \end{array}\right.$$

(1)

where x and y are the degrees of freedom of the first and second Duffing oscillator, respectively, ${\xi }_i\left(t\right)$ is a delta-correlated white noise, $ϵ$ is a noise amplitude, $U\left(x\right)=-\frac{\alpha }{2}{x}^2+\frac{1}{4}{x}^4+E_{0}xcos\left(\omega t\right)$ is a periodically modulated double well potential. The parameter $\alpha$ determines the positions of the potential minima, $x_{min}=\pm \sqrt{\alpha }$, and the bistability barrier, $\Delta U={\alpha }^2/4$, separating the minima. The coupling forces $G_i(x,y)$ are introduced in the form:

$$\left\{\begin{array}{c}{G}_1(x,y)=b(y-x)\hfill \\ G_2(x,y)=b(x-y)\hfill \end{array}\right.$$

(2)

where $b>0$ is the coupling constant. It is easy to see that the coupling tries to synchronize the dynamics of bistable systems by minimization of the difference between x and y.

In what follows, we will use numerical simulations of the equations (1) using, e.g., Heun’s method, to represent the synchronization effect. Specifically, we choose dimensionless parameters $\alpha =10$, ${ϵ}^2=18$ for which the Kramers rate $r_K\propto exp\left(-\frac{{\alpha }^2}{2{ϵ}^2}\right)$ is close to those obtained in molecular simulations, and we varied coupling constants b from $b=0-0.75$. Given potential parameters, the critical value of the external field, at which the energy injected by the field in one period is comparable to the barrier value, is $E_{cr}=\Delta U/x_{min}\approx 7.9$.

In the absence of coupling and periodic modulations ($b=0$ and $E=0$), the random nose activates jump-like transitions over the bistability barrier with with the mean life time of the states corresponding to Kramers rate.

In this case, the system of two Duffing oscillators has four attractors shown on the left panel of Fig. 1, which correspond to the independent SV of the Duffing oscillators. The same brightness of the attractors indicates that the system equally visits them. In the case of coupling (Figure 1, right panel), two attractors are visited more often, indicating that random transitions of two Duffing oscillators are synchronous.

A weak periodic field applied to the system (see (1)) cyclically tilts the bistable potential of each Duffing oscillator and induces periodic modulation of the Kramers rate of the transitions between the two states. It is well known that if the period of modulation coincides with the double inverse $r_K$, the SR occurs, and the random transitions between the two states become highly regular. We refer readers interested in fine-tuning the SR time-scale matching conditions to Gammaitoni et al. [14] and references therein.

It is clear that modulation of Duffing oscillators by a periodic field should make an additional contribution to the synchronization effect. However, even in SR mode, the improvement in synchronization depends much more on the coupling than on the modulating field (see Figure 2).

2.2. Pyridine-Furan Springs

As described in the previous section, the synchronization is achieved by an appropriate coupling of two bistable systems. In our previous paper [19] we have shown that oligo-PF springs a few nanometers in size subjected to critical tension exhibit bistable dynamics characteristic of Duffing oscillators. Accordingly, they were also chosen for this paper to study the synchronization phenomena.

Pyridine-furan copolymer (Figure 3a) is a conductive polymer consisting of 5- and 6-member heterocyclic pyridine and furan rings (respectively)[20]. These PF copolymers in the cis-configuration tend to form a helix-like shape, which is stabilized by the $\pi -\pi$ interactions of aromatic groups located on the adjacent turns of the spring [21]. In particular, the cis-configuration of oligo-PF of 5 monomer units, oligo-PF-5, was selected for study in this paper (see Figure 1a-b). Since the coupling between the two oligo-PF springs was supposed to be like a rigid oligomeric bar, the various chains for the bar were tested, and an optimal option was defined as an oligo-pyrrole chain of six monomer units in length (see Figure 3), which provided a smooth but rigid connection between both springs.

2.3. Simulation Details

Molecular dynamics is the best simulation method available to model effects like spontaneous synchronization where the composite behavior of the system depends on the physical motion of molecules. The PF-nanospring system was modeled in a fully atomistic representation in the canonical ensemble (number of particles/volume/temperature [NVT]). The springs were surrounded by enviornmental water, and the total simulation box size was $7.0*7.0*7.0$$n{m}^3)$. Simulations were conducted with a time step of 2 fs using the Gromacs 2019 package [22], OPLS-AA [23] force field parameters for the oligomer, and the SPC/E model [24] for water (for more details, see Parameters for Molecular Dynamics Simulation section of Supplementary Materials). Following characterization of the equilibrated state of PF-springs at $280\mathrm{K}$[21], this temperature was set by the velocity-rescale thermostat [25]. Each dynamic trajectory was 300–350 ns long and was repeated three times to obtain better statistics; therefore the effective length of the trajectories was about 1 $\mu s$ for each sample.

At first, two oligo-PF-5 springs with fixed bottom (conditionally) ends and movable upper ends connected by an oligomeric bar were equilibrated at $280\mathrm{K}$. For convenience and in accordance with the designations introduced in Section 2.1 we call (conditionally) left spring X and right spring Y. Then, springs were subjected to an external tensile forces F (in Figure 3c): at that, the bottom ends of the springs were kept fixed, and the upper ends were pulled by forces applied normally to the connecting bar. The end-to-end distance of each of the oligo-PF-5 springs, $R_e^X$ and $R_e^Y$, (see Figure 3b) was chosen as an order parameter describing the long-term dynamics of the spring. Vibrations of coupled oligo-PF-5 springs occur between two well-reproduced states of each spring with end-to-end distances equal to $R_e\approx 1.0nm$ and $R_e\approx 1.4nm$. These states are referred to as the squeezed and the stress-strain states, respectively[19].

The mutual movements of the springs were plotted into 2D phase diagrams, which were the long time statistics of the 2D-probability density for simulteneously visiting states $R_e^X$ and $R_e^Y$, respectively. In Figure 4 and Figure 6,brighter colors correspond to higher probability density and dark violet means zero probability density. The statistics of the two states map were extracted directly from the $R_e\left(t\right)$ series. The synchronization coefficient was defined as the height of the main peak of the normalized cross-correlation function of time trajectories of $R_e^X$ and $R_e^Y$.

3. Results

3.1. Synchronization of spontaneous vibrations of PF-5 springs

Next, the dynamics of the oligo-PF-5 springs connected by an oligomeric bridge were examined as they were subjected to pulling forces of different values. Under weak tension, the springs were stretched slightly in accordance with linear elasticity. However, as soon as the pulling force reached a critical value $F=220pN$, the coupled oligo-PF-5 springs became bistable and started to vibrate spontaneously similar to SV of a single oligo-PF-5 spring described in detail in our previous paper[19]. Above this critical load, the springs spontaneously vibrate between the squeezed and stress–strain states. SV are observed in a wide region of pulling forces from $220pN$ to $300pN$. In the regime of SV, the end-to-end distance fluctuations of each individual PF-spring have a mean amplitude of 0.2 nm, so the stress–strain states of both springs are clearly distinguished from the squeezed states.

Panel a of Figure 4 shows 2D-phase diagram of two copled PF-springs, i.e., 2D stationary probability density for simultaneous location of the ends of left and right springs at the positions $R_e^X$ and $R_e^Y$, respectively. Figure 4b shows a shortened part of typical SV-trajectories (for the full trajectory see the Supplementary Materials) of the coupled oligo-PF-5 springs in the symmetric bistability region ($F=250$) pN. In this regime, SV are very pronounced, because neither the squeezed nor stress–strain states dominate, so the mean lifetimes of the two states for both springs were approximately the same, and equal to $\tau =6.14$ ns. The bistability barrier of the oligo-PF-5 springs is circa 10 $k_{B}T$, as was shown in our previous paper[19]. Therefore, twith the except of value of the critical pulling force, all other characteristics of bistability of individual PF-springs in coupled system. e.g. the amplitude of vibrations, mean lifetimes of the squeezed state and the stress-stain state, the intensity of thermal noise activating the spontaneous vibrations, etc., are close to the same characteristics of a freestanding PF-5 spring.

Figure 4. SV of two coupled oligo-PF-5 springs at $F=250$ pN: (a) 2D stationary probability density for simultaneous location of the left and right springs at the positions $R_e^X$ and $R_e^Y$, respectively; (b) corresponding trajectories of the end-to-end distances $R_e^X$ and $R_e^Y$.

Figure 4. SV of two coupled oligo-PF-5 springs at $F=250$ pN: (a) 2D stationary probability density for simultaneous location of the left and right springs at the positions $R_e^X$ and $R_e^Y$, respectively; (b) corresponding trajectories of the end-to-end distances $R_e^X$ and $R_e^Y$.

Figure 5 a shows cross-correlation function between SV-trajectories of the left and right springs coupled by an oligomeric bar. The power spectrum of this function is given on Figure 5b. The cross-correlation function clearly demonstrates that spontaneous vibrations of coupled oligo-PF-5 nanosprings are significantly synchronized with the correlation coefficient $k\approx 0.5$. The power spectrum of the autocorrelation function indicates that the spectra of spontaneous vibrations of both nanosprings do not have natural frequencies, with the exception of a weak specific response near the so-called SR-frequency, the reciprocal of which is equal to twice the lifetime of the state when the system experiences spontaneous oscillations [13,14,15].

Figure 5. Correlation between SV of the PF-nanosprings coupled by an oligomeric bar: (a) cross-correlation function for the SV-trajectories of the left and right springs; (b) power spectrum of the cross-correlation function.

Figure 5. Correlation between SV of the PF-nanosprings coupled by an oligomeric bar: (a) cross-correlation function for the SV-trajectories of the left and right springs; (b) power spectrum of the cross-correlation function.

3.2. Synchronization of the spring vibrations in the stochastic resonance regime.

Next, we examined the SR mode of the coupled oligo-PF-5 springs by applying an additional oscillating forces that weakly modulated the upper ends of the coupled springs. The oscillating electrical field $E=E_{0}cos\left(2\pi \nu t\right)$, $E_0$ = 0.1 V/nm, $\nu =1/12.5$ GHz acting on unit charges present at the upper ends of the springs was used to model SR, while the compensating charges were at the bottom (fixed) ends of the springs (for more details, see Parameters of the periodic signal section of Supplementary Materials). Note, that the frequency $\nu =1/12.5$ GHz is the exactly the SR frequency equal to the half of the Kramers frequency of the transitions over the bistability barriers of individual PF-springs. Typical vibrations of coupled oligo-PF-5 springs in the SR mode are shown on Figure 6.

Figure 6. Vibrations of the oligo-PF-5 springs coupled by an oligomeric bar in the SR-regime: (a) 2D stationary probability density for simultaneous location of the left and right springs at the positions $R_e^X$ and $R_e^Y$, respectively; (b) typical SR-trajectories $R_e^X$ and $R_e^Y$ of individual PF-5 springs.

Figure 6. Vibrations of the oligo-PF-5 springs coupled by an oligomeric bar in the SR-regime: (a) 2D stationary probability density for simultaneous location of the left and right springs at the positions $R_e^X$ and $R_e^Y$, respectively; (b) typical SR-trajectories $R_e^X$ and $R_e^Y$ of individual PF-5 springs.

The cross-correlation function between the SR-trajectories of the left and right springs coupled by an oligomeric bar is shown in Figure 7a. The power spectrum of this function is shown on Figure 7b. The cross-correlation function clearly demonstrates that SR enhance synchronization effect: the correlation coefficient is raised up to $k\approx 0.7$. In its turn, the power spectrum of the correlation function indicates that close to the half of the Kramers frequency of the transitions over the bistability barriers the vibrations of the PF-springs have well-defined SR signal. Thus, the synchronization effect can be improved by applying a weak oscillating force with SR-frequency.

Figure 7. Correlation between vibrations of the PF-springs coupled by an oligomeric bar in the SR regime: (a) cross-correlation function of the vibrations the left and right springs; (b) Power spectrum of the cross-correlation function.

Figure 7. Correlation between vibrations of the PF-springs coupled by an oligomeric bar in the SR regime: (a) cross-correlation function of the vibrations the left and right springs; (b) Power spectrum of the cross-correlation function.

4. Discussion

Based on the results, the main finding is that coupled oligo-PF-5 springs exhibit synchronized dynamics in SV and SR regimes. Similar to a single spring, two coupled springs under the critical stretching show SV. The value of the critical tensile in the case of coupled springs is slightly lower than for a single oligo-PF-5 spring (as previously reported[19]). The difference in critical tensiles and dynamic behavior between single and coupled oligo-PF-5 springs arises since coupled springs are bridged by a connectivity bar, which causes additional tensiles as compared to single springs. This connectivity bar could form an additional $\pi$ bond with the springs and thus affect the bistability area and the SV scale. This means that the collective dynamics of any hybrid system comprised of a complex set of coupled springs should be checked independently since specific couplings could introduce new particular tensions and additional intramolecular interactions.

When spontaneous vibrations occur, the presence of a connectivity bar allows the non-correlated jumping behavior to transform into correlated transitions. This would represent a transition either from a squeezed state to a stress-strain state and vice-versa. The correlation coefficient associated with a such transition is circa $0.5$ (i.e. partially correlated regime). Notably, the correlation can be improved in the regime of SR by an introduction of a weak external force. The external force causes a well-defined frequency and amplitude to appear that characterizes the system dynamics. In this regime, the correlation coefficient grows up to $0.7$ (i.e. high correlation regime).

These findings are consistent with the mathematical model described in Section 2.1 in terms of synchronization upon introduction of the bond between the bistable systems. However, Figure 2 shows that at the SR regime the correlation coefficient between the bistable systems has a very weak dependence on the modulation amplitude E. In other words, there is no significant difference in correlation coefficient in regimes of SV (when $E=0$) and SR ($E>0$). Also, there are no shifts between the trajectories in the SV regime. On the other hand, in the simulated system the correlation coefficient differs significantly between the SV and SR regimes and at the same time the trajectories of SV for left and right springs were observed to be time-shifted by 3 ns (the peak of the cross-correlation function for $R_e^X$ and $R_e^Y$ isn’t located at zero). The shifted maximum of the cross-correlation function indicates a systematic time delay in the dynamics of the springs. The finite rigidity of the connectivity bar and the system’s geometrical freedom are the reasons why one spring doesn’t feel the displacement of another spring immediately. As a result, one of the springs becomes a driver for the synchronous behavior of both. Importantly, no shift is observed in the SR regime since the movement of each spring was as well consistent to the oscillations of the weak external electrical field. Such difference between the theory and simulations can be explained due to the mathematical model of the bond, which is in theory infinitely rigid and immediately transfers the displacements between the bistable systems.

Based on the results presented above, the following best practices can be recommended for designing a particular connecting bar considering the key variables of rigidity and length. In particular, to build a working bistable and controllable structure, it is important to first select a connectivity bar of high rigidity (i.e. high persistent length). This can be a broad class of $\pi$-conjugated structures which includes pyrrole, as used in this study. The more rigid connection is introduced, the stroger coupling of bistable systems is achieved and thus the synchronization effect increases.

The second crucial design factor is the connectivity bar’s length. A small bar obstructs the mutual motion of the coupled springs due to sterical hindrance. This hindrance occurs due to the complex nature of spring motion. The end of spring is moving not only up, but also a bit sideways (motion reminding a helical path). These tiny sideways motions are oppositely directed for left and right springs, thus causing slight rotation in the plane of the connectivity bar around its center. If the connectivity bar is short, the rotation angle is limited, and therefore sideways movements are limited. As a consequence vertical movements are locked as well.

A connectivity bar of longer length will avoid sterical obstruction. However, an increase of the length of the connectivity bar naturally causes a decrease in its rigidity and causes a delay between the displacements of left and right oligo-PF springs. A large connectivity bar will work as a big lever arm and force springs to move further from their equilibrium, causing an artificially large amplitude. In other words, interplay between the elastic properties of the oligo-PF springs and the torsional deformation of the crossbar are defining the overall system dynamics. Taking into account both these factors, choosing a semirigid bar of circa 2 nm in size would be preferred - one that is rigid enough to provide a bond but does not influence the motion of the springs. Future conjugated structures that have high persistence lengths could be considered as candidates for a best-performing coupling.

5. Conclusions

In this paper, we have presented a study of the mechanic-like behavior of a novel nanospring system composed of two bistable oligomeric pyridin-furan springs coupled by an oligomeric pyrrole bridge. Through molecular dynamics simulations in all-atom representation, we have investigated the long-term conformational dynamics of the system and found that the sponteneus vibrations of coupled springs become synchronized.

The system exhibits synchronized vibrations under certain conditions, which can be controlled, for example, by changing the tensile forces, the length and stiffness of the oligomeric bridge, and the elasticity of the springs. The simulations of two PF-nanospring coupled by an pyrrole bridge show that SV of the springs are synchronized with the correlation coefficient $k\approx 0.5$. In the case of spontaneous synchronization of oligo-PF-5 springs in the SR mode, the correlation coefficient reaches a relatively high value of $k\approx 0.7$. In summary, our study provides new insights into the bistable nanomechanics of coupled PF-nanosprings and highlights the fact that the synchronization can be constructively achieved in relatively simple ways.