<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Comput. Neurosci.</journal-id>
<journal-title>Frontiers in Computational Neuroscience</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Comput. Neurosci.</abbrev-journal-title>
<issn pub-type="epub">1662-5188</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fncom.2024.1347748</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Neuroscience</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Noise-induced synchrony of two-neuron motifs with asymmetric noise and uneven coupling</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Jagdev</surname> <given-names>Gurpreet</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2560232/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Yu</surname> <given-names>Na</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/754434/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Mathematics, Toronto Metropolitan University</institution>, <addr-line>Toronto, ON</addr-line>, <country>Canada</country></aff>
<aff id="aff2"><sup>2</sup><institution>Institute of Biomedical Engineering, Science and Technology (iBEST), Unity Health Toronto, and Toronto Metropolitan University</institution>, <addr-line>Toronto, ON</addr-line>, <country>Canada</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Lubica Benuskova, Comenius University, Slovakia</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Matja&#x0017E; Perc, University of Maribor, Slovenia</p>
<p>Ergin Yilmaz, Bulent Ecevit University, T&#x000FC;rkiye</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Na Yu <email>nayu&#x00040;torontomu.ca</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>23</day>
<month>02</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>18</volume>
<elocation-id>1347748</elocation-id>
<history>
<date date-type="received">
<day>01</day>
<month>12</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>06</day>
<month>02</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2024 Jagdev and Yu.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Jagdev and Yu</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>Synchronous dynamics play a pivotal role in various cognitive processes. Previous studies extensively investigate noise-induced synchrony in coupled neural oscillators, with a focus on scenarios featuring uniform noise and equal coupling strengths between neurons. However, real-world or experimental settings frequently exhibit heterogeneity, including deviations from uniformity in coupling and noise patterns. This study investigates noise-induced synchrony in a pair of coupled excitable neurons operating in a heterogeneous environment, where both noise intensity and coupling strength can vary independently. Each neuron is an excitable oscillator, represented by the normal form of Hopf bifurcation (HB). In the absence of stimulus, these neurons remain quiescent but can be triggered by perturbations, such as noise. Typically, noise and coupling exert opposing influences on neural dynamics, with noise diminishing coherence and coupling promoting synchrony. Our results illustrate the ability of asymmetric noise to induce synchronization in such coupled neural oscillators, with synchronization becoming increasingly pronounced as the system approaches the excitation threshold (i.e., HB). Additionally, we find that uneven coupling strengths and noise asymmetries are factors that can promote in-phase synchrony. Notably, we identify an optimal synchronization state when the absolute difference in coupling strengths is maximized, regardless of the specific coupling strengths chosen. Furthermore, we establish a robust relationship between coupling asymmetry and the noise intensity required to maximize synchronization. Specifically, when one oscillator (receiver neuron) receives a strong input from the other oscillator (source neuron) and the source neuron receives significantly weaker or no input from the receiver neuron, synchrony is maximized when the noise applied to the receiver neuron is much weaker than that applied to the source neuron. These findings reveal the significant connection between uneven coupling and asymmetric noise in coupled neuronal oscillators, shedding light on the enhanced propensity for in-phase synchronization in two-neuron motifs with one-way connections compared to those with two-way connections. This research contributes to a deeper understanding of the functional roles of network motifs that may serve within neuronal dynamics.</p></abstract>
<kwd-group>
<kwd>network motifs</kwd>
<kwd>coupled oscillators</kwd>
<kwd>synchrony</kwd>
<kwd>asymmetric noise</kwd>
<kwd>heterogeneity</kwd>
</kwd-group>
<counts>
<fig-count count="9"/>
<table-count count="0"/>
<equation-count count="9"/>
<ref-count count="35"/>
<page-count count="11"/>
<word-count count="8105"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>1 Introduction</title>
<p>The synchronization of neural oscillations is acknowledged for its role in facilitating communication between neurons (Varela et al., <xref ref-type="bibr" rid="B30">2001</xref>; Fell and Axmacher, <xref ref-type="bibr" rid="B6">2011</xref>), and is essential for higher order cognitive processes including memory formation, motor coordination, and sensory information processing (Varela et al., <xref ref-type="bibr" rid="B30">2001</xref>; Womelsdorf and Fries, <xref ref-type="bibr" rid="B31">2007</xref>; Kawasaki et al., <xref ref-type="bibr" rid="B10">2018</xref>). Noise, prevalent in real neural networks, plays a crucial role in giving rise to emergent dynamics that may serve important physiological functions (Varela et al., <xref ref-type="bibr" rid="B30">2001</xref>; Fell and Axmacher, <xref ref-type="bibr" rid="B6">2011</xref>). One intriguing constructive effect of noise is noise-induced synchrony, wherein noise alone can induce coordinated and synchronized dynamics. This phenomenon has been widely explored, especially in neural networks, as shown by studies such as (Perc and Marhl, <xref ref-type="bibr" rid="B18">2006</xref>; Perc, <xref ref-type="bibr" rid="B17">2007</xref>; Touboul et al., <xref ref-type="bibr" rid="B28">2020</xref>). Network motifs are frequently recurring structural patterns in neuronal networks, commonly regarded as the fundamental building blocks of complex networks (Milo et al., <xref ref-type="bibr" rid="B12">2002</xref>; Reigl et al., <xref ref-type="bibr" rid="B21">2004</xref>). Among network motifs, those featuring two neurons are notably more prevalent than other multi-neuron motifs (Reigl et al., <xref ref-type="bibr" rid="B21">2004</xref>; Song et al., <xref ref-type="bibr" rid="B26">2005</xref>), making them central to shaping the collective behavior of the network. Consequently, this study focuses on the noise-induced synchrony of the most over-represented motifs: two-neuron motifs.</p>
<p>In the context of two-neuron motifs or two coupled oscillators, prior research has extensively studied the noise-induced synchrony under homogeneous configurations such as common noise (where each neuron encounters identical noise), symmetric noise where each neuron has independent noise but with equal noise intensity), and couplings of exactly equal strengths between neurons. These investigations are well-documented in reviews and books (Rosenblum et al., <xref ref-type="bibr" rid="B22">2001</xref>; Pikovsky, <xref ref-type="bibr" rid="B19">2002</xref>; Boccaletti et al., <xref ref-type="bibr" rid="B3">2006</xref>). In real-world or experimental scenarios, however, coupling and noise patterns often deviate from uniformity (Song et al., <xref ref-type="bibr" rid="B26">2005</xref>; Morgan and Soltesz, <xref ref-type="bibr" rid="B13">2008</xref>). Our understanding of the broader implications associated with the effects of asymmetric noise (where each neuron receives independent noise sources with different noise intensities) and uneven coupling (e.g., unequal coupling strength) is limited, as only a few studies have explored the influence of heterogeneity in noise and/or coupling on the dynamics of two coupled neuronal oscillators. For example, the interplay between uneven coupling strengths and symmetric noise in two coupled oscillators has been shown to promote synchronization (Blasius, <xref ref-type="bibr" rid="B2">2005</xref>) and enhance the transmission of sub-threshold external signals (Masoliver and Masoller, <xref ref-type="bibr" rid="B11">2018</xref>). In a pair of coupled oscillators with uneven coupling and asymmetric additive noise (Amro et al., <xref ref-type="bibr" rid="B1">2015</xref>) finds that the phase coherence of one oscillator is a non-monotonic function of the additive noise applied to the other oscillator: as the phase coherence of one oscillator decreases, the other increases.</p>
<p>To address these less explored aspects, this paper investigates the synchronous dynamics of a pair of excitable neurons in the presence of two sources of heterogeneity: asymmetric noise and uneven coupling. Here the excitable neurons are modeled by the normal form of a Hopf bifurcation (HB), a deterministic framework underpinning critical transitions between quiescent and oscillatory states in complex systems. Therefore our findings provide a versatile framework for illustrating a diverse array of dynamical patterns near a HB within a heterogeneous configuration. Given that the two-neuron motif constitutes a fundamental element in a neuronal network, this study contributes valuable insights into the generation of diverse network behaviors in the context of heterogeneity in both noise and coupling. The structure of the remainder of this paper is as follows: Section 2 provides an introduction to the mathematical model and the methodologies used. In Section 3.1, we present the bifurcation diagram of two deterministic oscillators and Section 3.2 introduces noise-induced oscillations. Section 3.3 focuses on demonstrating noise-induced synchrony and how asymmetric noise affects this synchrony. In Section 3.4, we study the effects of the bifurcation parameter on synchronization, and in Section 3.5, we examine the dual effects of uneven coupling and asymmetric noise on synchrony. Finally, Section 4 offers a summary and engages in a discussion of the findings.</p>
</sec>
<sec sec-type="materials and methods" id="s2">
<title>2 Materials and methods</title>
<p>Synchrony in neural oscillators is not solely determined by their phase; rather, the statistical prevalence of amplitude dynamics influencing synchronization is noteworthy (Gambuzza et al., <xref ref-type="bibr" rid="B8">2016</xref>). This phenomenon, referred to as amplitude-sensitive synchrony, is a distinctive characteristic of oscillators located near an HB. Hence, our research employs the &#x003BB;&#x02212;&#x003C9; system, recognized as a minimal model (or normal form), capable of capturing both amplitude and phase dynamics of an oscillator in the proximity of a HB. Therefore, the insights gained from this study have broader applicability to other dynamical systems situated in the vicinity of a HB.</p>
<sec>
<title>2.1 Model</title>
<p>We examine a duo of coupled &#x003BB;&#x02212;&#x003C9; oscillators, selected with specific parameters that position the model near a supercritical HB. In the absence of noise, these oscillators remain quiescent, but they become excited upon the introduction of an intrinsic noise stimulus. Their coupling strength is uneven, and they are subjected to asymmetric intensities of additive noise. Moreover, both oscillators are represented by the set of stochastic differential equations (SDEs).</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mi>&#x003BB;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mi>&#x003C9;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>&#x003B7;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E2"><label>(2)</label><mml:math id="M2"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mi>&#x003C9;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:mi>&#x003BB;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E3"><label>(3)</label><mml:math id="M3"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msubsup><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:msubsup><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>i, j</italic> &#x0003D; 1, 2. <inline-formula><mml:math id="M4"><mml:msub><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:msubsup><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mrow></mml:msqrt></mml:math></inline-formula> represents the amplitude of the <italic>i</italic>th oscillator. <inline-formula><mml:math id="M5"><mml:mi>&#x003BB;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003BB;</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:mi>&#x003B1;</mml:mi><mml:msubsup><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:mi>&#x003B3;</mml:mi><mml:msubsup><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msubsup></mml:math></inline-formula> controls the increment and decrement of the amplitude of the <italic>i</italic>th oscillator. &#x003BB;<sub>0</sub> is the control parameter and a HB occurs at &#x003BB;<sub>0</sub> &#x0003D; 0 (see <xref ref-type="fig" rid="F1">Figure 1A</xref>). <inline-formula><mml:math id="M6"><mml:mi>&#x003C9;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C9;</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C9;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msubsup><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:math></inline-formula> determines the increment and decrement of the frequency of the <italic>i</italic>th oscillator. We consider a supercritical HB by setting the parameter values as &#x003B1; &#x0003D; &#x02212;0.2, &#x003B3; &#x0003D; &#x02212;0.2, &#x003C9;<sub>0</sub> &#x0003D; 2, and &#x003C9;<sub>1</sub> &#x0003D; 0. <italic>d</italic><sub><italic>i</italic></sub>(<italic>x</italic><sub><italic>j</italic></sub> &#x02212; <italic>x</italic><sub><italic>i</italic></sub>) and <italic>d</italic><sub><italic>i</italic></sub>(<italic>y</italic><sub><italic>j</italic></sub> &#x02212; <italic>y</italic><sub><italic>i</italic></sub>) represent diffusive coupling between oscillators <italic>i</italic> and <italic>j</italic> with oscillator specific coupling strength, <italic>d</italic><sub><italic>i</italic></sub>, <italic>i</italic> &#x0003D; 1, 2. &#x003B4;<sub><italic>i</italic></sub><italic>d&#x003B7;</italic><sub><italic>i</italic></sub>(<italic>t</italic>) represents an independent intrinsic noise applied to <italic>x</italic><sub><italic>i</italic></sub> (i.e., the noise is unique to each oscillator) where &#x003B7;<sub><italic>i</italic></sub>(<italic>t</italic>) is a Wiener process with zero mean and unity variance and &#x003B4;<sub><italic>i</italic></sub> is the noise intensity. We restrict our attention to the excitatory coupling in the range 0.01 &#x02264; <italic>d</italic><sub>1</sub>, <italic>d</italic><sub>2</sub> &#x02264; 0.3 (as in Yu et al., <xref ref-type="bibr" rid="B35">2008</xref>).</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Bifurcation diagrams of the deterministic system (&#x003B4;<sub>1</sub> &#x0003D; &#x003B4;<sub>2</sub> &#x0003D; 0) vs. control parameter &#x003BB;<sub>0</sub>. <bold>(A)</bold> single/uncoupled oscillator (<italic>x</italic> &#x0003D; <italic>x</italic><sub>1</sub> &#x0003D; <italic>x</italic><sub>2</sub>) with <italic>d</italic><sub>1</sub> &#x0003D; <italic>d</italic><sub>2</sub> &#x0003D; 0, a single Hopf bifurcation (HB) occurs at &#x003BB; &#x0003D; 0. <bold>(B)</bold> symmetrically coupled oscillators (<italic>x</italic><sub>1</sub>, <italic>x</italic><sub>2</sub>) with <italic>d</italic><sub>1</sub> &#x0003D; <italic>d</italic><sub>2</sub> &#x0003D; 0.05, <bold>(C, D)</bold> asymmetrically coupled oscillators (<italic>x</italic><sub>1</sub> and <italic>x</italic><sub>2</sub>, respectively) with <italic>d</italic><sub>1</sub> &#x0003D; 0.1 and <italic>d</italic><sub>2</sub> &#x0003D; 0.01. Two HB points, HB<sub>1</sub> and HB<sub>2</sub>, occur and are indicated by blue and orange dots, respectively. The value of &#x003BB;<sub>0</sub> corresponding to HB<sub>1</sub> is zero in <bold>(B&#x02013;D)</bold>. Stable solutions are marked by solid blue lines and solid black lines, and unstable solutions are marked by dashed blue lines and dashed black lines. Other parameters are: &#x003B1; &#x0003D; &#x02212;0.2; &#x003B3; &#x0003D; &#x02212;0.2; &#x003C9;<sub>0</sub> &#x0003D; 2; and &#x003C9;<sub>1</sub> &#x0003D; 0.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-18-1347748-g0001.tif"/>
</fig>
</sec>
<sec>
<title>2.2 Methods</title>
<p>To study the interplay of noise and coupling and its effect on the synchronization of our model we analyze the synchrony of both oscillators when subject to the additive noise, &#x003B4;<sub><italic>i</italic></sub><italic>d&#x003B7;</italic><sub><italic>i</italic></sub>, <italic>i</italic> &#x0003D; 1, 2. Since the oscillators rotate about the fixed point (<italic>x</italic><sub><italic>i</italic></sub>, <italic>y</italic><sub><italic>i</italic></sub>) &#x0003D; (0, 0), <italic>i</italic> &#x0003D; 1, 2, when driven by noise, the time-dependent phase of each oscillator is taken to be the natural phase (Rosenblum et al., <xref ref-type="bibr" rid="B22">2001</xref>),</p>
<disp-formula id="E4"><label>(4)</label><mml:math id="M7"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mo class="qopname">arctan</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>/</mml:mo><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p><italic>i</italic> &#x0003D; 1, 2. In the classical treatment of phase analysis, synchrony measures are often based on the distribution of the phase difference, &#x00394;&#x003D5;<sub><italic>n,m</italic></sub> &#x0003D; <italic>n&#x003D5;</italic><sub>2</sub> &#x02212; <italic>m&#x003D5;</italic><sub>1</sub>, where <italic>n, m</italic> &#x02208; &#x02115; characterize the order of locking (Rosenblum et al., <xref ref-type="bibr" rid="B22">2001</xref>, <xref ref-type="bibr" rid="B23">2022</xref>). However, in the presence of noise, the phase of the oscillators can exhibit random jumps of &#x000B1;2&#x003C0;, called phase slips, which can cause the phase difference, &#x00394;&#x003D5;<sub><italic>n,m</italic></sub>, to compound errors, and lead to erroneous results. Therefore, instead of considering the natural phase in <xref ref-type="disp-formula" rid="E4">Equation 4</xref>, we consider the cyclic relative phase (Mormann et al., <xref ref-type="bibr" rid="B14">2000</xref>; Rosenblum et al., <xref ref-type="bibr" rid="B22">2001</xref>),</p>
<disp-formula id="E5"><label>(5)</label><mml:math id="M8"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>&#x003C6;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mtext class="textrm" mathvariant="normal">mod</mml:mtext><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>which is the natural phase wrapped over the unit circle. This procedure ensures that errors in &#x00394;&#x003C6;<sub><italic>n,m</italic></sub> caused by phase slips do not compound which leads to more stable numerical results. Furthermore, for simplicity, we consider only 1 &#x02212; 1 synchronization: &#x00394;&#x003C6; &#x0003D; &#x00394;&#x003C6;<sub>1,1</sub>.</p>
<p>The bifurcation diagrams (i.e., <xref ref-type="fig" rid="F1">Figures 1</xref>, <xref ref-type="fig" rid="F2">2</xref>) are generated using XPPAUT software (Ermentrout, <xref ref-type="bibr" rid="B5">2012</xref>). All further analysis (i.e., <xref ref-type="fig" rid="F3">Figures 3</xref>&#x02013;<bold>9</bold>) is conducted using MATLAB. To simulate the SDEs in <xref ref-type="disp-formula" rid="E1">Equations 1</xref>, <xref ref-type="disp-formula" rid="E2">2</xref> we use the Euler-Maruyama method with time-step <italic>dt</italic> &#x0003D; 0.01 and arbitrary random initial conditions <inline-formula><mml:math id="M9"><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0007E;</mml:mo><mml:mi>N</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:mn>00</mml:mn><mml:msup><mml:mrow><mml:mn>8</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula>, <italic>i</italic> &#x0003D; 1, 2. Due to the nature of white noise, noise-induced oscillations contain high-frequency fluctuations at very low or very high noise intensities, which poses difficulties in the computation of phase. Hence, to achieve more consistent numerical results a low-pass filter is applied to remove high-frequency fluctuations. The signal-to-noise ratio, &#x003B2;, and synchronization measures |&#x00394;&#x003C6;|, <italic>R</italic>, and &#x003C1; in <xref ref-type="disp-formula" rid="E6">Equations 6</xref>&#x02013;<xref ref-type="disp-formula" rid="E9">9</xref>, respectively, are averaged over <italic>N</italic> &#x0003D; 200 trials. The XPPAUT and MATLAB source code can be found at: <ext-link ext-link-type="uri" xlink:href="https://github.com/TMUcode/CoupleNeurons">https://github.com/TMUcode/CoupleNeurons</ext-link>.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Two-parameter bifurcation diagrams (coupling strength <italic>d</italic><sub><italic>i</italic></sub> vs. &#x003BB;<sub>0</sub>) for three coupling cases: <bold>(A)</bold> symmetric coupling strengths, <italic>d</italic><sub>1</sub> &#x0003D; <italic>d</italic><sub>2</sub>, <bold>(B)</bold> asymmetric coupling strengths with varying <italic>d</italic><sub>1</sub> and fixed <italic>d</italic><sub>2</sub> &#x0003D; 0.05, and <bold>(C)</bold> asymmetric coupling strengths with varying <italic>d</italic><sub>2</sub> and fixed <italic>d</italic><sub>1</sub> &#x0003D; 0.05. The branches labeled HB<sub>1</sub> (solid blue) and HB<sub>2</sub> (dashed blue) represent the same two distinct HB points in <xref ref-type="fig" rid="F1">Figure 1</xref>. The value of &#x003BB;<sub>0</sub> corresponding to HB<sub>1</sub> is consistently zero across three cases, aligning with the observation in <xref ref-type="fig" rid="F1">Figure 1</xref>. Conversely, the value of &#x003BB;<sub>0</sub> corresponding to HB<sub>2</sub> linearly increases with <italic>d</italic><sub><italic>i</italic></sub>. Other parameters are: &#x003B4;<sub>1</sub> &#x0003D; &#x003B4;<sub>2</sub> &#x0003D; 0; &#x003B1; &#x0003D; &#x02212;0.2; &#x003B3; &#x0003D; &#x02212;0.2; &#x003C9;<sub>0</sub> &#x0003D; 2; and &#x003C9;<sub>1</sub> &#x0003D; 0.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-18-1347748-g0002.tif"/>
</fig>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Time series of <italic>x</italic><sub>1</sub>and<italic>x</italic><sub>2</sub> in the excitable regime (&#x003BB;<sub>0</sub> &#x0003C; 0) with various noise intensities: <bold>(A)</bold> &#x003B4;<sub>1</sub> &#x0003D; &#x003B4;<sub>2</sub> &#x0003D; 0, <bold>(B)</bold> &#x003B4;<sub>1</sub> &#x0003D; &#x003B4;<sub>2</sub> &#x0003D; 0.01, <bold>(C)</bold> &#x003B4;<sub>1</sub> &#x0003D; 0.01 and &#x003B4;<sub>2</sub> &#x0003D; 0.05, and <bold>(D)</bold> &#x003B4;<sub>1</sub> &#x0003D; 0.01 and &#x003B4;<sub>2</sub> &#x0003D; 3. The blue lines represent <italic>x</italic><sub>1</sub> and the orange lines represent <italic>x</italic><sub>2</sub>. <bold>(A)</bold> shows the deterministic system with a stable fixed point at <italic>x</italic><sub><italic>i</italic></sub> &#x0003D; 0, resulting in an overlap between <italic>x</italic><sub>1</sub> and <italic>x</italic><sub>2</sub>. <bold>(B&#x02013;D)</bold> illustrate the noise-induced oscillations, with <bold>(C)</bold> showing higher levels of regularity and synchrony in comparison to <bold>(B, D)</bold>. Other parameters are: &#x003B1; &#x0003D; &#x02212;0.2; &#x003B3; &#x0003D; &#x02212;0.2; &#x003C9;<sub>0</sub> &#x0003D; 2; &#x003C9;<sub>1</sub> &#x0003D; 0; &#x003BB;<sub>0</sub> &#x0003D; &#x02212;0.5; <italic>d</italic><sub>1</sub> &#x0003D; 0.3; and <italic>d</italic><sub>2</sub> &#x0003D; 0.01.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-18-1347748-g0003.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<sec>
<title>3.1 Bifurcation analysis</title>
<p>The dynamics of the deterministic system (&#x003B4;<sub>1</sub> &#x0003D; &#x003B4;<sub>2</sub> &#x0003D; 0) are examined in this section. We consider three cases and compare their bifurcation diagrams: single (i.e., uncoupled) oscillators with <italic>d</italic><sub>1</sub> &#x0003D; <italic>d</italic><sub>2</sub> &#x0003D; 0 (<xref ref-type="fig" rid="F1">Figure 1A</xref>), two symmetrically coupled oscillators with <italic>d</italic><sub>1</sub> &#x0003D; <italic>d</italic><sub>2</sub> &#x02260; 0 (<xref ref-type="fig" rid="F1">Figure 1B</xref>), and two asymmetrically coupled oscillators with <italic>d</italic><sub>1</sub> &#x02260; <italic>d</italic><sub>2</sub> (<xref ref-type="fig" rid="F1">Figures 1C</xref>, <xref ref-type="fig" rid="F1">D</xref>). When &#x003BB;<sub>0</sub> &#x0003C; 0, both oscillators in all three cases are quiescent, they rest at stable fixed points, as shown in <xref ref-type="fig" rid="F1">Figure 1</xref> (solid black line). Conversely, when &#x003BB;<sub>0</sub> &#x0003E; 0, the fixed points are destabilized (dashed black line), and stable periodic orbits emerge (solid blue line) in all three cases. Both isolated (<xref ref-type="fig" rid="F1">Figure 1A</xref>) and coupled (<xref ref-type="fig" rid="F1">Figures 1B</xref>&#x02013;<xref ref-type="fig" rid="F1">D</xref>) oscillators undergo a supercritical HB (denoted as HB<sub>1</sub>) at &#x003BB;<sub>0</sub> &#x0003D; 0. However, when the oscillators are coupled, the system exhibits a second HB (denoted as HB<sub>2</sub>) which leads to unstable periodic orbits (dashed blue line in <xref ref-type="fig" rid="F1">Figures 1B</xref>&#x02013;<xref ref-type="fig" rid="F1">D</xref>).</p>
<p>When the oscillators are symmetrically coupled (e.g., <italic>d</italic><sub>1</sub>, <italic>d</italic><sub>2</sub> &#x0003D; 0.05 in <xref ref-type="fig" rid="F1">Figure 1B</xref>), the amplitudes of the periodic orbits generated by both oscillators are identical for each &#x003BB;<sub>0</sub> value, however, the unstable periodic orbits have an amplitude that is slightly less than the amplitude of the stable periodic orbits. When the oscillators are coupled asymmetrically (e.g., <italic>d</italic><sub>1</sub> &#x0003D; 0.1 and <italic>d</italic><sub>2</sub> &#x0003D; 0.01 in <xref ref-type="fig" rid="F1">Figures 1C</xref>, <xref ref-type="fig" rid="F1">D</xref>), the bifurcation diagram for one oscillator (<xref ref-type="fig" rid="F1">Figure 1C</xref>) is identical to its counterpart in <xref ref-type="fig" rid="F1">Figure 1B</xref>, whereas, the bifurcation diagram of the other oscillator (<xref ref-type="fig" rid="F1">Figure 1D</xref>) shows a single difference: the amplitude of the unstable periodic orbit is near zero, which is much smaller counterpart in <xref ref-type="fig" rid="F1">Figure 1B</xref>.</p>
<p>To further explore the effects of coupling on the deterministic system, we calculate the two-parameter bifurcation diagrams of our model by taking both the coupling strength, <italic>d</italic><sub><italic>i</italic></sub>, and &#x003BB;<sub>0</sub> as control parameters. The trajectory of HB<sub>1</sub> (solid blue line) in the two-parameter bifurcation diagram in <xref ref-type="fig" rid="F2">Figure 2</xref> indicates the transition of the system between the quiescent and oscillatory states in the deterministic regime remains the same (&#x003BB;<sub>0</sub> &#x0003D; 0) no matter the coupling strength. Moreover, the deterministic system is quiescent when &#x003BB;<sub>0</sub> &#x0003C; 0, but can be excited by an external stimulus (e.g., noise) to produce oscillations. Thus, we refer to the left side of HB<sub>1</sub> as the excitable regime, and the right side of HB<sub>1</sub> as the oscillatory regime.</p>
<p>Although HB<sub>2</sub> is not a determinant factor of the excitability in the deterministic regime, we find that it has a relationship with the coupling strength: &#x003BB;<sub>0</sub> &#x02248; <italic>d</italic><sub>1</sub> &#x0002B; <italic>d</italic><sub>2</sub>. Specifically, when the coupling is symmetric (<italic>d</italic><sub>1</sub> &#x0003D; <italic>d</italic><sub>2</sub>), HB<sub>2</sub> occurs at &#x003BB;<sub>0</sub> &#x02248; 2<italic>d</italic><sub><italic>i</italic></sub>, <italic>i</italic> &#x0003D; 1or2, resulting in a line with a slope approximately equal to 1/2, as depicted by the dashed blue line in <xref ref-type="fig" rid="F2">Figure 2A</xref>. For asymmetric couplings, such those obtained by varying <italic>d</italic><sub>1</sub> while fixing <italic>d</italic><sub>2</sub> &#x0003D; 0.05, HB<sub>2</sub> appears as a straight line with a slope approximately equal to 1 (dashed blue line in <xref ref-type="fig" rid="F2">Figure 2B</xref>), since &#x003BB;<sub>0</sub> &#x02248; <italic>d</italic><sub>1</sub> &#x0002B; 0.05. Another example of asymmetric coupling is obtained by varying <italic>d</italic><sub>2</sub> while fixing <italic>d</italic><sub>1</sub> &#x0003D; 0.05, which yields a straight line with the equation &#x003BB;<sub>0</sub> &#x02248; <italic>d</italic><sub>2</sub> &#x0002B; 0.05 (dashed blue line in <xref ref-type="fig" rid="F2">Figure 2C</xref>).</p>
</sec>
<sec>
<title>3.2 Noise-induced oscillations</title>
<p>In this section, we study noise-induced oscillations. The control parameter &#x003BB;<sub>0</sub> is set to be in the excitable regime, for example, &#x003BB;<sub>0</sub> &#x0003D; &#x02212;0.5. As a result, the deterministic system, with &#x003B4;<sub>1</sub> &#x0003D; &#x003B4;<sub>2</sub> &#x0003D; 0, displays damped oscillations that eventually converge to the fixed point (<italic>x</italic><sub><italic>i</italic></sub> &#x0003D; <italic>y</italic><sub><italic>i</italic></sub> &#x0003D; 0, <italic>i</italic> &#x0003D; 1, 2), as illustrated in <xref ref-type="fig" rid="F3">Figure 3A</xref>. However, when the intrinsic noise, &#x003B4;<sub><italic>i</italic></sub><italic>d&#x003B7;</italic><sub><italic>i</italic></sub>, is introduced, it can trigger repeated excursions from the stable fixed point which leads to oscillatory motion; otherwise known as noise-induced oscillations. Examples are presented in <xref ref-type="fig" rid="F3">Figures 3B</xref>&#x02013;<xref ref-type="fig" rid="F3">D</xref> which display sample time series of <italic>x</italic><sub>1</sub> and <italic>x</italic><sub>2</sub> in the presence of the intrinsic noise stimulus with different intensities. When the noise intensity is small and symmetric (e.g. &#x003B4;<sub>1</sub> &#x0003D; &#x003B4;<sub>2</sub> &#x0003D; 0.01 in <xref ref-type="fig" rid="F3">Figure 3B</xref>) there are intermittent periods of phase drift and phase locking; both in-phase and anti-phase. When the noise intensities become asymmetric by increasing one of the noise intensities, for example, &#x003B4;<sub>1</sub> &#x0003D; 0.01 and &#x003B4;<sub>2</sub> &#x0003D; 0.05 as in <xref ref-type="fig" rid="F3">Figure 3C</xref>, the time series for <italic>x</italic><sub>1</sub> (blue line) and <italic>x</italic><sub>2</sub> (orange line) exhibit increased regularity and longer epochs of in-phase locking. Furthermore, if &#x003B4;<sub>2</sub> is increased further, for example, &#x003B4;<sub>2</sub> &#x0003D; 3 as in <xref ref-type="fig" rid="F3">Figure 3D</xref>, the oscillations of <italic>x</italic><sub>2</sub> become less regular and more chaotic. Overall, these examples demonstrate that synchrony can be optimized by tuning just one of the noise intensities to be an intermediate value.</p>
<p>To determine an appropriate range of noise intensities, &#x003B4;<sub>1</sub> and &#x003B4;<sub>2</sub>, we make use of the signal-to-noise ratio measure (Gang et al., <xref ref-type="bibr" rid="B9">1993</xref>; Pikovsky and Kurths, <xref ref-type="bibr" rid="B20">1997</xref>).</p>
<disp-formula id="E6"><label>(6)</label><mml:math id="M10"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>&#x003B2;</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x00394;</mml:mo><mml:mi>&#x003C9;</mml:mi><mml:mo>/</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C9;</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>h</italic><sub><italic>p</italic></sub> and &#x003C9;<sub><italic>p</italic></sub> denote the height and central frequency of the power spectrum density peak of <italic>x</italic>, respectively, and &#x00394;&#x003C9; denotes the width of the power spectrum density peak at half-maximal power, <inline-formula><mml:math id="M11"><mml:msup><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msub><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>. We have computed the relationship between &#x003B2; and &#x003B4; for a single &#x003BB;&#x02212;&#x003C9; oscillator, (<italic>x, y</italic>), and present the results in <xref ref-type="fig" rid="F4">Figure 4</xref>. When &#x003B4; &#x0003C; 0.01, we find &#x003B2; is significantly small, indicating that the noise level is too weak. However, as the level of noise gradually increases, the &#x003B2; curve displays a peak around &#x003B4; &#x0003D; 0.9, indicating that this noise intensity is optimal. For larger values of &#x003B4;, specifically, when &#x003B4; &#x0003E; 0.9, we observe a sharp decrease in &#x003B2;, suggesting that the noise intensity is now overpowering the regularity of oscillations. Therefore, based on these observations, we select the range 0.01 &#x02264; &#x003B4;<sub>1</sub>, &#x003B4;<sub>2</sub> &#x02264; 5 for our analysis.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Signal-to-noise ratio measure &#x003B2; vs. noise intensity &#x003B4; for a single/uncoupled neuron oscillator in the excitable regime (&#x003BB;<sub>0</sub> &#x0003C; 0). Greater regularity in noise-induced oscillations is indicated by larger values of &#x003B2;. The largest &#x003B2; occurs when noise intensity is &#x003B4; &#x0003D; 0.9. Other parameters are: &#x003B1; &#x0003D; &#x02212;0.2; &#x003B3; &#x0003D; &#x02212;0.2; &#x003C9;<sub>0</sub> &#x0003D; 2; &#x003C9;<sub>1</sub> &#x0003D; 0; and &#x003BB;<sub>0</sub> &#x0003D; &#x02212;0.5.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-18-1347748-g0004.tif"/>
</fig>
</sec>
<sec>
<title>3.3 Noise-induced synchronization and asymmetric noise intensity</title>
<p>The results presented in Section 3.2 suggest that the synchrony of our system with asymmetric noise and uneven coupling can be optimized by adjusting the noise intensities, &#x003B4;<sub>1</sub> and &#x003B4;<sub>2</sub>. To quantify the effects of noise on the synchrony of our model more rigorously, we introduce three measures of synchrony. The first measure is the absolute phase difference, denoted |&#x00394;&#x003C6;|, which is defined as</p>
<disp-formula id="E7"><label>(7)</label><mml:math id="M12"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mo>|</mml:mo><mml:mo>&#x00394;</mml:mo><mml:mi>&#x003C6;</mml:mi><mml:mo>|</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>T</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mo>|</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C6;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C6;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>|</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003C6;<sub><italic>i</italic></sub>(<italic>t</italic>) is the phase of the <italic>i</italic>th oscillator. Since our model considers excitatory coupling, we focus on the in-phase dynamics of the oscillators, such that smaller values of |&#x00394;&#x003C6;| correspond to a greater degree of synchrony. The second measure we use is the mean phase coherence, denoted <italic>R</italic>, which is defined as (Mormann et al., <xref ref-type="bibr" rid="B14">2000</xref>; Rosenblum et al., <xref ref-type="bibr" rid="B22">2001</xref>),</p>
<disp-formula id="E8"><label>(8)</label><mml:math id="M13"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>R</mml:mi><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>T</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mo class="qopname">sin</mml:mo><mml:mo>&#x00394;</mml:mo><mml:mi>&#x003C6;</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>&#x0002B;</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>T</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo class="qopname">&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mo class="qopname">cos</mml:mo><mml:mo>&#x00394;</mml:mo><mml:mi>&#x003C6;</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:msqrt></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x00394;&#x003C6; &#x0003D; &#x003C6;<sub>1</sub> &#x02212; &#x003C6;<sub>2</sub>. The value of <italic>R</italic> ranges from 0 to 1, and larger values of <italic>R</italic> indicate a greater degree of synchrony. The third synchronization measure we use is the normalized synchronization index, denoted &#x003C1;, which is defined as (Rosenblum et al., <xref ref-type="bibr" rid="B22">2001</xref>),</p>
<disp-formula id="E9"><label>(9)</label><mml:math id="M14"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>&#x003C1;</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <inline-formula><mml:math id="M15"><mml:mi>S</mml:mi><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>M</mml:mi></mml:mrow></mml:munderover><mml:msub><mml:mrow><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo class="qopname">ln</mml:mo><mml:msub><mml:mrow><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> is the Shannon entropy, <italic>S</italic><sub><italic>max</italic></sub> &#x0003D; ln<italic>M</italic> is the maximum entropy, <italic>p</italic><sub><italic>k</italic></sub> refers to the probability of finding &#x00394;&#x003C6; in the <italic>k</italic>th bin of the histogram, and <italic>M</italic> is the total number of bins. The quantity &#x003C1; is normalized to 0 &#x02264; &#x003C1; &#x02264; 1, and smaller values of &#x003C1; indicate a narrower distribution of &#x00394;&#x003C6; and thus a greater degree of synchrony.</p>
<p>To begin we consider the effect of a single noise intensity on the degree of synchrony between two unevenly coupled oscillators. To achieve this, we set <italic>d</italic><sub>1</sub> &#x02260; <italic>d</italic><sub>2</sub>, fix one noise intensity to be &#x003B4;<sub>1</sub> &#x0003D; 0.05, and vary the other noise intensity, &#x003B4;<sub>2</sub>, from 0.01 to 5 and compute the level of synchrony using the three measures listed above. We present the results in <xref ref-type="fig" rid="F5">Figure 5A</xref>, where the blue, orange, and black curves correspond to |&#x00394;&#x003C6;|, <italic>R</italic>, and &#x003C1;, respectively. For weak &#x003B4;<sub>2</sub>, both &#x003C1; and |&#x00394;&#x003C6;| show a rapid decrease while <italic>R</italic> displays a rapid increase. These consistent results suggest that the level of synchrony is increasing. However, for strong &#x003B4;<sub>2</sub>, we notice an opposite trend, where &#x003C1; and |&#x00394;&#x003C6;| increase while <italic>R</italic> decreases, indicating a decrease in the level of synchrony.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p><bold>(A)</bold> Three measures of phase synchronization: |&#x00394;&#x003C6;| (the absolute value of the phase difference), <italic>R</italic> (the mean phase coherence), and &#x003C1; (the synchronization index) vs. &#x003B4;<sub>2</sub>. The orange, blue, and black curves correspond to <italic>R</italic>, |&#x00394;&#x003C6;|, and &#x003C1;, respectively. The three vertical dashed lines labeled by a, b, and c correspond to &#x003B4;<sub>2</sub> &#x0003D; 0.05, 0.95, and3, respectively. Smaller values of |&#x00394;&#x003C6;| and &#x003C1;, and larger values of <italic>R</italic> indicate greater noise-induced synchrony. All three measures demonstrate consistent synchrony dynamics as &#x003B4;<sub>2</sub> varies, reaching maximum synchrony at &#x003B4;<sub>2</sub> &#x0003D; 0.95. <bold>(B)</bold> The empirical probability density functions of &#x00394;&#x003C6;, with the curves a, b, and c corresponding to noise intensities &#x003B4;<sub>2</sub> &#x0003D; 0.05, 0.95, and 3 as in <bold>(A)</bold>, respectively. The peaks of density functions corresponding to a, b, and c are &#x00394;&#x003C6; &#x0003D; &#x02212;0.12, 0.08, and 0.12, respectively. Curve b with &#x003B4;<sub>2</sub> &#x0003D; 0.95 exhibits the most concentrated distribution of &#x00394;&#x003C6;, indicating optimal noise-induced synchrony, consistent with the finding in <bold>(A)</bold>. Other parameters are: &#x003B4;<sub>1</sub> &#x0003D; 0.05; &#x003B1; &#x0003D; &#x02212;0.2; &#x003B3; &#x0003D; &#x02212;0.2; &#x003C9;<sub>0</sub> &#x0003D; 2; &#x003C9;<sub>1</sub> &#x0003D; 0; &#x003BB;<sub>0</sub> &#x0003D; &#x02212;0.5; <italic>d</italic><sub>2</sub> &#x0003D; 0.01; and <italic>d</italic><sub>1</sub> &#x0003D; 0.3.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-18-1347748-g0005.tif"/>
</fig>
<p>Furthermore, synchrony is maximized at an intermediate noise intensity, which we find to be &#x003B4;<sub>2</sub> &#x0003D; 0.95 for all three measures in <xref ref-type="fig" rid="F5">Figure 5A</xref>. Notice that it is close to the optimal noise intensity of a single oscillator as shown in <xref ref-type="fig" rid="F4">Figure 4</xref>. These changes in synchrony may also be observed from the probability density functions of &#x00394;&#x003C6; shown in <xref ref-type="fig" rid="F5">Figure 5B</xref>. The dotted blue curve (labeled a), solid blue curve (labeled b), and dashed blue curve (labeled c) correspond to noise intensities &#x003B4;<sub>2</sub> &#x0003D; 0.05, 0.95, and3, respectively (i.e., points a, b, and c in <xref ref-type="fig" rid="F5">Figure 5A</xref>). As &#x003B4;<sub>2</sub> increases from weak to intermediate levels, the center of the distribution of &#x00394;&#x003C6; shifts from small negative values toward the point &#x00394;&#x003C6; &#x0003D; 0, which represents in-phase synchronization. Moreover, the distribution of &#x00394;&#x003C6; is centered at 0.08 at the optimal noise intensity &#x003B4;<sub>2</sub> &#x0003D; 0.95 in <xref ref-type="fig" rid="F5">Figure 5A</xref>. Additionally, the width and peak of the distribution become narrower and higher, respectively, indicating a more concentrated distribution of phase differences between the two oscillators (i.e., a higher level of synchrony). As the noise intensity continues to increase, the center of the &#x00394;&#x003C6; distribution remains almost unchanged, but its width and peak increase and decrease, respectively, which is indicative of a lower degree of synchrony.</p>
<p>Next, we investigate how both noise intensities, &#x003B4;<sub>1</sub> and &#x003B4;<sub>2</sub>, affect the level of synchrony between two unevenly coupled oscillators. To systematically examine their influence, we generate heat maps, which represent the level of synchrony using the measures |&#x00394;&#x003C6;| and <italic>R</italic> over the parameter space 0.01 &#x02264; &#x003B4;<sub>1</sub>, &#x003B4;<sub>2</sub> &#x02264; 5, as illustrated in <xref ref-type="fig" rid="F6">Figures 6A</xref>, <xref ref-type="fig" rid="F6">B</xref>, respectively. In <xref ref-type="fig" rid="F6">Figure 6A</xref>, warmer colors indicate larger values of |&#x00394;&#x003C6;| (i.e., lower levels of synchrony), while cooler colors correspond to smaller values of |&#x00394;&#x003C6;| (i.e., higher levels of synchrony). Note, a three-dimensional representation of <xref ref-type="fig" rid="F6">Figure 6A</xref> is also presented in <xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1A</xref>. The heat map of <italic>R</italic> in <xref ref-type="fig" rid="F6">Figure 6B</xref> corroborates the results obtained from <xref ref-type="fig" rid="F6">Figure 6A</xref>, but with a different color scheme; warmer colors indicate higher levels of synchrony (i.e., larger <italic>R</italic> values), while cooler colors denote lower levels of synchrony (i.e., smaller <italic>R</italic> values). As such, we will primarily focus on <xref ref-type="fig" rid="F6">Figure 6A</xref>. The highest level of synchrony (i.e., the absolute minimum of |&#x00394;&#x003C6;|) is observed in the dark blue triangular region [where &#x003B4;<sub>1</sub> &#x02208; [0.01, 0.06], &#x003B4;<sub>2</sub> &#x02208; [0.1, 1], and &#x003B4;<sub>1</sub>/&#x003B4;<sub>2</sub> &#x02248; 0.2]. A local minimum of |&#x00394;&#x003C6;| is also observed in the lower right part of the heat map [where &#x003B4;<sub>1</sub> &#x02208; [0.8, 3] and &#x003B4;<sub>2</sub> &#x02208; [0.01, 0.015]], with a relatively high level of synchrony. Due to uneven coupling strengths (<italic>d</italic><sub>1</sub> = 0.3 and <italic>d</italic><sub>2</sub> = 0.01) used in <xref ref-type="fig" rid="F6">Figure 6</xref>, the locations of the absolute and local minima are not symmetrical. Our numerical results further indicate a relationship between the coupling strength and the optimal noise intensity required to achieve the maximum level of synchrony (i.e., the absolute maximum of |&#x00394;&#x003C6;|). For instance, as shown in <xref ref-type="fig" rid="F6">Figure 6</xref>, oscillator 1 has a larger coupling strength than oscillator 2 (i.e., <italic>d</italic><sub>1</sub> &#x0003E; <italic>d</italic><sub>2</sub>), thus requiring &#x003B4;<sub>1</sub> &#x0003C; &#x003B4;<sub>2</sub> to achieve the largest level of synchrony. In Section 3.5 we study the interplay of coupling and noise intensity and their effects on the synchrony in greater detail.</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p><bold>(A)</bold> Heat map of the mean absolute phase difference |&#x00394;&#x003C6;|, as a function of noise intensities &#x003B4;<sub>1</sub> and &#x003B4;<sub>2</sub>. Colder colors correspond to smaller values of |&#x00394;&#x003C6;| and a decrease in |&#x00394;&#x003C6;| indicates greater synchrony. The maximum synchrony is located within the dark blue triangle on the left side. The three-dimensional visualization of (A) is depicted in <xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1A</xref>. <bold>(B)</bold> Heat map of mean phase coherence, <italic>R</italic>, as a function of &#x003B4;<sub>1</sub> and &#x003B4;<sub>2</sub>. Warmer colors correspond to larger values of <italic>R</italic> and an increase in <italic>R</italic> indicates better noise-induced synchrony. The maximum synchrony is within the yellow triangle on the left side, consistent with the observation in <bold>(A)</bold>. Other parameters are: &#x003B1; &#x0003D; &#x02212;0.2; &#x003B3; &#x0003D; &#x02212;0.2; &#x003C9;<sub>0</sub> &#x0003D; 2; &#x003C9;<sub>1</sub> &#x0003D; 0; &#x003BB;<sub>0</sub> &#x0003D; &#x02212;0.5; <italic>d</italic><sub>1</sub> &#x0003D; 0.3; and <italic>d</italic><sub>2</sub> &#x0003D; 0.01.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-18-1347748-g0006.tif"/>
</fig>
</sec>
<sec>
<title>3.4 The effects of &#x003BB;<sub>0</sub> on synchronization</title>
<p>As explained in Section 2.1, our system operates in the excitable regime, &#x003BB;<sub>0</sub> &#x0003C; 0, where HB<sub>1</sub> occurs at &#x003BB;<sub>0</sub> &#x0003D; 0. Previous research (e.g., Yu et al., <xref ref-type="bibr" rid="B32">2009</xref>, <xref ref-type="bibr" rid="B33">2021</xref>) has demonstrated that in excitable networks, the distance of the control parameter from a critical point (or excitation threshold) is negatively correlated with synchronization. Moreover, in this section, we investigate the effects (if any) of the control parameter &#x003BB;<sub>0</sub> on the level of synchrony. We only use |&#x00394;&#x003C6;| as a measure of synchrony since it shows results consistent with &#x003C1; and <italic>R</italic>, as demonstrated in <xref ref-type="fig" rid="F5">Figures 5</xref>, <xref ref-type="fig" rid="F6">6</xref>.</p>
<p>We begin by evaluating |&#x00394;&#x003C6;| for various &#x003BB;<sub>0</sub> values within the excitable regime, &#x003BB;<sub>0</sub> &#x0003C; 0, as illustrated in <xref ref-type="fig" rid="F7">Figure 7</xref> (|&#x00394;&#x003C6;| vs. &#x003B4;<sub>2</sub> with fixed &#x003B4;<sub>1</sub> &#x0003D; 0.1). The solid blue, dashed orange, and dotted black lines correspond to &#x003BB;<sub>0</sub> &#x0003D; &#x02212;0.03, &#x02212;0.5, and &#x02212;1, respectively. As observed for the blue curve in <xref ref-type="fig" rid="F5">Figure 5A</xref>, all three |&#x00394;&#x003C6;| curves in <xref ref-type="fig" rid="F7">Figure 7A</xref> exhibit a decreasing trend until they reach a minimum value, after which they start increasing again. This suggests that the level of synchrony can be optimized by adjusting the intrinsic noise intensity &#x003B4;<sub>2</sub>. Furthermore, for a fixed &#x003B4;<sub>2</sub> within the weak intensity range (e.g., &#x003B4;<sub>2</sub> &#x0003C; 0.5 in <xref ref-type="fig" rid="F7">Figure 7A</xref>), |&#x00394;&#x003C6;| is smaller when &#x003BB;<sub>0</sub> is closer to zero. This indicates that in the weak noise regime synchronization is enhanced as &#x003BB;<sub>0</sub> approaches the excitation threshold &#x003BB;<sub>0</sub> &#x0003D; 0 (i.e., HB<sub>1</sub> in <xref ref-type="fig" rid="F1">Figures 1</xref>, <xref ref-type="fig" rid="F2">2</xref>). Conversely, when &#x003B4;<sub>2</sub> is fixed within the strong intensity range (e.g., &#x003B4;<sub>2</sub> &#x0003E; 1.5 in <xref ref-type="fig" rid="F7">Figure 7A</xref>), |&#x00394;&#x003C6;| increases as &#x003BB;<sub>0</sub> approaches the excitation threshold, resulting in lower levels of synchrony.</p>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p><bold>(A)</bold> Absolute value of the phase difference, |&#x00394;&#x003C6;|, vs. noise intensity &#x003B4;<sub>2</sub> for &#x003BB;<sub>0</sub> &#x0003D; &#x02212;0.03, &#x02212;0.5, and &#x02212; 1 with &#x003B4;<sub>1</sub> &#x0003D; 0.1. The blue solid, orange dashed, and black dotted lines correspond to &#x003BB;<sub>0</sub> &#x0003D; &#x02212;0.03, &#x02212;0.5, and &#x02212; 1, respectively. Increasing &#x003BB;<sub>0</sub> leads to a decrease in the minimum value of |&#x00394;&#x003C6;|, implying enhanced synchrony as &#x003BB;<sub>0</sub> approaches zero. <bold>(B)</bold> Minimum |&#x00394;&#x003C6;|, denoted as min{|&#x00394;&#x003C6;|}, as a function of &#x003BB;<sub>0</sub> by varying both &#x003B4;<sub>1</sub> and &#x003B4;<sub>2</sub> over from 0.01 to 5. min{|&#x00394;&#x003C6;|} significantly decreases when &#x003BB;<sub>0</sub> &#x0003E; &#x02212;0.3, suggesting a rapid improvement in synchrony as &#x003BB;<sub>0</sub> approaches zero. <bold>(C)</bold> The sum of the optimal noise intensities, <inline-formula><mml:math id="M19"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>, decreases as &#x003BB;<sub>0</sub> increases, indicating that smaller noise intensities are needed to achieve min{|&#x00394;&#x003C6;|} as &#x003BB;<sub>0</sub> approaches zero. Other parameters are: &#x003B1; &#x0003D; &#x02212;0.2; &#x003B3; &#x0003D; &#x02212;0.2; &#x003C9;<sub>0</sub> &#x0003D; 2; &#x003C9;<sub>1</sub> &#x0003D; 0; <italic>d</italic><sub>1</sub> &#x0003D; 0.3; and <italic>d</italic><sub>2</sub> &#x0003D; 0.01.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-18-1347748-g0007.tif"/>
</fig>
<p>Furthermore, we observe that the minimum values of |&#x00394;&#x003C6;| vary for different &#x003BB;<sub>0</sub>, and the corresponding optimal noise intensities differ as well. The minimum |&#x00394;&#x003C6;| values, denoted by min{|&#x00394;&#x003C6;|}, are marked by black dots in <xref ref-type="fig" rid="F7">Figure 7A</xref>. When &#x003BB;<sub>0</sub> is closer to the excitation threshold, our model achieves a greater degree of synchrony, and a smaller optimal &#x003B4;<sub>2</sub> is required to produce it (for fixed &#x003B4;<sub>1</sub> &#x0003D; 0.1). To explore this further, we plot the values of min{|&#x00394;&#x003C6;|} and the sum of the optimal noise intensities, denoted <inline-formula><mml:math id="M16"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>, for a range of &#x003BB;<sub>0</sub> values in <xref ref-type="fig" rid="F7">Figures 7B</xref>, <xref ref-type="fig" rid="F7">C</xref>, respectively. For &#x003BB;<sub>0</sub> values that are relatively far from the excitation threshold (e.g., &#x003BB;<sub>0</sub> &#x0003C; &#x02212;0.3), the curve in <xref ref-type="fig" rid="F7">Figure 7B</xref> is relatively flat and the curve in <xref ref-type="fig" rid="F7">Figure 7C</xref> has a relatively low slope (&#x02248;&#x02212;0.332). This indicates that min{|&#x00394;&#x003C6;|} and the corresponding optimal intensities remain stable over this region of &#x003BB;<sub>0</sub>. Hence, it follows that small changes in &#x003BB;<sub>0</sub> do not significantly affect the synchronization of our model when &#x003BB;<sub>0</sub> is far from the threshold &#x003BB;<sub>0</sub> &#x0003D; 0. However, when &#x003BB;<sub>0</sub> approaches the excitation threshold (e.g., &#x02212;0.5 &#x0003C; &#x003BB;<sub>0</sub> &#x0003C; 0 in <xref ref-type="fig" rid="F7">Figures 7B</xref>, <xref ref-type="fig" rid="F7">C</xref>), min{|&#x00394;&#x003C6;|} decreases exponentially, and <inline-formula><mml:math id="M17"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula> decreases with a steeper slope (&#x02248; &#x02212;1.84) as <inline-formula><mml:math id="M18"><mml:msub><mml:mrow><mml:mi>&#x003BB;</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>&#x02192;</mml:mo><mml:msup><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mo>-</mml:mo></mml:mrow></mml:msup></mml:math></inline-formula>. These findings suggest that by shifting &#x003BB;<sub>0</sub> closer to the excitation threshold, the noise-induced synchrony of our system can be enhanced, with smaller intensities of noise required to produce this effect.</p>
</sec>
<sec>
<title>3.5 The dual effects of uneven coupling and asymmetric noise on synchrony</title>
<p>In Section 3.3, we briefly described the relationship between the necessary coupling strength and the optimal noise intensity for attaining the maximum level of synchrony. <xref ref-type="fig" rid="F6">Figure 6</xref> illustrates that when the coupling strength <italic>d</italic><sub>1</sub> of oscillator 1 is greater than that of oscillator 2 (<italic>d</italic><sub>2</sub>), achieving the greatest level of synchrony requires &#x003B4;<sub>1</sub> &#x0003C; &#x003B4;<sub>2</sub>. In this section, we investigate how uneven coupling and asymmetric noise intensities interact with one another to influence synchrony. For the sake of simplicity, we set &#x003BB;<sub>0</sub> &#x0003D; &#x02212;0.5 and use |&#x00394;&#x003C6;| to quantify synchrony.</p>
<p>We begin by examining the change in |&#x00394;&#x003C6;| as a function of &#x003B4;<sub>1</sub> and &#x003B4;<sub>2</sub> in four disparate coupling cases: <italic>d</italic><sub>1</sub> &#x0226B; <italic>d</italic><sub>2</sub> (e.g., <italic>d</italic><sub>1</sub> &#x0003D; 0.3 and <italic>d</italic><sub>2</sub> &#x0003D; 0.01 in <xref ref-type="fig" rid="F6">Figure 6A</xref>); <italic>d</italic><sub>1</sub> &#x0003E; <italic>d</italic><sub>2</sub> (e.g., <italic>d</italic><sub>1</sub> &#x0003D; 0.2 and <italic>d</italic><sub>2</sub> &#x0003D; 0.1 in <xref ref-type="fig" rid="F8">Figure 8A</xref>); <italic>d</italic><sub>1</sub> &#x0003C; <italic>d</italic><sub>2</sub> (e.g., <italic>d</italic><sub>1</sub> &#x0003D; 0.1 and <italic>d</italic><sub>2</sub> &#x0003D; 0.2 in <xref ref-type="fig" rid="F8">Figure 8B</xref>), and <italic>d</italic><sub>1</sub> &#x0226A; <italic>d</italic><sub>2</sub> (e.g., <italic>d</italic><sub>1</sub> &#x0003D; 0.01 and <italic>d</italic><sub>2</sub> &#x0003D; 0.3 in <xref ref-type="fig" rid="F8">Figure 8C</xref>). The four cases are also visualized in three dimensions using surface plots and presented in <xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1</xref>. Both heat maps and surface plots show that |&#x00394;&#x003C6;| has both a unique absolute minimum and unique local minimum which correspond to points of increased synchrony. Furthermore, each panel in <xref ref-type="fig" rid="F8">Figure 8</xref> exhibits a relatively low degree of synchrony along the line &#x003B4;<sub>1</sub> &#x0003D; &#x003B4;<sub>2</sub> (i.e., a relatively larger |&#x00394;&#x003C6;|), which is in agreement with our results in <xref ref-type="fig" rid="F6">Figure 6</xref>. These observations suggest that noise-induced synchrony tends to be maximized by asymmetric noise intensities (i.e., &#x003B4;<sub>1</sub>/&#x003B4;<sub>2</sub> &#x02260; 1) as opposed to symmetric noise intensities (i.e., &#x003B4;<sub>1</sub>/&#x003B4;<sub>2</sub> &#x0003D; 1).</p>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p>Heat maps of |&#x00394;&#x003C6;|, as a function of &#x003B4;<sub>1</sub> and &#x003B4;<sub>2</sub> for three cases of coupling strengths: <bold>(A)</bold> <italic>d</italic><sub>1</sub> &#x0003E; <italic>d</italic><sub>2</sub> (with <italic>d</italic><sub>1</sub> &#x0003D; 0.2, <italic>d</italic><sub>2</sub> &#x0003D; 0.1); <bold>(B)</bold> <italic>d</italic><sub>1</sub> &#x0003C; <italic>d</italic><sub>2</sub> (with <italic>d</italic><sub>1</sub> &#x0003D; 0.1, <italic>d</italic><sub>2</sub> &#x0003D; 0.2); and <bold>(C)</bold> <italic>d</italic><sub>1</sub> &#x0226A; <italic>d</italic><sub>2</sub> (with <italic>d</italic><sub>1</sub> &#x0003D; 0.01, <italic>d</italic><sub>2</sub> &#x0003D; 0.3). When <italic>d</italic><sub>1</sub> &#x0226B; <italic>d</italic><sub>2</sub> (<xref ref-type="fig" rid="F6">Figure 6A</xref>) and <italic>d</italic><sub>1</sub> &#x0003E; <italic>d</italic><sub>2</sub> [<bold>(A)</bold> here], the absolute minimum of |&#x00394;&#x003C6;| is located on the upper left side of the heat map, where &#x003B4;<sub>1</sub> &#x0003C; &#x003B4;<sub>2</sub>. When <italic>d</italic><sub>1</sub> &#x0003C; <italic>d</italic><sub>2</sub> <bold>(B)</bold> and <italic>d</italic><sub>1</sub> &#x0226A; <italic>d</italic><sub>2</sub> <bold>(C)</bold>, the absolute minimum of |&#x00394;&#x003C6;| lies on the lower right side, where &#x003B4;<sub>1</sub> &#x0003E; &#x003B4;<sub>2</sub>. The maximum synchrony varies with both &#x003B4;<sub><italic>i</italic></sub> and <italic>d</italic><sub><italic>i</italic></sub>, but it is never achieved when &#x003B4;<sub>1</sub> &#x0003D; &#x003B4;<sub>2</sub>, regardless of the chosen coupling strengths. Other parameters are: &#x003B1; &#x0003D; &#x02212;0.2; &#x003B3; &#x0003D; &#x02212;0.2; &#x003C9;<sub>0</sub> &#x0003D; 2; &#x003C9;<sub>1</sub> &#x0003D; 0; and &#x003BB;<sub>0</sub> &#x0003D; &#x02212;0.5. The three-dimensional visualization of <xref ref-type="fig" rid="F8">Figure 8</xref> (surface plots) is illustrated in <xref ref-type="supplementary-material" rid="SM1">Supplementary Figures 1B</xref>&#x02013;<xref ref-type="supplementary-material" rid="SM1">D</xref>.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-18-1347748-g0008.tif"/>
</fig>
<p><xref ref-type="fig" rid="F6">Figures 6A</xref>, <xref ref-type="fig" rid="F8">8</xref> further indicate that the ratio of coupling strengths, <italic>d</italic><sub>1</sub>/<italic>d</italic><sub>2</sub>, determines whether the minimum point of |&#x00394;&#x003C6;| is local or absolute. Let abs.min{|&#x00394;&#x003C6;|} denote the absolute minimum of {|&#x00394;&#x003C6;|} and <inline-formula><mml:math id="M22"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula> denote the corresponding ratio of optimal noise intensities. We find that when <italic>d</italic><sub>1</sub> &#x0226B; <italic>d</italic><sub>2</sub> or <italic>d</italic><sub>1</sub> &#x0003E; <italic>d</italic><sub>2</sub>, abs.min{|&#x00394;&#x003C6;|} is located on the upper left (i.e., &#x003B4;<sub>2</sub> &#x0003E; &#x003B4;<sub>1</sub>) of the heat maps in <xref ref-type="fig" rid="F6">Figures 6A</xref>, <xref ref-type="fig" rid="F8">8A</xref>, with <inline-formula><mml:math id="M23"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>&#x0226A;</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula>. Conversely, when <italic>d</italic><sub>1</sub> &#x0003C; <italic>d</italic><sub>2</sub> or <italic>d</italic><sub>1</sub> &#x0226A; <italic>d</italic><sub>2</sub>, abs.min{|&#x00394;&#x003C6;|} is located on the lower right (i.e., &#x003B4;<sub>1</sub> &#x0003E; &#x003B4;<sub>2</sub>) of the heat maps in <xref ref-type="fig" rid="F8">Figures 8B</xref>, <xref ref-type="fig" rid="F8">C</xref>, with <inline-formula><mml:math id="M24"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>&#x0226B;</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula>. <xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1</xref> provides a clear depiction of the gradual transition of abs.min{|&#x00394;&#x003C6;|} from the left side (with &#x003B4;<sub>1</sub> &#x0003C; &#x003B4;<sub>2</sub>) to the right side (with &#x003B4;<sub>1</sub> &#x0003E; &#x003B4;<sub>2</sub>) as <italic>d</italic><sub>1</sub>/<italic>d</italic><sub>2</sub> changes from larger to smaller values.</p>
<p>To systematically explore the interaction and effects of the noise intensities and coupling strengths on synchrony, we consider the parameter space <italic>d</italic><sub>1</sub>, <italic>d</italic><sub>2</sub> &#x02208; [0.01, 0.3] and &#x003B4;<sub>1</sub>, &#x003B4;<sub>2</sub> &#x02208; [0.01, 0.3], and compute abs.min{|&#x00394;&#x003C6;|} and <inline-formula><mml:math id="M25"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>. We present the results in <xref ref-type="fig" rid="F9">Figures 9A</xref>, <xref ref-type="fig" rid="F9">B</xref> in the form of heat maps. The upper-left &#x0201C;triangle&#x0201D; of <xref ref-type="fig" rid="F9">Figures 9A</xref>, <xref ref-type="fig" rid="F9">B</xref>, where <italic>d</italic><sub>1</sub> &#x0003C; <italic>d</italic><sub>2</sub>, corresponds to the absolute minimum points on the lower right of the heat maps in <xref ref-type="fig" rid="F8">Figures 8B</xref>, <xref ref-type="fig" rid="F8">C</xref>. And the lower-right &#x0201C;triangle&#x0201D; of <xref ref-type="fig" rid="F9">Figures 9A</xref>, <xref ref-type="fig" rid="F9">B</xref>, with <italic>d</italic><sub>1</sub> &#x0003E; <italic>d</italic><sub>2</sub>, correspond to the absolute minimum points on the upper left of the heat maps in <xref ref-type="fig" rid="F6">Figures 6A</xref>, <xref ref-type="fig" rid="F8">8A</xref>.</p>
<fig id="F9" position="float">
<label>Figure 9</label>
<caption><p><bold>(A)</bold> Heat map illustrates the ratio of the optimal noise intensities, <inline-formula><mml:math id="M20"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>, required to minimize |&#x00394;&#x003C6;|, in relation to varying both coupling strengths, <italic>d</italic><sub>1</sub> and <italic>d</italic><sub>2</sub>. <bold>(B)</bold> Heat map of the absolute minimum of |&#x00394;&#x003C6;|, denoted as abs.min{|&#x00394;&#x003C6;|}, as a function of <italic>d</italic><sub>1</sub> and <italic>d</italic><sub>2</sub>. The maximal noise-induced synchrony is not achieveed when <inline-formula><mml:math id="M21"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula> or when <italic>d</italic><sub>1</sub> &#x0003D; <italic>d</italic><sub>2</sub>. <bold>(C)</bold> Schematic illustration of a two-neuron motif used to enhance noise-induced synchrony. One oscillator (receiver neuron) receives robust coupling from the other oscillator (source oscillator) along with minimal noise input. Conversely, the source neuron experiences minimal coupling from the receiver neuron but substantial noise input. <bold>(D)</bold> Two-neuron motif featuring a uni-directional connection, optimizing noise-induced synchrony. Other parameters in <bold>(A, B)</bold> are: &#x003B1; &#x0003D; &#x02212;0.2; &#x003B3; &#x0003D; &#x02212;0.2; &#x003C9;<sub>0</sub> &#x0003D; 2; &#x003C9;<sub>1</sub> &#x0003D; 0; and &#x003BB;<sub>0</sub> &#x0003D; &#x02212;0.5.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fncom-18-1347748-g0009.tif"/>
</fig>
<p><xref ref-type="fig" rid="F9">Figure 9A</xref> illustrates how the choice of coupling affects <inline-formula><mml:math id="M26"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>. Previously, we have shown that the optimal <inline-formula><mml:math id="M27"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula> does not equal 1 for some example pairs of <italic>d</italic><sub>1</sub> and <italic>d</italic><sub>2</sub> in <xref ref-type="fig" rid="F6">Figures 6</xref>, <xref ref-type="fig" rid="F8">8</xref>. Here our result is more robust: noise-induced synchrony is never maximized at &#x003B4;<sub>1</sub>/&#x003B4;<sub>2</sub> &#x0003D; 1 no matter the choice of coupling strengths <italic>d</italic><sub>1</sub> and <italic>d</italic><sub>2</sub>. In particular, even when <italic>d</italic><sub>1</sub> &#x0003D; <italic>d</italic><sub>2</sub>, <inline-formula><mml:math id="M28"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>&#x02260;</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula>. <xref ref-type="fig" rid="F9">Figure 9A</xref> also indicates that there is a strong relationship between <italic>d</italic><sub>1</sub>/<italic>d</italic><sub>2</sub> and <inline-formula><mml:math id="M29"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>. Specifically, when <italic>d</italic><sub>1</sub>/<italic>d</italic><sub>2</sub> &#x0003C; 1 (upper-left &#x0201C;triangle&#x0201D;), <inline-formula><mml:math id="M30"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>&#x02208;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mn>15</mml:mn><mml:mo>,</mml:mo><mml:mn>35</mml:mn></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>&#x0226B;</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula>; whereas, when <italic>d</italic><sub>1</sub>/<italic>d</italic><sub>1</sub> &#x0003E; 1 (lower-right &#x0201C;triangle&#x0201D;), <inline-formula><mml:math id="M31"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>&#x02208;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:mn>025</mml:mn><mml:mo>,</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:mn>04</mml:mn></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>&#x0226A;</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula>. We also compute <inline-formula><mml:math id="M32"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula> and <inline-formula><mml:math id="M33"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula> for both <italic>d</italic><sub>1</sub>/<italic>d</italic><sub>2</sub> cases. When <italic>d</italic><sub>1</sub>/<italic>d</italic><sub>2</sub> &#x0003C; 1, the mean optimal noise intensities are <inline-formula><mml:math id="M34"><mml:mover accent="false" class="mml-overline"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:mrow><mml:mo accent="true">&#x000AF;</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:mn>3462</mml:mn></mml:math></inline-formula> and <inline-formula><mml:math id="M35"><mml:mover accent="false" class="mml-overline"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:mrow><mml:mo accent="true">&#x000AF;</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:mn>015</mml:mn></mml:math></inline-formula>, which results in an average ratio of 0.3462/0.015 &#x02248; 23.08. And conversely, when <italic>d</italic><sub>1</sub>/<italic>d</italic><sub>2</sub> &#x0003E; 1, <inline-formula><mml:math id="M36"><mml:mover accent="false" class="mml-overline"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:mrow><mml:mo accent="true">&#x000AF;</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:mn>015</mml:mn></mml:math></inline-formula> and <inline-formula><mml:math id="M37"><mml:mover accent="false" class="mml-overline"><mml:mrow><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup></mml:mrow><mml:mo accent="true">&#x000AF;</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:mn>3462</mml:mn></mml:math></inline-formula>, which results in an average ratio of 0.015/0.3462 &#x02248; 0.043 (note that the latter ratio is the reciprocal the former).</p>
<p>Next, we examine how the coupling strengths influence the degree of synchrony of our system, as measured by abs.min{|&#x00394;&#x003C6;|}. As in <xref ref-type="fig" rid="F9">Figure 9B</xref>, abs.min{|&#x00394;&#x003C6;|} exhibits symmetry with respect to the diagonal line <italic>d</italic><sub>1</sub> &#x0003D; <italic>d</italic><sub>2</sub>. Along this line, abs.min{|&#x00394;&#x003C6;|} decreases as both <italic>d</italic><sub>1</sub> and <italic>d</italic><sub>2</sub> increase. However, the smallest value of abs.min{|&#x00394;&#x003C6;|} is not located along this line. This suggests that when the oscillators are symmetrically coupled, the degree of synchrony is positively correlated with the coupling strength <italic>d</italic><sub>1</sub> &#x0003D; <italic>d</italic><sub>2</sub>, but does not lead to maximal synchronization.</p>
<p>In contrast, when the coupling strengths are asymmetric, abs.min{|&#x00394;&#x003C6;|} decreases as the absolute difference between two coupling strength, |<italic>d</italic><sub>1</sub> &#x02212; <italic>d</italic><sub>2</sub>|, becomes larger, as seen in <xref ref-type="fig" rid="F9">Figure 9B</xref>. Furthermore, the lowest values of abs.min{|&#x00394;&#x003C6;|} occur at the upper-left and lower-right corners of <xref ref-type="fig" rid="F9">Figure 9B</xref>, where one coupling strength is nearly zero and the other is the largest in the considered range. This indicates that a one-way coupling maximizes noise-induced synchronization between two coupled oscillators. Overall, synchrony is maximized when <italic>d</italic><sub>1</sub> &#x0226A; <italic>d</italic><sub>2</sub> and &#x003B4;<sub>1</sub>/&#x003B4;<sub>2</sub> &#x02248; 15 (the upper-left corner in <xref ref-type="fig" rid="F9">Figures 9A</xref>, <xref ref-type="fig" rid="F9">B</xref>), or <italic>d</italic><sub>1</sub> &#x0226B; <italic>d</italic><sub>2</sub> and &#x003B4;<sub>1</sub>/&#x003B4;<sub>2</sub> &#x02248; 1/15 (lower-right corner in <xref ref-type="fig" rid="F9">Figures 9A</xref>, <xref ref-type="fig" rid="F9">B</xref>).</p>
<p>Our findings indicate that the noise-induced synchrony between two coupled oscillators is enhanced in a specific scenario, as illustrated in <xref ref-type="fig" rid="F9">Figure 9C</xref>. Specifically, one oscillator (referred to as the source neuron) experiences strong intrinsic noise and receives weak coupling input from the other oscillator (the receiver neuron). Conversely, the receiver neuron experiences weak intrinsic noise but receives substantial coupling input from the source neuron. The optimum noise-induced synchrony occurs in the extreme case, i.e., a uni-directional connection where the input from the receiver neuron to the source neuron is absent, as demonstrated in <xref ref-type="fig" rid="F9">Figure 9D</xref>.</p>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>4 Discussion</title>
<p>In many investigations exploring the noise-inuced synchronization of coupled excitable neurons, assumptions of homogeneity or near-homogeneity (e.g., common noise, symmetric noise, and equal coupling strength) are undertaken. To move beyond these seemingly unrealistic setups, our study considers a heterogeneous configuration through the introduction of independent noise sources with varying intensities and distinct coupling strengths for each neuron. We investigate the synchronous dynamics of the most representative neuron motif, characterized mathematically as a pair of coupled &#x003BB;&#x02212;&#x003C9; oscillators in the heterogeneous context. These neural oscillators remain are quiescent in the absence of noise but may be excited by adding an intrinsic noise stimulus. Our results show that noise can induce synchronization in coupled oscillators, and this synchronization is enhanced by positioning the model closer to the bifurcation point (i.e., excitation threshold), which is consistent with previous studies (e.g., Yu et al., <xref ref-type="bibr" rid="B34">2006</xref>, <xref ref-type="bibr" rid="B35">2008</xref>, <xref ref-type="bibr" rid="B33">2021</xref>; Thompson et al., <xref ref-type="bibr" rid="B27">2012</xref>). In previous research, the noise-induced synchrony of excitable systems has been studied extensively with a focus on symmetrical interactions between oscillators (e.g., Neiman et al., <xref ref-type="bibr" rid="B15">1999</xref>; Rosenblum et al., <xref ref-type="bibr" rid="B22">2001</xref>; Freund et al., <xref ref-type="bibr" rid="B7">2003</xref>). However, in biological systems, interactions between oscillators are often asymmetric, and thus, the assumption of symmetric interactions may be overly restrictive (Cimponeriu et al., <xref ref-type="bibr" rid="B4">2003</xref>; Sheeba et al., <xref ref-type="bibr" rid="B24">2009</xref>).</p>
<p>We investigate the impacts of two sources of heterogeneity (uneven coupling and asymmetric intrinsic noise) on the noise-induced synchrony of two coupled oscillators. Our results indicate that when the two noise intensities, &#x003B4;<sub>1</sub> and &#x003B4;<sub>2</sub>, are asymmetric (i.e., &#x003B4;<sub>1</sub>/&#x003B4;<sub>2</sub> &#x02260; 1), synchrony is promoted. We further find that the synchronization between two coupled oscillators is optimized when the absolute difference between the coupling strengths, |<italic>d</italic><sub>1</sub> &#x02212; <italic>d</italic><sub>2</sub>|, is as large as possible, irrespective of the coupling strengths selected. Furthermore, we have identified a correlation between the asymmetry of coupling and the intensity of noise required to maximize it. Specifically, when <italic>d</italic><sub>1</sub>/<italic>d</italic><sub>2</sub> &#x0003C; 1, it leads to <inline-formula><mml:math id="M38"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>&#x0003E;</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula>, whereas, when <italic>d</italic><sub>1</sub>/<italic>d</italic><sub>2</sub> &#x0003E; 1, it leads to <inline-formula><mml:math id="M39"><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>/</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msubsup><mml:mo>&#x0003C;</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula>.</p>
<p>Our results suggest that there is a strong relationship between uneven coupling and asymmetric noise in coupled oscillators, which have potential applications in various real-world problems that exhibit asymmetry in the interactions between oscillators, such as cardio-respiratory electroencephalogram (EEG) interactions (Palu&#x00161; and Stefanovska, <xref ref-type="bibr" rid="B16">2003</xref>; Sheeba et al., <xref ref-type="bibr" rid="B24">2009</xref>), optical communication systems, the detection of radar signals in presence of channel noise (Tsang and Lindsey, <xref ref-type="bibr" rid="B29">1986</xref>), and neuronal dynamics (Singer, <xref ref-type="bibr" rid="B25">1999</xref>; Sheeba et al., <xref ref-type="bibr" rid="B24">2009</xref>). Moreover, our results indicate that two-neuron motifs with uni-direction (i.e., one-way) connections may exhibit a greater propensity for in-phase synchronization than two-neuron motifs with bi-directional (i.e., two-way) connections, which reinforces our understanding of the functional significance of network motifs in neuronal dynamics.</p>
<p>However, further investigations are necessary to explore the relationship between the ratios <italic>d</italic><sub>1</sub>/<italic>d</italic><sub>2</sub> and &#x003B4;<sub>1</sub>/&#x003B4;<sub>2</sub> more extensively. An extension of our current work may include investigating: (1) the role of asymmetric noise and coupling in the anti-phase synchronization of coupled oscillators by considering inhibitory coupling (i.e., <italic>d</italic><sub><italic>i</italic></sub> &#x0003C; 0 for <italic>i</italic> &#x0003D; 1or2); and (2) investigation of the role of noise and coupling asymmetries in other common network motifs, such as three-neuron feed-forward-loops.</p>
<p>The HB is one of the most common bifurcation schemes as it describes the appearance or disappearance of regular orbits with a slight change of a parameter. In this study, we employ the normal form of HB to investigate a specific scenario: noise-induced synchrony near a supercritical Hopf point. By focusing on this particular region, we aim to uncover the underlying mechanisms that lead to synchronized oscillations in physical or biological systems. Our results not only shed light on the phenomena of interest but also offer potential explanations or predictions for various systems that undergo similar shifts from a quiescent state to sustained oscillations. However, it is important to acknowledge that different systems may exhibit other bifurcation structures, necessitating further investigation to comprehensively understand their dynamic behavior.</p>
</sec>
<sec sec-type="data-availability" id="s5">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="sec" rid="s9">Supplementary material</xref>, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec sec-type="author-contributions" id="s6">
<title>Author contributions</title>
<p>GJ: Formal analysis, Investigation, Methodology, Software, Visualization, Writing - original draft, Writing - review &#x00026; editing. NY: Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Supervision, Visualization, Writing - original draft, Writing - review &#x00026; editing.</p>
</sec>
</body>
<back>
<sec sec-type="funding-information" id="s7">
<title>Funding</title>
<p>GJ was supported by the Queen Elizabeth II Graduate Scholarship in Science and Technology and the Ontario Graduate Scholarship in Canada. NY received startup funding from the Faculty of Science at Toronto Metropolitan University and the Discovery Grant from the Natural Sciences and Engineering Research Council of Canada (NSERC).</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s8">
<title>Publisher&#x00027;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec sec-type="supplementary-material" id="s9">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fncom.2024.1347748/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fncom.2024.1347748/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Data_Sheet_1.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Amro</surname> <given-names>R. M.</given-names></name> <name><surname>Lindner</surname> <given-names>B.</given-names></name> <name><surname>Neiman</surname> <given-names>A. B.</given-names></name></person-group> (<year>2015</year>). <article-title>Phase diffusion in unequally noisy coupled oscillators</article-title>. <source>Phys. Rev. Lett.</source> <volume>115</volume>:<fpage>e034101</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.115.034101</pub-id><pub-id pub-id-type="pmid">26230796</pub-id></citation></ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Blasius</surname> <given-names>B.</given-names></name></person-group> (<year>2005</year>). <article-title>Anomalous phase synchronization in two asymmetrically coupled oscillators in the presence of noise</article-title>. <source>Phys. Rev. E</source> <volume>72</volume>:<fpage>e066216</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevE.72.066216</pub-id><pub-id pub-id-type="pmid">16486049</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Boccaletti</surname> <given-names>S.</given-names></name> <name><surname>Latora</surname> <given-names>V.</given-names></name> <name><surname>Moreno</surname> <given-names>Y.</given-names></name> <name><surname>Chavez</surname> <given-names>M.</given-names></name> <name><surname>Hwang</surname> <given-names>D. U.</given-names></name></person-group> (<year>2006</year>). <article-title>Complex networks: structure and dynamics</article-title>. <source>Phys. Rep.</source> <volume>424</volume>, <fpage>175</fpage>&#x02013;<lpage>308</lpage>. <pub-id pub-id-type="doi">10.1016/j.physrep.2005.10.009</pub-id></citation>
</ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cimponeriu</surname> <given-names>L.</given-names></name> <name><surname>Rosenblum</surname> <given-names>M. G.</given-names></name> <name><surname>Fieseler</surname> <given-names>T.</given-names></name> <name><surname>Dammers</surname> <given-names>J.</given-names></name> <name><surname>Schiek</surname> <given-names>M.</given-names></name> <name><surname>Majtanik</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2003</year>). <article-title>Inferring asymmetric relations between interacting neuronal oscillators</article-title>. <source>Progr. Theoret. Phys. Suppl.</source> <volume>150</volume>, <fpage>22</fpage>&#x02013;<lpage>36</lpage>. <pub-id pub-id-type="doi">10.1143/PTPS.150.22</pub-id></citation>
</ref>
<ref id="B5">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Ermentrout</surname> <given-names>B.</given-names></name></person-group> (<year>2012</year>). <source>XPPAUT</source> (<publisher-loc>Dordrecht</publisher-loc>: <publisher-name>Springer Netherlands</publisher-name>), <fpage>519</fpage>&#x02013;<lpage>531</lpage>.</citation>
</ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fell</surname> <given-names>J.</given-names></name> <name><surname>Axmacher</surname> <given-names>N.</given-names></name></person-group> (<year>2011</year>). <article-title>The role of phase synchronization in memory processes</article-title>. <source>Nat. Rev. Neurosci.</source> <volume>12</volume>, <fpage>105</fpage>&#x02013;<lpage>118</lpage>. <pub-id pub-id-type="doi">10.1038/nrn2979</pub-id><pub-id pub-id-type="pmid">21248789</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Freund</surname> <given-names>J. A.</given-names></name> <name><surname>Schimansky-Geier</surname> <given-names>L.</given-names></name> <name><surname>H&#x000E4;nggi</surname> <given-names>P.</given-names></name></person-group> (<year>2003</year>). <article-title>Frequency and phase synchronization in stochastic systems</article-title>. <source>Chaos</source> <volume>13</volume>, <fpage>225</fpage>&#x02013;<lpage>238</lpage>. <pub-id pub-id-type="doi">10.1063/1.1500497</pub-id><pub-id pub-id-type="pmid">12675429</pub-id></citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gambuzza</surname> <given-names>L.</given-names></name> <name><surname>G&#x000F3;mez-Garde&#x000F1;es</surname> <given-names>J.</given-names></name> <name><surname>Frasca</surname> <given-names>M.</given-names></name></person-group> (<year>2016</year>). <article-title>Amplitude dynamics favors synchronization in complex networks</article-title>. <source>Sci. Rep.</source> <volume>6</volume>:<fpage>24915</fpage>. <pub-id pub-id-type="doi">10.1038/srep24915</pub-id><pub-id pub-id-type="pmid">27108847</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gang</surname> <given-names>H.</given-names></name> <name><surname>Ditzinger</surname> <given-names>T.</given-names></name> <name><surname>Ning</surname> <given-names>C. Z.</given-names></name> <name><surname>Haken</surname> <given-names>H.</given-names></name></person-group> (<year>1993</year>). <article-title>Stochastic resonance without external periodic force</article-title>. <source>Phys. Rev. Lett.</source> <volume>71</volume>, <fpage>807</fpage>&#x02013;<lpage>810</lpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.71.807</pub-id><pub-id pub-id-type="pmid">10055373</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kawasaki</surname> <given-names>M.</given-names></name> <name><surname>Kitajo</surname> <given-names>K.</given-names></name> <name><surname>Yamaguchi</surname> <given-names>Y.</given-names></name></person-group> (<year>2018</year>). <article-title>Sensory-motor synchronization in the brain corresponds to behavioral synchronization between individuals</article-title>. <source>Neuropsychologia</source> <volume>119</volume>, <fpage>59</fpage>&#x02013;<lpage>67</lpage>. <pub-id pub-id-type="doi">10.1016/j.neuropsychologia.2018.07.026</pub-id><pub-id pub-id-type="pmid">30055179</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Masoliver</surname> <given-names>M.</given-names></name> <name><surname>Masoller</surname> <given-names>C.</given-names></name></person-group> (<year>2018</year>). <article-title>Sub-threshold signal encoding in coupled fitzhugh-nagumo neurons</article-title>. <source>Sci. Rep.</source> <volume>27</volume>:<fpage>8276</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-018-26618-8</pub-id><pub-id pub-id-type="pmid">29844354</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Milo</surname> <given-names>R.</given-names></name> <name><surname>Shen-Orr</surname> <given-names>S.</given-names></name> <name><surname>Itzkovitz</surname> <given-names>S.</given-names></name> <name><surname>Kashtan</surname> <given-names>N.</given-names></name> <name><surname>Chklovskii</surname> <given-names>D.</given-names></name></person-group> (<year>2002</year>). <article-title>Network motifs: simple building blocks of complex networks</article-title>. <source>Science</source> <volume>298</volume>, <fpage>824</fpage>&#x02013;<lpage>827</lpage>. <pub-id pub-id-type="doi">10.1126/science.298.5594.824</pub-id><pub-id pub-id-type="pmid">12399590</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Morgan</surname> <given-names>R. J.</given-names></name> <name><surname>Soltesz</surname> <given-names>I.</given-names></name></person-group> (<year>2008</year>). <article-title>Nonrandom connectivity of the epileptic dentate gyrus predicts a major role for neuronal hubs in seizures</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>105</volume>, <fpage>6179</fpage>&#x02013;<lpage>6184</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.0801372105</pub-id><pub-id pub-id-type="pmid">18375756</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mormann</surname> <given-names>F.</given-names></name> <name><surname>Lehnertz</surname> <given-names>K.</given-names></name> <name><surname>David</surname> <given-names>P.</given-names></name> <name><surname>E. Elger</surname> <given-names>C.</given-names></name></person-group> (<year>2000</year>). <article-title>Mean phase coherence as a measure for phase synchronization and its application to the eeg of epilepsy patients</article-title>. <source>Phys. D</source> <volume>144</volume>, <fpage>358</fpage>&#x02013;<lpage>369</lpage>. <pub-id pub-id-type="doi">10.1016/S0167-2789(00)00087-7</pub-id></citation>
</ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Neiman</surname> <given-names>A.</given-names></name> <name><surname>Schimansky-Geier</surname> <given-names>L.</given-names></name> <name><surname>Cornell-Bell</surname> <given-names>A.</given-names></name> <name><surname>Moss</surname> <given-names>F.</given-names></name></person-group> (<year>1999</year>). <article-title>Noise-enhanced phase synchronization in excitable media</article-title>. <source>Phys. Rev. Lett.</source> <volume>83</volume>, <fpage>4896</fpage>&#x02013;<lpage>4899</lpage>.</citation>
</ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Palu&#x00161;</surname> <given-names>M.</given-names></name> <name><surname>Stefanovska</surname> <given-names>A.</given-names></name></person-group> (<year>2003</year>). <article-title>Direction of coupling from phases of interacting oscillators: an information-theoretic approach</article-title>. <source>Phys. Rev. E</source> <volume>67</volume>:<fpage>055201</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevE.67.055201</pub-id><pub-id pub-id-type="pmid">12786211</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Perc</surname> <given-names>M.</given-names></name></person-group> (<year>2007</year>). <article-title>Effects of small-world connectivity on noise-induced temporal and spatial order in neural media</article-title>. <source>Chaos Solit. Fract.</source> <volume>31</volume>, <fpage>280</fpage>&#x02013;<lpage>291</lpage>. <pub-id pub-id-type="doi">10.1016/j.chaos.2005.10.018</pub-id></citation>
</ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Perc</surname> <given-names>M.</given-names></name> <name><surname>Marhl</surname> <given-names>M.</given-names></name></person-group> (<year>2006</year>). <article-title>Pacemaker enhanced noise-induced synchrony in cellular arrays</article-title>. <source>Phys. Lett. A</source> <volume>353</volume>, <fpage>372</fpage>&#x02013;<lpage>377</lpage>. <pub-id pub-id-type="doi">10.1016/j.physleta.2005.12.107</pub-id></citation>
</ref>
<ref id="B19">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Pikovsky</surname> <given-names>A.</given-names></name></person-group> (<year>2002</year>). <source>Synchronization: A Universal Concept in Nonlinear Sciences</source>. <publisher-loc>Cambridge</publisher-loc>: <publisher-name>Cambridge University Press</publisher-name>.</citation>
</ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pikovsky</surname> <given-names>A. S.</given-names></name> <name><surname>Kurths</surname> <given-names>J.</given-names></name></person-group> (<year>1997</year>). <article-title>Coherence resonance in a noise-driven excitable system</article-title>. <source>Phys. Rev. Lett.</source> <volume>78</volume>, <fpage>775</fpage>&#x02013;<lpage>778</lpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.78.775</pub-id></citation>
</ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Reigl</surname> <given-names>M.</given-names></name> <name><surname>Alon</surname> <given-names>U.</given-names></name> <name><surname>Chklovskii</surname> <given-names>D.</given-names></name></person-group> (<year>2004</year>). <article-title>Search for computational modules in the <italic>C. elegans</italic> brain</article-title>. <source>BMC Biol.</source> <volume>2</volume>:<fpage>25</fpage>. <pub-id pub-id-type="doi">10.1186/1741-7007-2-25</pub-id><pub-id pub-id-type="pmid">15574204</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Rosenblum</surname> <given-names>M.</given-names></name> <name><surname>Pikovsky</surname> <given-names>A.</given-names></name> <name><surname>Kurths</surname> <given-names>J.</given-names></name> <name><surname>Sch&#x000E4;fer</surname> <given-names>C.</given-names></name> <name><surname>Tass</surname> <given-names>P.</given-names></name></person-group> (<year>2001</year>). <article-title>Phase synchronization: from theory to data analysis</article-title>, in <source>Neuro-Informatics and Neural Modelling, Volume 4 of Handbook of Biological Physics</source> (<publisher-loc>Amsterdam</publisher-loc>: <publisher-name>Elsevier Science</publisher-name>), <fpage>279</fpage>&#x02013;<lpage>321</lpage>.</citation>
</ref>
<ref id="B23">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Rosenblum</surname> <given-names>M. G.</given-names></name> <name><surname>Pikovsky</surname> <given-names>A. S.</given-names></name> <name><surname>Kurths</surname> <given-names>J.</given-names></name></person-group> (<year>2022</year>). <article-title>Synchronization approach to analysis of biological systems</article-title>, in <source>The Random and Fluctuating World: Celebrating Two Decades of Fluctuation and Noise Letters</source> (<publisher-loc>Singapore</publisher-loc>: <publisher-name>World Scientific</publisher-name>), <fpage>335</fpage>&#x02013;<lpage>344</lpage>.</citation>
</ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sheeba</surname> <given-names>J. H.</given-names></name> <name><surname>Chandrasekar</surname> <given-names>V. K.</given-names></name> <name><surname>Stefanovska</surname> <given-names>A.</given-names></name> <name><surname>McClintock</surname> <given-names>P. V. E.</given-names></name></person-group> (<year>2009</year>). <article-title>Asymmetry-induced effects in coupled phase-oscillator ensembles: Routes to synchronization</article-title>. <source>Phys. Rev. E</source> <volume>79</volume>:<fpage>e046210</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevE.79.046210</pub-id><pub-id pub-id-type="pmid">19518315</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Singer</surname> <given-names>W.</given-names></name></person-group> (<year>1999</year>). <article-title>Striving for coherence</article-title>. <source>Nature</source> <volume>397</volume>, <fpage>391</fpage>&#x02013;<lpage>393</lpage>. <pub-id pub-id-type="doi">10.1038/17021</pub-id></citation>
</ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Song</surname> <given-names>S.</given-names></name> <name><surname>Sj&#x000F6;str&#x000F6;m</surname> <given-names>P.</given-names></name> <name><surname>Reigl</surname> <given-names>M.</given-names></name> <name><surname>Nelson</surname> <given-names>S.</given-names></name> <name><surname>Chklovskii</surname> <given-names>D.</given-names></name></person-group> (<year>2005</year>). <article-title>Highly nonrandom features of synaptic connectivity in local cortical circuits</article-title>. <source>PLoS Biol.</source> <volume>3</volume>:<fpage>e68</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pbio.0030068</pub-id><pub-id pub-id-type="pmid">15737062</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thompson</surname> <given-names>W. F.</given-names></name> <name><surname>Kuske</surname> <given-names>R.</given-names></name> <name><surname>Li</surname> <given-names>Y.-X.</given-names></name></person-group> (<year>2012</year>). <article-title>Stochastic phase dynamics of noise driven synchronization of two conditional coherent oscillators</article-title>. <source>Discr. Contin. Dyn. Syst.</source> <volume>32</volume>:<fpage>2971</fpage>. <pub-id pub-id-type="doi">10.3934/dcds.2012.32.2971</pub-id></citation>
</ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Touboul</surname> <given-names>J. D.</given-names></name> <name><surname>Piette</surname> <given-names>C.</given-names></name> <name><surname>Venance</surname> <given-names>L.</given-names></name> <name><surname>Ermentrout</surname> <given-names>G. B.</given-names></name></person-group> (<year>2020</year>). <article-title>Noise-induced synchronization and antiresonance in interacting excitable systems: applications to deep brain stimulation in parkinson&#x00027;s disease</article-title>. <source>Phys. Rev. X</source> <volume>10</volume>:<fpage>e011073</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevX.10.011073</pub-id></citation>
</ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tsang</surname> <given-names>C.-S.</given-names></name> <name><surname>Lindsey</surname> <given-names>W.</given-names></name></person-group> (<year>1986</year>). <article-title>Bit synchronization in the presence of asymmetric channel noise</article-title>. <source>IEEE Trans. Commun.</source> <volume>34</volume>, <fpage>528</fpage>&#x02013;<lpage>537</lpage>. <pub-id pub-id-type="doi">10.1109/TCOM.1986.1096577</pub-id></citation>
</ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Varela</surname> <given-names>F.</given-names></name> <name><surname>Lachaux</surname> <given-names>J. P.</given-names></name> <name><surname>Rodriguez</surname> <given-names>E.</given-names></name> <name><surname>Martinerie</surname> <given-names>J.</given-names></name></person-group> (<year>2001</year>). <article-title>The brainweb: phase synchronization and large-scale integration</article-title>. <source>Nat. Rev. Neurosci.</source> <volume>2</volume>, <fpage>229</fpage>&#x02013;<lpage>239</lpage>. <pub-id pub-id-type="doi">10.1038/35067550</pub-id><pub-id pub-id-type="pmid">11283746</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Womelsdorf</surname> <given-names>T.</given-names></name> <name><surname>Fries</surname> <given-names>P.</given-names></name></person-group> (<year>2007</year>). <article-title>The role of neuronal synchronization in selective attention</article-title>. <source>Curr. Opin. Neurobiol.</source> <volume>17</volume>, <fpage>154</fpage>&#x02013;<lpage>160</lpage>. <pub-id pub-id-type="doi">10.1016/j.conb.2007.02.002</pub-id><pub-id pub-id-type="pmid">17306527</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname> <given-names>G.</given-names></name> <name><surname>Yi</surname> <given-names>M.</given-names></name> <name><surname>Jia</surname> <given-names>Y.</given-names></name> <name><surname>Tang</surname> <given-names>J.</given-names></name></person-group> (<year>2009</year>). <article-title>A constructive role of internal noise on coherence resonance induced by external noise in a calcium oscillation system</article-title>. <source>Chaos Solit. Fract.</source> <volume>41</volume>, <fpage>273</fpage>&#x02013;<lpage>283</lpage>. <pub-id pub-id-type="doi">10.1016/j.chaos.2007.12.001</pub-id></citation>
</ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname> <given-names>N.</given-names></name> <name><surname>Jagdev</surname> <given-names>G.</given-names></name> <name><surname>Morgovsky</surname> <given-names>M.</given-names></name></person-group> (<year>2021</year>). <article-title>Noise-induced network bursts and coherence in a calcium-mediated neural network</article-title>. <source>Heliyon</source> <volume>7</volume>:<fpage>e08612</fpage>. <pub-id pub-id-type="doi">10.1016/j.heliyon.2021.e08612</pub-id><pub-id pub-id-type="pmid">35024481</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname> <given-names>N.</given-names></name> <name><surname>Kuske</surname> <given-names>R.</given-names></name> <name><surname>Li</surname> <given-names>Y.</given-names></name></person-group> (<year>2006</year>). <article-title>Stochastic phase dynamics: multiscale behavior and coherence measures</article-title>. <source>Phys. Rev. E</source> <volume>73</volume>:<fpage>e056205</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevE.73.056205</pub-id><pub-id pub-id-type="pmid">16803024</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname> <given-names>N.</given-names></name> <name><surname>Kuske</surname> <given-names>R.</given-names></name> <name><surname>Li</surname> <given-names>Y. X.</given-names></name></person-group> (<year>2008</year>). <article-title>Stochastic phase dynamics and noise-induced mixed-mode oscillations in coupled oscillators</article-title>. <source>Chaos</source> <volume>18</volume>:<fpage>e015112</fpage>. <pub-id pub-id-type="doi">10.1063/1.2790369</pub-id><pub-id pub-id-type="pmid">18377093</pub-id></citation></ref>
</ref-list>
</back>
</article>