Abstract
The classical three-body harmonic system in \({\mathbb {R}}^d\) (\(d>1\)) with finite rest lengths and zero total angular momentum \(L=0\) is considered. This model describes the dynamics of the \(L=0\) near-equilibrium configurations of three point masses \((m_1,m_2,m_3)\) with arbitrary pairwise potential \(V(r_{ij})\) that solely depends on the relative distances between bodies. It exhibits an interesting mixed regular and chaotic dynamics as a function of the energy and the system parameters. The corresponding harmonic quantum system plays a fundamental role in atomic and molecular physics. In this work we report on a novel electronic experimental realization of the model as a complementary tool to analyze the rich dynamics of the classical system. Our setup allows us to experimentally explore different regions of behavior due to the fact that the intrinsic parameters and initial states of the system are independently set by voltage inputs. Chaotic and periodic motions are characterized employing time series, phase planes, and the largest Lyapunov exponents as a function of the energy and system parameters. The results show an excellent qualitative as well as quantitative agreement between theory and experiment.
Similar content being viewed by others
Introduction
The three-body harmonic oscillator (TBHO) in \({\mathbb {R}}^d\) (\(d>1\)) with finite rest lengths is a 9-parametric system depending on three arbitrary masses, three rest lengths and three spring constants. For generic pairwise potentials, \(V=U(r_{12}) +U(r_{13}) +U(r_{23})\), that only depend on the three relative mutual distances \(r_{ij}=\mid {\mathbf{r}}_i-{\mathbf{r}}_j\mid\) between bodies, this model describes the dynamics near the equilibrium configurations and as such it possesses a wide range of relevant problems it can be applied to both in classical1,2 and quantum mechanics3. In the simplest planar case (\(d=2\)) with equal masses and equal spring constants the classical system already displays a very rich dynamics for different values of the energy where a power-law statistics that fits the Levi-walk model4 emerges5,6. Even on the invariant manifold of zero total angular momentum, regions of regular and chaotic dynamics coexist in the space of parameters due to the intrinsic nonlinearities induced by non-zero rest lengths. If these rest lengths are set equal to zero the system becomes integrable7. Not surprisingly, at zero rest lengths the corresponding quantum system, restricted to the subspace of zero total angular momentum, turns out to be an exactly-solvable problem with spectra linear in quantum numbers8. In variables \(\rho _{ij}=r_{ij}^2\), the relevant Hamiltonian operator possesses a hidden algebra \(s\ell (4)\) and it admits a quasi-exactly solvable extension. However, for non zero rest lengths a single exact solution to the Schrodinger equation has not been found so far. Hence, the TBHO also represents a realistic model to test theoretical and numerical tools to elucidate the connection between classical and quantum mechanics in chaotic systems9.
In this Letter, we report the physical realization of the classical TBHO using analog electrical components where the basic building blocks are non linearly coupled electrical inductance-capacitance oscillators. In other words, the dynamics of the electronic device is governed by the same equations as the TBHO. Therefore, it allows us to experimentally investigate the rich dynamics of the system in the whole space of parameters. Also, using computer simulations the average Lyapunov exponent, regular and chaotic regions can be identified. We observe an excellent agreement between the numerical simulation and the experiment.
Results
Let us consider in \({\mathbb {R}}^d\) (\(d>1\)) the TBHO at zero total angular momentum with the center of mass chosen as the origin of the reference frame. The corresponding Hamiltonian of the collision- and collinear configuration-free TBHO is of the form10,11
here \(m_{ij}\equiv \frac{m_i\,m_j}{m_i+m_j}\) are reduced masses (\(i\ne j=1,2,3\)), \(r_{ij}={\mid {\mathbf{r}}_i-{\mathbf{r}}_j\mid }\) denote the relative distances, \(P_{ij}\) the associated canonical momentum variables, and
where \(\omega > 0\) is the frequency, \(R_{ij}\ge 0\) are the rest lengths of the system and \(\nu _{ij}\) define spring constants, see Fig. 1.
The Hamiltonian (Eq. 1) is the symmetric reduction of the original 3d-dimensional problem. It governs all the trajectories of the original problem possessing zero total angular momentum. This is valid for any pairwise potential \(V=V(r_{12},r_{13},r_{23})\). Interestingly, it can also be interpreted as the Hamiltonian for a three-dimensional moving body in a curved space10.
Alternatively, one can promote the mutual distances (squared), namely \(\rho _{ij}=r_{ij}^2\), as new dynamical coordinates complemented by their corresponding canonical momentum variables. The benefits of this representation yield in the fact that (1) it encodes the discrete symmetries of the original 3d-dimensional kinetic energy, (2) the kinetic energy in the reduced Hamiltonian is simply a polynomial expression and (3) it nicely establishes an interesting link between the dynamics of the TBHO and the geometrical quantities of the so called triangle of interaction, the triangle formed by the vector positions of the three particles.
The corresponding co-metric in Eq. (1) is given by
and its determinant
(\(M=m_1+m_2+m_3\)) factorizes. The quantities
possess a clear geometrical interpretation. The term \(S_{\triangle }^2\) is the area squared of the triangle of interaction whilst \(\mathscr{I}\) is the moment of inertia of the system. In the case of 3 equal masses, one can further introduce the 3 lowest elementary symmetric polynomials \(\sigma _i \ (i=1,2,3)\) in \(\rho\)-coordinates as generalized coordinates. They are invariant under the action of \({\mathbb {Z}}_2^3 \oplus {\mathscr{S}}_3^2\), reflections (\(r_{ij} \rightarrow -r_{ij}\)) and permutations of the \(\rho\)-variables. In this representation, the associated reduced kinetic energy is a polynomial again. In particular, at \(d=3\) the discriminant of the 4th degree polynomial equation, solution of which corresponds to the physically relevant 3-body Newtonian potential in terms of \(\sigma -\)coordinates, do coincide with one of the factors of the determinant associated to the corresponding co-metric10.
Chaotic and periodic regimes
We perform simulations by numerically integrating the equations of motion obtained from Hamiltonian (Eq. (1)), in order to characterize periodic and chaotic regions for different values of the energy. We study the nature of the trajectories by computing the largest Lyapunov exponents. The average largest Lyapunov exponents as a function of the energy H and the mass \(m_1\) are shown in Fig. 2(i). We compute the average Lyapunov exponent for each pair (H,\(m_1\)) by integrating the motion equations for 2500 initial conditions in the state space defined by H, \(P_{13}=0\), \(P_{23}=0\) and \(r_{23}=1\). Because of Eq. (3), the dynamics associated with configurations for which the area of the triangle of interaction, during time evolution, vanishes are not considered to calculate the average largest Lyapunov exponents. Here, the masses and spring constants are set to \(m_2=m_3=1\), \(\nu _{12}=\nu _{31}=\nu _{23}=1\), \(\omega =1\) and \(R_{12}=R_{13}=R_{23}\equiv R= 1\). We selected three different energies, whose values are indicated on the phase map by blue points: (a) \(H=0.2\), (b) \(H=0.4\) and (c) \(H=0.6\). Detailed maps of the largest Lyapunov exponents in the plane (\(P_{12}\), \(r_{12}\)), corresponding to the selected energies are shown on the lower panel of Fig. 2. For low energies, regular dynamics can be found only (\(H=0.2\)). For energies close to \(H =0.4\), mixed zones where chaos and regularity coexist acquire relevance in the dynamics. When the energy increases up to \(H = 0.6\), chaotic dynamics predominate.
To further illustrate the rich dynamic of the system, in Fig. 3 the nontrivial dependence of the largest Lyapunov exponent as a function of the energy H and the equilibrium distance R is displayed, whereas in the Fig. 4 the ground state energy \(E_0\) of the quantum system calculated with the Lagrange-mesh method12 is presented. In Fig. 3, at fixed energy \(H=E\) the largest Lyapunov exponent is not a monotonic function of the rest length R. Interestingly, in the quantum case we observe from Fig.4 that a local shallow minimum of \(E_0(R)\) appears.
Importantly, the three-body chain of harmonic oscillators can be electronically implemented using the working principle of an analog computer. In this device, the physical quantities from the problem to be solved such as fluid, pressure, or velocity, are represented by voltage signals, so the temporal evolution of the variables of interest is directly measured on the analog circuit. This has enabled the implementation of experimental systems that offer certain advantages over the original physical systems13,14,15,16,17. Especially, analog computers are suited to simulate dynamical systems18,19,20,21,22,23, and even, stochastic dynamical systems24,25. A potential advantage of analog computers is the simulation of continuous physical quantities by performing real-time simultaneous operations. Particularly, when it is sought to integrate dynamical systems that include high-frequency components, digital computers demand high computational resources because they require small steps during numeral integration whereas analog computers carry out a continuous integration. Here, we exploit this fact to build an electronic version of the three-body chain of harmonic oscillators and to experimentally investigate the presence of regular and chaotic behaviors on the analog model. The projections of the trajectories in the planes (\(P_{12}\), \(r_{12}\)), for the numerical simulations and experimental results are shown in Fig. 2(ii)(a–c). Additionally, we plot their respective time series in Fig. 2(iii)(a–c). Here, the purple line denotes the temporal evolution of \(r_{12}\) and blue one of \(P_{12}\). For energy \(H \sim\) 0.2, the trajectories are periodic [see Fig 2(ii)(a)], characterized by a Lyapunov exponent close to zero and a completely regular oscillation as displayed in Fig 2(iii)(a). Note that as energy increases up to \(H \sim\) 0.4, a chaotic trajectory is observed in Fig 2(ii)(b) and confirmed for an irregular time series in Fig. 2(iii)(b). As Fig 2(ii)(c),(iii)(c) show the phase plane and time series exhibit chaotic dynamic for \(H \sim\) 0.7. Note that the experimental results are in quantitative and qualitative agreement with the numerical simulations.
Discussion
In summary, we have presented a novel electronic experimental realization of the classical three-body closed chain of harmonic oscillators. Our setup allowed us to explore different regions of behavior due to the fact that the intrinsic parameters and initial states of the system are independently set by voltage inputs. Complementary phase planes and Lyapunov exponents were computed to characterize periodic and chaotic regions for different values of the energy. Importantly, the experimental results are in quantitative and qualitative agreement with the numerical simulations.
Methods
We design and implement an experimental setup based on active electrical networks and passive linear electrical components to reproduce the dynamics described by the Hamiltonian (Eq. (1)). The active networks consist of interconnecting passive electronic components, and operational amplifiers (OPAMPs) to satisfy basic configurations whose transfer functions are equivalent to mathematical operations, namely, addition, subtraction, integration, and amplification. These electronic circuits are typically employed to implement linear operations; however, it is possible to carry out nonlinear functions, such as products and divisions between signals from logarithm and antilogarithm configurations connected in cascade or parallel. Both functions comprise linear circuits in which the output voltage is proportional to the natural logarithm and antilogarithm of the input, respectively. Interestingly, nonlinear operations are typically performed through integrated analog multipliers. The main reason for this is that log and anti-log amplifiers include diodes and transistors in their feedback, which leads to a high dependence on the temperature. In addition, the bandwidth is limited by the signal amplitude. It is worth mentioning that through connecting linear and nonlinear basic configurations in a structured way, it is possible to build complex operations where the involved voltage signals represent physical variables of the system of interest. We perform nonlinear functions between signals such as multiplication/division and squaring by employing analog multipliers AD633JN. In the analog multipliers, the manufacturer includes a factor that divides the output voltage over ten, in order to avoid saturation of the device.
Traditionally, the parameters of the physical system are mapped in passive components on the electronic circuit, specifically, in the relationship of two resistors or resistor and capacitor. Under this condition, the values that can take the parameters are limited to the commercial values of the components. To overcome this limitation, we introduce the system parameters (\(m_1\), \(m_2\), \(m_3\), \(\nu _{12}\), \(\nu _{23}\), \(\nu _{31}\), \(\omega\), R) through voltage signals generated by digital-to-analog converters, which are manufactured by Microchip with series MCP4921. These devices can take discrete values between 0V and 5V with a resolution of 1.22 mV (12 bits) by allowing to select parameters from a broad range of possible values. Note that we have transferred the issue from replacing passive electronic components to carrying out the product of voltage signals with analog multipliers.
Because of their importance, the initial conditions [\(P_{12}(0)\), \(P_{23}(0)\), \(P_{13}(0)\), \(\rho _{12}(0)\), \(\rho _{23}(0)\), \(\rho _{31}(0)\)] are also set via voltage values. In this sense, each of the voltage signals is individually controlled via software. Using a serial peripheral interface (SPI) protocol, digital-to-analog converters are communicated to a master system board, which incorporates a 32-bits STM32F401CC microcontroller from the ARM Cortex-M4 family. Figure 5a shows the scheme of the designed electronic circuit. Because canonical variables (\(P_{ij},r_{ij}\)) of the Hamiltonian (Eq. (1)) are described by motion equations with a similar structure, we show solely the electronic circuit for the variables \(\dot{P}_{12}\) and \(\dot{r}_{12}\). There, \(R_j\) , \(C_j\), \(U_j\), and \(M_j\) stand for resistors and capacitors, general-purpose operational amplifiers LF353, and analog multiplier AD633JN, respectively. Here the input and output signals are indicated with red and blue nodes, individually.
Figure 5b shows the printed circuit board of the classical three-body closed chain of harmonic oscillators. Electrical contact failures were avoided by mounting and soldering all electronic components on a printed circuit board (PCB, \(20 \times 20\) cm). In Fig. 5b, the initial conditions and control parameter modules are denoted with magenta and red squares, respectively. Note that the motion equations occupy a great part of the PCB part (green square). A stabilized DC power supply (KEITHLEY Triple channel, 2231A-30-3) was used to energize our setup, whereas the acquisition of the output signals (orange square), corresponding to the physical variables of the system, was performed with Rohde-Schwarz oscilloscope. The experimental setup mounted in the laboratory is shown in Fig. 6 (“Supplementary information”).
Data availability
The data that support the findings of this study are available from the corresponding authors upon reasonable request.
References
Musielak, Z. E. & Quarles, B. The three-body problem. Rep. Progress Phys. 77, 065901 (2014).
Moore, C. Braids in classical dynamics. Phys. Rev. Lett. 70, 3675–3679. https://doi.org/10.1103/PhysRevLett.70.3675 (1993).
Estienne, C. S., Busuttil, M., Moini, A. & Drake, G. W. F. Critical nuclear charge for two-electron atoms. Phys. Rev. Lett. 112, 173001. https://doi.org/10.1103/PhysRevLett.112.173001 (2014).
Zaburdaev, V., Denisov, S. & Klafter, J. Lévy walks. Rev. Mod. Phys. 87, 483–530. https://doi.org/10.1103/RevModPhys.87.483 (2015).
Saporta Katz, O. & Efrati, E. Self-driven fractional rotational diffusion of the harmonic three-mass system. Phys. Rev. Lett. 122, 024102. https://doi.org/10.1103/PhysRevLett.122.024102 (2019).
Saporta Katz, O. & Efrati, E. Regular regimes of the harmonic three-mass system. Phys. Rev. E 101, 032211. https://doi.org/10.1103/PhysRevE.101.032211 (2020).
de Castro, A. S. & Sugaya, M. Exact solution for a three-dimensional three-body problem with harmonic interactions. Eur. J. Phys. 14, 259–261 (1993).
Turbiner, A. V., Miller, W. & Escobar-Ruiz, M. A. Three-body closed chain of interactive (an)harmonic oscillators and the algebra \(s{\ell} (4)\). J. Phys. A 53, 055302. https://doi.org/10.1088/1751-8121/ab5f39 (2020).
Gutzwiller, M. C. Periodic orbits and classical quantization conditions. J. Math. Phys. 12, 343–358 (1971).
Escobar-Ruiz, A. M., Linares, R., Turbiner, A. V. & Miller, W. Jr. Classical n-body system in geometrical and volume variables: I. Three-body case. Int. J. Mod. Phys. A 36, 2150140. https://doi.org/10.1142/S0217751X21501402 (2021).
Murnaghan, F. D. A symmetric reduction of the planar three-body problem. Am. J. Math. 58, 829–32. https://doi.org/10.2307/2371252 (1936).
Baye, D. The lagrange-mesh method. Phys. Rep. 565, 1–107 (2015).
León-Montiel, R. D. J., Svozilík, J. & Torres, J. P. Generation of a tunable environment for electrical oscillator systems. Phys. Rev. E. 90, 012108 (2014).
Lamata, L., Parra-Rodriguez, A., Sanz, M. & Solano, E. Digital–analog quantum simulations with superconducting circuits. Adv. Phys. X 3, 1457981 (2018).
Quiroz-Juárez, M. A., Chávez-Carlos, J., Aragón, J. L., Hirsch, J. G. & León-Montiel, R. D. J. Experimental realization of the classical Dicke model. Phys. Rev. Res. 2, 033169 (2020).
Quiroz-Juárez, M. A. et al. Reconfigurable network for quantum transport simulations. Phys. Rev. Res. 3, 013010 (2021).
Jiménez-Ramírez, O., Cruz-Domínguez, E., Quiroz-Juárez, M., Aragon, J. & Vázquez-Medina, R. Experimental detection of hopf bifurcation in two-dimensional dynamical systems. Chaos Solitons Fractals X 6, 100058 (2021).
Vazquez-Medina, R., Jimenez-Ramirez, O., Quiroz-Juarez, M. & Aragon, J. Arbitrary waveform generator biologically inspired. Chaos Solitons Fractals 51, 36–51 (2013).
Noordergraaf, A., Verdouw, P. D. & Boom, H. B. The use of an analog computer in a circulation model. Progress Cardiovasc. Diseases 5, 419–439 (1963).
Moore, C. Dynamical recognizers: Real-time language recognition by analog computers. Theoret. Comput. Sci. 201, 99–136 (1998).
León-Montiel, R. D. J. et al. Observation of slowly decaying eigenmodes without exceptional points in Floquet dissipative synthetic circuits. Commun. Phys. 1, 1–11 (2018).
Quiroz-Juárez, M. et al. Generation of ECG signals from a reaction-diffusion model spatially discretized. Sci. Rep. 9, 1–10 (2019).
Quiroz-Juárez, M. A. et al. On-demand parity-time symmetry in a lone oscillator through complex, synthetic gauge fields. arXiv preprint. arXiv:2109.03846 (2021).
León-Montiel, R. D. J. et al. Noise-assisted energy transport in electrical oscillator networks with off-diagonal dynamical disorder. Sci. Rep. 5, 1–5 (2015).
Quiroz-Juárez, M. et al. Emergence of a negative resistance in noisy coupled linear oscillators. EPL (Europhys. Lett.) 116, 50004 (2017).
Acknowledgements
The Authors are grateful to L. Jiménez-Lara for their interest in the work and useful discussions. This work is partially supported by Programa Especial de Apoyo a la Investigación 2021, UAM-I.
Author information
Authors and Affiliations
Contributions
A.M.E.R., M.A.Q.J., J.L.R.C and N.A. provided the theoretical and numerical analysis. M.A.Q.J. and A.M.E.R. conceived the experimental setup. M.A.Q.J. designed and implemented the experimental setup. All authors contributed extensively to the planning, discussion and writing up of this work.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Escobar-Ruiz, A.M., Quiroz-Juarez, M.A., Del Rio-Correa, J.L. et al. Classical harmonic three-body system: an experimental electronic realization. Sci Rep 12, 13346 (2022). https://doi.org/10.1038/s41598-022-17541-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-17541-0
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.