Tuning Logical Phi-Bit State Vectors in an Externally Driven Nonlinear Array of Acoustic Waveguides via Drivers’ Phase

: We experimentally navigate the Hilbert space of two logical phi-bits supported by an externally driven nonlinear array of coupled acoustic waveguides by parametrically changing the relative phase of the drivers. We observe sharp phase jumps of approximately 180 ◦ in the individual phi-bit states as a result of the phase tuning of the drivers. The occurrence of these sharp phase jumps varies from phi-bit to phi-bit. All phi-bit phases also possess a common background dependency on the drivers’ phase. Within the context of multiple time scale perturbation theory, we develop a simple model of the nonlinear array of externally driven coupled acoustic waveguides to shed light on the possible mechanisms for the experimentally observed behavior of the logical phi-bit phase. Finally, we illustrate the ability to experimentally initialize the state of single-and multiple-phi-bit systems by exploiting the drivers’ phase as a tuning parameter. We also show that the nonlinear correlation between phi-bits enables parallelism in the manipulation of two-and multi-phi-bit superpositions of states.


Introduction
A "phi-bit" [1] is a two-state degree of freedom in an acoustic wave (an acoustic spin), which can be in a coherent superposition of states with complex amplitude coefficients.Phi-bits are, therefore, analogs of qubits.So far, phi-bits have been theoretically, computationally, and experimentally studied in the externally driven nonlinear arrays of coupled acoustic waveguides [2].Acoustic quantum analogs readily extend beyond the qubit paradigm into qutrit and even qudit analogs [3].Subsequently, we theoretically, computationally, and experimentally (demonstrated) introduced the notion of a logical phi-bit [4].A logical phi-bit is a two-state degree of freedom in the spectral domain of nonlinear acoustic modes supported by the externally driven arrays of coupled waveguides.We demonstrated the exponentially complex and scalable Hilbert space of the states of multiple logical phi-bits (≤16 with 2 16 -dimensional space) and the nonseparability of coherent superpositions [4].The applicability of logical phi-bits to the development of a classical, acoustic-based quantum-inspired information processing platform necessitates that phi-bits and their supporting array of waveguides satisfy the DiVincenzo five criteria for the physical construction of a quantum computer [5].
Criterion 1: "A scalable physical system with well-characterized qubit".Logical phi-bits live within the same driven nonlinear array of coupled acoustic waveguides.They occupy the same real estate-the phi-bit physical system does not need to scale physically.Since phi-bits are acoustic modes in the spectral domain with a bandwidth of a few tens of Hz, there is a lot of room in the ultrasonic domain for numbers of phi-bits exceeding 50, which is considered to be the threshold for quantum advantage.
Criterion 2: "The ability to initialize the state of the qubits to a simple fiducial state".A logical phi-bit can be prepared in states spanning the Bloch sphere [4].Multiple phi-bits can be Quantum Rep. 2023, 5 326 prepared in classical entangled states (i.e., nonseparable superpositions of states) [4].That is, a multipartite system composed of N two-level correlated logical phi-bits can support coherent superpositions of nonlinear acoustic modes spanning exponentially complex Hilbert spaces of dimension 2 N [4].The complete Hilbert space of multi-phi-bit systems can be tiled through their representation (i.e., the choice of the basis) and by tuning the system's driving conditions, such as the drivers' frequency [6].
Criterion 3: "Long relevant decoherence times".Coherent superpositions of phi-bit states are complex acoustic amplitudes.These states are not subjected to decoherence, in contrast to the probability amplitude of quantum systems, which are composed of superconducting or trapped ion qubits.
Criterion 4: "A 'universal' set of quantum gates".We experimentally demonstrated the phase, Hadamard [7], and C-NOT gates [8] that form the components of a universal set.We further showed a three-phi-bit gate that would challenge qubit-based computing platforms [9].
Criterion 5: "A qubit-specific measurement capability".Measurements in conventional quantum computing approaches, including photonic platforms [10], entail the collapse of the wavefunction.This is the destruction of the coherent superposition of states that is at the heart of the exponential scaling that leads to quantum advantage.In contrast, measurements in phi-bit systems do not alter the wave function, which persists as long as the drivers are applied.
An understanding of the nonlinear phenomena that underlie the manipulation of and correlations between phi-bits is critical to the use of phi-bits as qubit analogs for the development of quantum-analog computing platforms that can support and navigate scalable, exponentially complex Hilbert spaces.The navigation of phi-bit states in multiple correlated phi-bit Hilbert spaces is necessary to initialize (criterion 2) and operate (criterion 4) on the state of phi-bits in ways that are analogous to quantum initialization, operations, and gates.This is carried out by exploiting the correspondence between the state of correlated logical phi-bits represented in a low-dimensional linearly scaling physical space and their state representation in a high-dimensional, exponentially scaling Hilbert space.The navigation necessitates control of the state of logical phi-bits in response to parametric changes in the physical system that supports them.In Ref. [6], we utilized the frequency of the drivers to parametrically explore the phi-bit states.Here, we experimentally demonstrate the manipulation of logical phi-bits states by tuning the relative phase of the external drivers applied to the nonlinear array of coupled acoustic waveguides.Within simple models of nonlinear arrays of externally driven coupled acoustic waveguides, we also propose the possible mechanisms for the experimentally observed behaviors.

Materials and Methods
Logical phi-bits are supported by a nonlinear metamaterial composed of coupled elastic waveguides [4].The waveguides are aluminum rods arranged in a linear array with epoxy-filled lateral gaps (Figure 1).Ultrasonic longitudinal contact transducers drive and detect the acoustic field at the rod ends.Separate waveform generators excite the top and bottom driving transducers with sinusoidal signals of the same magnitude at frequency f 1 = f 3 .Another signal generator drives the middle transducer at the frequency f 2 .With the same magnitude as the outer transducers.The relative phase, ∆θ, between the signal generators' exciting transducers 1 and 3 is controlled and can be tuned over the entire range [0-360 • ].Three detecting transducers at the opposite ends collect data on the displacement field, providing temporal and spatial information across the array.Details of the experimental system are given in Ref. [6].Temporal Fourier transforms of output signals generate spectral information.The displacement field measured at the waveguide's detection ends is the Fourier sum of modes with the primary frequencies f 1 and f 2 as well as secondary nonlinear modes, the frequencies of which are a linear combination of the driving frequencies: p f 1 + q f 2 , where p and q are integers.The measured displacement field of a nonlinear mode can be expressed in the compact form: where ω i = 2π f i ; i = 1, 2 is the angular frequency and C j is the magnitude of the displacement of the waveguide j; j = 1, 2, 3.The phases ϕ 1 , ϕ 2 , and ϕ 3 arise from the complex nature of the resonant amplitudes of the driven system as well as any nonlinear effect combining these resonant amplitudes and the drivers' phases.Equation ( 1) can be reformulated by normalizing the amplitude to that of the first waveguides and expressing the field in terms of phase differences only: where Ĉ2 and Ĉ3 are normalized to C 1 and ϕ 12 = ϕ 2 − ϕ 1 and ϕ 13 = ϕ 3 − ϕ 1 .We reduce this expression by dropping the direct reference to waveguide 1.We now define the displacement field at the end of the waveguides by the renormalized 2 × 1 vector: = ĉ2 e iϕ 12 ĉ3 e iϕ 13 e i(pω We redefine the state of the logical phi-bit {p, q} in terms of phase difference of the displacement field only by constructing the nontemporal part of the field as the normalized complex amplitude state vectors: This state vectors live in single phi-bit 2D Hilbert space, with a basis of 1 0 , 0 The quantities e iϕ 12 and e iϕ 13 represent the complex coefficient in a linear combination of the 2D basis vectors.
We can now define a logical phi-bit as a two-level, secondary nonlinear mode of vibration, the state of which is characterized by the set of frequency coefficients {p, q} and the spatial mode associated with two independent relative phases of the displacement between the waveguides.
The mixed frequency, pω 1 + qω 2 , serves as a good quantum number for defining a logical phi-bit.The amplitude and phases at the waveguide ends can be measured unambiguously.A single phi-bit state lives in a 2D Hilbert space h (j) .Employing Dirac's ket notation to represent the basis vectors of that space, Equation (4) reduces to A single phi-bit state represented in this form spans the Bloch sphere and is analogous to a quantum bit (qubit).Equation ( 5) is effectively a coherent superposition of |0 and |1 states with complex amplitudes, propose the possible mechanisms for the experimentally observed behaviors.

Materials and Methods
Logical phi-bits are supported by a nonlinear metamaterial composed of couple elastic waveguides [4].The waveguides are aluminum rods arranged in a linear array wi epoxy-filled lateral gaps (Figure 1).In (b), we also show the schematic of the experimental system for generating and detecting logical phi-bits.This includes separate signal generators and amplifiers that are used to drive piezoelectric transducers, driving and detecting the transducers attached to the opposite ends of the waveguides by the pressure of three independent rubber bands.A thin layer of honey is used as an ultrasonic coupling agent between the transducers and the waveguide ends.The detected signals enter an oscilloscope via independent input channels for analysis.The waveguides are suspended by thin threads for isolation.

Results
We have conducted an experiment with the system in Figure 1 by driving waveguides 1 and 3 at the same frequency f 1 = f 3 = 62 kHz and the middle waveguide 2 at a frequency of f 2 = 66 kHz.The reason for choosing these driving frequencies is that, at these frequencies, the experimental aluminum rod's (McMaster-Carr 1615T172: diameter = 1/2 inch, length = 0.6096 m, and density = 2660 kg/m 3 ) longitudinal waves have a wavelength of around 10 cm, making propagation along rod-like waveguides almost one-dimensional [2].Additionally, the transducers provide appropriate driving and detecting amplitudes at these frequencies, and the linear longitudinal modes of finite-length waveguides are welldefined [2].Initially, drivers 1 and 3 are in phase, and their relative phase, ∆θ, is increased by increments of 2.5 • up to 360 • .The displacement field at the detection ends of the waveguides is Fourier transformed.The Fourier spectrum includes primary peaks at the two frequencies f 1 and f 2 and secondary peaks at multiple frequencies, p f 1 + q f 2 , corresponding to nonlinear phi-bit modes.The complex amplitude of each phi-bit mode at the ends of the three waveguides is used to calculate the phase differences ϕ 12 and ϕ 13 .In Figure 2, we illustrate these phase differences for two phi-bit modes, namely a phi-bit (a) with nonlinear frequency f a = 4 f 1 − 2 f 2 ({p = 4, q = −2}) and a phi-bit (b) with frequency f b = 4 f 1 − f 2 ({p = 4, q = −1}).In addition to the response to ∆θ of the phase differences ϕ 12 and ϕ 13 , from phi-bits (a) and (b), we have also measured the phases ϕ 12 ( f 1,2 ) and ϕ 13 ( f 1,2 ) for the primary modes observed at the frequencies f 1 and f 2 .We have also calculated the quantities ϕ 0 12 (a) = 4ϕ 12 ( f 1 ) − 2ϕ 12 ( f 2 ) and ϕ 0 13 (a) = 5ϕ 13 ( f 1 ) − 2ϕ 13 ( f 2 ), as well as . These quantities (ϕ 0 12 and ϕ 0 13 , with the superscript 0) would represent the phase differences of the phi-bits if these were simple linear combinations of ϕ 12 and ϕ 13 in the primary linear modes.The response of the two phi-bits is composed of two separate sets of features.ϕ 12 and ϕ 13 of both phi-bits (a) and (b) follow the trend of ϕ 0 12 (a), ϕ 0 12 (b), ϕ 0 13 (a), and ϕ 0 13 (b) as functions of ∆θ.The variations following the linear combinations of primary mode phases will be subsequently called backgrounds.Sharp phase jumps that happen within a few degrees overlap with the backgrounds.These phase jump amount to less than 180 • .Similar behaviors are observed for other phi-bit nonlinear modes.
These experimental results broach the challenging questions of the origin of the experimentally observed trends and jumps in the phases ϕ 12 and ϕ 13 ; this is in the context of enabling manipulation of phi-bit states (Equation ( 5)) within the Bloch sphere.The subsequent section develops a perturbative model of a nonlinear array of externally driven acoustic waveguides to offer possible answers to this question.

Model of Nonlinear Logical Phi-Bit and Effect of Drivers' Phase on Phi-Bit State Vector
We consider a model of the acoustic metamaterial composed of three one-dimensional elastic waveguides coupled elastically along their length (Figure 3).Each waveguide is driven externally at its end at the position  0.  represents the position along the waveguides.The nonlinear elastic wave equation in the long wavelength limit is written as The parameter β is proportional to the speed of sound along the waveguides.The parameter  represents damping. ⃡ is the identity matrix. measures the elastic coupling strength between waveguides due to epoxy. ⃖ ⃗ is the matrix characterizing the elastic coupling between the three waveguides.In the case of our planar array of waveguides, the coupling matrix takes the form ⃑ ,  ⃑ , and  ⃑ are 3 × 1 vectors representing the external driving harmonic forces for the three different driving angular frequencies  2 ,  2 , and  2 . ,  , and  are the phases of the driving forces.
The displacement in waveguides 1, 2, and 3 is represented by the 3 × 1 vector  ⃑  ,  ,  . ⃑  ⃑ is a nonlinear term with strength .

Model of Nonlinear Logical Phi-Bit and Effect of Drivers' Phase on Phi-Bit State Vector
We consider a model of the acoustic metamaterial composed of three one-dimensional elastic waveguides coupled elastically along their length (Figure 3).Each waveguide is driven externally at its end at the position x = 0. x represents the position along the waveguides.

Model of Nonlinear Logical Phi-Bit and Effect of Drivers' Phase on Phi-Bit State Vector
We consider a model of the acoustic metamaterial composed of three one-dimensional elastic waveguides coupled elastically along their length (Figure 3).Each waveguide is driven externally at its end at the position  0.  represents the position along the waveguides.The nonlinear elastic wave equation in the long wavelength limit is written as The parameter β is proportional to the speed of sound along the waveguides.The parameter  represents damping. ⃡ is the identity matrix. measures the elastic coupling strength between waveguides due to epoxy. ⃖ ⃗ is the matrix characterizing the elastic coupling between the three waveguides.In the case of our planar array of waveguides, the coupling matrix takes the form ⃑ ,  ⃑ , and  ⃑ are 3 × 1 vectors representing the external driving harmonic forces for the three different driving angular frequencies  2 ,  2 , and  2 . ,  , and  are the phases of the driving forces.
The displacement in waveguides 1, 2, and 3 is represented by the 3 × 1 vector  ⃑  ,  ,  . ⃑  ⃑ is a nonlinear term with strength .The nonlinear elastic wave equation in the long wavelength limit is written as The parameter β is proportional to the speed of sound along the waveguides.The parameter µ represents damping.
↔ I is the identity matrix.α measures the elastic coupling strength between waveguides due to epoxy.
↔ M is the matrix characterizing the elastic coupling between the three waveguides.In the case of our planar array of waveguides, the coupling matrix takes the form F 1 , F 2 , and F 3 are 3 × 1 vectors representing the external driving harmonic forces for the three different driving angular frequencies ω 1 = 2π f 1 , ω 2 = 2π f 2 , and ω 3 = 2π f 3 .θ 1 , θ 2 , and θ 3 are the phases of the driving forces.
The displacement in waveguides 1, 2, and 3 is represented by the 3 × 1 vector We seek analytical approximations to the nonlinear Equation ( 6) with three frequency excitations.For this, we consider a variation on that equation that enables us to use frequency detuning with multiple time scale perturbation theory [11,12].
Note the additional parameter δ and the dependency of damping on ε.A factor of 2 for multiplying the damping parameter is introduced in Equations ( 6) and ( 8) to simplify the derivation that follows.The function, G( U), models the nonlinearities that act on the waveguides along their length.It may, therefore, qualitatively represent the effect of the nonlinear elasticity of epoxy.
We introduce two time scales: τ 0 = t and τ 1 = εt.We also expand the displacement field as the sum of the zero-order (linear) and first-order (nonlinear) terms, as follows: (τ 0 , τ 1 ).The first-order and second-order time derivatives take the form The wave equation to the 0th order in ε is effectively the linear equation: (10)   Note that there is no damping coefficient in Equation (10) as the effect of damping is now included in the first-order equation.We will see later that the phase associated with damping (included now in the first-order equation) will come back as a correction to the complex amplitude of the zero-order solution.
We can solve this equation by defining λ n and E n with n = 1, 2, 3: the eigenvalues and eigen vectors of the ↔ M matrix, where E n represents the spatial eigenmodes across the waveguides with components E n,j , j = 1, 2, 3. We write The eigenmodes have eigenvalues λ 1 = 0, λ 2 = 1 and λ 3 = 3, and are given by We now expand the displacement vector on the complete orthonormal basis, { E n }: The 3 × 1 vectors, F l , are also expanded on the basis { E n }: Inserting Equations ( 11)-(13) into Equation ( 10) yields a set of three equations, giving the form of By employing plane wave expansions, the solutions to the homogeneous equation take the form of with the dispersion relation: k jn is the jth wave number for mode n.The star in Equation (15) stands for the complex conjugate, and the superscript h stands for the homogeneous solution.The summations in Equation (15) cover a discrete set of wave numbers since the waveguides have a finite length.
The particular solutions of Equation ( 14) that follow the temporal dependency of the driving forces are found in the form are real resonant amplitudes.In Equation ( 16), the superscript p stands for a particular solution.In order to obtain Equation (16), we have expanded the cosine functions into complex exponentials.
The complete solution, therefore, takes the form To obtain a first order in perturbation, the wave equation is given by ) Note that the external driving force does not appear in Equation ( 18) as its full contribution was accounted for in the zeroth order equation.Here, the zeroth order solution now serves as a driver for the first-order equation.
We now choose a form for the nonlinear term G( U (0) ) that enables us to proceed analytically and illustrate the effect of the phases θ 1 , θ 2 , and θ 3 on nonlinear modes: This form assumes that the spatial modes, E n , do not interact with each other.However, for each spatial mode, the plane wave modes, e ik jn x , may interact with each other.We use a third order nonlinearity for the sake of analytical tractability.

Defining U
(1) = ∑ 3 n=1 E n u 1,n and using Equation ( 14), we can rewrite Equation (18) as a set of three equation, each one corresponding to a different spatial mode: The terms on the right-hand side of Equation ( 20) can lead to secular behavior.In order to calculate (u 0,n ) 3 , we rewrite equation (17) as with ξ n = A n e iω 0,n τ 0 + Λ n,1 e i(ω 1 τ 0 +θ 1 ) + Λ n,2 e i(ω 2 τ 0 +θ 2 ) + Λ n,3 e i(ω 3 τ 0 +θ 3 ) .With this, we rewrite the cubic term in Equation ( 20) as In this expression, since we have the product of three summations over the wave numbers, we introduced three different wave numbers, k jn , k jn , k jn .However, to simplify the problem, we investigate the self-interaction of one single mode k jn .We will more carefully examine two types of terms from Equation ( 22), those that contain the resonant frequencies of the system, ω 0,n , and those corresponding to the driving frequencies, ω n .Beginning with those involving resonant frequencies, we seek terms in Equation ( 22) with a time dependency given by e iω 0,n τ 0 .These terms will contribute to the secular behavior of Equation (20) and are grouped as n A * n e ik jn x e iω 0,n τ 0 + 3A n e ik jn x e iω 0,n τ 0 2 Furthermore, we seek terms in e ik jn x but without a time dependency of the form e iω 0,n τ 0 , which correspond to a mixing of the driving frequencies and phases.We focus on terms mixing the frequency and phase of driver 2 with either drivers 1 or 3, which we group as n,2 Λ n,1 e ik jn x e i((2ω 2 −ω 1 )τ 0 +2θ 2 −θ 1 ) + Λ 2 n,2 Λ n,3 e ik jn x e i((2ω 2 −ω 3 )τ 0 +2θ 2 −θ 3 ) ) In order to make contact with the experimental conditions, we now take ω 1 = ω 3 and θ 1 = −∆θ/2, θ 3 = +∆θ/2, where ∆θ is the phase difference between drivers 1 and 3.
3.1.1.Case I: n = 1 or 3 In this case, we rewrite Equation (24) in the simpler form: We now introduce a detuning parameter: σ = 1 ε (2ω 2 − ω 1 − ω 0,n ) and redefine T 2 as Λ n,1 e ik jn x e iω 0,n τ 0 e i(στ By regrouping all the terms contributing to secular behavior and setting them to zero, we obtain the relation Using u 0,n = u h 0,n + u p 0,n , Equations ( 15) and ( 16) for the wave number k n , and Equations ( 23) and ( 26), after some algebraic manipulations, we obtain where In Equation ( 28), the primed quantity is a derivative with respect to τ 1 .We have also chosen θ 2 = 0 without consequence on the rest of the derivation.We now take A n (τ 1 ) = 1 2 a(τ 1 )e ib(τ 1 ) .We also define: η = στ 1 − b such that b = σ − η .By inserting these into Equation (28) and considering only steady state behavior (that is, a = 0 and η = 0), the regrouping of the real and imaginary terms into two separate equations results in We can get the amplitude-frequency response of the specific nonlinear mode, 2ω 2 − ω 1 due to self-interaction by eliminating the phase η using the trigonometric relation sin 2 η + cos 2 η = 1.That response takes the form of Recall that a is the correction to the amplitude A n k jn solution of the homogenous zeroth order wave equation as a result of the nonlinear perturbation.In some frequency ranges, σ, the amplitude frequency response may not be a single-valued function.The amplitude may show overhangs to a high frequency or low frequency depending on the sign of δ.It is centered on the backbone curve given by σ = δΓ 1 + 3δ 8ω 0,n a 2 , where the amplitude decays to zero when σ = δΓ 1 .The maximum value of the amplitude is obtained by setting the square root term to zero, . The maximum amplitude is, therefore, dependent on the phase difference between drivers 1 and 3.
In Figure 4, we note that the amplitude of the nonlinear mode depends on the phase difference between the two drivers, 1 and 3. We note that for ∆θ = 180 • , the amplitude a = 0 and the spatial modes E 1 and E 3 do not contribute to the nonlinear correction of the solution of Equation (15).In particular, we note that, at a fixed frequency, σ~3.1 (vertical line in the figure), changing the drivers' phase difference ∆θ from 0 • to 90 • reduces the number of solutions for the amplitude from 3 to 1.This dependency will impact the behavior of the phase, η.We obtain the phase from Equation (29a) as sin η = sgn(δΓ 2 ) a a m .In Figure 4, we illustrate the effect of a change in the drivers' phase on η as a function of frequency.Let us suppose that we prepared the nonlinear mode shown by the closed circle in Figure 5.This is carried out by setting ∆θ = 0 and by slowly increasing the frequency σ from −2 to 3.1, i.e., by moving along the upper part of the solid line amplitude-frequency response of  Under this condition, the amplitude-frequency response is single value and corresponds to the state illustrated as an open circle in Figure 5.This state has a single value ~0.The change in the driver phase difference leads to a near 130° change in the phase of the nonlinear mode.This jump in phase carries to  ⃑ through  and  , .Indeed, we can now rewrite  in the form Here, we are using indices 1 and 3 to specify that the amplitude and phase correspond to the spatial modes  ⃑ and  ⃑ .Inserting the definition of:  , 2   , , ,   and  , we rewrite Equation (31) as For illustrative purposes, we have taken 3δ 8ω 0,n = 0.5, (δΓ 2 ) 2 = 0.008, and µ 2 = 0.00425.We consider three values for ∆θ = 0, 90, 150   Under this condition, the amplitude-frequency response is single value and corresponds to the state illustrated as an open circle in Figure 5.This state has a single value ~0.The change in the driver phase difference leads to a near 130° change in the phase of the nonlinear mode.This jump in phase carries to  ⃑ through  and  , .Indeed, we can now rewrite  in the form Here, we are using indices 1 and 3 to specify that the amplitude and phase correspond to the spatial modes  ⃑ and  ⃑ .Inserting the definition of:  , 2   , , ,   and  , we rewrite Equation (31) as Under this condition, the amplitude-frequency response is single value and corresponds to the state illustrated as an open circle in Figure 5.This state has a single value η ∼ 0. The change in the driver phase difference leads to a near 130 • change in the phase of the nonlinear mode.This jump in phase carries to U (0) through A n and u 0,n .Indeed, we can now rewrite A n in the form Here, we are using indices 1 and 3 to specify that the amplitude and phase correspond to the spatial modes E 1 and E 3 .Inserting the definition of: ), τ 0 = t and τ 1 = εt, we rewrite Equation (31) as Considering a single k jn and θ 2 = 0 in Equation ( 21) and focusing our attention on the complex term (keeping aside the complex conjugate in that equation), we obtain u 0,n=1,3 = A n=1,3 e iω 0,n=1,3 τ 0 + 2Λ n=1,3,1 e iω 1 τ 0 2 cos ∆θ 2 + Λ n=1,3,2 e iω 2 τ 0 e ik jn=1,3 x (33) Inserting Equation (32) into Equation (33) yields u 0,n=1,3 = a 1,3 (σ 1,3 )e −iη 1,3 e i(2ω 2 −ω 1 )τ 0 + 2Λ n=1,3,1 e iω 1 τ 0 2 cos ∆θ 2 + Λ n=1,3,2 e iω 2 τ 0 e ik jn=1,3 x (34) The last two terms correspond to the primary linear modes observed experimentally.The first term corresponds to one of the possible nonlinear phi-bit modes, namely {p = 2, q = −1}.

Case II: n = 2
In this case, Equation (24) becomes Λ n,1 e ik jn x e iω 0,n τ 0 e i(στ 1 +2θ 2 ) 2i sin The secular Equation ( 27) yields where again and 2 a(τ 1 )e ib(τ 1 ) and defining η = στ 1 − b such that b = σ − η , Equation (36) becomes at steady state (i.e., a = 0 and η = 0) Equations (37a) and (37b) result from the imaginary and parts of Equation (36), respectively.Note that in case II, the effect of the phase difference ∆θ occurs through a sine function instead of a cosine function, as was the case in case I. Additionally, note the swapping between sin η and cos η compared to case I.
We can obtain the amplitude-frequency response of the specific nonlinear mode, 2ω 2 − ω 1 , for mode n = 2 due to self-interaction by eliminating phase η using the trigonometric relation: sin 2 η + cos 2 η = 1.That response takes the form of This equation has the same backbone curve as that of case 1.The only difference resides in the way the amplitude varies with the phase difference ∆θ through sin ∆θ 2 instead of cos ∆θ 2 .Now, for ∆θ = 0, that is sin ∆θ 2 = 0, Equations (37a) and (37b) give a = 0. Spatial mode 2 does not contribute to the nonlinear displacement field.Case II mirrors case I with a π phase lag in ∆θ.The second spatial mode contributes to phase jumps in η with, again, a π phase lag in ∆θ.Considering a single k jn in Equation ( 21) for spatial mode E 2 , we obtain u 0,n=2 = a 2 (σ 2 )e −iη 2 e i(2ω 2 −ω 1 )τ 0 + 2Λ n=2,1 e iω 1 τ 0 2i sin ∆θ 2 + Λ n=2,2 e iω 2 τ 0 e ik jn=2 x (39) Again, in Equation (39), we apply a subscript of 2 to the amplitude a and phase η to highlight the fact that, in case II, one deals with the spatial mode n=2.
Let us now combine cases I and II.The first-order correction (due to the self-interaction of the k jn wave number to the complete displacement field of Equation (12), In the Fourier spectrum of the displacement field, the displacement associated with the nonlinear logical phi-bit mode with a frequency of 2ω 2 − ω 1 measured at x = 0 is given by The driver phase difference, ∆θ, can, therefore, be used through its impact on η 1 , η 2 , and η 3 to tune the complex amplitudes 1   √ 2 e iϕ 12 and 1 √ 2 e iϕ 13 of a logical phi-bit.These amplitudes may exhibit large phase jumps, which are associated with the phase jumps in η 1 , η 2 , and η 3 , which is a behavior that is consistent with the experimental observation of the sharp phase jumps reported in Section 3.
We now address the possible origin of the phi-bit phase background as a function of ∆θ.Let us consider the particular solution part of Equation (40) evaluated at x = 0: In the Fourier spectrum of the total displacement field, the two terms in Equation (42) correspond to the primary modes.We rewrite Equation (42) by introducing In Equation ( 43), we have replaced τ 0 with t.As seen in Figure 1, the array of coupled acoustic waveguides is only part of a system that includes transducers and their ultrasonic coupling agent at the ends of the waveguides.These system components may possess nonlinear as well as damping characteristics, which will be able to mix the frequencies associated with the linear displacement of Equation (43).Indeed, let us consider a model of the transducer/coupling agent component of the system in the form of extrinsic nonlinear damped oscillators at the waveguide ends driven by the displacement U (0) p (ω 1 , ω 2 ): where ω R and µ R are the characteristic frequency and damping coefficient, respectively.
G is the nonlinear function with a strength ε R .
V 3 ) where the V i ; i = 1, 2, 3, are displacement degrees of freedom associated with the extrinsic oscillators at the ends of rods, i = 1, 2, 3. K is some proportionality constant converting displacement into force.Note that since the transducers/coupling agent at the ends of the waveguides are independent of each other, the components, V i , in Equation (44) are also independent of each other.So we can rewrite Equation (44) in the form of three independent equations: In Equation (45), we have chosen the nonlinear function, G i , in the form of a power V Q i .We can solve Equation (45) using perturbation theory.We seek solutions in the form: i .Regarding the zeroth order in ε R , the equation reduces to a linear form: By seeking zeroth order harmonic solutions with the same frequencies as the driving displacement, we obtain Regarding the first order in ε R , Equation (45) becomes Here, the zeroth order solution drives through the nonlinear term to the first-order solution.Inserting Equation (47) into the nonlinear Equation (48) leads to a series of terms with the general form of , where D is a proportionality constant and p + q = Q.By focusing on one of these terms, corresponding to one of the logical phi-bits with {p, q}, Equation (48) yields A particular solution of Equation (49) may be written as The fraction prefactors in Equation ( 50) are the same for each component "i" of and add the same general phase of each component.Therefore, we have Note that the coefficients (C i (ω 1 , ∆θ)) p (C i (ω 2 , ∆θ)) q are real.The Fourier spectrum of the displacement detected by the transducers will include the contribution of Equation (51) to the phase of the nonlinear phi-bit mode, {p,q}.These phases are linear combinations of the phase of primary linear modes at frequencies ω 1 and ω 2 with linear coefficients p and q.They will, therefore, contribute the background of the phases of a phi-bit that can be determined from the phases of the primary modes, as was shown it Section 3. Here, the background in the logical phi-bit phase clearly depends on the parameter ∆θ in accordance with the experimental observation.However, to gain full insight into the effect of ∆θ on the phi-bit phase, one needs to go beyond the particular solution of Equation (42).Indeed, this solution does not include the effect of damping as it was lumped into the perturbation of Equation ( 8) in order to use frequency detuning with multiple time scale perturbation theory and shed light on the nonlinear amplitude-frequency response of a phi-bit.In order to illuminate the dependency of the phase of the primary modes, that is, the background of the phi-bit phase, on ∆θ, we now reconsider the linear version of Equation ( 6), with a single driving frequency, ω, applied to the three waveguides with different phases, namely We simplify this equation by working with a complex exponential instead of a cosine function for the driving force.That is, we reformulate the driving force on the right-hand side of Equation (52) as where we have expressed the phase difference between drivers 2 and 3 relative to driver 1 as θ 12 and θ 13 .The coefficients A, B, and C are the magnitudes of the forces.We define C = rA.Note that our experiment is conducted with r = 1.However, because we use different transducers for each waveguide, that condition may not be exactly realized.
Expanding → F on the basis of spatial modes E 1 , E 2 , and E 3 gives the three components: We also expand the left-hand side of Equation (52) on the basis of spatial modes; that is, we define U = ∑ n u n E n .Equation (52) takes the form of three equations, each one corresponding to a spatial mode.We seek solutions for u n in the form of superpositions of plane waves.The complex resonant amplitudes are obtained for each spatial mode as , and , where ω 2 0,n k jn = β 2 k 2 jn + α 2 λ n for n = 1, 2, 3, wave number, k jn , and spatial mode eigenvalue, λ n , as was defined earlier.Inserting Equations (54a)-(54c) into the resonant amplitudes and determining their corresponding phase, as defined by r sin θ 13 1 + r cos θ 13 (55a) These phases are independent of θ 12 .We can now redefine θ 13 = ∆θ.We note that as ∆θ passes through 180 • for r = 1, the phases φ (1) = φ (3) undergo a π jump.There is no such jump in φ (2) .In order to see how such a jump would carry to a phase jump in the phase of the primary linear mode and subsequently to the background of a phi-bit, we write the complete solution for U: reflect only the sharp phase jumps associated with changes in ∆θ.For the sake of illustration, we reduce the complex representation of the phi-bit states to their real part, namely In Figure 6, we report the components of ψ r (a) and ψ r (b) as functions of the drivers' phase tuning parameter ∆θ.
Figure 6 illustrates the ability to initialize the system composed of the two phi-bits (a) and (b) into all possible two-phi-bit states, which are the tensor products of selected two single phi-bit states.The single phi-bit states lie in a two-dimensional Hilbert space.The Hilbert space of the two phi-bit states has exponentially scaled dimensions: 2 2 .By driving the physical system with ∆θ ∼ 122 o , one can realize the two-phi-bit state labeled (1): . Similarly, at ∆θ ∼ 157 o , 187 o , and 207 o , one can initialize the system into the states: (2): The two phi-bit states are constructed as the tensor products of single-phi-bit states and are, therefore, separable (i.e., not classically entangled).However, entanglement is not necessary to correlate the two phi-bits, as would be the case with two qubits, in order to enable parallelism in two-bit operations.Here, the two phi-bits (a) and (b) are correlated via the nonlinearity of the elasticity of the physical system.For phi-bits, classical entanglement would only be needed to initialize two different phi-bit states.This can be achieved with the current data by using a different representation (i.e., basis) of the two phi-bit states.We can see that through a single action on the nonlinear array of coupled acoustic waveguides, that is, tuning the drivers' phase difference ∆θ, one can parametrically change the states of phi-bits (a) and (b) and, therefore, simultaneously change the four components of the two phi-bit state vector.This action is actually a controllable operation on two phi-bit states.
For instance, changing ∆θ from 122 o to 157 o results in the following unitary transformation acting on state (1) and producing state (2) to within a global phase of π: This unitary transformation is a phase gate that adds π to the second and fourth components of the two-phi-bit state vector.Changing ∆θ from 157 o to 187 o adds a phase of π to the second and third components of state vector (2) and produces state vector (3).Changing ∆θ from 187 o to 207 o adds a phase of π to the first and third components of state vector (3) to produce state vector (4).This is only a simple illustration of how one can exploit the sharp phase jumps in logical phi-bit phases associated with the parametric variation of the relative phase between drivers in a nonlinear array of acoustic waveguides.Here, the difference in the phase response of the two phi-bits provides the complexity necessary to realize a complete set of two phi-bit states.This approach is scalable, and it can be applied readily to more than many logical phi-bits (which exist effectively in the spectral domain) without increasing the footprint of the physical system.
For the sake of illustration, we reduce the complex representation of the phi-bit states to their real part, namely In Figure 6, we report the components of   and   as functions of the drivers' phase tuning parameter Δ. Figure 6 illustrates the ability to initialize the system composed of the two phi-bits (a) and (b) into all possible two-phi-bit states, which are the tensor products of selected two single phi-bit states.The single phi-bit states lie in a two-dimensional Hilbert space.The Hilbert space of the two phi-bit states has exponentially scaled dimensions: 2 2 .By driving the physical system with Δ~122 , one can realize the two-phi-bit state labeled (1): . Similarly, at Δ~157 , 187 , and 207 , one can initialize the system into the states: (2): and (4): , respectively.The

Example of Initialization Using Phi-Bit Background Phase
As seen in Figure 7 and observed in Section 2, a background in the phi-bit phase response is common to all phi-bits.This background arises from the phases of the primary mode at the frequency f 1 = 2π/ω 1 : ϕ 12 ( f 1 ) and ϕ 13 ( f 1 ).Since ∆θ only applies to drivers 1 and 3, which are driven at f 1 , the phases of the primary modes at f 2 = 2π/ω 2 , ϕ 12 ( f 2 ) and ϕ 13 ( f 2 ) remain constant.Correlations between phi-bits make their background phases react simultaneously to a change in ∆θ in the same way.
two phi-bit states are constructed as the tensor products of single-phi-bit states and are, therefore, separable (i.e., not classically entangled).However, entanglement is not necessary to correlate the two phi-bits, as would be the case with two qubits, in order to enable parallelism in two-bit operations.Here, the two phi-bits (a) and (b) are correlated via the nonlinearity of the elasticity of the physical system.For phi-bits, classical entanglement would only be needed to initialize two different phi-bit states.This can be achieved with the current data by using a different representation (i.e., basis) of the two phi-bit states.We can see that through a single action on the nonlinear array of coupled acoustic waveguides, that is, tuning the drivers' phase difference Δ, one can parametrically change the states of phi-bits (a) and (b) and, therefore, simultaneously change the four components of the two phi-bit state vector.This action is actually a controllable operation on two phibit states.For instance, changing Δ from 122 to 157 results in the following unitary transformation acting on state (1) and producing state (2) to within a global phase of π: This unitary transformation is a phase gate that adds π to the second and fourth components of the two-phi-bit state vector.Changing Δ from 157 to 187 adds a phase of π to the second and third components of state vector (2) and produces state vector (3).Changing Δθ from 187 to 207 adds a phase of π to the first and third components of state vector (3) to produce state vector (4).This is only a simple illustration of how one can exploit the sharp phase jumps in logical phi-bit phases associated with the parametric variation of the relative phase between drivers in a nonlinear array of acoustic waveguides.Here, the difference in the phase response of the two phi-bits provides the complexity necessary to realize a complete set of two phi-bit states.This approach is scalable, and it can be applied readily to more than many logical phi-bits (which exist effectively in the spectral domain) without increasing the footprint of the physical system.

Example of Initialization Using Phi-Bit Background Phase
As seen in Figure 7 and observed in Section 2, a background in the phi-bit phase response is common to all phi-bits.This background arises from the phases of the primary mode at the frequency  2  ⁄ :   and   .Since Δ only applies to drivers 1 and 3, which are driven at  , the phases of the primary modes at  2  ⁄ ,   and   remain constant.Correlations between phi-bits make their background phases react simultaneously to a change in Δ in the same way.When considering (for illustrative purposes) N phi-bits with the same background phases, (we have experimentally verified that in addition to phi-bits (a) and (b) presented ⎝ 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 ⎠ This is only an example of the use of representation as a means of initializing multiple phi-bit states.While we are here, we have exploited the background phase of the phi-bit, and we can also use the nonlinear phase jumps in combination with representations to add complexity to the multi-phi-bit states that can be realized.Note that since the phi-bits are nonlinearly correlated, an initializing operation, such as  ⃡ , can be achieved experimentally with a single action of tuning the drivers' phase difference.

Conclusions
We demonstrated how one could navigate the Hilbert space of two logical phi-bits by parametrically changing the physical system that supports them.We experimentally showed the manipulation of logical phi-bits state vectors by tuning the relative phase of the external drivers applied to a nonlinear array of coupled acoustic waveguides.We observed sharp phase jumps in the phi-bit states as a result of tuning the relative phase between drivers.These sharp phase jumps amount to approximately 180°.While the occurrence of the sharp phase jumps varies from phi-bit to phi-bit, all phi-bit phases possess a common background that is dependent on the drivers' phase.We introduced simple models of the nonlinear arrays of externally driven coupled acoustic waveguides to shed light on the possible mechanisms that may underlie the experimentally observed behavior.Finally, by exploiting the sharp phase jumps, we illustrated the ability to initialize a system composed of two logical phi-bits into all possible two phi-bit states, which are tensor products of a selected subset of single phi-bit states.We also showed an approach to initialize exponentially complex states of multiple phi-bits by introducing multi-phi-bit representations and by tuning the phase background common to all phi-bits.The correlation between phi-bits resulting from the nonlinear elasticity of the physical system is critical to the parallelism of unitary operations on multiple phi-bit states.
This work provides only a simple demonstration of how phi-bits, which are acoustic, classical analogs of qubits, can bring us closer to realizing a classical acoustic physical One may notice the simple relations between some of these states, which may be viewed as resulting from nontrivial unitary transformations.
For instance, 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 This is only an example of the use of representation as a means of initializing multiple phi-bit states.While we are here, we have exploited the background phase of the phi-bit, and we can also use the nonlinear phase jumps in combination with representations to add complexity to the multi-phi-bit states that can be realized.Note that since the phi-bits are nonlinearly correlated, an initializing operation, such as ↔ T, can be achieved experimentally with a single action of tuning the drivers' phase difference.

Conclusions
We demonstrated how one could navigate the Hilbert space of two logical phi-bits by parametrically changing the physical system that supports them.We experimentally showed the manipulation of logical phi-bits state vectors by tuning the relative phase of the external drivers applied to a nonlinear array of coupled acoustic waveguides.We observed sharp phase jumps in the phi-bit states as a result of tuning the relative phase between drivers.These sharp phase jumps amount to approximately 180 • .While the occurrence of the sharp phase jumps varies from phi-bit to phi-bit, all phi-bit phases possess a common background that is dependent on the drivers' phase.We introduced simple models of the nonlinear arrays of externally driven coupled acoustic waveguides to shed light on the possible mechanisms that may underlie the experimentally observed behavior.Finally, by exploiting the sharp phase jumps, we illustrated the ability to initialize a system composed of two logical phi-bits into all possible two phi-bit states, which are tensor products of a selected subset of single phi-bit states.We also showed an approach to initialize exponentially complex states of multiple phi-bits by introducing multi-phi-bit representations and by tuning the phase background common to all phi-bits.The correlation between phi-bits resulting from the nonlinear elasticity of the physical system is critical to the parallelism of unitary operations on multiple phi-bit states.

Figure 1 .
Figure 1.Picture (a) and exploded view (b) of the array of acoustic waveguides coupled with epo resin.In (b), we also show the schematic of the experimental system for generating and detecti

Figure 1 .
Figure 1.Picture (a) and exploded view (b) of the array of acoustic waveguides coupled with epoxy resin.In (b), we also show the schematic of the experimental system for generating and detecting logical phi-bits.This includes separate signal generators and amplifiers that are used to drive piezoelectric transducers, driving and detecting the transducers attached to the opposite ends of the waveguides by the pressure of three independent rubber bands.A thin layer of honey is used as an ultrasonic coupling agent between the transducers and the waveguide ends.The detected signals enter an oscilloscope via independent input channels for analysis.The waveguides are suspended by thin threads for isolation.

Figure 2 .
Figure 2. Measured  and  , of phi-bits (a) (left) and (b) (right) as functions of the drivers' phase difference Δ.The open circles are the experimental results, with the thin line serving as a guide for the eyes.The solid lines represent the quantities   4  2  ,   5  2  ,   4    , and   5    calculated from the primary modes.

Figure 3 .
Figure 3. Schematic of the metamaterial composed of a parallel array of three coupled waveguides.

Figure 2 .
Figure 2. Measured  and  , of phi-bits (a) (left) and (b) (right) as functions of the drivers' phase difference Δ.The open circles are the experimental results, with the thin line serving as a guide for the eyes.The solid lines represent the quantities   4  2  ,   5  2  ,   4    , and   5    calculated from the primary modes.

Figure 3 .
Figure 3. Schematic of the metamaterial composed of a parallel array of three coupled waveguides.

Figure 3 .
Figure 3. Schematic of the metamaterial composed of a parallel array of three coupled waveguides.
according to Equation (13), we need to consider two cases.Case I corresponds to the spatial modes n = 1, 3 for which the F (l) n are the same, leading to Λ n,1 = Λ n,3 .Case II corresponds to the spatial mode n = 2, for which Λ n,3 = −Λ n,1 .

Figure 4 10 Figure 4 .
Figure 4 up to an amplitude a~2.0.In that state, η ∼ 130 o .Increasing the phase difference between drivers 1 and 3 to ∆θ = 90 • reduces the maximum amplitude of the nonlinear mode (long dashed curve) in Figure 4, which now occurs at a frequency σ < 3. Quantum Rep. 2023, 5, FOR PEER REVIEW 10

Figure 5 .
Figure 5. Schematic illustration of the phase-frequency response of the forced array of coupled waveguides due to self-interaction.The phase is in degrees.At the fixed frequency ~3.1, a change in driver phase Δ of 90° (horizontal arrow) may take the nonlinear mode (closed circle) into the single-valued state shown as the open circle; that is, a change in  of approximately 130° (vertical arrow).

Figure 4 .
Figure 4. Schematic illustration of the amplitude-frequency response of the forced array of coupled waveguides due to self-interaction.The frequency  2   , is normalized to Γ .For illustrative purposes, we have taken , 0.5, Γ 0.008, and  0.00425.We con-

Figure 5 .
Figure 5. Schematic illustration of the phase-frequency response of the forced array of coupled waveguides due to self-interaction.The phase is in degrees.At the fixed frequency ~3.1, a change in driver phase Δ of 90° (horizontal arrow) may take the nonlinear mode (closed circle) into the single-valued state shown as the open circle; that is, a change in  of approximately 130° (vertical arrow).

Figure 5 .
Figure 5. Schematic illustration of the phase-frequency response of the forced array of coupled waveguides due to self-interaction.The phase is in degrees.At the fixed frequency σ~3.1, a change in driver phase ∆θ of 90 • (horizontal arrow) may take the nonlinear mode (closed circle) into the single-valued state shown as the open circle; that is, a change in η of approximately 130 • (vertical arrow).

Figure 8 .
Figure 8. Components of a three-phi-bit representation (see text for details) as functions of the drivers' phase difference.

Figure 8 .
Figure 8. Components of a three-phi-bit representation (see text for details) as functions of the drivers' phase difference.