Separability and Nonseparability of Elastic States in Arrays of One-Dimensional Elastic Waveguides Separability and Nonseparability of Elastic States in Arrays of One-Dimensional Elastic Waveguides

We show that the directional projection of longitudinal waves propagating in a parallel array of N elastically coupled waveguides can be described by a nonlinear Dirac-like equation in a 2 N dimensional exponential space. This space spans the tensor product Hilbert space of the two-dimensional subspaces of N uncoupled waveguides grounded elastically to a rigid substrate (called φ -bits). The superposition of directional states of a φ -bit is analogous to that of a quantum spin. We can construct tensor product states of the elastically coupled system that are nonseparable on the basis of tensor product states of N φ -bits. We propose a system of coupled waveguides in a ring configuration that supports these nonseparable states.


Introduction
Quantum bit-based computing platforms can capitalize on exponentially complex entangled states which allow a quantum computer to simultaneously process calculations well beyond what is achievable with serially interconnected transistor-based processors. Ironically, a pair of classical transistors can emulate some of the functions of a qubit. While current manufacturing can fabricate billions of transistors on a chip, it is inconceivable to connect them in the exponentially complex way that would be required to achieve nonseparable quantum superposition analogues. In contrast, quantum systems possess such complexity through the nature of the quantum world. Outside the quantum world, the notion of classical nonseparability [1][2][3] has been receiving a lot of attention from the theoretical and experimental point of views in the field of optics. Degrees of freedom of photon states that span different Hilbert spaces can be made to interact in a way that leads to local correlations. Correlation has been achieved between degrees of freedom that include spin angular momentum and orbital angular momentum (OAM) [4][5][6][7][8][9], OAM, polarization and radial degrees of freedom of a beam of light [10] as well as propagation direction [11,12]. Recently, we have extended this notion to correlation between directional and OAM degrees of freedom in elastic systems composed of arrays of elastic waveguides [13]. This classical nonseparability lies only in the tensor product Hilbert space of the subspaces associated with these degrees of freedom. This Hilbert space does not possess the exponential complexity of a multiqubit Hilbert space, for instance. It has been suggested theoretically and experimentally that classical systems coupled via nonlinear interactions may have computational capabilities approaching that of quantum computers [14][15][16].
We demonstrated in Ref. [17] that nonlinear elastic media can be used to produce phonons that can be correlated simultaneously in time and frequency. We have also shown an analogy between the propagation of elastic waves on elastically coupled one-dimensional (1D) wave guides and quantum phenomena [18][19][20][21]. More specifically, the projection on the direction of propagation of elastic waves in an elastic system composed of a 1D waveguide grounded to a rigid substrate (denoted φ-bit) is isomorphic to the spin of a quantum particle. The pseudospin states of elastic waves in these systems can be described via a Dirac-like equation and possess 2 Â 1 spinor amplitudes. Unlike the quantum systems, these amplitudes are, however, measurable through the measurement of transmission coefficients. The notion of measurement is an important one as it has been realized that separability is relative to the choice of the partitioning of a multipartite system. Indeed, it is known that given a multipartite physical system, whether quantum or classical, the way to subdivide it into subsystems is not unique [22,23]. For instance, the states of a quantum system may not appear entangled relative to some decomposition but may appear entangled relative to another partitioning. The criterion for that choice may be the ability to perform observations and measurements of some degrees of freedom of the subsystems [23].
The objective of this paper is to investigate the notion of separability and nonseparability of multipartite classical mechanical systems supporting elastic waves. These systems are composed of 1D elastic waveguides that are elastically coupled along their length to each other and/or to some rigid substrate. The 1D waveguides support spinor-like amplitudes in the twodimensional (2D) subspace of directional degrees of freedom. The amplitudes of N coupled waveguides span an N-dimensional subspace. Subsequently, the Hilbert space spanned by the elastic modes is a 2N-dimensional space, comprised of the tensor product of the directional and waveguides subspaces. This representation is isomorphic to the degrees of freedom of photon states in a beam of light. While beams of light cannot be decomposed into subsystems, an elastic system composed of coupled 1D waveguides can. Indeed, the elastic system considered here forms a multipartite system composed of N 1D waveguide subsystems. We show that, since each waveguide possesses two directional degrees of freedom, one can represent the elastic states of the N-waveguide system in the 2 N dimensional tensor product Hilbert space of N 2D spinor subspaces associated with individual waveguides. The elastic modes in this representation obey a 2 N dimensional nonlinear Dirac-like equation. These modes span the same space as that of uncoupled waveguides grounded to a rigid substrate, i.e., N φ-bits. However, the modes' solutions of the nonlinear Dirac equation cannot be expressed as tensor products of the states of N uncoupled grounded waveguides, i.e., φ-bit states.
In Section 2 of this chapter, we introduce the mathematical formalism that is needed to demonstrate the nonseparability of elastic states of coupled elastic waveguides in an exponentially complex space. Throughout this section, we use illustrations of the concepts in the case of systems composed of small numbers of waveguides. However, the approach is fully scalable and can be generalized to any large number of coupled waveguides. In Section 3, we draw conclusions concerning the applicability of this approach to solve complex problems.

Models and methods
We have previously considered systems constituted of N one-dimensional (1D) waveguides coupled elastically along their length [13]. In this section, we summarize the results of these previous investigations to develop a formalism to address our current considerations. The parallelly coupled waveguides can be arranged in any desired way. The propagation of elastic modes is limited to longitudinal modes along the waveguides in the long wavelength limit, i.e., the continuum limit. We consider the representations of the modes of the coupled waveguide systems in two spaces. The first space scales linearly with N. The second space scales as 2 N and leads to a description of the elastic system with exponential complexity. The linear representation enables us to operate easily on the states in the exponential space.

Representation of elastic states in a space scaling linearly with N
A compact form for the equations of motion of the N coupled waveguides is: Here, the propagation of elastic waves in the direction x along the waveguides is modeled by the dynamical differential operator, H ¼ ∂ 2 ∂t 2 À β 2 ∂ 2 ∂x 2 . The parameter β is proportional to the speed of sound in the medium constituting the waveguides and the parameter α 2 characterizes the strength of the elastic coupling between them (here, we consider that the strength is the same for all coupled waveguides). u NÂ1 is a vector with components, u i , i ¼ 1; N ð Þ , representing the displacement of the ith waveguide. The coupling matrix operator M NÂN describes the elastic coupling between waveguides which, in the case of N = 3 parallel waveguides in a closed ring arrangement with first neighbor coupling, takes the form: Eq. (1) takes the form of a generalized Klein-Gordon (KG) equation and its Dirac factorization introduces the notion of the square root of the operator H: In this factorization, the dynamics of the system are represented in terms of first derivatives with respect to time, t, and position along the waveguides, x. There are two possible Dirac equations: In Eq. (3), U NÂN and U 2NÂ2N are antidiagonal matrices with unit elements. σ x ¼ 0 1 1 0 and σ y ¼ 0 i Ài 0 are two of the Pauli matrices. Ψ 2NÂ1 is a 2 N dimensional vector which represents the modes of vibration of the N waveguides projected in the two possible directions of propagation (forward and backward) and ffiffiffiffiffiffiffiffiffiffiffiffiffi M NÂN p is the square root of the coupling matrix. The square root of a matrix is not unique but we will show later that we can pick any form without loss of generality.
We choose components of the Ψ 2NÂ1 vector in the form of plane waves ψ I ¼ a I e ikx e iωt with I ¼ 1, …, 2N and k and ω being the wave number and angular frequency, respectively, Eq. (3) becomes: In Eqs. (4) and (5), σ z ¼ 1 0 0 À1 is the third Pauli matrix, I NÂN is the identity matrix of order N and a 2Nx1 is a 2 N dimensional vector whose components are the amplitudes a I . In obtaining Eq. (4), we have multiplied all terms in Eq. (3) on the left by U 2NÂ2N .
Writing Eq. (4) as a linear combination of tensor products of N Â N and 2 Â 2 matrix operators: we seek solutions in the form of tensor products: .
While the degrees of freedom associated with E NÂ1 span an N dimensional Hilbert subspace, the degrees of freedom associated with s 2Â1 span a 2D space.
Replacing a 2NÂ1 from Eq. (7) in Eq. (6) yields: Choosing E NÂ1 to be an eigenvector, e n , of the matrix ffiffiffiffiffiffiffiffiffiffiffiffiffi M NÂN p with eigen value λ n Eq. (8) reduces to: For nontrivial eigenvectors e n , the problem in the space of the directions of propagation reduces to finding solutions of In obtaining Eq. (9), we have also used the fact that e n is an eigen vector of I NÂN with eigen value 1 and we note that Eq. (9) is the 1D Dirac equation for an elastic system which solutions, s 2Â1 , have the properties of Dirac spinors [18][19][20][21]. The components of the spinor represent the amplitude of the elastic waves in the positive and negative directions along the waveguides, respectively.
Eq. (10), now written in the matrix form, can now be solved for a given λ n ; ω n À βk AEαλ n AEαλ n ω n þ βk s 1 This eigen equation gives the dispersion relation ω 2 n ¼ βk À Á 2 þ αλ n ð Þ 2 (vide infra) and the following eigen vectors projected into the space of directions of propagation: To determine the eigen vectors of ffiffiffiffiffiffiffiffiffiffiffiffiffi M NÂN p , we note that they are identical to the eigen vectors of the coupling matrix M NÂN and the eigen values of M NÂN are also λ 2 n . These properties indicate that we do not have to determine the square root of the coupling matrix to find the solutions a 2NÂ1 . All that is required is to calculate the eigen vectors and the eigen values of the coupling matrix. Hence, the nonuniqueness of ffiffiffiffiffiffiffiffiffiffiffiffiffi M NÂN p does not introduce difficulties in determining the elastic modes of the coupled system in the Dirac representation.
In the case of the coupling matrix, M 3Â3 , presented in Eq. (2), the eigen values and real eigen vectors are obtained as λ 2 Eq. (3) being linear, its solutions can be written as linear combinations of elastic wave functions in the form: In Eq. (14), we have expressed the dependencies on the wave number k and the number of waveguides N. The eigen vectors e n N ð Þ depend on the connectivity of the N waveguides. The space spanned by these solutions scales linearly with the number of waveguides, i.e., as 2N.

Representation of elastic states in a space scaling as 2 N
We first illustrate the notion of exponential space in the case of three waveguides. Each guide is connected to a rigid substrate and therefore constitutes a φ-bit. The waveguides are not coupled to each other. The dynamics of the system can be described by a single equation which is constructed as follows: In Eq. (15), we are now defining a positional variable for each waveguide, namely, The quantity α is a measure of the strength of the elastic coupling to the rigid substrate. The . When seeking solutions in the form of tensor products of spinor solutions for the three waveguides (as indicated by the upper scripts) it is straightforward to show that one recovers from Eq.
The Hilbert space spanned by the solutions of Eq. (15) is the product space of the three 2D subspaces associated with each waveguide. The states of a system composed of N φ-bits span a space when dimension is 2 N .
The question that arises then concerns the possibility of writing an equation in the exponential Hilbert space for N waveguides coupled to each other. For instance, we wish to obtain the states of the system composed of three waveguides coupled in a ring arrangement from an equation of the form: The matrix αε 8Â8 represents the coupling between the waveguides in the 2 N¼3 space. We are still seeking solutions in the form of tensor products (Eq. (16)). After a lengthy algebraic manipulation, we find that we can reproduce Eq. (3) with the coupling matrix of Eq.
(2) if one chooses ε ij ¼ 0 excepting The Dirac equation of the three coupled waveguides in the exponential space is therefore nonlinear. Generalization to N coupled chains will result in the following nonlinear equation: The solutions of the nonlinear Dirac equation for the coupled waveguides span the same space as that of the system of φ-bits, i.e., uncoupled waveguides connected to rigid substrates. The next subsection addresses the question of separability of the coupled waveguide system into a system of uncoupled φ-bits.

Elastic states in the exponential space
For a system of waveguides that are not coupled, the elastic states, solutions of linear equations of the form of Eq. (15), are tensor products but also linear combinations of tensor products of spinor solution for individual waveguides (see Eq. (17)). It is therefore possible to construct nonseparable states in the exponential space for systems of uncoupled waveguides. For example, if we consider a system of two uncoupled waveguides, a possible state of the system in the 2 2 space can be constructed in the form of the following linear combination of tensor products: Choosing s 0 ¼ s 0 0 and writing Eq. (21) at the location x ¼ 0, one gets: The bracket takes the form: The vector is not separable into a tensor product of two 2 Â 1 vectors. Considering on the basis 0 j i ¼ 1 0 and 1 j i ¼ 0 1 , one can write the state given in Eq. (22) in the form of the nonseparable Bell state: Since the waveguides are not coupled, it is, however, not possible to manipulate the state of one of the waveguides by manipulating the state of the other one. Simultaneous manipulation of the state of waveguides in the exponential space requires coupling. We now address elastic states in the coupled waveguides system.
For a system of N coupled waveguides, we construct a solution of Eq. (3) that takes the form of a linear combination of solutions given in Eq. (14): The n and n 0 correspond to two different nonzero eigen values, λ n and λ n 0 , i.e., they correspond to two different dispersion relations ω n k ð Þ and ω n 0 k ð Þ. We also choose the wave number k 0 such that ω n 0 k 0 ð Þ ¼ ω n k ð Þ ¼ ω 0 . These modes are illustrated in Figure 1 in the case of an N = 9 waveguide system. χ n and χ n 0 are the coefficients of the linear combination.
With e n N ð Þ ¼ where the specific values of the components A I and A 0 I are determined by the connectivity and coupling of the waveguides, the state of Eq. (25) can be rewritten as: Here, we have chosen, for the sake of simplicity, the þ of the AE in the s 2Â1 terms. The first two terms in Eq. (26) form a 2 Â 1 spinor, φ 1 ð Þ ¼ φ , which corresponds to the first waveguide, the next two terms form a spinor φ 2 ð Þ for the second waveguide, etc. We can then construct a solution of the nonlinear Dirac Eq. (20) in the exponential space as the tensor product: Since Eq. (20)  . If one desires to express Φ 2 N Â1 as a nonseparable state, one has to define a new basis in which this wave function cannot be expressed as a tensor product. This is done in Section 2.5. However, prior to demonstrating this, we illustrate in the next subsection how one can manipulate states of the form Φ 2 N Â1 in the exponential space.

Operating on exponentially-complex tensor product elastic states
In this subsection, we expand tensor product states of the form given in Eq. (27) in linear combinations of tensor products of pure states in the exponential space. We illustrate this expansion in the case of three parallel waveguides elastically coupled to each other. Each waveguide is also coupled elastically to a rigid substrate. We treat the case where the strength of all the couplings is the same. In that case, the coupling matrix is: This matrix has three nonzero eigen values λ 2 0 ¼ 1, and λ 2 1 ¼ λ 2 2 ¼ 4 corresponding to two dispersion relations ω 2 n ¼ βk À Á 2 þ αλ n ð Þ 2 with cutoff frequencies. The second band is doubly degenerate. The eigen vectors are also given in Eq. (13). We now consider an elastic mode in the linear space that is a linear combination of these eigen modes (see Eq. (26)): In Eq. (28), the A I 's can be the components of the eigen vector e 0 and the A 0 I 's can be linear combinations of the components of the eigen vectors e 1 and e 2 . We can calculate the tensor product of the spinor components in the form of Eq. (27) Eq. (29) can be rewritten after some algebraic manipulations in the form of the linear combination: In Eq. (30), we have defined The tensor product of Eq. (30) then reduces to The spinors s 2Â1 and s 0 2Â1 can be expressed on the basis 0 j i ¼ 1 0 and 1 j i ¼ 0 1 : . Inserting Eqs. (33a) and (33b) into Eq. (32), we can express the tensor product Φ 2 3 Â1 on the basis 0 In defining the basis vectors for the exponential space, we have omitted the symbols ⊗ . It is also implicit that the left, middle, and right elements in the tensor product a j i b j i c j i correspond to the first, second, and third waveguides, respectively. We find In a true quantum system composed of three spins for instance, states can be created in the form of linear combinations like m 1 0 interdependent. While somewhat restrictive compared to true quantum systems, the coefficients T I depend on an extraordinary number of degrees of freedom which allows exploration of a large volume of the exponential tensor product space. In the case of the three waveguides, these degrees of freedom include (a) the components (or linear combinations) of the eigen vectors of the coupling matrix through the choice of the eigen modes or the application of a rotational operation that creates cyclic permutations of the eigen vector components, (b) the linear coefficients χ n and χ n 0 used to form the multiband linear superposition of states in the linear space, (c) the frequency and therefore wave number which affect the spinor states and the phase factors e ikx and e ik 0 x , and (d) a phase added to the terms e ikx and e ik 0 x .
In the case of N > 3, Eq. (21) can be extended to linear combinations of more than two modes with the same frequency, leading to additional freedom in the control of the T I . Furthermore, the elastic coefficients β of the waveguides and the coupling elastic coefficient α could also be modified by using constitutive materials with tunable elastic properties via, for instance, the piezoelectric, magneto-elastic or photoelastic effects [24][25][26]. Also note that in all the examples we considered, the coupling of the waveguides had the same strength. Tunability of the coupling elastic medium would lead to the ability to modify the connectivity of the waveguides and therefore the coupling matrix. Exploration of the elastic modes given in Eq. (34) can be realized by varying any number of these variables. We illustrate in Figure 2 an example of operation in a very simple case. Figure 2b shows that by varying a single parameter one may achieve a wide variety of states. For instance, one can obtain states with T 1 > 0 and T 2 > 0 or T 1 > 0 and T 2 ¼ 0 or T 1 > 0 and T 2 < 0 or T 1 ¼ 0 and T 2 ¼ 0. Another interesting example occurs at χ n 0˜0:6, there only T 5 and T 8 are different from zero. Then, with q 11 ¼ T1 0:25 ð Þ T1 χ n 0 ð Þ and q 44 ¼ T4 0:25 ð Þ T4 χ n 0 ð Þ . Another interesting state occurs at χ n 0 ¼ 0:5. Here, we have T 1 ¼ T 5 , T 2 ¼ T 3 , T 4 ¼ T 8 and T 6 ¼ T 7 . We also have T 3 ¼ ÀT 5 and T 4 ¼ ÀT 6 This state can be written as the tensor product T 1 0:5 ð Þ 0 j i À T 6 0:5 ð Þ 1 j iÞ 0 j i À 1 j iÞ 0 j i À 1 j iÞ ð ð ð .
This simple example indicates the large variability in T I 's (i.e., of states) that we can achieve with a single variable. The large number of available variables will lead to even more flexibility in defining states and operators in the exponential tensor product space.

Nonseparability of states in exponentially complex space
States given in Eq. (27) are tensor products on the basis 0 They are therefore always separable in that basis. Consequently these states cannot be written as nonseparable Bell states. However, we might be able to identify a basis in which Eq. (27) is not separable.
Since the Dirac equation (Eq. (15)) for the uncoupled waveguides is linear, its solutions can be a tensor product of a linear combination of N different individual spinors. For instance, we can write: r I and μ I , I ¼ 1, …, N are linear coefficients.
The state of the coupled waveguide system will be separable in the exponential space into the state of φ-bits if A necessary condition for satisfying Eq. (38) is that Furthermore, the first two terms in the column vector Φ 2 N Â1 of the coupled system are φ Similarly, the ratio of the first two terms in the vector Ψ 2 N Â1 of the φ-bit system is given by Anecessary condition for Eq. (38) to be satisfied is that This leads to which can be rewritten as: This condition takes the more compact form: P, Q, R, S are real. Eq. (40) is true for all values of position x. At x ¼ 0, we obtain the relation P þ Q þ R þ S ¼ 0. Inserting that relation into Eq. (27) and eliminating Q yields: For this condition to be satisfied, one needs the real part of the right-hand side of the equation to be equal to zero. This can be achieved for all x's by setting k 0 ¼ k. In this case, equating the imaginary parts leads to R ¼ ÀP. However, when k 0 6 ¼ k, Eq. (39) and therefore Eq. (38) are not satisfied. When k 0 6 ¼ k which corresponds to considering a linear combination of multiband states, Φ 2 N Â1 is not separable into the tensor product of individually uncoupled φ-bit waveguides. Therefore, we conclude that there are a large number of solutions of the nonlinear Dirac equation (Eq. (20)) representing states of arrangements of elastically coupled 1-D waveguides that are not separable in the 2 N dimensional tensor product Hilbert space of individual φ-bits.
We illustrate the notion of nonseparability of exponentially complex states of a coupled system composed of N ¼ 2 waveguides on a basis in the exponential Hilbert space of two individual φ-bits. The waveguides are coupled to each other but also to a rigid substrate such that the coupling matrix, M NÂN , takes the form: The eigen values and real eigen vectors of this coupling matrix are λ 2 0 ¼ 1, and λ 2 1 ¼ 3 and Following the procedure of Section 2.4, we construct a tensor product state in the 2 2 exponential space: Eq. (42) is equivalent to Eq. (32) but for two coupled waveguides.
We can find the coefficients c ij by exploiting the orthogonality of the τ i 0 s. For instance, we can multiply Eq. (49a) to the left by the Hermitian conjugate τ † 1 , leading to The matrix C 4Â4 can be diagonalized. Let d 1 , d 2 , d 3 , and d 4 be the four eigen values of C 4Â4 with their associated eigen vectors We can construct the following matrix out of the four eigen vectors: On the new basisτ 1 ;τ 2 ;τ 3 ;τ 4 f gconstructed by using the relationτ 4Â1 ¼ V À1 4Â4 τ 4Â1 V 4Â4 , the matrix that couples the η basis and the τ basis takes the form:C 4Â4 ¼ On theτ basis, Eq. (43) can be rewritten as Then on the basisτ, we can investigate the separability or nonseparability of Φ 2 2 Â1 by calculating the determinant of the linear coefficients in Eq. (53), that is det Only in the unlikely event of degenerate eigen values, d 1 , d 2 , d 3 , and d 4 , would this determinant be equal to zero. A nonzero determinant given in Eq. (54) indicates that the state Φ 2 2 Â1 is nonseparable on the basisτ 1 ;τ 2 ;τ 3 ;τ 4 f g .
The existence of nonseparable solutions to the nonlinear Dirac equation raises the possibility of exploiting these solutions for storing and manipulating data within the 2 N dimensional tensor product Hilbert space. The exploration of algorithms for exploiting these solutions is beyond the scope of this chapter; however, we note that these solutions may well be observed in physical systems including elastic waveguides which are embedded in a coupling matrix. The manipulation of the system could be achieved either by externally altering the parameters of the system, i.e., the elastic properties of the waveguides, or by changing the frequency and wavenumber of input waves. These possibilities are illustrated for a five-waveguide system driven by transducers in Section 2.6. Figure 3 illustrates a possible realization of a five waveguide system. The parallel elastic waveguides are embedded in an elastic medium which couples them elastically. The waveguides are arranged in a ring pattern.

Physical realization and actuation
Modes of the form given in Eq. (21) can be excited with N transducers attached to the input ends of the N waveguides and connected to N phase-locked signal generators to excite the appropriate eigen vectors e n and e n 0 . These modes can be excited by applying a superposition of signals on the transducers with the appropriate phase, amplitude and frequency relations. The frequencies ω n k ð Þ and ω n 0 k 0 ð Þ are used to control the spinor parts of the wave functions s 2Â1 k ð Þ and s 2Â1 k 0 ð Þ . The spinor components which represent a quasistanding wave can be quantified by measuring the transmission coefficient (normalized transmitted amplitude) along any one of the waveguides. It is then possible to operate on the eigen vectors e n and e n 0 without affecting the spinor states or vice versa. For instance, one could apply a rotation that permutes cyclically the components of e n by changing the phase of the signal generators. Such an operation could be quantified by measuring the phase of the transmission amplitude at the output end of the waveguides.

Conclusions
We have shown that the directional projection of elastic waves supported by a parallel array of N elastically coupled waveguides can be described by a nonlinear Dirac-like equation in a 2 N dimensional exponential space. This space spans the tensor product Hilbert space of the twodimensional subspaces of N uncoupled waveguides grounded elastically to a rigid substrate (which we called φ-bits). We demonstrate that we can construct tensor product states of the elastically coupled system that are nonseparable on the basis of tensor product states of N uncoupled φ-bits. A φ-bit exhibits superpositions of directional states that are analogous to those of a quantum spin, hence it acts as a pseudospin. Since parallel arrays of coupled waveguides span the same exponentially complex space as that of uncoupled pseudospins, the type of elastic systems described here may serve as a simulator of interacting spin networks. The possibility of tuning the elastic coefficients and the elastic coupling constants of the waveguides would allow us to explore the properties of spin networks with variable connectivity and coupling strength. The mapping between the 2N dimensional and the 2 N dimensional representations of the elastic system leads to the capacity for exploring an exponentially scaling space by handling a linearly growing number of waveguides (i.e., preparation, manipulation, and measurement of these states). The scalability of the elastic system, the coherence of elastic waves at room temperature, and the ability to measure classical superpositions of states may offer an attractive way for addressing exponentially complex problem through the analogy with quantum systems.