Light-Fueled Synchronization of Two Coupled Liquid Crystal Elastomer Self-Oscillators

The synchronization and group behaviors of self-excited coupled oscillators are common in nature and deserve to be explored, for self-excited motions have the advantages of actively collecting energy from the environment, being autonomous, making equipment portable, and so on. Based on light-powered self-excited oscillators composed of liquid crystal elastomer (LCE) bars, the synchronization of two self-excited coupled oscillators is theoretically studied. Numerical calculations show that self-excited oscillations of the system have two synchronization modes, in-phase mode and anti-phase mode, which are mainly determined by their interaction. The time histories of various quantities are calculated to elucidate the mechanism of self-excited oscillation and synchronization. For strong interactions, the system always develops into in-phase synchronization mode, while for weak interaction, the system will evolve into anti-phase synchronization mode. Furthermore, the effects of initial conditions, contraction coefficient, light intensity, and damping coefficient on the two synchronization modes of the self-excited oscillation are investigated extensively. The initial condition generally does not affect the synchronization mode and its amplitude. The amplitude of self-oscillation always increases with increasing contraction coefficient, gravitational acceleration, and light intensity, while it decreases with the increasing damping coefficient. This work will deepen people’s understanding of the synchronization behaviors of self-excited coupled oscillators, and the theoretical framework could be extended to scenarios involving large-scale synchronization of the systems with numerous interacting oscillators.


Introduction
Self-excited oscillation is a phenomenon of periodic state change within a system under constant external stimulation [1][2][3][4][5][6][7].Because of its unique advantages, self-oscillation has wide application prospects in many fields, such as energy harvesting [8,9], signal sensors [10], soft robotics [11][12][13][14][15], mechano-logistic devices [16,17], biomimetic designs [18], and so on.First, self-oscillation can directly harvest energy from the constant external environment in order to maintain its own periodic motion, which is similar to the effect of biological active feeding.Second, the periodic motion of self-oscillation does not need periodic external stimulation, only constant external stimulation.This feature greatly reduces the requirements of system motion control and does not need a complex controller.Third, the above characteristics of no controller and battery greatly reduce the complexity of the active machine and make it more portable, which is expected to achieve high power [19][20][21].In many conventional and active material systems, different self-excited oscillations have recently been constructed [22][23][24].
Although a great deal of work has been carried out on single self-excited oscillators, the interaction and collective motion of multiple self-excited oscillators need to be further explored.Synchronization and collective motion are ubiquitous in nature, such as the circadian rhythm and cardiac pacemaker cells [56][57][58][59][60].The first work on synchronization can be traced back to Huygens' observations in 1665 on the synchronization of coupled pendulums [61].He observed that two identical clocks oscillate synchronously with two pendulums in opposite directions.Recent studies have confirmed that the synchronization between two pendulums results from the coupling caused by small mechanical vibrations that propagate in the wooden structure connecting the clocks [62].Furthermore, the synchronous movement of a large number of metronomes with more degrees of freedom on a free moving base is also demonstrated experimentally [63].Recently, based on lightresponsive LCE, Ghislaine et al. experimentally studied the synchronized oscillations of thin plastic actuators fueled by light and found two kinds of synchronous oscillation phenomena of in-phase and anti-phase in steady state [64,65].Their numerical simulations qualitatively explain the origin of the synchronized motion and found that the motion can be tuned by the mechanical properties of the coupling joint.
At present, there are few studies on the interaction and group phenomenon of selfexcited coupled oscillators based on active materials [64], and the synchronization mechanism and its behavior need to be further explored.In this paper, based on a self-excited oscillator made of photoresponsive LCE previously proposed by us [66], we investigate the synchronous behavior of two identical self-excited oscillators powered by steady illumination.This paper is as follows.First, based on dynamic LCE model [67], the dynamic governing equation for two identical self-excited oscillators under steady illumination is formulated in Section 2. Second, two kinds of synchronization mode of the self-excited oscillations are presented in Section 3. In Section 4 and 5, the detailed self-excited mechanism of in-phase and anti-phase modes are elucidated, respectively.Meanwhile, the influences of initial conditions and spring constant on the synchronization mode, amplitude, and period of the self-excited oscillations are investigated.Finally, the conclusion is given in Section 6.

Model and Formulation
In this section, a theoretical model is proposed for two identical LCE oscillators connected by torsional springs under steady illumination.It mainly includes the dynamic model of two LCE oscillators, evolution law of number fraction in the two LCE oscillators, non-dimensionalization and solution method of the differential governing equations with variable coefficients.can oscillate, because the light-driven deformation periodically changes the center of gravity and reverses the resultant moment of the system due to the self-shadowing effect.In the following, we will investigate the synchronization of the two self-excited oscillations.According to the theorem of moment of momentum, the differential equations for the dynamics of the two LCE bars rotating around the fixed z axis are

Dynamic Model of the two LCE Oscillators
where the angular momenta of where the moments of inertia of about z axis are, respectively, where is the moment of inertia of about z axis, and According to the theorem of moment of momentum, the differential equations for the dynamics of the two LCE bars rotating around the fixed z axis are where the angular momenta of A 1 O 1 B 1 and A 2 O 2 B 2 are, respectively, where the moments of inertia of A 1 O 1 B 1 and A 2 O 2 B 2 about z axis are, respectively, where It is worth noting that the length depends on the light-driven contraction strain.Then, the length l 1 of the O 1 A 1 and the length l 2 of the O 2 A 2 can be expressed as where ε 1 (t) and ε 2 (t) are the light-driven contraction strains of O 1 A 1 and O 2 A 2 , respec- tively.For simplicity, we assume that the light-driven contraction strain of the material is proportional to the number fraction ϕ(t), where C 0 is the contraction coefficient.The number fraction ϕ(t) is given in Section 2.2.In Equation ( 1), M z1 and M z2 are the resultant moments of all external forces of bar where the driving moments of A 1 O 1 B 1 and A 2 O 2 B 2 can easily be calculated as where g is the gravitational acceleration.
The damping force is assumed to be proportional to the velocity, and then the damping moments of A 1 O 1 B 1 and A 2 O 2 B 2 can be easily calculated as where β is the damping coefficient, dθ 1 dt = .
θ 1 is the angular velocity of The moment exerted by the torsion spring on the two bars is assumed to be linear with the angle difference, where α is the spring coefficient of the torsion spring.

Evolution Law of Number Fraction in the Two LCE Oscillators
In order to calculate the light-driven contraction strain and the lengths of the LCE bars, we need to obtain the number fractions in the LCE bars.According to the research of Yu et al., the trans-to-cis isomerization of LCE can be induced by UV or laser with wavelength less than 400 nm [68].Generally, the cis-to-trans isomerization driven by UV light and the thermal trans-to-cis excitation can be neglected.Therefore, the number fraction of cisisomers depends on the thermal excitation from trans to cis, the thermally driven relaxation from cis to trans and the light-driven trans-to-cis isomerization.Then, the number fraction of bent cis-isomers in LCE can be governed by the following equation [67,68], where T 0 is the thermal relaxation time of cis state to trans state, I 0 is the light intensity, and η 0 is the light absorption constant.By solving Equation (10), the number fraction of cis-isomers can be expressed as: where ϕ 0 is the number fraction of cis-isomers at t = 0.In the light zone, for initially zero number fraction of cis-isomers, i.e., ϕ 0 = 0, Equation ( 11) can be simplified as In the dark zone, namely I 0 = 0, ϕ 0 can be set as the maximum value of ϕ(t) in Equation (12), namely, ϕ 0 = η 0 T 0 I 0 η 0 T 0 I 0 +1 , and Equation ( 11) can be simplified as: Polymers 2023, 15, 2886 5 of 23

Solution Method
Equations ( 16)-( 19) are ordinary differential equations with variable coefficients, and there are no analytic solutions.Hereon, the classical fourth-order Runge-Kutta method is utilized to solve ordinary differential equations by MATLAB software.Extensively used Runge-Kutta explicit iterative formulas are detailed as follows, where H is the time step, and k i (i = 1 to 4) are listed as below, Each of Equations ( 16)- (19), second-order ordinary differential equations (ODEs), needs to be transformed into two first-order ODEs in accordance with the numerical rules  21) takes a two-dimensional vector form by zoning as below, in the light zone: in the dark zone: Equations ( 22)-( 25) in combination with the initial conditions θ 0 1 , .θ 0 1 , and θ 0 2 , .θ 0 2 constitute the closed system of dynamic equations describing the two coupled LCE oscillators connected by a torsion spring under steady illumination.The self-sustained motion of the LCE oscillators, i.e., the variation of angle and angular velocity with time, can be obtained by iteration.

Two Synchronization Modes
In this section, based on solving the governing Equations ( 22)- (25), we present two synchronization modes, namely in-phase mode and anti-phase mode.In order to study the synchronization mode of two LCE oscillators under steady illumination, it is necessary to determine the typical values of the dimensionless parameters.The typical material properties and geometric parameters obtained from the existing tests [37,69] are shown in Table 1, and the corresponding dimensionless parameters are shown in Table 2.In the following, these values of parameters are used to study the synchronization modes of two LCE oscillators under steady illumination.Figure 2 shows two typical synchronization modes of self-excited motion of LCE oscillators: in-phase mode and anti-phase mode.In the computation, we fix I = 0.33, The time history curve and domain of attraction of the in-phase mode for α = 0.1 are given in Figure 2a,c, respectively.The results show that the two LCE bars oscillate in in-phase mode.Figure 2b,d present the time history curve and domain of attraction of the anti-phase mode for α = 0.03, respectively.The results show that the two LCE bars oscillate in anti-phase mode.Through careful calculation, it is found that there exists a critical spring constant α crit = 0.075 for the two synchronization modes.The two LCE oscillators will oscillate in in-phase mode for α > α crit , while they oscillate in anti-phase mode for α < α crit .
In the following, we will discuss the two synchronization modes in turn.

In-Phase Synchronization Mode
In the above mechanical model of two coupled LCE oscillators, there are five dimensionless parameters, including I , g , α , β , and 0 C .In this section, the mechanism of self-oscillation in in-phase mode, and the influences of dimensionless parameters on the in-phase synchronous mode are studied in details.

Mechanisms of the Self-Excited Oscillation in In-Phase Mode
To investigate the mechanism of the self-excited oscillation in in-phase mode of the two LCE oscillators under steady illumination, Figure 3 presents time-histories of some key physical quantities of the in-phase mode in Figure 2b,d.Figure 3a plots the time histories of 1 θ or 2 θ , which shows that the two LCE bars oscillate periodically in in-phase mode.Figure 3b plots time histories of the number fractions of cis-isomers in the two LCE bars.The number fractions of cis-isomers increase in the light zone, while they decrease in the dark zone.Therefore, the contraction strains in the two LCE bars increase in the light zone while decreasing in the dark zone, as shown in Figure 3c.In Figure 3d, the moment of the torsion spring on the two bars is zero, because the angle difference between the two bars in in-phase mode is zero.In Figure 3e, the driving moments of the two oscillators also change periodically in in-phase mode.Figure 3f

In-Phase Synchronization Mode
In the above mechanical model of two coupled LCE oscillators, there are five dimensionless parameters, including I, g, α, β, and C 0 .In this section, the mechanism of self-oscillation in in-phase mode, and the influences of dimensionless parameters on the in-phase synchronous mode are studied in details.

Mechanisms of the Self-Excited Oscillation in In-Phase Mode
To investigate the mechanism of the self-excited oscillation in in-phase mode of the two LCE oscillators under steady illumination, Figure 3 presents time-histories of some key physical quantities of the in-phase mode in Figure 2b,d.Figure 3a plots the time histories of θ 1 or θ 2 , which shows that the two LCE bars oscillate periodically in in-phase mode.Figure 3b plots time histories of the number fractions of cis-isomers in the two LCE bars.The number fractions of cis-isomers increase in the light zone, while they decrease in the dark zone.Therefore, the contraction strains in the two LCE bars increase in the light zone while decreasing in the dark zone, as shown in Figure 3c.In Figure 3d, the moment of the torsion spring on the two bars is zero, because the angle difference between the two bars in in-phase mode is zero.In Figure 3e, the driving moments of the two oscillators also change periodically in in-phase mode.Figure 3f delineates the dependence of the driving moment on θ 1 or θ 2 .In Figure 3f, the area surrounded by the closed curve represents the net work performed by the steady illumination during one cycle of the self-excited oscillation, which compensates for the energy dissipation caused by the damping to maintain the oscillation of the LCE oscillators.

Effect of the Spring Constant on the In-Phase Mode
Figure 4a,b plot the time histories of the two oscillators for different spring constants.In this computation, we set 4a,b, for the two different spring constants, the two curves of the two LCE oscillators coincide after a period of time, which means that the two bars oscillate synchronously in in-phase mode.Figure 4c,d plot the phase diagrams for the two spring constants, in which 1 θ and 2 θ evolve from the initial disorder into a final attraction domain.Figure 4e further plots the domains of attraction of 1 θ and 2 θ for different spring constants.The calculation shows that the two domains of attraction are identical.The results show that the spring constant has no effect on the in-phase mode.
Furthermore, Figure 4f describes the limit cycles of 1 θ and 1 θ  for different spring constants.Similarly, the two limit cycles are also the same.This implies that the spring constant has no effect on the amplitude and frequency of the LCE oscillators.More calculations show that the influence of the spring constant is the same for In this computation, we set 4a,b, for the two different spring constants, the two curves of the two LCE oscillators coincide after a period of time, which means that the two bars oscillate synchronously in in-phase mode.Figure 4c,d plot the phase diagrams for the two spring constants, in which θ 1 and θ 2 evolve from the initial disorder into a final attraction domain.Figure 4e further plots the domains of attraction of θ 1 and θ 2 for different spring constants.The calculation shows that the two domains of attraction are identical.The results show that the spring constant has no effect on the in-phase mode.spring constant has no effect on its amplitude and frequency.In in-phase mode, the system is equivalent to the single oscillator, as shown in Figure 5a.In in-phase mode, the system is equivalent to the single oscillator.In anti-phase mode, the system is equivalent to the single oscillator with half-length torsion spring.Furthermore, Figure 4f describes the limit cycles of θ 1 and .θ 1 for different spring constants.Similarly, the two limit cycles are also the same.This implies that the spring constant has no effect on the amplitude and frequency of the LCE oscillators.More calculations show that the influence of the spring constant is the same for α > α crit = 0.075.This is because for α > α crit , the LCE oscillators are in in-phase mode, and both the angle difference between the two bars and the moment of the spring are zero.Therefore, the spring constant has no effect on its amplitude and frequency.In in-phase mode, the system is equivalent to the single oscillator, as shown in Figure 5a.In in-phase mode, the system is equivalent to the single oscillator.In anti-phase mode, the system is equivalent to the single oscillator with half-length torsion spring.Furthermore, Figure 6f describes the limit cycles of 1 θ and 1 θ  for different initial conditions.Similarly, the two limit cycles are also identical.This implies that the amplitude and frequency of the LCE oscillators are independent on the initial condition, which Furthermore, Figure 6f describes the limit cycles of θ 1 and .θ 1 for different initial conditions.Similarly, the two limit cycles are also identical.This implies that the amplitude and frequency of the LCE oscillators are independent on the initial condition, which is further validated by more calculations.It is noted that the effect of initial condition on amplitude and frequency is similar to that of the single LCE oscillator [66].Similar to the effect of the contraction coefficient, its amplitude increase with the increase of the light intensity.This is because with the rise of light intensity I, the light-driven deformation of the LCE bar increases and then the center of gravity of the LCE bar changes greatly, which leads to the augmentation of the driving moment and the amplitude of the self-excited oscillation.for different damping coefficients.It is found that its amplitude decreases with the increase of the damping coefficient, which is consistent with the physical intuition [70].As the β increases, the greater the energy consumed and the smaller the amplitude of the oscillation.θ 1 for different damping coefficients.It is found that its amplitude decreases with the increase of the damping coefficient, which is consistent with the physical intuition [70].As the β increases, the greater the energy consumed and the smaller the amplitude of the oscillation.

Effect of the Gravitational Acceleration on the In-Phase Mode
Figure 10 plots the effect of gravitational acceleration g on the in-phase mode for and θ 0 2 = 10 • .Figure 10a plots the domains of attraction of θ 1 and θ 2 for different gravitational accelerations.Furthermore, Figure 10b describes the limit cycles of θ 1 and .θ 1 for different gravitational accelerations.The results show that the amplitude increases with the increase of gravitational acceleration.According to the physical meaning of g = T 0 / l 0 /g 2 , the larger g is, the slower the cis-to-trans conversion is and therefore, the smaller the total moment M D 1 , M D 2 is, as shown in Equation ( 7), and the larger the amplitude are.

Anti-Phase Synchronization Mode
In this section, we further study the mechanism of self-oscillation in anti-phase mode and the influences of dimensionless parameters on the anti-phase synchronous mode in details.

Mechanisms of the Self-Excited Oscillation in Anti-Phase Mode
To investigate the mechanism of the self-excited oscillation in anti-phase mode of the two LCE oscillators under steady illumination, Figure 11 presents time histories of some key physical quantities of the anti-phase mode in Figure 2a,c.Figure 11a plots the time histories of 1  and 2  , which show that the two LCE bars oscillate periodically in anti- phase mode.Figure 11b plots time histories of the number fractions of cis-isomers in the two LCE bars.Similarly, the number fractions of cis-isomers also increase in the light zone while decreasing in the dark zone.Therefore, the contraction strains in the two LCE bars increase in the light zone while decreasing in the dark zone, as shown in Figure 11c.In Figure 11d, the moment of the torsion spring on the two bars changes periodically, because the angle difference between the two bars in anti-phase mode varies periodically.In Figure 11e, the driving moments of the two oscillators in anti-phase mode also fluctuates periodically.Figure 11f delineates the dependence of the driving moment on In Figure 11f, the area surrounded by the closed curve represents the net work performed

Anti-Phase Synchronization Mode
In this section, we further study the mechanism of self-oscillation in anti-phase mode and the influences of dimensionless parameters on the anti-phase synchronous mode in details.

Mechanisms of the Self-Excited Oscillation in Anti-Phase Mode
To investigate the mechanism of the self-excited oscillation in anti-phase mode of the two LCE oscillators under steady illumination, Figure 11 presents time histories of some key physical quantities of the anti-phase mode in Figure 2a,c.Figure 11a plots the time histories of θ 1 and θ 2 , which show that the two LCE bars oscillate periodically in anti-phase mode.Figure 11b plots time histories of the number fractions of cis-isomers in the two LCE bars.Similarly, the number fractions of cis-isomers also increase in the light zone while decreasing in the dark zone.Therefore, the contraction strains in the two LCE bars increase in the light zone while decreasing in the dark zone, as shown in Figure 11c.In Figure 11d, the moment of the torsion spring on the two bars changes periodically, because the angle difference between the two bars in anti-phase mode varies periodically.In Figure 11e, the driving moments of the two oscillators in anti-phase mode also fluctuates periodically.Figure 11f delineates the dependence of the driving moment on θ 1 or θ 2 .In Figure 11f, the area surrounded by the closed curve represents the net work performed by the steady illumination during one cycle of the self-excited oscillation, which compensates for the energy loss caused by damping to maintain the oscillation of the LCE oscillators.

Effect of the Spring Constant on the Anti-Phase Mode
and θ 0 2 = 10 • .For α = 0, the phase difference between the two bars after a period of time is a constant value that is generally not 180 • or 0 • , as shown in Figure 12a.Figure 12d plots the phase diagrams for α = 0, in which θ 1 and θ 2 evolve from the initial disorder into a final attraction domain.The domains of attraction of θ 1 and θ 2 and limit cycles of θ 1 and .θ 1 are further plotted in Figure 12g,h, respectively.Through calculation, we find that for α = 0, the final phase difference of two bars depends on initial conditions, while the limit cycles are independent on the initial conditions.It can be understood that for α = 0, the system is equivalent to the single oscillator [66], and the phase difference is determined by their independent self-excited oscillations.
by the steady illumination during one cycle of the self-excited oscillation, which compensates for the energy loss caused by damping to maintain the oscillation of the LCE oscillators.Through calculation, we find that for 0 = α , the final phase difference of two bars depends on initial conditions, while the limit cycles are independent on the initial conditions.It can be understood that for 0 = α , the system is equivalent to the single oscillator [66], and the phase difference is determined by their independent self-excited oscillations.In Figure 12b,c, for α = 0.03 and α = 0.05, the self-excited oscillations of the two LCE oscillators have a phase difference of half a cycle after a period of time, which means that the two bars oscillate synchronously in anti-phase mode.Figure 11e,f plot the phase diagrams for α = 0.03 and α = 0.05, in which θ 1 and θ 2 evolve from the initial disorder into a final attraction domain.Furthermore, Figure 12g,h plot the domains of attraction of θ 1 and θ 2 and the limit cycles of θ 1 and .θ 1 .Obviously, the domains of attraction and the limit cycles are also different.It can be seen that its amplitude decreases with the increase of the spring constant.Careful calculation shows that the period also decreases with the increase of the spring constant.In fact, the system in anti-phase mode for α < α crit is equivalent to a single self-excited oscillator constrained by a fixed torsion spring with half the original length, as shown in Figure 5b.The greater the spring constant, the smaller the amplitude and the period.This result is consistent with the physical intuition [70].mode.Figure 13c,d plot the phase diagrams for the two initial conditions, in which θ 1 and θ 2 evolve from the initial disorder into a final attraction domain.Figure 13e further plots the domains of attraction of θ 1 and θ 2 for different initial conditions.The calculation shows that the two domains of attraction are the same.The results also show that the initial condition has no effect on the synchronous mode.θ and 2 θ and the limit cycles of 1 θ and 1 θ  .Obviously, the domains of attraction and the limit cycles are also different.It can be seen that its amplitude decreases with the increase of the spring constant.Careful calculation shows that the period also decreases with the increase of the spring constant.In fact, the system in anti-phase mode for crit α α < is equivalent to a single self-excited oscillator constrained by a fixed torsion spring with half the original length, as shown in Figure 5b.The greater the spring constant, Furthermore, Figure 13f describes the limit cycles of θ 1 and .θ 1 for different initial conditions.Similarly, the two limit cycles are also identical.This implies that the amplitude and frequency of the LCE oscillators are independent on the initial condition.More calculations also show that the influence of initial conditions is the same.It is noted that the effect of initial condition on amplitude and frequency of the anti-phase mode is similar to that of the single LCE oscillator [66].

Effect of Initial Conditions on the Anti-Phase Mode
Furthermore, Figure 13f describes the limit cycles of 1  and 1   for different initial conditions.Similarly, the two limit cycles are also identical.This implies that the amplitude and frequency of the LCE oscillators are independent on the initial condition.More calculations also show that the influence of initial conditions is the same.It is noted that the effect of initial condition on amplitude and frequency of the anti-phase mode is similar to that of the single LCE oscillator [66].14a plots the domains of attraction of θ 1 and θ 2 for different contraction coefficients.Furthermore, Figure 14b describes the limit cycles of θ 1 and .θ 1 for different contraction coefficients.Similar to the effect of contraction coefficient on the in-phase mode, the amplitude in anti-phase mode increases with the increase of the contraction coefficient.With the rise of contraction coefficient C 0 , the strain of the LCE bar increases and then the center of gravity of the LCE bar changes greatly.Therefore, the driving moment increases and the amplitude of the self-excited oscillation ascends.  for different damping coefficients.Similarly, its amplitude also decreases with the rise of the damping coefficient.This result can be understood as follows: the more  increases, the greater the energy consumed and the smaller the amplitude.for different damping coefficients.θ  for different gravitational accelerations.Similarly, its amplitude also increases with the rise of the gravitational acceleration.This result can be understood as follows: the more g increases, the slower the cis-to-trans conversion, the smaller the total θ 1 for different gravitational accelerations.Similarly, its amplitude also increases with the rise of the gravitational acceleration.This result can be understood as follows: the more g increases, the slower the cis-to-trans conversion, the smaller the total moment M D 1 , M D 2 in Equation ( 7), and thus the smaller the amplitude.

Conclusions
The synchronization and group behaviors of self-excited coupled oscillators are common in nature and deserve to be explored, for self-excited motions have the advantages of actively collecting energy from the environment, being autonomous, making equipment portable, and so on.Based on self-excited oscillators composed of LCE bars, the synchronization of two identical self-excited oscillators is theoretically investigated in this paper.By combing the well-established dynamic LCE model, a theoretical model for the self-excited motion of the two self-excited oscillators connected by a torsion spring and synchronization of the self-excited motion is numerically calculated by MATLAB software.It is found that self-excited oscillation of the system has two synchronization modes: in-phase mode and anti-phase mode.By plotting the time histories of various quantities, we elucidate the mechanism of self-excited oscillation and the two synchronization modes.
The numerical calculations show that the synchronization mode is mainly determined by the interaction between the two identical self-excited oscillators.For strong interactions with the spring constant, the system always develops into in-phase synchronization mode, and both its amplitude and period are independent on the spring constant.In this case, the system is equivalent to two identical self-excited oscillators without interaction.For weak interaction with a small spring constant, the system will evolve into anti-

Conclusions
The synchronization and group behaviors of self-excited coupled oscillators are common in nature and deserve to be explored, for self-excited motions have the advantages of actively collecting energy from the environment, being autonomous, making equipment portable, and so on.Based on self-excited oscillators composed of LCE bars, the synchronization of two identical self-excited oscillators is theoretically investigated in this paper.By combing the well-established dynamic LCE model, a theoretical model for the self-excited motion of the two self-excited oscillators connected by a torsion spring and synchronization of the self-excited motion is numerically calculated by MATLAB software.It is found that self-excited oscillation of the system has two synchronization modes: in-phase mode and anti-phase mode.By plotting the time histories of various quantities, we elucidate the mechanism of self-excited oscillation and the two synchronization modes.
The numerical calculations show that the synchronization mode is mainly determined by the interaction between the two identical self-excited oscillators.For strong interactions with the spring constant, the system always develops into in-phase synchronization mode, and both its amplitude and period are independent on the spring constant.In this case, the system is equivalent to two identical self-excited oscillators without interaction.For weak interaction with a small spring constant, the system will evolve into anti-phase synchronization mode, and its amplitude decreases with the increase of the interaction.In this case, it is equivalent to single self-excited oscillator constrained by fixed torsion spring with half original length.
Furthermore, the effects of initial conditions, contraction coefficient, light intensity and damping coefficient on the two synchronization modes of the self-excited oscillation are investigated systematically, by plotting their domain of attraction of the system.The initial condition generally does not affect the synchronization mode or its amplitude.The amplitude of the self-oscillation always increases with the increasing contraction coefficient, gravitational acceleration, and light intensity, while decreasing with the increasing damping coefficient.
In the future, the proposed design in this work is worth experimental implementation to validate the predicted phenomena.It is recommended that the LCE sample with a large contraction coefficient and small thickness is fabricated and a high light intensity is used in experiments.This study will deepen people's understanding of the synchronization of self-excited oscillators and provide promising applications in energy acquisition, power generation, monitoring, soft robotics, medical equipment, and micro/nanodevices.Meanwhile, the theoretical framework could be extended to scenarios involving large-scale synchronization of systems with numerous interacting oscillators and suggests possible avenues for future research in this direction.

Figure
Figure 1a sketches the dynamic model of two identical LCE oscillators connected by a torsion spring under steady illumination.O 1 A 1 and O 2 A 2 are light-responsive LCE bars, while O 1 B 1 and O 2 B 2 are passive material bars.The four bars have the same mass m.Bars A 1 O 1 B 1 and A 2 O 2 B 2 can rotate around the z axis.The original length of O 1 A 1 and O 2 A 2 is l 0 and the length of O 1 B 1 and O 2 B 2 is kl 0 , with k being the length ratio.The connecting torsion spring between the two bars allows them to apply corresponding torque to one

Figure 1 .
Figure 1.Schematics of dynamic model of two identical LCE oscillators connected by a torsion spring under steady illumination.The LCE oscillators are made up of a light-responsive LCE bar and a passive material bar.The region from d θ to u θ is steadily illuminated.(a) Oblique view.(b) Side view.

Figure 1 .
Figure 1.Schematics of dynamic model of two identical LCE oscillators connected by a torsion spring under steady illumination.The LCE oscillators are made up of a light-responsive LCE bar and a passive material bar.The region from θ d to θ u is steadily illuminated.(a) Oblique view.(b) Side view.

Polymers 2023 , 25 Figure 2 .
Figure 2. Two types of synchronization modes of the light-powered self-excited LCE oscillators.The parameters are 33 .0 = I , 4 .0 0 = C , 8 .9 = g , 2 .9 = k , delineates the dependence of the driving moment on 1 θ or 2 θ .In Figure, 3f, the area surrounded by the closed curve represents the net work performed by the steady illumination during one cycle of the selfexcited oscillation, which compensates for the energy dissipation caused by the damping

Figure 4 .
Figure 4. Effect of spring constant on the in-phase mode.The parameters are 33 .0 = I , 4 .0 0 = C ,

Figure 5 .
Figure 5. Equivalent systems of (a) in-phase synchronization mode, and (b) anti-phase synchronization mode.In in-phase mode, the system is equivalent to the single oscillator.In anti-phase mode, the system is equivalent to the single oscillator with half-length torsion spring.

Figure 4 .
Figure 4. Effect of spring constant on the in-phase mode.The parameters are 33 .0 = I , 4 .0 0 = C ,

Figure 5 .
Figure 5. Equivalent systems of (a) in-phase synchronization mode, and (b) anti-phase synchronization mode.In in-phase mode, the system is equivalent to the single oscillator.In anti-phase mode, the system is equivalent to the single oscillator with half-length torsion spring.

Figure 5 . 25 4. 3 ..
Figure 5. Equivalent systems of (a) in-phase synchronization mode, and (b) anti-phase synchronization mode.In in-phase mode, the system is equivalent to the single oscillator.In anti-phase mode, the system is equivalent to the single oscillator with half-length torsion spring.

Figure 8 .
Figure 8 plots the effect of light intensity I on the in-phase mode for 4 .0 0 = C , 1 .0 = α

Polymers 2023 ,
15,  x FOR PEER REVIEW 13 of 25 sity I , the light-driven deformation of the LCE bar increases and then the center of grav- ity of the LCE bar changes greatly, which leads to the augmentation of the driving moment and the amplitude of the self-excited oscillation.

Figure 9 .
Figure 9 plots the effect of damping coefficient β on the in-phase mode for 33 .0 = I

cles of 1  and 1 
 for different contraction coefficients.Similar to the effect of contraction coefficient on the in-phase mode, the amplitude in anti-phase mode increases with the increase of the contraction coefficient.With the rise of contraction coefficient 0 C , the strain of the LCE bar increases and then the center of gravity of the LCE bar changes greatly.Therefore, the driving moment increases and the amplitude of the self-excited oscillation ascends.

Figure 15 Figure 14 .. 1  and 2  1  and 1 
Figure 15 plots the effect of light intensity I on the anti-phase mode for 4 .0 0 = C ,

Table 1 .
Material properties and geometric parameters.