The State of the Art in Constraining Axion-to-Nucleon Coupling and Non-Newtonian Gravity from Laboratory Experiments

Constraints on the Yukawa-type corrections to Newton's gravitational law and on the coupling constant of axionlike particles to nucleons obtained from different laboratory experiments are reviewed and compared. The constraints on non-Newtonian gravity under discussion cover the wide interaction range from nanometers to millimeters and follow from the experiments on neutron scattering, measuring the Casimir force and Cavendish-type experiments. The constraints on the axion-to-nucleon coupling constant following from the magnetometer measurements, Cavendish-type experiments, Casimir physics, and experiments with beams of molecular hydrogen are considered which refer to the region of axion masses from $10^{-10}$ eV to 200 eV. Particular attention is given to the recent constraints obtained from measuring the Casimir force at nanometer separation distance between the test bodies. Several proposed experiments focused on constraining the non-Newtonian gravity, axionlike particles and other hypothetical weakly interacting particles, such as chameleons and symmetrons, are discussed.


Introduction
In recent years, more and more experiments are devoted to a search for some uncharged weakly interacting elementary particles which are not accessible for the accelerator techniques used in high energy physics. Many particles of this kind were predicted in different extensions of the Standard Model, supergravity and supersymmetry [1,2], and their discovery would have a profound effect on several branches of physics and change our picture of the world. Thus, an exchange of light scalars (for instance, dilatons [3]) between atoms of neighboring macrobodies results in the Yukawa-type correction to the Newtonian gravitational law and opens up a new area to the experimenter: the search for non-Newtonian gravity [4]. The predicted light pseudoscalar particles called axions [5,6] not only provide a plausible explanation for the problem of strong CP violation and large electric dipole moment of a neutron in QCD, but suggest themselves as possible constituents of dark matter [7].
The Yukawa-type interaction at short separations between the test masses has long been searched for as the so-called fifth force in the gravitational experiments of Eötvos and Cavendish-type [4]. According to the results of [8,9], both the Yukawa-and power-type corrections to Newtonian gravity can be recorded also in experiments on measuring the van der Waals and Casimir forces caused by the zero-point and thermal fluctuations of the electromagnetic field. The development of extradimensional theories, which anticipate the compactification scale at energies of the order of 1 TeV [10,11], made one more prediction of the Yukawa-type corrections to Newton's law at short separations [12,13] and stimulated their search.
An exchange of one pseudoscalar particle between two fermions results in the spin-dependent interaction potential [14,15] which turns into zero after an averaging over the volumes of two unpolarized macroscopic bodies. Because of this, the process of one-axion exchange does not lead to any fifth force. An effective potential due to a simultaneous exchange of two axions is more complicated and depends on the details of an axion-to-fermion interaction. If this interaction is described by the pseudoscalar Lagrangian, the spin-independent effective potential leading to some fifth force between two macroscopic bodies is obtained [14,15]. If, however, an axion-to-fermion interaction is described by the pseudovector Lagrangian, as it holds for QCD axions which are the pseudo-Nambu-Goldstone bosons, the character of an effective potential due to the exchange of two axions remains a mystery [16].
The constraints on the constant of Yukawa-type interaction α are usually obtained as functions of the interaction range λ = M −1 , where M is a mass of the exchange scalar particle (here and below we puth = c = 1). For non-Newtonian gravity originating from extra dimensions λ has the meaning of the compactification scale. Depending on the interaction range, the strongest constraints on α follow from different experiments. Over a long period of time, experiments on measuring the Casimir force were used for obtaining the strongest constraints on α over the wide interaction range below a few micrometers [17]. For larger λ up to a few centimeters the strongest constraints on α were found from checking the inverse square law of Newtonian gravity (the Cavendish-type experiments [18][19][20]) and for even greater λ from tests of the equivalence principle (the Eötvos-type experiments [21,22]). At the moment, the strongest constraints on α in the interaction range below approximately 10 nm follow from the experiments on neutron scattering [23][24][25]. The Casimir force measurements (and their modification where the Casimir force is nullified [26]) allow obtaining the strongest constraints on α within the interaction range from approximately λ = 10 nm to λ = 8 µm. For larger λ the gravitational experiments preserve their leading role in constraining the Yukawa-type corrections to Newtonian gravity.
Constraining the parameters of axions and axionlike particles is the subject of extensive studies. There are many types of axions in the literature grouped in the models of hadronic (QCD) axions [27,28] and Grand Unified Theory (GUT) axions [29,30] which in turn are divided into many submodels. Axionlike particles can interact with photons, electrons, and nucleons, and each interaction channel is used in the laboratory experiments and astrophysical observations dedicated to the axion quest (see reviews [1,7,[31][32][33]). Our interest here is focussed on the interaction of axions with nucleons. The constraints on the strength of an axion-to-nucleon interaction constant g an are usually obtained as functions of the nucleon mass m a . In fact the strongest constraints on g an follow from astrophysical arguments on stellar cooling accounted for by the axion emission [34]. It has been known, however, that the theory of dense nuclear matter suffers from significant uncertainties which make the obtained constraints not fully reliable [34]. At the moment, the strongest laboratory constraints on g an in the range of of m a < 1 µeV follow from the magnetometer measurements [35] and for heavier axions from the Cavendish-type [19,36] and Casimir-less [26,37] experiments. For m a > 0.5 eV the strongest laboratory limits on g an were obtained from the experiment with beams of molecular hydrogen [38,39] and for m a > 200 eV from the experiment on nuclear magnetic resonance [39] (see also review [40] of the laboratory constraints on both the axionlike particles and Yukawa-type interactions with a focus on neutron experiments).
In this article, we review constraints on the Yukawa-type interaction and interaction of axions with nucleons obtained during the last few years from measurements of the Casimir force. These constraints are compared with the strongest limits following from the gravitational and some other laboratory experiments in the same interaction ranges. We also discuss several proposals focussed on further strengthening the constraints on both non-Newtonian gravity and axionlike particles.
The article structure is the following In Section 2, we consider constraints on the strength of Yukawa-type correction to Newtonian gravity obtained from different measurements of the Casimir force and Casimir-less experiments. In Section 3, the same is done with respect to the constant of axion-to-nucleon interaction. Section 4 presents the constraints on both non-Newtonian gravity and axionlike particles obtained recently from measuring the Casimir force at nanometer separation range. These constraints are compared with those outlined in Sections 2 and 3 and with the most strong results following from some other laboratory experiments. In Section 5, several proposed experiments are discussed allowing further strengthening of the obtained constraints. Section 6 contains our discussion. Finally, in Section 7 the reader will find the conclusions.

The Yukawa-Type Correction to Newtonian Gravity and Constraints on It from the Casimir Effect
The Yukawa-type correction to the Newtonian gravitational potential between two particles with masses m 1 and m 2 situated at the points r 1 and r 2 is conventionally written in the form [4] where r 12 = |r 1 − r 2 |, α is the dimensionless interaction constant of the Yukawa interaction, and G is the Newtonian gravitational constant. As mentioned in Section 1, the interaction range λ has a meaning of either the Compton wavelength of an exchange scalar particle of mass M, λ = M −1 , or the compactification scale of extra spatial dimensions. For two test bodies V 1 and V 2 at the closest separation a the Yukawa interaction energy is obtained by interacting (1) over their volumes where ρ 1 and ρ 2 are the respective mass densities. The Yukawa-type force acting between the two bodies can be calculated as the negative derivative of (2) It should be noted that although the Yukawa-type force is often called a "correction" to Newtonian gravity, at short separations below a micrometer such an experimentally allowed "correction" may exceed the gravitational force by many orders of magnitude. For this reason, when obtaining the constraints, it is customary to omit the contribution of the genuinely Newtonian gravitation. This does not make any impact on the obtained results.
We begin with constraints found at the shortest interaction range from measuring the lateral Casimir force acting between the sinusoidally corrugated surfaces. It is common knowledge that the Casimir force between two parallel plates or between a plate and a sphere acts in perpendicular direction to the surfaces. However, for surfaces covered with the sinusoidal corrugations of a common period but having some nonzero phase shift Φ (see Figure 1) the Casimir free energy F C (a, Φ, T) depends not only on a and on temperature T, but also on Φ. This leads to initiation of the lateral Casimir force The sinusoidally corrugated surfaces are subjected to both the normal Casimir force acting perpendicular to the surface and the lateral Casimir force acting in a tangential direction. Due to the presence of corrugations, the closest separation between the two surfaces is smaller in this case than for the smooth bodies. Because of this, by measuring both the lateral and normal forces between corrugated surfaces it becomes possible to constrain hypothetical interactions within the shortest interaction range. Here, we discuss the constraints obtained from experiments on measuring both the lateral and normal Casimir forces starting from the lateral one. The constraints found from measurements of the normal Casimir force acting between the smooth surfaces are also considered.
The lateral Casimir force (4) was first predicted theoretically in [41] and observed experimentally in [42,43]. In [44,45] it was measured more precisely in the configuration of a corrugated sphere in front of a corrugated plate with increased amplitudes at room temperature. The schematic of the experimental setup is shown in Figure 2 [44]. A sinusoidally corrugated grating was vertically mounted on the piezoelectric tube. A smooth sphere was placed at the end of the cantilever of an atomic force microscope. A mica sheet was attached to the bottom of this sphere. The second sphere was attached to the bottom free end of the mica sheet. Then, the sinusoidal corrugations were imprinted on the second sphere from the grating as a template using a special procedure [44]. This ensures a parallelism of the corrugation axes on both test bodies. The interacting surfaces were covered with an Au layer.
The lateral Casimir force F C, lat was measured at different values of the phase shift Φ i at different separations from 120 to 190 nm, and the total experimental errors ∆ i F C, lat have been determined at the 95% confidence level. The lateral Casimir force between the sinusoidally corrugated surfaces was also calculated using the generalization of the Lifshitz theory based on the Rayleigh scattering approach, and good agreement with the measurement data in the limits of the errors ∆ i F C, lat was observed. This means that the lateral force of Yukawa-type, which also should arise due to the presence of corrugations, satisfies the inequality The hypothetical force F Yu,lat in the experimental configuration of Figure 2 was calculated in [46] using (1) and (2), where a dependence on Φ arises for the corrugated plates, and subsequent negative differentiation with respect to Φ similar to (4). This gave the possibility to determine the value of α and λ which satisfy (5). The obtained results are shown in Figure 3 by the line 1, where here and below the region of the (λ, α) plane above the line is excluded and the region below it is allowed by the experimental results. The constraints of the line 1 are stronger by up to a factor of 2.4 × 10 7 than the previously known constraints obtained at a high confidence level [46] (although the strength of constraints obtained in [47] from the measurement data of [48] is slightly stronger, their confidence level remains indefinite [17,49]).
At larger interaction range λ > 11.6 nm the stronger constraints on the strength of Yukawa interaction follow from measurements of the normal Casimir force between sinusoidally corrugated surfaces of a sphere and a plate at different orientation angles of corrugations [50,51]. The schematic of the experimental setup is shown in Figure 4 [50]. In this case, unlike Figure 2, a sinusoidally corrugated grating was mounted on the piezoelectric tube laterally. An initially smooth sphere was attached to the end of the cantilever of an atomic force microscope. The sinusoidal corrugations on it were imprinted from the grating as a template. Both interacting surfaces were coated with an Au layer.
The normal Casimir force F C, nor was measured at different separations from 127 to 200 nm for different angles θ i between the axes of corrugations from 0 to 2.4 • at T = 300 K. The total experimental errors ∆ i F C, nor have been determined at the 67% confidence level. However, we have recalculated them to the 95% confidence level in order to obtain constraints consistent with other results discussed here. In the limits of these errors the measurement data were found in good agreement with the theoretical Casimir force calculated at T = 300 K using the method of derivative expansion [52,53]. Thus, the perpendicular to the surface Yukawa force should satisfy the inequality The Yukawa force in the configuration of Figure 4 was calculated in [54]. The resulting constraints on α and λ are shown by the line 2 in Figure 3. They were considered as the strongest ones up to 17.2 nm [55].  At λ > 17.2 nm the stronger constraints have been obtained from measuring the effective Casimir pressure between two parallel plates by means of a micromechanical torsional oscillator [56,57]. The schematic diagram of this measurement is shown in Figure 5 [17]. The oscillator consisted of an Au-coated plate suspended at two opposite points on the midplane. Under the plate there were two electrically contacted electrodes used to measure the capacitance and to induce oscillation in the plate at the resonant frequency of the oscillator. Above the plate there was an Au-coated sphere of radius R attached to the side of an Au-coated optical fiber. The separation between the plate and the sphere was determined as where the meanings of z meas , D 1 , D 2 , θ, and b are clear from Figure 5. In the dynamic regime used in this experiment, an immediately measured quantity was a shift of the resonant frequency of the oscillator under an impact of the Casimir force F C,sp acting between a sphere and a plate. The solution for the equation of motion of this oscillator allows recalculation of the frequency shift into the gradient of the Casimir force ∂F C,sp /∂a. Using the proximity force approximation [17], which is valid in this case up to a fraction of a percent [58,59], one obtains the effective Casimir pressure between two parallel plates The effective Casimir pressure (8) was measured within the separation range from 162 to 746 nm at T = 300 K, and the total experimental errors ∆ i P C were determined at the 95% confidence level. A very good agreement was found between the measurement data and theoretical predictions of the Lifshitz theory using the optical data for the complex index of refraction of Au extrapolated down to zero frequency by the dissipationless plasma model [56,57]. An alternative extrapolation using the dissipative Drude model was excluded by the data. A discussion of both extrapolations can be found in [17,49,60]. Later an exclusion of the Drude extrapolation was conclusively confirmed by the experiment [61] where the theoretical predictions from two alternative extrapolations differ by up to a factor of 1000. We Shematic setup Figure 5. Schematic diagram of the experimental setup using a micromechanical oscillator for determination of the effective Casimir pressure between two parallel plates (see the text for further discussion). emphasize that a difference between predictions of the excluded and confirmed versions of the Lifshitz theory cannot be modeled by the Yukawa force [62]. This means that the effective Yukawa pressure between two parallel plates should satisfy the inequality The Yukawa-type pressure between two parallel Au-coated plates was calculated in [56,57] and used for constraining the values of α and λ in accordance to (9). The obtained constraints are shown by line 3 in Figure 3. For some time these constraints were considered as the strongest ones up to λ = 89 nm.
At larger λ the strong constraints on α were obtained from the so-called Casimir-less experiment where an impact of the Casimir force was nullified [63]. This experiment was performed using the micromechanical oscillator shown in Figure 5. The plate of this oscillator was, however, composite. It consisted of two halves coated by Au and Ge films with a common Au overlayer. The oscillator was moved back and forth below the sphere, so that only a difference of forces from two halves of the plate has been measured. In doing so a common Au overlayer was sufficiently thick to make the Casimir forces from two halves of the plate equal. As to the Yukawa-type forces between the sphere and the plate halves coated with Au and Ge films, they were predicted to be different. The sensitivities ∆ i F diff of this setup to the differential force determined at the 95% confidence level were of the order of 1 fN at separations a i varying from 150 to 500 nm. This means that the difference of Yukawa forces from the Au and Ge films should satisfy the inequality Both these Yukawa forces were calculated and used for constraining of α and λ in [63]. The constraints following from (10) are shown by line 5 in Figure 3 (line 4 is discussed in the end of this section). They were the best ones within the range of λ from 89 to 465 nm.
At larger λ stronger constraints were obtained by sensing the difference of lateral forces between an Au sphere attached to the end of soft cantilever and a density modulation source mass [64]. To make a cantilever sensitive to the lateral force, it was placed normal to the source-mass surface. The obtained results are shown by line 6 in Figure 3. At λ = 2.09 µm they are replaced by the stronger ones found from measurement of the Casimir force between the Au-coated surfaces of a plate and a spherical lens of centimeter-size radius by means of the torsion pendulum [65]. These constraints are shown by line 7 in Figure 3. They extend from λ = 2.09 µm to λ = 3.16 µm. For larger λ the strongest constraints on the Yukawa-type interaction follow from the gravitational experiments (see Section 4).
In the end of this section, we present the results of the second Casimir-less experiment [26] which improved the results of the first one [63] significantly and resulted in stronger constraints over a wider interaction range than those ones obtained from several experiments listed above. As compared to the first Casimir-less experiment, the experimental setup for measuring the difference of forces was considerably upgraded. A differential force was measured between an Au-coated sphere and either an Au or a Si sectors of the rotating disc coated with an Au overlayer of sufficient thickness. In this way, the difference of Casimir forces between the sphere and the two neighboring sectors was again nullified, and the measurement result was determined by only the difference of Yukawa forces between the same sphere and either an Au or a Si sectors. Furthermore, the force sensitivity of the setup was improved, so that the minimum detectable force ∆ i F diff varied from 0.1 to 0.2 fN when separation distance increased from 200 to 1000 nm.
The constraints on the Yukawa-type interaction have been again derived from inequality (10) where the lower index Ge should be replaced with Si. The obtained results [26] are shown by line 4 in Figure 3. It is seen that they extend the previously known constraints over a wide interaction range and strengthen them by up to a factor of 10 3 . The constraints of line 4 are the strongest ones in the interaction region from λ = 40 nm to λ = 8 µm. Thus, for λ > 40 nm they improve the constraints following from measurements of the Casimir pressure and all other constraints obtained from measuring the Casimir force at larger λ.
A summary of the strongest constraints on the Yukawa-type interaction over the interaction range from 1 nm to 1 mm obtained from measurements of the Casimir force and other laboratory experiments is presented in Section 4.

Constraints on the Coupling of Axions to Nucleons from the Casimir Effect
As was mentioned in Section 1, the effective potential due to an exchange of one axion between two nucleons depends on their spins σ 1 and σ 2 . An explicit form of this potential is given by [14,15,66] V (1) e −m a r 12 . (11) Here, we assume that the dimensionless coupling constants of an axion to a proton and a neutron are approximately equal g an = g ap . The nucleon mass is notated m, the axion mass is m a , and k = r 12 /|r 12 | is the unit vector. Equation (11) is valid for both the QCD and GUT axions. The effective potential (11) averages to zero after an integration over the volumes of two macrobodies. Because of this, the process of one-axion exchange does not lead to any additional force which could be searched for in the Casimir experiments (an experiment on measuring the Casimir force between two polarized test bodies has been proposed [67] but is not realized yet). If, however, the axion-to-nucleon interaction is described by the pseudoscalar Lagrangian characteristic for GUT axions, a simultaneous exchange by two axions results in the spin-independent effective potential [14,15,68] V (2) an (r 12 ) = − where K 1 (z) is the modified Bessel function of the second kind. This result, which is valid under the condition r 12 ≫ m −1 , resembles the Yukawa-type potential (1). As to the pseudovector Lagrangian, which can be also used to describe the interaction of axions to nucleons, it is characteristic for the QCD axion models. For this Lagrangian, the form of effective potential describing the two-axion exchange remains unknown [16]. The results presented below in this section are based on the use of (12) and, thus, are applicable to the GUT axions and axionlike particles. Similar to the case of Yukawa-type forces, the interaction energy of two test bodies V 1 and V 2 due to an exchange of two GUT axions between their protons and neutrons is given by an m a 32π 3 m 2 n 1 n 2 where n 1 and n 2 are the numbers of nucleons per unit volume of the first and second test bodies, respectively. Then, the additional force arising due to two-axion exchange is given by Different experiments on measuring the Casimir force were used to constrain the axion-to-nucleon interaction in the same way as described in Section 2 for the interaction of Yukawa type. We begin with the second Casimir-less experiment [26] (see Section 2). Its results were used for constraining axions in [37]. For this purpose, the difference of additional forces arising in the experimental configuration due to two-axion exchange was found using (13) and (14). Taking into account that no force was registered up to the sevsitivities ∆ i F diff varying from 0.1 to 0.2 fN, this difference should satisfy the inequality The constraints on an axion-to-nucleon interaction are usually displayed graphically as g 2 an /(4π) depending on the axion mass m a . Similar to the Yukawa-type interaction, the region of the [m a , g 2 an /(4π)] plane above the line is excluded by the results of respective experiment and below the line is allowed. By the line 1 in Figure 6 we show the constraints on the axion mass and interaction constant following from the Casimir-less experiment. Below they are compared with similar constraints following from several other experiments.
Next, we consider an experiment on measuring the effective Casimir pressure by means of a micromechanical oscillator shown schematically in Figure 5 [56,57]. This experiment was used for constraining the parameters of axions in [69]. For this purpose, an additional pressure P (2) an was calculated arising in the experimental configuration due to the two-axion exchange. The constraints were obtained from the inequality |P where the meaning of ∆ i P C is explained in Section 2. These constraints are shown by the line 2 in Figure 6. As is seen in Figure 6, the constraints on an axion following from measuring the Casimir pressure become stronger that those found from the Casimir-less experiment for axion masses m a > 0.9 eV. For axion masses m a > 8 eV the stronger constraints on g an follow from measurements of the lateral Casimir force between corrugated surfaces (see Figure 2 for a schematic of the experimental setup and [44,45] for a description of the measurement procedure). To obtain these constraints, the additional force F (2) an,lat was calculated in [70]. Then the constraints were obtained from the inequality where ∆ i F C,lat is the total experimental error in measuring the lateral Casimir force (see Section 2 for details). These constraints are shown by the line 3 in Figure 6. In addition to the constraints shown by the lines 1-3, which are the strongest ones over the respective regions of axion masses, in Figure 3 we also present the constraints obtained from three more Casimir experiments. Although these experiments lead to somewhat weaker constraints, they are useful by providing a qualitatively similar picture based on the independent measurement data.
The line 4 in Figure 6 was obtained in [71] from the measure of agreement between experiment and theory in measuring the difference in Casimir forces between a Ni-coated sphere and either an Au or a Ni sectors of a rotating disc coated with a thin Au overlayer [61]. As opposed to the Casimir-less experiments [26,63], the Au overlayer was sufficiently thin in this case and did not make equal the Casimir forces when a sphere was spaced above an Au and a Ni sectors. It is seen that the line 4 is sandwiched between the lines 1 and 2 discussed above. Thus, the constraints following from measuring the difference in Casimir forces are in agreement with those obtained from the Casimir-less experiment and experiment on measuring the Casimir pressure.
The line 5 in Figure 6 was obtained in [72] from the dynamic experiment measuring the gradient of the Casimir force between an Au-coated sphere and an Au-coated plate performed by means of the dynamic force microscope [73]. Finally, the line 6 in Figure 6 follows [74] from the experiment on measuring the Casimir-Polder force between 87 Rb ultracold atoms belonging to the Bose-Einstein condensate and a silica glass plate [75]. As is seen in Figure 6, the lines 5 and 6 indicate the weaker constraints on the coupling constant of axions to nucleons than those shown by the line 2 following from measuring the effective Casimir pressure.
In the next section the constraints of Figure 6 are compared with some other constraints obtained from the laboratory experiments unrelated to measurements of the Casimir force.

Constraints from Measuring the Casimir Force in Nanometer Separation Range and Other Laboratory Experiments
Here, we consider the constraints on the Yukawa-type hypothetical interaction and on the coupling constant of axions to nucleons obtained recently [76] from measuring the Casimir force between an Au-coated sphere and a silicon carbide (SiC) plate at T = 300 K in N 2 atmosphere [77]. This experiment has the advantage that the surfaces of both interacting bodies are very smooth. This allowed to measure the Casimir force by means of an atomic force microscope over the separation region from 10 to 200 nm. The experimental results for the Casimir force were compared with theoretical predictions of the Lifshitz theory. For this purpose, the optical properties of SiC plate have been determined by means of ellipsometers [78], as well as the optical properties of Au coating on the sphere which were characterized earlier [79]. Both the experimental and theoretical errors were determined at the 67% confidence level. In [76] the total measure of agreement between experiment and theory ∆ i F C was found at the 95% confidence level. Within these errors no additional force was registered originating either from the non-Newtonian gravity or from the axion-to-nucleon interaction.
We begin with the Yukawa-type interaction (1). From (1)-(3) one obtains the following expression for the Yukawa force between an Au-coated sphere and a plate [76]: where ρ 1 and ρ 2 are the densities of Au and SiC, respectively, d is the thickness of an Au coating on the sphere, and the function Ψ is given by In the derivation of (18) it was taken into account that the sphere core material (borosilicate glass) gives only a minor contribution to the Yukawa force which does not influence the obtained results. The constraints on the Yukawa-type correction to Newtonian gravity have been found from the inequality They are shown by the line labeled "new" in Figure 7. In the same figure we reproduce from Figure 3 the lines 1, 2, 3, and 4 obtained from measuring the lateral and normal Casimir forces between corrugated surfaces, the Casimir pressure, and from the second Casimir-less experiment, respectively. As is seen in Figure 7, the new constraints are stronger than those of line 1 within the interaction range from 1 to 3.7 nm. The largest strengthening by the factor of 5 × 10 5 is reached at λ = 1 nm. Now we compare the constraints on the Yukawa interactions obtained from the Casimir effect and from some other laboratory experiments. The strongest constraints on α in the nanometer interaction range have been found from the experiments on neutron scattering [23,24]. They are indicated in Figure 7 by the two lines labeled "n". As is seen in Figure 7, the neutron scattering results in the strongest constraints on α up to λ = 10 nm.
Over a wide interaction range from λ = 10 nm to λ = 8 µm the strongest constraints on α follow at the moment from the Casimir physics. In doing so measurements of the lateral and normal Casimir forces between corrugated surfaces lead to the most strong results only within the very narrow ranges from 10 to 11,6 nm and from 11.6 to 17.2 nm, respectively (lines 1 and 2 in Figure 7). The measurement of the Casimir pressure results in the strongest constraints in the range from λ = 17.2 nm to λ = 40 nm (see line 3 in Figure 7). Finally, the second Casimir-less experiment leads to the most strong constraints in the wide range from λ = 40 nm to λ = 8 µm (line 4 in Figure 7).
For larger λ the strongest constraints on the Yukawa-type corrections to Newtonian gravitational law are currently obtained from the gravitational experiments. Specifically, within the interaction range from λ = 8 µm to λ = 1 cm the strongest constraints follow from the Cavendish-type experiments [18][19][20]. Thus, in the range from 8 to 9 µm they were obtained in [18], in the range from 9 µm to 4 mm in [19], and for larger λ up to 1 cm in [20]. The initial part of the line representing these constraints in the region from 8 to 34 µm is labeled "gr" in Figure 7. Recently the constraints of [19] were strengthened by up to a factor of 3 in the region of λ from 40 to 350 µm by means of the improved Cavendish-type experiment [80]. Now we deal with the constraints on an axion which follow [76] from measuring the Casimir force in nanometer separation region [77]. An expression for the additional force acting between an Au-coated sphere and a SiC plate due to two-axion exchange is obtained from (13) and (14) [76] where the function χ is defined as and n 1 , n 2 are the numbers of nucleons per unit volume of Au and SiC [4]. The constraints on the coupling constant of axions to nucleons were obtained from the inequality which is similar to (20) used to constrain the Yukawa-type interactions from the results of the same experiment. This constraints are shown by line "new" in Figure 8. In this figure the solid lines 1, 2, and 3 following from the second Casimir-less experiment, measurements of the Casimir pressure, and the lateral Casimir force between corrugated surfaces, respectively, are reproduced from Figure 6. As is seen in Figure 8, the new constraints are stronger than those found from measuring the lateral Casimir force for axion masses m a > 17.8 eV.
We are coming now to a comparison of the above constraints on an axion-to-nucleon coupling constant obtained from the Casimir effect with those following from some other laboratory experiments. Strong constraints on the coupling constant of axions to neutrons over the wide range of axion masses from m a = 10 −4 µeV to m a = 1 µeV were obtained from the magnetometer measurements employing spin-polarized K and 3 He atoms [35]. These constraints in the range of m a from 10 −2 µeV to 1 µeV are shown by line labeled "m" in Figure 8. We note that the constraints found by means of a magnetometer are derived from the one-axion exchange between two neutrons which is described by the effective potential (11). Because of this they are equally applicable to both the QCD and GUT axions. Within the wide range of m a from 1 µeV to 1.7 meV the strongest constraints on the axion-to-nucleon coupling constant valid for the GUT axions follow from the gravitational experiments of Cavendish type [19]. These constraints derived in [36] are shown by line labeled "gr" in Figure 8. For m a > 1.7 meV the strongest constraints on g an are given by line 1 reproduced from Figure 6 which was obtained from the Casimir-less experiment.
In the case of axions, this experiment leads to competitive constraints up to m a = 0.5 eV. For m a > 0.5 eV the constraints following from the Casimir physics (lines 1-3 and "new" in Figure 8) become much weaker than the constraints obtained from experiment on measuring the forces between protons in the beam of molecular hydrogen [38,39]. These constraints are shown by line labeled "H 2 " in Figure 8. For axions with m a > 200 eV the constraints of line "H 2 " were further strengthened by comparing with theory the experimental results on the nuclear magnetic resonance for nucleons in deuterated molecular hydrogen [39].

Proposed Experiments
In this section we consider several experiments suggested in the literature which allow further strengthening the constraints on non-Newtonian gravity and axionlike particles but are not performed yet.
All the experiments on measuring the Casimir interaction considered above have been planned and performed for the purposes which were not directly connected with investigation of hypothetical interactions. The obtained constraints can be considered as some by-product. In [81] the parameters of these experiments were optimized in order the stronger constraints could be obtained on both the Yukawa-type interaction and the interaction due to two-axion exchange. It was shown [81] that the prospective constraints are up to an order of magnitude stronger than those ones already obtained and outlined above.
All precise measurements of the Casimir interaction performed during the last few years employed the sphere-plate geometry and the separation distances below a micrometer. For this reason much attention was attracted to the proposed CANNEX (Casimir And Non-Newtonian force EXperiment) test designed to measure the Casimir force between two parallel plates spaced at separations up to a few micrometers [82,83]. In [84] the parameters of this experiment were optimized for obtaining the strongest constraints on both the non-Newtonian gravity and axionlike particles. Here, we illustrate the potentialities of the proposed setup.
In Figure 9, the constraints on the Yukawa interaction constant α, which could be obtained from the CANNEX test, are shown by the line labeled "proposed". For comparison purposes, the lines 4-7 and "gr" following from experiments discussed in Sections 2 and 4 are reproduced from Figure 7. As is seen in Figure 9, the proposed line provides an order of magnitude stronger constraints on α than those given by line 4 following from the Casimir-less experiment. It is seen also that the proposed constraints encroach on the province where the Cavendish-type experiments presently lead to the strongest constraints shown by the line labeled "gr".
A similar state of affairs holds for the additional interaction between two parallel plates due to two-axion exchange. The constraints on g an , which could be obtained from the CANNEX test, are shown by the line labeled "proposed" in Figure 10. For comparison purposes, the lines 1, 2, and labeled "gr" are reproduced from Figure 8. From Figure 10 it is seen that the proposed line suggests stronger constraints than those indicated by line 1 following from the Casimir-less experiment. Some additional potentialities are also suggested by the CANNEX setup in out of thermal equilibrium conditions [85].
There are several more proposed laboratory experiments aiming to strengthen the constraints on non-Newtonian gravity and interaction of axionlike particles to nucleons in the ranges of parameters considered above. Thus, it was suggested to measure the Casimir-Polder force between a Rb atom and  a movable Si plate separated with a metallic shield in the form of an Au film [86]. It is expected that the experiment of this kind will help to obtain the stronger constraints on non-Newtonian gravity in the interaction range of λ ∼ 1 µm [86]. The respective strengthening of the constraints on g an is also expected.
According to [87], the weakly bound molecules can be used as sensitive sensors for the non-Newtonian gravity. Calculations show that the precision spectroscopy of these systems will allow to strengthen constraints on the interaction of Yukawa type obtained from the experiments on neutron scattering (lines labeled "n" in Figure 7) at least by the two orders of magnitude. Note that a few constraints obtained from atomic spectroscopy in the range of extremely small λ are discussed in [88].
One more proposed detection of non-Newtonian gravity at short range exploits the levitated sensors [89]. It is based on the effect of normal mode splitting in the optomechanical cavity and aims the detection of large extra dimensions proposed in [10,11] and the resulting Yukawa-type corrections to Newtonian gravity [12,13].
In addition to the Yukawa-type interaction and forces arising due to two-axion exchange between two nucleons, some other types of the fifth force are predicted which are suppressed in the regions with high density of matter. These forces are connected with such hypothetical scalar particles of variable mass as chameleons and symmetrons. It was already pointed out [81][82][83]90] that experiments on measuring the Casimir force at large separations can be used for a search of chameleons as possible candidates for the constituents of dark energy. Recently the Casimir experiment was proposed [91] which should be capable to place stronger constraints on symmetrons.

Discussion
In this article, we have reviewed the present state of the art in constraining the Yukawa-type corrections to Newton's law of gravitation and the axion-to-nucleon interaction from laboratory experiments. In was shown that measurements of the Casimir force and Casimir-less experiments lead to strong constraints on the Yukawa interaction constant α and on the coupling constant of axions to nucleons g an over the wide ranges of λ and m a , respectively. The very recent results were presented on constraining both α and g an from the experiment on measuring the Casimir force between an Au-coated sphere and a SiC plate performed at short separations down to 10 nm. These constraints are much stronger that those found previously from measuring the lateral Casimir force between corrugated surfaces.
It should be noted, however, that within the interaction range of λ below 10 nm and above 8 µm the strongest constraints on the Yukawa-type interaction follow not from the Casimir experiments, but from the laboratory experiments of some other kind. Thus, for λ < 10 nm the most strong constraints on α were obtained from experiments on neutron scattering whereas for λ > 8 µm from the Cavendish-type experiments. In a similar way, for axion masses m a below 1.7 meV the strongest laboratory constraints on the coupling constant of axions to nucleons follow not from the Casimir physics, but from the magnetometer measurements and experiments of Cavendish type. At m a > 0.5 eV the strongest constraints on g an were obtained from measuring the forces between protons in the beam of molecular hydrogen.
Thus, the future trends in this field of research are connected with different experimental approaches to constraining the Yukawa-type forces and the interaction of axions with nucleons.

Conclusions
To conclude, recent advances in experimental technique made possible performing more precise laboratory experiments in different branches of physics. The results of these experiments and the measure of their agreement with theory were used for obtaining stringent laboratory constraints on the strength of Yukawa-type interactions and the coupling constant of axions to nucleons. In this article, we summarized the strongest constraints of this kind following from the diversified experiments each of which was used for obtaining the best result in some specific range of parameters. One may hope that by improving the characteristics of these experiments and by implementing some of the recent proposals discussed above it will be possible to obtain even stronger constraints on the non-Newtonian gravity, axionlike particles and on some other related objects, such as chameleons and symmetrons.