On quantum gravity tests with composite particles.

Models of quantum gravity imply a fundamental revision of our description of position and momentum that manifests in modifications of the canonical commutation relations. Experimental tests of such modifications remain an outstanding challenge. These corrections scale with the mass of test particles, which motivates experiments using macroscopic composite particles. Here we consider a challenge to such tests, namely that quantum gravity corrections of canonical commutation relations are expected to be suppressed with increasing number of constituent particles. Since the precise scaling of this suppression is unknown, it needs to be bounded experimentally and explicitly incorporated into rigorous analyses of quantum gravity tests. We analyse this scaling based on data from past experiments involving macroscopic pendula, and provide tight bounds that exceed those of current experiments based on quantum mechanical oscillators. Furthermore, we discuss possible experiments that promise even stronger bounds thus bringing rigorous and well-controlled tests of quantum gravity closer to reality.

Introduction

A minimum length scale of the order of Planck length is a feature of many models of quantum gravity (QG) that seek to unify quantum mechanics and gravitation[1]. But, the existence of such a length scale has so far eluded experimental verification. Direct detection of the Planck length, 10−35 m, is infeasible with current and foreseeable technology because the effects of quantum gravity are expected to become directly relevant only at energies of the order of Planck energy, 1019 GeV, which is 15 orders of magnitude larger than the energy scales achievable in the Large Hadron Collider today. So in order to experimentally probe quantum gravity, we must rely on indirect tests of signatures of Planck length.

One class of indirect tests of a minimal length scale is based on observations of photon arrival times from gamma-ray bursts in distant galaxies[2]. Such experiments, however, are hard to control as they rely on a wide range of hard to verify model assumptions. These include, but are not limited to, the properties of the models of the origins of gamma-ray bursts and the perturbations due to billions of light years of interstellar medium traversed by the gamma-rays. Furthermore, the precision of such experiments is difficult to extend due to limitations on the distance to observable gamma-ray bursts and the maximal energy of the gamma-rays. This motivates the quest for highly controlled table-top experiments to test for Planck scale physics[3-8].

The underlying concept on which several such table-top experiments rely upon is the deformations of the canonical commutation relations of position and momentum as a consequence of a variety of formulations of quantum gravity[1,9-12]. These deformations are a phenomenological approach to modelling the existence of a minimum length scale in quantum gravity and different models of quantum gravity involve different forms of modification of the commutator[13-15]. Here we focus on a paradigmatic model[13]where Mp = 2.176435 × 10−8 kg and c = 299,792,458 m s−1 are the Planck mass and speed of light, respectively and the variables x and p are the x – components of the position and momentum, respectively. β0 is a dimensionless parameter, which is expected to be of the order of unity if the minimal length scale is of the order of Planck length, but it is not fixed by theory. Measuring or placing bounds on this parameter is thus an open experimental challenge.

One approach to bounding the value of β0 is to use single particle systems, for example, using measurements of Landau levels, of the Lamb shift, or of electron tunnelling through a potential barrier[16]. However, the best bound obtained with these methods is β0 < 1020, which is far from the expected β0 ~ 1. To improve the bounds significantly, recent experimental proposals suggest using massive composite systems rather than elementary particles. These experiments aim to exploit the fact that the quantum gravity signal is enhanced with larger momenta, which result from larger system mass. Experiments and proposals in this direction include those based on the change in resonant frequency of a harmonic oscillator[5,8,17], the change in broadening times of large molecular wave-packets[18], and optomechanical schemes[3,6,7].

However, the implications of using multi-particle systems to probe quantum gravity are not clear. This is because the deformations of the canonical commutation relations like Eq. (1) have been derived for point particles and not for centre of mass (COM) modes of multi-particle objects[19]. The deformations for the COM modes are expected to decrease with the number of constituent particles in the test object[19,20], but the exact expression for this suppression is not known and therefore needs to be bounded by experiment. Studies of the soccer-ball problem[21-23], which arises in a different framework of quantum gravity, also point toward a suppression of the Planck-scale corrections with number of particles. Even if the scaling with particle number of this suppression was known, the question of what constitutes a fundamental particle remains open.

To address these challenges, we define a new parameter that accounts for the suppression of quantum gravity corrections with the number of constituent particles, irrespective of how these particles are defined. Such a parameter may be obtained from theory once a full quantum description of gravity is available, and meanwhile it can be estimated or bounded from experimental observations. We study past experiments and propose new experiments to obtain bounds on this and established parameters of quantum gravity deformations.

Results

Overview

To account for the unknown scaling law that governs the suppression of corrections to the canonical commutation relations with the number of particles, we define the parameter α such that the deformation is suppressed by N, i.e.,where N is the number of constituent particles in the test object. The reason for considering such a polynomial scaling with number of particles N (with exponent α) is the expected polynomial dependence in quasi-rigid macroscopic bodies, i.e., those whose centre-of-mass momenta are the sum of constituents’ individual momenta[19]. Once phenomenological evidence of deformed commutator models begins to accumulate, more refined models can be considered accounting for differences in the exact nature of commutator deformation experienced by bodies with different nature of interaction. Furthermore, we note that other models of quantum gravity could be considered, for instance, those involving non-commuting spacetime coordinates. However, our focus is on deformed position and momentum commutators because these underlie proposed tabletop tests of quantum gravity.

In this work, we argue that the inclusion and hence estimation of α is essential for the rigorous interpretation of any experiment that uses composite test objects to measure β0. We propose to assess any such experiment by the exclusion area in a two-dimensional parameter space spanned by α and β and carry out such an analysis for three experiments. While the precise value of α is unknown, it is commonly accepted that it needs to be positive[19,20]. We shall see that the best bounds that can be calculated from recent experiments based on micro- and nano-scale quantum harmonic oscillators[5,8] are in fact negative for β0 = 1 (and in fact any β0 < 106).

Here we show that measured data of a macroscopic pendulum reveals the first positive bound on α for any value of β0 > 10−2 based on a careful analytical examination of the effect of deformations of the canonical commutation relations (Eq. (2)) on the time period of a pendulum. Specifically, we obtain α > 0.07 for a value of β0 = 1, which is expected in various models of quantum gravity. The reason for this significant enhancement of the bound over those obtained from micro- and nano-scale quantum harmonic oscillators can be traced back to the low achievable momenta in those experiments which in turn lead to very weak deformations of the canonical commutation relations, for example, in the second term in Eq. (2). We argue that our bound on α can be improved further by moving a pendulum to a vacuum setup with optimised low damping suspension or by moving to diamagnetically levitated systems which exhibit extremely low damping rates[24,25] on earth and promise even better values when located in a space probe.

We begin with a study of the effects of deformed canonical commutation relations proposed in theories of quantum gravity on the time period of a macroscopic pendulum and use these analytical results to place a bound on the parameter α. Since a pendulum is not strictly a harmonic oscillator, we extend the calculations of corrections to the time period of a harmonic oscillator to that of a pendulum. The time period of a pendulum as a function of its amplitude can be determined experimentally with considerable precision. Using data from one such experiment[26] we provide the first positive bound on α and moreover, suggest refined experiments for substantially improved bounds.

Correction to time period of pendulum

Here we calculate the corrections to the time period of a pendulum due to quantum gravity deformations of the canonical commutation relations so that it can then be compared against experimental data. We consider a pendulum of mass m and length L. Its Hamiltonian iswhere p is the generalised momentum conjugate to x and θ is the angle between the radius vector and the vertical, i.e., . The time period can be obtained from the Hamiltonian in two ways. Here, we take the approach that defines a new momentum operator which satisfies the standard canonical commutation relations and hence the standard Heisenberg equations of motion. In this approach the Hamiltonian is modified. In an alternative, equivalent, approach the Hamiltonian can be left unchanged while the equations of motion are modified due to the deformed commutator[27]. Here we follow the first method due to ease of calculation.

For illustrative purposes, we perform the calculations using classical mechanics but later show that the results hold even by performing quantum calculations with deformed commutators. Here, we deform the standard Poisson brackets in analogy to the deformation of the canonical commutation relation due to quantum gravity[28-30], i.e.,where

To ensure that the equations of motion are unchanged, a new momentum operator is defined such that we recover the standard Poisson bracket, i.e.,

Without the small momentum assumption of refs. [5,8], we find that is related to p as

Writing the Hamiltonian in terms of the new momentum operator , we note that it is modified to

This modification of the Hamiltonian compared with Eq. (3) ensures that the equations of motion remain unchanged with respect to x and , i.e.,

Separating the variables and integrating over half a time period as detailed in the Methods, we obtain the time period for small amplitudes A and small β as approximately

Here, we have calculated the correction to the time period to first order in the quantum gravity parameter β. We argue that if β were not, in fact, small, then a much larger deviation in the time period would be observed. In such a case, more detailed calculations would be required to obtain accurate bounds. We show that the time period of a pendulum obtained from experimental data in the next section corroborates this assumption.

In the “Methods” section, we show that the Poisson equation approach and a fully quantum mechanical calculation yield the same result for the time period of a harmonic oscillator. This shows that the calculated expression for the time period (Eq. (10)) is not merely a consequence of using classical dynamics but extends to the quantum regime and thus provides further evidence for the correctness of our approach. The detailed calculations are presented in the “Methods” section, but we summarise the arguments here. To perform these calculations, we start with the eigenvalues and eigenfunctions of the quantum harmonic oscillator derived in ref. [13]. Since ladder operators are defined differently due to commutator deformation, we use a generalised Heisenberg algebra[31] to find the action of the ladder operators on the eigenstates. Using this algebra, we derive the expressions for the position and momentum operators in terms of the ladder operators. With this, the operators are well defined. In order to choose the most classical pure state in our calculations, we choose a definition of coherent states, the Gazeau-Klauder states[32], such that the states remain coherent states during the evolution under this Hamiltonian. Calculating the expectation value of the position operator with respect to these Gazeau-Klauder coherent states, we recover the classical calculations.

As a side remark we note that these matching results also connects two different approaches to studying deformed commutators, namely modifying the Poisson bracket[28-30] and modifying the commutator[13,33,34]. These two approaches have so far been thought to be separate[35].

Bounds on QG parameters from experimental data

Having calculated the correction to the time period of a pendulum, we use experimental data from precise measurements of the time period of a pendulum to place bounds on the quantum gravity parameters α and β0 using Eq. (10).

To this end, we consider an experiment[26] which measures the time-period of a pendulum as a function of its amplitude. The experimental data is represented in Fig. (3) of ref. [26], which plots the measured time-period as a function of the square of the amplitude. The figure contains data from two experiments: one using a conventional suspension of the pendulum and another using a cycloidal suspension. In this manuscript, we consider only the conventional suspension because such an experiment can be modelled by the calculations shown in the previous section.

We extracted the data of the experiment with the conventional suspension from Fig. (3)[26] to calculate the slope and intercept. From Eq. (10), we see that the intercept T0 of this line is and the slope is Since the expression for the slope of this plot includes corrections from quantum gravity, we calculate the slope from the extracted data from Fig. (3)[26] to compare with theory. These data are reported and the method of extraction of the data is detailed in Section.

We use the data extracted to perform a linear fit. The reported error of measurements is 2% in amplitude measurement and 3 × 10−5 s in time-period measurement[26]. However, additional error arising from the datapoint extraction have to be accounted for. For this, we use the size of the markers as a conservative estimate. This leads to an error of 5 × 10−3 m2 in the square of the amplitude and 10−4 s in the time period measurements.

We use the extracted data and account for both the sources of errors to perform a linear fit using the orthogonal distance regression method[36]. This fit yields a reduced χ2 value of 0.07, which indicates that the extracted data points do indeed agree with the linear fit. We note that the current work is only a blueprint that uses experimental data from 1964 and such an analysis can be repeated with improved experiments to determine tighter and more precise bounds.

From this fit we obtain the value of the intercept

Using this value and following the analysis of ref. [26], we can precisely infer the effective length of the pendulum to be

In this calculation, the local value of acceleration due to gravity g = 9.80393 m s−2 has been used[26]. The pendulum is an iron cylinder of radius 2.54 cm and height 5.08 cm and a mass of ~1.2 kg.

Using the obtained value of the length L and mass m, the slope is numerically evaluated to be

The uncertainty in the obtained slope due to uncertainty in the derived length of the pendulum is 9 × 10−5, which is small enough to be ignored in the following calculations.

From the linear fit of the data, the slope of the line is obtained to be 0.0232 ± 0.0012 (95% confidence interval). We see that for the fit and Eq. (13) to be consistent, we obtain . For β0 positive,

For the determination of N we assume that the nucleons form the elementary particles which leads to N = 7.32 × 1026 and therefore we obtain α > 0.07 for β0 = 1. Note that the bound on α is quite insensitive to the precise number of nucleons.

We note that such tight bounds were possible in our work as compared with those of previous works because the nonlinearities in a pendulum can be computed precisely. In other systems, the intrinsic nonlinearities are unknown and therefore contribute to larger bounds on quantum gravity parameters.

These bounds can therefore be made tighter by including the effect of dissipation on the change in frequency of a pendulum. We propose that future experiments measure the amount of dissipation and include its effect on the slope and intercept of the T versus A2 plot in order for the bounds to be more precise. We note that dissipation is not a feature of classical systems alone and quantum systems such as those proposed for other tests also couple to the environment, and this dissipation should be considered as well.

Recent experiment using oscillators in the quantum regime have been used to provide bounds on β0 under the assumption that α = 0. We argue however that any test of consequences of deformed canonical commutation relations due to quantum gravity need to account for both α and β0. The best bound on α from the experiments of ref. [5] is α > −0.33 for β0 = 1. Similarly, from ref. [8] we obtain α > −0.25. Note that these bounds are significantly worse than those obtained in the present work using the data from[26]. In Fig. 1 we present the parameter ranges that have been excluded in the α, β0-plane in the three experiments discussed here.

Fig. 1

Excluded regions of parameter space from different experiments.

This figure depicts the values of the two quantum gravity parameters β0, α that are excluded based on experimental observations. Solid lines represent bounds obtained from experimental data and dashed lines represent expected bounds from proposed experiments. The shaded areas represent the region excluded by these experiments. The present work based on[26] provides the largest excluded region of parameters which, in particular, excludes the key point β0 = 1, α = 0, thereby showing that suppression of quantum gravity deformations should be accounted for if β0 ~ 1 as expected from quantum gravity models. The proposal to use massive levitated diamagnetic objects described in the text promises significant improvement in bounds if the optimistic parameters required for such an experiment can be obtained.

Diamagnetic levitation for enhanced tests of QG

In order to explore how far we can possibly bound the α-parameter, here we propose an experiment that relies on the precise measurement of the oscillation frequency of a diamagnetic levitated particle to obtain enhanced bounds on the quantum gravity parameters. In an experiment on a space probe, such as LISA pathfinder, one could imagine to levitate a particle in a uniform magnetic field gradient. In this case the frequency of oscillation would be[37]where ρ is the mass density of the object, μ0 is the vacuum permeability, χv is the magnetic volume susceptibility of the material and we assume a constant magnetic field gradient. The change in frequency resulting from deformed commutators can be obtained from Eq. (10) and is given by

The oscillation frequency of levitated objects can be measured very precisely due to very low damping rates. Here we assume, optimistically, that the damping rate is the only source of error in frequency measurement. At low pressures of 266 × 10−10 Pa, the damping rate is expected to be γ = 1.2 × 10−7 Hz[38,39]. If, in the experiment, no deviation from the expected frequency is observed, then where Nm is the number of measurements taken.

We calculate the bounds that one would obtain if such an experiment can be performed. We consider optimistic parameters of a gold sphere of diameter 10 cm that is levitated in a uniform magnetic field gradient of 103 T/m that is initially displaced with an amplitude of 10 cm. The density and magnetic volume susceptibility of gold are ρ = 19,300 kg/m3 and χv = 3.287 × 10−5[40], respectively. This leads to the frequency of oscillations ω = 36.71 Hz, from which we estimate that for β0 = 1, we obtain α > 0.35 for a single measurement Nm = 1. However, we realise that such estimates are based on optimistic assessments of all the involved parameters. More realistically, the error in the frequency might not limited by damping alone and the different parameters feasible in a single experiment might be somewhat suboptimal once all components are integrated together. Thus, we can conservatively assume that the actual precision in frequency is three orders of magnitude worse than the optimistic value suggested above, i.e., Δω ≲ 10−4 Hz, we obtain the bounds α > 0.24 for β0 = 1.

The region of parameter space that can be excluded from such an experiment is shown in Fig. 1. This bound may be improved by performing the experiment in space, where the pressure is about 2000 times lower which leads to a further reduction in γ.

Using optomechanical tests to bound α

Experiments to place bounds on quantum gravity parameters are not restricted to the framework of measuring the change in frequency of oscillators. Here we show that other schemes can also be used to placed stringent bounds on α if appropriate initial states are considered.

Specifically, we consider the optomechanical scheme of a light pulse striking a mechanical resonator repeatedly and acquiring a phase that depends on quantum gravity corrections to the canonical commutation relations[3,6,7]. In this scheme, the mechanical resonator is initially in a thermal state very close to the ground state. After interacting with the resonator, the field of the light pulse is[3]where ξ is the amplitude of the coherent state of light, N is the mean photon number in the light pulse and λ is the optomechanical interaction strength. As before, m and ω are the mass and frequency of the oscillator. We see that the phase acquired by light has a term dependent on β, which originates from the additional contribution from quantum gravity. Using the experimental parameters of ref. [3] and the error analysis of ref. [7], and assuming that experiment returns a null result, we obtain the bound β0N− < 106 which leads to α > −0.3 for β0 = 1. The excluded parameter range is presented in Fig. 1.

This bound can be improved if, instead of the ground state, the resonator is initially in a coherent state with a large enough momentum. In this case, the output field is given byfor the initial state of the oscillator in a coherent state with mean initial momentum .

To understand better the mass dependence of the quantum gravity terms, we consider the mass dependence of λ using an established model[3]where is the finesse of the cavity and λL is the wavelength of light used. With this scaling of λ with the mass m, we see that the first negative term scales as while the second one scales as m. Hence, we see that the extra term arising from a non-zero initial momentum can, in principle, be made large by choosing a more massive oscillator. Making the optimistic assumption that all other parameters remain the same but we increase the mass to 10−3 kg a positive bound on α would be achieved. While such a large mass is likely to reduce the optomechanical coupling it nevertheless suggests that the measurement of the phase of the output light for macroscopic systems may provide another method to obtain a good bound on the deformation parameters.

Discussion

Quantum gravity suggests corrections to the canonical commutation relations that are proportional to a parameter β0. This parameter is expected to be of order of unity if physics exhibits a minimum length of order of the Planck length but is also expected to scale as N− where N is the number of constituent particles of the test mass—a consequence of the soccer ball problem of quantum gravity. We strongly argue that any test of such physics needs to account for both parameters αandβ0 in its analysis. We perform an analysis of several quantum regime experiments in those terms to show that they cannot provide positive bounds on α while we find that a macroscopic pendulum can provide the first positive bound on α assuming β0 = 1. This shows that the suppressions with the number of particles cannot be ignored in tests of quantum gravity and that entering the deep quantum regime is not essential for the observation of quantum gravity corrections to physical dynamics.

Our results pertaining to the time period are derived using the method of deformed Poisson brackets[28-30]. We also put this method on a more rigorous footing by connecting it to the calculations based on deformed commutators[13,33,34] and show that the two results match. Hence our work connects these two approaches that have thus far been considered independent[35]. Finally, we show that the suppression of quantum gravity deformations is not just restricted to this one framework of oscillator frequency measurement. For instance, we consider the optomechanical system of refs. [3,6,7] and verify that broadly analogous considerations hold. We discuss possible advanced experimental designs and the parameter requirements to allow for entering tests in the α > 1 regime that is suggested by various models of quantum gravity.

Methods

Time period of pendulum: detailed calculations

In this section, we obtain the expression for the time period of a pendulum using classical mechanics. This is done by solving the equation of motion for the modified Hamiltonian obtained in the calculation of the correction to time period of pendulum in the Results section.

From the expression for the energy in Eq. (8), the total energy in the system is , where ϕ is the angular amplitude. The redefined momentum (7) can be expressed in terms of the angular displacement θ and amplitude ϕ aswhere

The equation of motion Eq. (9) can be rewritten in terms of θ asand simplified to

Separating variables and integrating over half a cycle, we obtain the time period for half an oscillation

Since β ≪ 1, as can be numerically verified from the above equation, the expression can be simplified to

Furthermore, for small amplitudes, we can approximate the time-period toand in terms of the amplitude A, where A = ϕL, it can be expressed aswhich is the expression that has been used to compare with experiments.

Rigorous calculations using deformed commutators

In this section, we perform detailed calculations in a quantum mechanical framework to find the time-dependent position of a harmonic oscillator modified by quantum gravity. To do so, we start by using the energy eigenvalues and eigenstates derived quantum-mechanically in refs. [13,41] for deformed commutators. We then choose appropriate definitions of the ladder operators and derive the position operator in terms of these ladder operators. The definition of a coherent state is chosen such that the state remains invariant under evolution of the Hamiltonian.

In the following calculations, we work in the momentum basis where the operators and are defined by their action on the momentum wave-functions as

Solving the Schrödinger equation with Hamiltonianthe energy eigenvalues are found to be[13]for 1/r = (2βmℏω)2. The eigenfunctions in the momentum basis are[41]where are Gegenbauer polynomials and

The phase (−i) has been introduced here so that in the limit β → 0, we recover results from quantum mechanics, namely .

In the rest of the calculations, we work in the Fock basis , using the eigenvalues given by Eq. (31) and the eigenfunctions given by Eq. (32). The number operator is defined such that .

Next we describe the Generalised Heisenberg algebra that is used in the subsequent calculations of the position and momentum operators. The algebra of the ladder operators should be modified to account for the deformed commutators. Using the version of generalised Heisenberg algebra used in ref. [31], we find the action of the annihilation operator on an energy eigenstate to befor ν = βmℏω/2. The number operator is related to the ladder operators asAlso, the commutator is derived to beto first order in β.

Now we are ready to derive the expression for the position and momentum operators in terms of the ladder operators. These calculations will eventually enable the calculation of the trajectory of the oscillator. The relationship between and is obtained by starting with the Eq. (35) in the momentum basis, i.e,and using recursion relations of the Gegenbauer polynomials to express ψ(p) in terms of ψ(p). This gives us the operator in the momentum basis which can be expressed in terms of the operators and . This relation is then inverted to obtain in terms of the ladder operators.

Using the following recursion relationsof Gegenbauer polynomials, we obtain the action of the annihilation operators on the wavefunction to be

Since the ψ form a complete basis[13], Eq. (41) can be written in operator form in the limit of β ≪ 1 using the definitions of the position and momentum operators, (28) and (29), to obtain

We invert this relation to find and similarly in terms of the ladder operators.

Next we move to the definition of the generalised coherent state. To define the most classical pure state, a natural choice is a coherent state, but due to the Hamiltonian being modified from deformed commutators, our usual definition of coherent states no longer hold because the coherent state does not remain one after evolution under this Hamiltonian. Hence, here we introduce generalised coherent states that are suited for this modified Hamiltonian. These states are the Gazeau-Klauder states which were introduced in ref. [32]. Here, we describe them in detail, closely following the details in ref. [32].

Since a coherent state is parametrised by one complex number, we generalise it slightly by considering states parametrised by two real parameters . To ensure that behaves like a classical state, we demand that it satisfies the following conditions with respect to a given Hamiltonian H:

The continuity condition

Resolution of identity

Temporal stability such that the time-evolved state is always a generalised coherent state

The energy of the state only depends on J

The conditions in Eqs. (47) and (48) are defined with respect to a Hamiltonian, and so coherent states do not satisfy them with respect to a modified Hamiltonian and we need these generalised states.

If the eigenvalues and eigenstates of the Hamiltonian are defined such thatwe can verify that the definition of the generalised coherent statewheresatisfies all the above properties.

As an example, we consider the harmonic oscillator Hamiltonian. Here, the eigenvalues are given by (after ignoring the constants) and therefore in this case, e = n. Therefore, from the definitions of ρ and N(J), we see that the Gazeau-Klauder state iswhich is exactly the definition of a coherent state if we define . Thus we see that the Gazeau-Klauder state reduces to the coherent state when the Hamiltonian is the harmonic oscillator Hamiltonian.

For the Hamiltonian modified with deformed commutators, we see that the eigenvalues can be written as (31)and therefore, ignoring the constants,

Note that here the states for different t are not necessarily the same as , because ω is not related directly to the measured frequency but is merely a parameter in the Hamiltonian of Eq. (31).

Finally, we can calculate the trajectory of the oscillator. Using the rules of the generalised Heisenberg algebra, and the expressions for the position and momentum operators derived in the above section, we calculate the expectation value of position and momentum in this state and obtainand

Since these states satisfy the relation , the time-evolved expectation values are easily obtained by replacing γ with γ + ωt. In these calculations, we have assumed not only that β ≪ 1 but also βγ ≪ 1 and βωt ≪ 1.

We choose the initial state such that the oscillator starts at rest with non-zero amplitude, i.e., and . This condition is satisfied when γ = 0 and , as can be seen from Eqs. (56) and (55). Therefore, for this initial state, the expectation value of position is

The classical limit of Eq. (57) (ℏ → 0) satisfies the low amplitude limit of the differential equation (9) obtained classically, thus showing that we obtain identical results with both the methods used.

Potential future directions include connecting our formalism, which is based on the Gazeau-Klauder coherent states, with that of ref. [42], which starts with a somewhat different definition of creation and annihilation operators. Another important open problem is to go beyond coherent states as the initial states to thermal states. Our analysis can straightforwardly be applied to thermal states in principle, but the definition of a thermal state under the deformed commutators is unclear.

Data from experiment of ref. [26]

Table 1 contains the data extracted from the experimental results obtained in ref. [26]. This data was extracted by magnifying from Fig. 3 of the paper, setting the figure as a plot background and plotting markers over the figure such that the plot points coincide exactly with the points in the background. The coordinates of the plotted points were recorded as the extracted data.

Table 1

Measured data of the time-period of a pendulum as a function of its amplitude extracted from ref. [26].

A2 (cm2)	T2π (s)	A2 (cm2)	T2π (s)
43	3.47315	709	3.47468
52	3.47308	837	3.47498
132	3.47341	1020	3.47538
168	3.47342	1228	3.47583
204	3.47351	1404	3.47633
244	3.47363	1760	3.47705
293	3.47373	1850	3.47736
360	3.47396	2115	3.47801
387	3.47394	2160	3.47798
443	3.47409	2295	3.47847
578	3.47438