Prospects for Searching Thermal Effects, Non-Newtonian Gravity and Axion-Like Particles: Cannex Test of the Quantum Vacuum

We consider the CANNEX (Casimir And Non-Newtonian force EXperiment) test of the quantum vacuum intended for measuring the gradient of the Casimir pressure between two flat parallel plates at large separations and constraining parameters of the chameleon model of dark energy in cosmology. A modification of the measurement scheme is proposed that allows simultaneous measurements of both the Casimir pressure and its gradient in one experiment. It is shown that with several improvements the CANNEX test will be capable to strengthen the constraints on the parameters of the Yukawa-type interaction by up to an order of magnitude over a wide interaction range. The constraints on the coupling constants between nucleons and axion-like particles, which are considered as the most probable constituents of dark matter, could also be strengthened over a region of axion masses from 1 to 100 meV.


I Introduction
Since the development of quantum field theory, it has been appreciated that the quantum vacuum is a fundamental type of physical reality which potentially contains all varieties of elementary particles and their interactions. Although an infinitely large energy density of zero-point oscillations of quantum fields (the so-called virtual particles) could be considered to be catastrophic [1], convenient self-consistent procedures have been elaborated on how to make it equal to zero in the empty Minkowski space and take it into account when calculating the probabilities of arbitrary processes in the framework of the Standard Model. In doing so, the quantum vacuum is usually responsible for some part of the measured quantity (for instance, the Lamb shift), whereas the rest of it is determined by the ordinary (real) particles.
There is, however, one physical phenomenon, where the measured quantity is determined entirely by the quantum vacuum. This is the Casimir effect arising in the quantization volumes restricted by some material boundaries or in cosmological models with nontrivial topology [2][3][4][5]. An essence of this effect is that although the vacuum energy in restricted or topologically nontrivial volumes remains infinitely large, it becomes finite when subtracting the vacuum energy of empty topologically trivial Minkowski space. The negative derivative of the obtained finite vacuum energy with respect to the length parameter (either a separation between the boundary surfaces or a scale of the topologically nontrivial manifold) results in the Casimir force, which generalizes the familiar van der Waals force in the case of larger separations when one should take into consideration the effects of relativistic retardation [5].
Recently it was understood that the quantum vacuum may be responsible for an impressive phenomenon in nature, i.e., for acceleration of expansion of the Universe [6]. This can be explained by an impact of the energy density of quantum vacuum (which is often referred to as dark energy) corresponding to a nonzero renormalized value of the cosmological constant in the Einstein equations of general relativity theory [7].
Another subject is that at sufficiently short separations between the boundary surfaces, vacuum forces are stronger than Newtonian gravitation. In this case, they form a background for testing new physics, such as Yukawa-type corrections to Newton's law of gravitation arising due to exchange of light hypothetical scalar particles [8] or due to spontaneous compactification of extra spatial dimensions at the low-energy compactification scale [9]. Forces of this kind would alter the energy eigenstates of a neutron in the gravity potential of the Earth, and are searched for by a technique called Gravity Resonance Spectroscopy [10][11][12] by the qBounce

II Experimental Setup with Improved Precision
We recently demonstrated the feasibility of Casimir pressure gradient measurements with the parallel plate Cannex setup [33]. In the present article, we propose several improvements to the experiment that will lead to a significant increase in sensitivity and reduce the influence of systematic effects.
In our setup shown in Figure 1, the parallel plate geometry is implemented by a rigid vertical SiO 2 cylinder (lower plate) of radius R = 5.742 mm separated by a vacuum gap of width a from a movable (upper) sensor plate that can be described as a lumped mass-spring system [32]. This sensor with an interacting area A, elastic constant k, and effective mass m eff has a free resonance frequency ω 0 = k/m eff . In [33] an interferometric detection scheme was used to measure the gradient of the Casimir pressure. Here, we propose to implement the same scheme to measure both pressure gradients and the pressure acting between the two flat plates. Gradients P (a) ≡ ∂P/∂a are detected using a phase-locked loop (PLL) that senses shifts of the sensor's resonance frequency-a technique widely used in the literature [34][35][36][37][38]. In the absence of thermal drifts, the pressure P between our flat surfaces can be measured by monitoring the extension ∆a = P A/k of the spring-mass system using the same interferometer. For the determination and maintenance of parallelism, however, we recently used a feedback mechanism based on the capacitance between the two plates. While this scheme has been demonstrated to work in principle, practice has shown that vibrations in combination with the high Q-factor of our sensor lead to unacceptably long integration times, and associated susceptibility to thermal drifts [33]. To improve the performance and reach the full potential of the setup, we propose to replace the capacitive scheme by an interferometric one. As shown in Figure 1, three Fabry Pérot interferometers below the sensor plate (created by the end faces of optical fibers and the reflective surface of the sensor plate) monitor the distance a at different positions around the lower plate. (Please note that in the two-dimensional scheme in the figure, only two of these cavities are depicted.) The three fiber ends can be polished together with the lower plate to exactly match the surface of the latter. In such a system, the three lower interferometers could measure the tilt and frequency shift at all times. Synchronously, the upper (fourth) interferometer, being mechanically connected to the sensor frame, can be used to measure the extension ∆a of the sensor and, hence, the pressure acting on it. Besides this conceptual change several technical improvements are planned, which will result in the following: First, we aim to improve the seismic attenuation by a factor of 10 around f 0 = ω 0 /(2π) with respect to the present performance by means of active control techniques [39]. This would not only reduce the direct influence on measurements via non-linear effects and rms noise, but would also improve the stability of various feedback circuits in a nontrivial way. Second, a newly designed distributed thermal control concept will guarantee mK stability throughout the setup, thereby eliminating drift and uncertainty in the sensor characteristics. Third, electrostatic patch effects that create a systematic pressure background could be reduced using in situ Ar ion bombardment [38]. In the following, we analyze the major sources of experimental uncertainty and show the predicted effect of the mentioned improvements.

II.1 Sources of Error
Here, we consider experimental errors in the pressure and its gradient due to vibrations, the determination of displacement, and frequency shift of our sensor, variations in the temperature, tilt of the plates, and discuss the role of electrostatic patch effects.
Previously, we identified vibrations as a major source of error and gave respective tolerable limits for pressure measurements [33]. Using a five-axis seismic attenuation system, we achieved a damping factor of 60 in vertical direction around the sensor resonance f 0 = 10.24 Hz. The residual vibrations caused at times a peak displacement noise δa p ≈ 300 nm of the sensor relative to the lower plate, which proved to be a severe nuisance during the measurements. For pressure gradients, the main influence is via non-linear effects [34,40]. For a function φ(a) representing either the pressure (φ = P ) or its gradient (φ = ∂P/∂a), we expand for small sensor movements δa. Recognizing that δa ≡ 0, as usual for stochastic fluctuations, the main contribution to the error comes from the second order term. Thus, for the non-linear shifts in the pressure and its gradient, we therefore have respectively. Similarly, the rms pressure and pressure gradient (entering measurements mainly as a nuisance within the mechanical bandwidth of the sensor 2δf BW = 2.8 mHz around f 0 ) are computed from the spectral sensor movement δa(f ) via Preprint -3 We note that vibrational noise also hampers the convergence of various feedback circuits and thereby influences the achievable sensitivity in a nontrivial way. Experience has shown that such effects are negligible below a peak amplitude of ∼ 20 nm.
Another fundamental source of error is the uncertainty in various parameters involved in the evaluation of the frequency shift. We consider the following calibration procedure. At large separation a cal ≈ 80 µm we apply the AC electrostatic excitation voltages V ex and V AC driving the sensor resonance and the surface potential compensation circuit, respectively [33]. At this separation, the Casimir pressure P is negligibly small with respect to the electrostatic pressure P ES . The free resonance frequency is determined from a measurement of the resonance frequency ω 0P under the influence of the well-known electrostatic and Casimir forces. Here, ε 0 is the dielectric permittivity of vacuum. The effective mass m eff is determined during a separate sweep recording ω 0P as a function of an applied DC electrostatic potential. Eventually, we evaluate the resonance frequency shift measured at different separations to obtain the Casimir pressure gradient from Equation (1) after subtraction of the electrostatic pressure gradient ∂P ES /∂a. For pressure measurements at the same separations, we evaluate the detected extension ∆a of the sensor that is related to the total pressure P tot = P + P ES = ∆a ω 2 0 m/A. The achievable sensitivity for both types of measurement is limited by the uncertainties in all experimental quantities entering the respective evaluation. Numerical values for these and some other uncertainties considered below are given in Table 1.
Variations in the temperature contribute in two ways to the experimental error. First, different material expansion coefficients influence the separation between the two interacting surfaces by roughly 63 nm/K-an effect most influential at smaller a. Second, the temperature influences the Youngs modulus of our sensor, which leads to an additional error in a, but also offsets the resonance frequency and, thereby, mimics a pressure gradient.
We also consider errors due to tilt of the plates with respect to each other. For small angles δα between the plates, we may estimate the influence on the pressure and its gradient by averaging over the sensor area A with the local separation a(x, y, δα) deviating from its nominal value a due to the tilt. Numerical calculations show that the respective relative corrections to the pressure and its gradient are both of the order (δαR/a) 2 , in agreement with the literature [5]. For the achievable values of δa given in Table 1, these corrections are negligible. As has been discussed in Ref. [33], electrostatic effects can have a significant influence. While we recently compared our measurements with the model [41], we now use the model [42] that has been shown to describe the observed forces more realistically. The average patch size = 0.82 µm and the value for 2 mentioned in Table 1 were derived from the auto-correlation of actual Kelvin probe data for our surfaces. As is much smaller than the plate separation, we can use the approximation Here, V rms = 0.634 mV, and ζ(z) is the Riemann zeta function. Please note that P patch represents a systematic effect that can be characterized and removed from experimental data.

II.2 Sensitivity Estimation
Based on the models described in Section II.1, we have calculated the expected level of experimental uncertainty in measurements of the pressure and its gradient for both the present situation and presuming a successful implementation of the proposed improvements. For these calculations, we have assumed the parameters listed in Table 1. As can be seen in Figure 2a, the largest uncertainty in force measurements comes from the determination of the sensor extension, which includes contributions from the calibration and interferometry. Vibrations play, under the assumption that the proposed measures work as expected, a minor role. For force gradients ( Figure 2b) the determination of the frequency shift is the limiting factor at separations larger than ∼ 4 µm, while at smaller separations sensitivity is limited by vibrations. Temperature variations are influential on force gradient measurements for all separations, as they modify the Youngs modulus of our sensor. For comparison, we also plot the achievable uncertainty for the present version of the setup (dashed upper red lines), which are determined by the same factors as for the improved version. Please note that the patch pressure can   be characterized separately and removed from the data. The results for the improved uncertainties in Figure 2 are based on very conservative estimates. In the experiment, especially at larger separations, the uncertainty in the pressure could be reduced statistically by longer measurements. For further calculations in Section IV, we therefore assume a pressure sensitivity of 1 nN/m 2 (corresponding to 0.1 pN) at a ≥ 4 µm, which is slightly higher than the value corresponding to the bottom red line in Figure 2a. The uncertainty in the pressure gradient is mainly determined by the resolution of the PLL's frequency measurement. Here, the possibility for a statistical reduction of the uncertainty may be more problematic and therefore we use the pressure gradient sensitivity given by the bottom red line in Figure 2b for calculations in Section IV.

III Possibilities to Measure Thermal Effects in the Casimir Force
As mentioned in Section I, most of the already performed experiments on measuring the Casimir interaction exploited the sphere-plate geometry. In doing so, the measured quantity was either the Casimir force acting between a sphere and a plate (in the static measurement scheme) or its gradient (in the dynamic measurement scheme). Due to the proximity force approximation, the latter quantity can be recalculated into the effective Casimir pressure between two parallel plates [5,21]. There is only one modern experiment on the direct measurement of the Casimir pressure between two parallel plates [43], but it is not of sufficient precision to observe the thermal effects (an attempt to measure the Casimir effect between two parallel Al-coated plates at separations larger than a few micrometers was unsuccessful due to the presence of large background forces [44]). The distinctive feature of the Cannex test of the quantum vacuum is that it can be adapted for simultaneous measurements of the Casimir pressure between two parallel plates and its gradient (see Section II).
The Lifshitz formula for the Casimir pressure between two material plates spaced at a separation a at temperature T is given by [5,21] Here, it is assumed that the parallel plates made of a nonmagnetic material described by the dielectric permittivity ε(ω) are in thermal equilibrium with the environment at temperature T and the following notations are introduced. The Boltzmann constant is k B , the prime on the summation sign divides the term with l = 0 by 2, k ⊥ is the magnitude of the projection of the wave vector on the plane of the plates, ξ l = 2πk B T l/ with l = 0, 1, 2, . . . are the Matsubara frequencies , and q l = (k 2 ⊥ + ξ 2 l /c 2 ) 1/2 . The reflection coefficients r α are defined for two independent polarizations of the electromagnetic field, transverse magnetic (α = TM) and transverse electric (α = TE). Explicitly they are given by where By differentiating Equation (8) with respect to separation between the plates a, one obtains the gradient of the Casimir pressure Equations (8) and (11) take an exact account of the effects of finite conductivity of the plate metal. As to corrections due to surface roughness, at separations exceeding 3 µm they are much smaller than an error in the pressure measurements [5].
For numerical computations it is convenient to introduce the dimensionless variables In terms of these variables Equations (8) and (11) take the form and r TM (iζ l , y) = ε l y − y 2 + (ε l − 1)ζ 2 l ε l y + y 2 + (ε l − 1)ζ 2 l , r TE (iζ l , y) = y − y 2 + (ε l − 1)ζ 2 For application to the experimental setup of Cannex described in Section II, computations should be made for Au plates. The dielectric permittivity of Au along the imaginary frequency axis is obtained by means of the Kramers-Kronig relations using the available tabulated optical data for the complex index of refraction extrapolated down to zero frequency [5,21]. According to Section I, this extrapolation can be made either by means of the Drude model taking into account the relaxation properties of conduction electrons or the plasma model disregarding these relaxation properties. In dimensionless variables, the dielectric permittivity of the Drude model along the imaginary frequency axis is given by where the dimensionless plasma frequencyω p and relaxation parameterγ are connected with the dimensional ones byω p = 2aω p /c andγ = 2aγ/c. For Au the standard values ω p = 9.0 eV and γ = 35 meV are used here. The dielectric permittivity of the plasma model ε (p) l is obtained from Equation (16) by puttingγ = 0. Please note that for computations at separations a ≥ 3 µm performed below the optical data contribute negligibly small, so that the obtained gradients of the Casimir pressure are mostly determined by the extrapolations to lower frequencies.
Before presenting the computational results, we note that in the high-temperature limit T T cr = c/(2k B a) all the terms in Equations (13) and (14) with l ≥ 1 are exponentially small and both the Casimir pressure and its gradient are given predominantly by the terms with l = 0. The zero-frequency term of the Lifshitz formula takes different forms depending on the extrapolation used. If the plasma model ε where in accordance to Equation (15) r TE (0, y) = y − y 2 +ω 2 p y + y 2 +ω 2 p .
Calculating the first integrals on the right-hand side of Equations (17) and (18) and expanding the second ones in the powers of a small parameter 1/ω p , we find [5] If, however, the Drude model ε (D) l is used for extrapolation, one arrives at the result and Please note that at T = 300 K Equations (20) Taking into account, however, that the powers in all exponentially small terms with l ≥ 1 depend on −2πT /T cr , one can see that at T = 300 K Equations (20)- (22) lead to rather precise results at all separations exceeding 6 or 7 µm.
The computational results for the magnitude of the Casimir pressure and its gradient as functions of separation between the plates are presented in Figure 3 by the red and blue solid lines computed at T = 300 K using extrapolations of the optical data of Au by means of the plasma and Drude models, respectively. For comparison purposes, the computational results at T = 0 are presented by the dashed lines. They are obtained by Equations (13) and (14) where summation over the discrete Matsubara frequencies is replaced with a continuous integration. It is taken into account that with decreasing temperature the relaxation parameter γ quickly decreases towards a very small residual value at T = 0 that is determined by the defects of the crystal lattice. As a result, the values of P (a) and P (a) computed at T = 0 using the plasma and Drude models coincide at high accuracy.
For a better visualization, in the insets to Figure 3a,b the computational results in the separation region from 5 to 10 µm are presented using a uniform scale on the vertical axis. It is clearly seen that the theoretical predictions from using the plasma and Drude model extrapolations of the optical data can be discriminated if to take into account the errors in measuring P (a) and P (a) in the improved version of Cannex discussed in Section II. Please note that if the Drude model extrapolation is used the thermal effect in the Casimir pressure and its gradient vanishes at approximately 6.4 and 8.7 µm separations, respectively, where the blue and dashed lines intersect.
To determine the specific role of thermal effects in the Casimir pressure and its gradient, we have also computed the relative thermal corrections defined as where on the right-hand sides we indicated the dependence on temperature explicitly. Computational results for δ T P and δ T P are presented in Figure 4a,b as functions of separation by the red and blue lines for the cases when extrapolation of the optical data for Au to lower frequencies is made by means of the plasma and Drude models, respectively. As is seen in these figures, within the separation range from 3 to 15 µm the thermal effects make a considerable contribution to the Casimir pressure and its gradient. If the plasma model extrapolation is used, it reaches 440% and 240% of the zero-temperature pressure and its gradient, respectively, at a = 15 µm. When using the extrapolation by means of the Drude model, the relative thermal effect in the Casimir pressure varies from approximately −34% at a = 3 µm to 130% at a = 15 µm and in the pressure gradient from approximately −40% at a = 4 µm to 60% at a = 15 µm. In this case, the thermal effect in the Casimir pressure vanishes at a ≈ 6.4 µm and in its gradient at a ≈ 8.7 µm in agreement with Figure 3a,b, respectively .

IV Prospective Constraints on Non-Newtonian Gravity and Axion-Like Particles
As mentioned in Section I, the Casimir forces originating from the quantum vacuum form a background for testing the Yukawa-type corrections to Newton's gravitational law and for searching the axion-like particles. It is convenient to parametrize the Yukawa-type potential between two point masses m 1 and m 2 situated at the points r 1 and r 2 as [8] Here, α and λ are the interaction constant and the range of Yukawa interaction, and G = 6.674 × 10 −11 m 3 kg −1 s −2 is the Newtonian gravitational constant (we note that in the experimental configuration under consideration one can neglect by the Newtonian gravitational pressure because it is less than an error in measurements of the Casimir pressure).
The Yukawa-type pressure between the top and bottom plates in the experimental setup of Section II should be calculated taking into account the layer structure of both plates shown in Figure 1. The top plate is made of high-resistivity Si of density ρ Si = 2.33 g cm −3 coated with a layer of Cr of density ρ Cr = 7.15 g cm −3 and thickness d Cr of SiO 2 quartz crystal with the density ρ q = 2.64 g cm −3 coated with a layer of Cr of thickness d Cr = 5 nm followed by a layer of Au of thickness d b Au = 1 µm. The thickness of SiO 2 substrate (6 mm) again allows to consider it as a semispace. Now we assume that one mass ρdV 1 belongs to the top plate and another one ρdV 2 to the bottom one and integrate Equation (25) over the volumes of both parallel plates separated by a distance a taking into account their layer structure. Calculating the negative derivative of the obtained interacting energy with respect to a, we find the Yukawa force and finally the pressure [45] P Yu (a) = −2πGαλ 2 e −a/λ Φ(λ), (26) where the function Φ is defined as The gradient of the Yukawa pressure is obtained by differentiating Equation (26) with respect to a P Yu (a) = 2πGαλe −a/λ Φ(λ).
Now the constraints on the parameters α and λ of Yukawa-type interactions can be obtained from the inequalities |P Yu (a)| < ∆P (a), where ∆P and ∆P are the total experimental errors in the measured Casimir pressure and its gradient estimated in Section II. The meaning of Equation (29) is that the experimental data are found in agreement with theoretical predictions for the Casimir pressure and its gradient and no extra contribution of unknown origin was observed.
To estimate the strength of prospective constraints, which can be obtained from the improved version of Cannex, we use the thickness of the top and bottom Au layers d t Au = 200 nm and d b Au = 1 µm and the total experimental errors ∆P (a) = 1 nN/m 2 and ∆P (a) corresponding to the bottom red line in Figure 2b. The computational results for α and λ obtained by using the first and second inequalities in Equation (29) are shown by the red lines in Figure 5a,b, respectively. In so doing, the regions of (α, λ)-planes above each line are excluded and below each line are allowed. For comparison purposes, we present in Figure 5a,b the strongest constraints obtained in the same interaction range from other experiments. The top and bottom blue lines demonstrate the constraints found from the results of the Casimir-less experiment [46] and its improved version [47], respectively. Both experiments were Preprint For citations please use the journal reference: Symmetry, 11(3) 407 (2019).
Prospects for Searching Thermal Effects, Non-Newtonian Gravity and Axion-Like Particles: Cannex Test of the Quantum Vacuum G. L. Klimchitskaya, V. M. Mostepanenko, R. I. P. Sedmik, and H. Abele performed by means of a micromechanical torsional oscillator. The green line shows the constraints obtained from measuring the difference in lateral forces [48]. The constraints indicated by the black line are found in [49] from the torsion pendulum experiment. Finally, the constraints of the gray line follow from the Cavendish-type experiments performed at short separations [50][51][52].
As is seen in Figure 5a,b, the strongest constraints obtained up to date follow from the improved Casimir-less experiment [47] (the bottom blue line). It is seen also that the largest strengthening of these constraints, which could be reached from the Cannex test of quantum vacuum, follows from measurements of the Casimir pressure (this is because the pressure gradient is linear in λ whereas the pressure is quadratic in λ). The prospective constraints are stronger by up to a factor of 10 over a wide interaction range.
Next we consider the prospective constraints on the axion-to-nucleon coupling constants which could be obtained from the Cannex test of the quantum vacuum. Taking into account that the exchange of one axion between two nucleons results in the spin-dependent interaction potential [53] and the test bodies in this experiment are not polarized, any additional force of the axion origin could arise due to two-axion exchange.
Below we deal with axion-like particles coupled to nucleons by means of pseudo-scalar interaction Lagrangian [54]. In this case, the effective potential between two nucleons, spaced at the points r 1 and r 2 of the top and bottom plates, arising due to exchange of two axions, takes the form [53,55] V an (r) ≡ V an (|r 1 Here, g an is the dimensionless coupling constant of an axion to a nucleon (we assume that the coupling constants to a neutron and a proton are equal [53]), the mean of the proton and neutron masses is denoted as m, the axion mass is m a , K 1 (z) is the modified Bessel function of the second kind and it is assumed that r /(mc). Similar to the case of a Yukawa-type potential, the additional pressure between the test bodies due to two-axion exchange, is obtained by integrating Equation (30) over the volumes of both plates with account of their layer structure, calculate the negative derivative of the obtained result with respect to a and finally find the pressure (see [15] for details) Here, m H is the mass of atomic hydrogen and the function Ψ(x) is defined as where the coefficients C for each material are given by In Equation (33), ρ is the density of the respective material already indicated above, Z and N are the number of protons and the mean number of neutrons in the atoms of materials with mean mass M and µ = M/m H . For materials under consideration Z/µ = 0.40422, 0.46518, 0.50238, and 0.503205 for Au, Cr, Si, and SiO 2 , respectively [8]. The values of N/µ are 0.60378, 0.54379, 0.50628, and 0.505179 for the same respective materials [8].
The gradient of the pressure due to two-axion exchange is obtained from Equation (31) The constraints on the parameters of axion-nucleon interaction, g an , m a , are obtained from the inequalities |P an (a)| < ∆P (a), P an (a) < ∆P (a), similar to Equation (29). To find the strongest prospective constraints on the axion-to-nucleon interaction, we have used the same parameters of the improved setup, as listed above when considering the Yukawa interaction. The computational results for g 2 an /(2π) as a function of m a c 2 obtained by using the first and second inequalities in Equation (35) Figure 6 : Constraints on the axion-to-nucleon coupling constant obtained from the improved Casimir-less experiment, the Cavendish-type experiment, measurements of the effective Casimir pressure, and the lateral Casimir force are shown by the blue, gray, black, and green lines, respectively. The proposed constraints which could be obtained from the Cannex test of the quantum vacuum when measuring (a) the Casimir pressure and (b) its gradient are indicated by the red lines.

V Discussion
In the foregoing, we have considered the Cannex test of the quantum vacuum intended to measure the gradient of the Casimir pressure between two parallel plates at separations exceeding a few micrometers. As already discussed in the literature [32,33], this experiment would be capable to place stronger constraints on the chameleon model of dark energy and, thus, bring new information concerning the origin of vacuum energy and the value of the cosmological constant. Here, we proposed a modification of the Cannex setup, which allows for simultaneous measurements of both the pressure and its gradient. We also considered several improvements in the already existing setup which will allow for a more precise determination of several parameters. According to our results, with these improvements the Cannex test of the quantum vacuum could be used to directly measure for the first time the thermal effects in the Casimir pressure and its gradient at separations from 5 to 10 micrometers. We have also shown that this experiment could strengthen the presently known best constraints on the parameters of non-Newtonian gravity by up to a factor of 10 over a wide interaction range. The constraints on the axion-to-nucleon coupling constant could be strengthened by up to a factor of 3 in the region of axion masses from 1 to 100 meV.

VI Conclusions
To conclude, Cannex can be considered as a promising new laboratory experiment for investigating unusual features of the quantum vacuum, thermal effects in the Casimir forces, non-Newtonian gravity and properties of axion-like particles being hypothetical constituents of dark matter. Together with other laboratory experiments, such as the demonstration of the dynamical Casimir effect in superconducting circuits [62,63], the large-scale accelerator projects, and astrophysical data, it may serve as a useful tool in resolving fundamental problems of the quantum vacuum.