The microscopic Ampère formulation for the electromagnetic force density in linear dielectrics

We present a detailed derivation of the electromagnetic force density and pressure in linear dielectric media according to the so-called microscopic Ampère formulation, which considers the classical dipolar sources in matter along with the hidden momentum contribution. It is seen that, among the other formulations existing in the literature, our proposal is the only one universally compatible with the experimental works reported to date. A new radiation pressure equation for non-magnetic dielectrics under oblique illumination from p-polarized beams is also derived.

The force density on charged matter is unambiguously given by the Lorentz force law, where ρ is the electric charge density, J is the electric current density, E is the electric field and B is the magnetic induction field.The different force densities in each formulation existing in the literature can then be attributed to different modeling of the electromagnetic sources ρ and J.These sources are, of course, not arbitrary -they must be consistent with Maxwell's equations and with the conservation of electric charge.Historically, besides the two formulations that name the Abraham-Minkowski dilemma, the specialized literature discusses three additional main formulations for the electromagnetic force density generated in dielectrics due to the application of external fields, namely, the formulations from Einstein and Laub [24], Ampère and Chu [25].It can be seen that all of them present their own problems [5].First, Minkowski's formulation totally neglects the bound charges inside the material [7].In fact, had he used the correct charge density, he would have obtained Ampère's conventional formulation without the hidden momentum contribution (see Sec. II).As Abraham's formulation only differs from Minkowski's in the momentum density, it should share this same problem.Ampère's formulation does adopt the appropriate classical microscopic model, i.e., the electric sources are given as ideal dipoles and the magnetic sources as tiny current loops [26].This formulation works properly if we are interested only in the movement of the center of mass of rigid bodies.However, as only the macroscopic effects of polarizations and magnetizations are considered, Ampère's formulation is not expected to correctly describe the microscopic force distribution inside materials [27].The Einstein-Laub formulation, on its turn, has been built considering the usual dipolar approximation for the electromagnetic fields inside matter, but adopts an incorrect model for the microscopic magnetization mechanism, as it assumes the existence of magnetic monopoles, which have never been observed and showed disagreement with experiments [26,28,29].For nonmagnetic materials, its related force density should be correct -however, there would still remain strong theoretical issues, like the absence of Lorentz invariance of the electromagnetic stress-energy tensor [30].Although attributing different contributions to light and matter, Chu's formulation follows the same microscopic model of Einstein and Laub.At last, every formulation except Einstein-Laub's and Chu's do not account naturally for electro-and magnetostriction effects, which are quadratic effects on the fields tending to compress the medium towards regions of higher field intensity [31].A table summarizing the different force densities from these formulations can be found in Appendix A. For a more detailed discussion on the different electromagnetic formulations, see for example Refs.[3,5,6,32].
In the context just presented, it seems natural to consider a formulation that arises from the classical charge and current distributions related to microscopic electric and magnetic dipoles.This has been attempted earlier in Refs.[33,34], but in the absence of the hidden momentum contribution (discussed in detail in next section).An axiomatic approach was also reported [35], where it is argued that no specific model for the electromagnetic sources in matter has to be adopted for a consistent de-scription of electromagnetism.In this work, we employ the classical dipolar approximation for electromagnetic sources in matter to derive an electromagnetic force density that will be shown to present all the characteristics necessary to explain the existing experiments.In doing so, the enigma of the hidden momentum is also clarified and properly added to the force density.
This work is organized as follows.In Sec.II we discuss the essential concept of hidden momentum in the context of the Abraham-Minkowski problem.In Sec.III we describe the electromagnetic sources associated with ideal dipoles for both static and time-dependent cases.In Sec.IV we derive the electromagnetic force density related to these sources, as well as the radiation pressure.In Sec.V we compare the obtained results with many experimental investigations reported so far, followed by the discussions in Sec.VI.Lastly, in Sec.VII we summarize our results.

II. HIDDEN MOMENTUM
In 1967, Shockley and James [36] identified a previously unrecognized source of linear momentum that should arise when a magnetic dipole moment m interacts with an electric field E -even if both did not vary in time.This source became known as hidden momentum, and is given by where c is the speed of light in vacuum.Since then, many works tried to properly interpret this puzzling term, which is inevitably tied to the Abraham-Minkowski controversy.Some authors claimed that it occurs as a relativistic effect in systems that are macroscopically at rest, but contains internally moving parts, such as a common electric circuit [37][38][39][40] (see also Ref. [41] and references therein).It was also shown that hidden momentum is necessary to keep the correct relativistic properties of energy, momentum and rest mass of a charge and current carrying body [42].Indeed, it was recently suggested that hidden momentum is a general relativistic concept, not exclusive to electromagnetic systems [43].
Apart from these historical conceptual issues, it was formally shown [44] that starting from the conventional Quantum Electrodynamics (QED) Lagrangian for a point, spinless charged particle in relativistic motion and properly applying the center of mass-energy theorem there must be an extra momentum given by Eq. (2).Classically, the correct interpretation of hidden momentum is actually quite simple: a moving electric dipole develops a magnetic dipole.More specifically, this occurs when the electric dipole moment p of a particle moving with velocity v, both measured in the laboratory frame, is Lorentz-transformed to the particle's rest frame [45].To first order in |v|/c, the new electric dipole moment is [46] where p ′ and m ′ are the particle's electric and magnetic dipole moments, respectively, in its rest frame.The hidden momentum contribution then comes exactly from this last term.
In 1984 Aharonov and Casher showed in a seminal work [47] that the hidden momentum arises as a topological quantum effect when describing the interaction between a charged particle and a magnetic moment, where they obtained Eq. ( 2) as a non-relativistic limit of the Dirac equation -in this context, the hidden momentum is also known in the literature as Aharonov-Casher interaction.In analogy with the Aharonov-Bohm effect [48], this interaction does not necessarily generate a force, but introduces a phase shift in the wave function of the system, which has already been observed -see Refs.[49,50] for example.
Notice that the symmetry inherent to Maxwell's equations requires the existence of an effect analogous to hidden momentum for magnetic dipoles, i.e., an effect due to moving magnetic dipoles generating electric dipoles.This indeed takes place and is known in the literature as Röntgen interaction.Its momentum is given by p R = −p × B and can also be rigorously obtained from the QED framework [44,51,52].This interaction can also generate a topological phase [53], but, to our knowledge, such effect has not been observed yet.
Although hidden momentum has certainly been subject of more intense discussions in the literature, both interactions presented here are of equal importance, and they are expected to take essential part in the eventual resolution of the Abraham-Minkowski problem.They are known to arise when the center of mass-energy of the system is regarded as a dynamic variable.However, the Röntgen term appears naturally even in non-relativistic derivations (see Refs. [52,54] for example), while the hidden momentum necessarily requires a relativistic treatment, as shown in Ref. [44].The Röntgen interaction and hidden momentum contribute to the electromagnetic force density as, respectively, f R = d(P × B)/ dt and f h = − d(M × E/c 2 )/ dt, where P and M are the medium's polarization and magnetization, respectively, and the minus signs added to both equations stem from the fact that the force densities are generated due to the fields losing their momentum.The latter contribution to the force density is important even in systems with non-relativistic velocities [42], and so it must be properly added ad hoc in the results from non-relativistic derivations.

III. DIPOLAR SOURCES
In this section we present the microscopic electromagnetic sources associated to a classical point electric and magnetic dipole.The dipole is composed of two point charges of opposite value ±q separated by a distance |d|.The chosen inertial reference frame is the rest frame of the dipole -consequently, the velocity of the dipole's cen-ter of mass is taken as zero.The dipole moments will first be considered as static and later to be time-dependent.

A. Static dipole
The charge density associated to the classical dipole just described at a point r is given by [27,55,56] where p = qd is the electric dipole moment and δ 3 (r) is the three-dimensional Dirac delta function, i.e., δ 3 (r) = δ(x)δ(y)δ(z).On its turn, the current density in this model is where m = (1/2) r × J d 3 r is the magnetic dipole moment.
Notice that a formal calculation would require the consideration of dynamic dipoles (i.e., both p and m generally time-dependent) from the start because their acceleration can generate radiation-related terms.Besides, the concept of retarded time must also be included to assure physical causality.These requirements are addressed in the next section.

B. Time-dependent dipole
It is well-known that in the Lorenz gauge the electromagnetic potentials φ and A are described by nonhomogeneous wave equations whose formal solutions are [56] and where t r = t − |r − r ′ |/c is the retarded time.Under the dipolar approximation, the electromagnetic potentials of dynamic point dipoles are given by [57] and where p(t) and m(t) are the now time-dependent electric and magnetic dipole moments, respectively, and t 0 = t − r/c is the retarded time at the origin.The electromagnetic fields are given in terms of the potentials as E = −∇φ − ∂ t A and B = ∇ × A. We can use Gauss' Law and Ampère-Maxwell's law to obtain the time-dependent charge and current densities as ρ(r, t) = −ε 0 ∇ 2 φ−ε 0 ∇•∂ t A and J(r, t) = µ −1 0 ∇×∇×A.This is, however, not a very convenient procedure because the implicit dependence of t 0 in r makes the calculation very difficult.Alternatively, it is sufficient to show that the sources adopted in the last section, when promoted to time-dependent, generate the correct electromagnetic potentials given in Eqs.(7) and (8) when calculated using Eqs.( 5) and (6).We start by calculating the electric potential due to ρ(r, t) = −(p(t)•∇)δ 3 (r) as Using the property Due to the implicit dependence of t r on r, we have , where ∇′ denotes the nabla operator acting only on the spatial coordinates.Specifically, in the point dipole approximation we have ∇′ • p(t 0 ) = 0 as p does not depend on r, so that the electric potential in Eq. ( 10) is reduced to Eq. ( 7), as expected.
To calculate the vector magnetic potential A, we first notice that when p is time-dependent there is an extra term ṗ(t)δ 3 (r) in J originating from the continuity equation, ∂ t ρ = −∇ • J, where the over-dot denotes time derivative.Therefore, J(r, t) = ṗ(t)δ 3 (r) − (m(t) × ∇)δ 3 (r), and the vector magnetic potential is (11) The first term is trivially integrated to (µ 0 /4πr) ṗ(t 0 ).The second term is analogous to Eq. (10), with p → m and divergence operator → curl operator, resulting in Eq. ( 8) and completing our calculation.Thus, we have shown that if the dipolar approximation can be suitably applied, the dynamical microscopic electromagnetic sources can be described as and Consequently, by employing the Lorentz force density to these sources we should be able to obtain the appropriate force density distribution inside matter within the adopted approximations.Lastly, we notice that the retarded time t r is known to be related to general scattering phenomena inside matter according to the Ewald-Oseen extinction theorem [58,59].As we adopt here the common effective, continuum description of dielectrics through the use of the macroscopic parameters ε and µ, the effects of the retarded time are already implicitly contained in the resultant fields.Thus, in our model the fields and forces can be evaluated at the regular time t.

IV. ELECTROMAGNETIC FORCE DENSITY AND PRESSURE
We consider electromagnetic fields inside dielectric materials within the optical bandwidth.To calculate the electromagnetic force density and pressure, we suppose that at this optical length scale the microscopic sources are well described by the dipolar approximation, as discussed in the last section.

A. Force density
The force acting on charged matter is unambiguously given by the continuous version of the Lorentz force law, which after some algebra (see calculations in Appendix C) leads to the force density in the Microscopic Ampère (MA) formulation, Here, by "force density" we mean the electromagnetic force acting on a small volume of the dielectric, δV , which is microscopically large (i.e., encompasses a large number of dipoles), but is still much smaller than the dielectric's macroscopic volume.Also, ε r = ε/ε 0 is the relative permittivity and µ r = µ/µ 0 the relative permeability, while H = µ −1 B is the magnetic field.The last term in Eq. ( 15) is the Röntgen interaction, which naturally appeared in our non-relativistic derivation for a dipole at rest, as anticipated in Sec.II.If we employed a relativistic derivation for a moving dipole from the very beginning, there would be an extra contribution -the hidden momentum -as shown in Ref. [44].An alternative non-relativistic derivation for the force density (in non-magnetic media) where the dipole is moving can be found in Ref. [54].
As discussed in Sec.II, when switching to the laboratory frame it is necessary to add the hidden momentum contribution to Eq. ( 15) as f h ≈ ∂ t (E × M)/c 2 .The time derivative approximation d/ dt ≈ ∂/∂t takes place because the dipole's velocity (as measured in the laboratory frame) is much smaller than c.The force density is then where n = √ ε r µ r is the refractive index and it is assumed that ε and µ do not depend on time.A simple rearrangement of the gradients as products yields, at last, the main result of our work, This force density equation is given in the laboratory frame.It is valid for linear, isotropic inhomogeneous media, with ε and µ independent of time, i.e., no dispersion.The presence of free sources would generate the extra terms ρ f E and J f × B in Eq. ( 14), and can be included if necessary.The first and second terms can be assigned as the electrostriction and magnetostriction force densities, respectively.The third and fourth terms are the usual Abraham-Minkowski force, which occur in non-homogeneous regions, and the last term is the wellknown Abraham force.This equation contemplates almost every aspect of the reported experiments (as will be discussed in Sec.V), and arises naturally from a clear and simple microscopic model, with no need for phenomenological approaches.An alternative derivation of Eq. ( 17) for homogeneous media using the Lagrangian approach is found in Appendix B.

B. Radiation pressure
The electromagnetic pressure, widely known in the literature as radiation pressure, is here denoted by P rad and calculated at a flat dielectric interface assuming a laser beam with azimuthal symmetry about the propagation axis.The deformations induced by radiation pressure in dielectric liquids are typically bulges of height of order 10 nm [60][61][62][63][64], rapidly decreasing over a length of about one beam waist, which is usually of order 100 µm.Thus, considering the interface flat even when the fields are acting on it is certainly a good approximation.
As an isotropic quantity, the radiation pressure has two contributions in this case: one due to the discontinuity of refractive index, P and M in the direction normal to the interface and another one due to the difference in the radial forces (electro-and magnetostriction effects) in each medium.The first contribution can be obtained from Eq. ( 17) by properly integrating the normal component of the force density.For example, for a monochromatic beam propagating in z direction incident from a nonmagnetic, linear and isotropic dielectric medium into a flat interface at z 0 = 0, as illustrated in Fig. 1 (a), we have We assumed in Eq. ( 18) there are no free charges at the dielectric interface.Here, |E| 2 avg is the average of the squared electrical field magnitude across the interface and the last term of Eq. ( 17) has been averaged out.
The second term of the radiation pressure at the interface is related to the radial direction, and for normal incidence is given by where the same assumptions were used again and the function inside the brackets was implicitly calculated at fixed r.
Summing the two contributions, we have the radiation pressure at the flat dielectric interface as Applying Maxwell's equations boundary conditions, we have then Here E z,t is the transmitted field component normal to the interface.As this component is not continuous across the interface, we averaged its squared magnitude with a simple arithmetic mean.The tangential components E x and E y are continuous across the interface, and therefore do not need the subscript indicating the current medium.Notice a very important subtlety introduced by this equation: the permittivities ε 1 and ε 2 are related to the direction of z, which in our convention always points from medium 1 to medium 2. The gradient ∇ε is calculated accordingly, resulting in the term outside the brackets in Eq. ( 21).On the other hand, the permittivities ε i and ε t are related to the beam propagation direction, i.e., they refer to incident and transmitted components, respectively.Thus, if the beam is propagating in the z direction, we have ε i = ε 1 and ε t = ε 2 , as in Fig. 1 (a); if the propagation direction is reversed, we have ε t = ε 1 and ε i = ε 2 , as in Fig. 1 (b).By using the conventions just described, Eq. ( 21) is then valid for any beam polarization, incidence angle and propagation direction.Specifically, for normal incidences we have E 2 z ≪ E 2 x +E 2 y , even for typical focused Gaussian beams.This leads to which is the widely known Abraham-Minkowski radiation pressure.In fact, as the electrostriction effects are seen to cancel at the interface, Eq. ( 20) could also be obtained through Abraham's force density, f Ab = −|E| 2 ∇ε/2however, electrostriction plays an important role in the stability of fluids, as will be discussed in Sec.V A. We see that keeping the incidence as normal and reversing beam propagation direction would generate the same radiation pressure equation.For oblique beam incidences, we must properly account the different reflection and transmission coefficients for each polarization.We will consider the beam is locally a plane wave -which is a good approximation, since even for focused Gaussian beams the field components in directions other than the polarization one are typically negligible.In this condition, we can apply Fresnel equations to describe the reflected and transmitted field amplitudes.First, for s polarization, the field is by definition perpendicular to the plane of incidence.In this case, we have in our convention E z = 0, which generates where E 0 is the field amplitude, t s is the transmission coefficient for s polarization and θ i is the incident angle relative to the interface's normal direction.
For p polarized beams, we have a non-zero normal component, so that the radiation pressure becomes where t p and r p are the transmission and reflection coefficients for p polarization, respectively, and θ t = sin −1 ((n 1 /n 2 ) sin θ i ) is the transmitted (refracted) angle.
The terms divided by two within the brackets correspond exactly to the spatial average of E 2 z across the interface.

V. COMPARISON TO EXPERIMENTS
As we have seen in Sec.IV, the MA formulation naturally accounts for the electro-and magnetostriction effects, presents a radiation pressure of the Abraham-Minkowski form and an Abraham-type momentum density.To the best of our knowledge, there is no force density in the literature valid for optical regime that simultaneously presents all these characteristics without including phenomenological approaches.These fundamental properties will be used in this section to analyze the main existing experiments related to electromagnetic force density.We emphasize that the great majority of these experiments were quantitatively described in terms of one of the previously known electromagnetic formulations, either in the original or in subsequent works.Our objective here is to show if and how MA formulation can also be applied to interpret them.For better organization, these experiments are grouped in four categories: radiation pressure experiments, photon momentum experiments, bulk force experiments and total force experiments.

A. Radiation pressure experiments
The surface deformation of water under normal laser incidence was successfully explained using the radiation pressure given in Eq. ( 21), both in old and recent measurements [60,65,66] -specifically, in Ref. [65] the reversed beam propagation direction was also considered.An interface of different fluids close to the critical point was studied in Ref. [67] and the observed surface deformations were also well described by Eq. (21).
The radiation pressure for oblique incidence adopted in the literature is [68][69][70][71][72] where n i is the incidence medium's refractive index, I is the beam intensity and R and T are the interface's reflectance and transmittance, respectively.This FIG. 2. Radiation pressure for oblique incidence and ppolarized beams according to the literature, Eq. ( 24), and to our suggestion, Eq. ( 23).In (a) we have an water-to-air incidence and the radiation pressure is calculated for angles smaller than the critical angle for total internal reflection, θc ≈ 48.77 • .In (b) we have an air-to-water incidence.The result for s polarization is also shown for completeness.Beam intensity is 1.0 W/m².The relative permittivities used for water and air were 1.769 and 1.0, respectively.
equation contemplates both polarizations in a single equation by using the appropriate R and T , and has been applied to explain the experiments reported in Refs.[63,64,68,69,72,73].However, this equation does not account properly for the discontinuity of the normal field at the interface for p polarization, as shown in Appendix D.
For an air-water interface, Fig. (2) shows the radiation pressures from Eqs. ( 23) and (24) for both incidence directions.The behavior is qualitatively the same, but the magnitude of the corrected version is about 5-10% larger.At first sight, our result seems compatible with the deformations observed in Refs.[63,64,73], but further investigations are necessary.Specifically, we notice that in Refs.[63,64] the ellipsoidal character of the incident beam's cross section at the interface under different incidence angles was apparently not taken into account -a fact that can significantly alter the employed theory.
Notice that when considering a fluid incompressible we are assuming any information about deformations in the fluid propagates instantaneously -thus, no transient response is present.Indeed, at equilibrium, applying the divergence operator to the Navier-Stokes equation for an incompressible fluid at rest yields [74] where f em is the electromagnetic body force and P is the fluid's pressure.This is an elliptic partial differential equation for P, known as Poisson's equation.It is well-known to possess unique solutions (up to an additive constant) for a very broad class of boundary conditions.Thus, in this situation, the Abraham-Minkowski pressure, Eq. ( 21) with E z = 0, arises naturally as a boundary condition uniquely related to the divergence of the body force density from Eq. ( 17).This is in accordance with the fact that striction effects play a significant role in the stability of fluids [68,[75][76][77].On the other hand, if the fluid develops a position-dependent velocity field, the pressure at the surface can not be uniquely related to the body force anymore, as there will be another source term in Eq. ( 25).In fact, the pressure (and consequently the surface deformation) can even change signs, as shown in Ref. [78].This is a possible explanation to the Abraham-type deformation of a free fluid surface reported in Ref. [79].

B. Photon momentum experiments
The Abraham force (last term of Eq. ( 17)) is known to be related to a photon momentum proportional to n −1 [5,7,80], which is usually referred to as Abraham-type momentum in the context of the Abraham-Minkowski controversy, while the Minkowski-type momentum is proportional to n.The recoil of a submerged mirror due to radiation pressure was measured twice [81,82] and the results were directly proportional to the refractive index of the background dielectric media -i.e., of Minkowski's form.This can effectively be explained by the Dopplershifted recoil of the mirror, while the field retains Abraham's form of momentum [7]: the incident photon has a momentum proportional to n −1 , but, due to Doppler's effect, the mirror's recoil is described by a momentum transfer linear in n if the recoil velocity is non-relativistic.Therefore, by considering this interaction of field and matter, the results from Refs.[81,82] are in accordance with MA force density.Alternatively, this can also be described by the mass-polariton quasi-particle model [83], where a light-induced atomic mass density wave carries the difference of the Minkowski and Abraham momenta as discussed in Sec.VI.This explanation also applies to Ref. [84], where the recoil of ultra cold atoms in a Bose-Einstein condensate due to radiation pressure was observed to be compatible with Minkowski momentum.
There is an old measurement of the photon drag effect in semiconductors that agreed with Minkowski's momentum [85], and the reasoning from last paragraph is valid again.This also applies to most cases reported in Ref. [86], where the same effect was measured in thin metal films.However, it is important to notice that one specific measurement in this last reference showed a negative dependence on Minkowski's momentum -a result that, according to the authors, still lacks theoretical explanation regarding the optical transduction and micro-scopic momentum exchange mechanisms.

C. Bulk force experiments
This section deals with the experiments where only the body forces from the striction effects, the first and second terms of Eq. ( 17), are relevant.The electrostriction effect was measured inside a fluid in Refs.[87,88] using high-intensity static fields (the latter one in a microgravity environment) and the results agreed with the so-called Helmholtz force density [3,89,90], given by f H = ε 0 ∇ (ρ m ∂ε r /∂ρ m ) T |E| 2 /2 .This equation is phenomenological and valid for media in thermodynamic equilibrium with external fields, where ρ m denotes mass density and T temperature.We notice that the results from Refs.[87,88] can also be explained by the MA formulation if we consider the local-field correction in the form of the Clausius-Mossotti relation [56,91], where we have ρ m ∂ε r /∂ρ m = (ε r − 1)(ε r + 2)/3.However, it should be emphasized that the equivalence between MA and Helmholtz force densities is limited to static and quasi-static conditions -although the mathematical form would be identical, the two force density equations are built under very different assumptions.In particular, it is certainly not expected that the medium is in thermodynamic equilibrium under optical excitation.Indeed, for optical excitation the local-field correction is typically given as the Lorentz-Lorenz relation [92].Nonetheless, recent measurements of the electrostriction effect in water for laser excitation at optical frequency were very well described by the MA formulation without the correction to the local-field [61,93].This can be justified by the argumentation presented in Ref. [94], where this correction is absent due to the optical electrostriction effect causing energetically non-conservative changes in the dipole moments through the variation of the material's local mass density.
A quite intricate measurement of the electromagnetic force inside optical fibers was reported in Ref. [95].It was concluded that the force density has a different symmetry than the expected from MA formulation for this case.It should be mentioned that the irregular positiondependent refractive index in the optical fiber due to its fabrication process may play a significant role in the force density symmetry through the terms proportional to ∇ε and ∇µ.
In Ref. [96] the observation of the Abraham force was reported in a liquid-filled hollow optical fiber, where the Abraham-Minkowski pressure at the free liquid surface was claimed to be carefully suppressed by the geometry of the waveguide.However, striction forces were not considered in the analysis of the results, where it is expected they would generate deformations with contrary direction to the observed one.Furthermore, there would be an additional Abraham-Minkowski force at the ring core/liquid interface, and adhesion effects are also expected to be important -indeed, these last two forces should partially cancel striction forces.If this cancellation is significant, the remaining force term according to MA formulation would be Abraham's one, in agreement with the authors' conclusions.Alternatively, all the effects but Abraham's force could also be relatively very small.The presence of more than one mode in the excitation wave can also be relevant.

D. Total force experiments
As electro-and magnetostriction are gradient forces, they always produce zero total force when integrated over the material's volume, and so they do not contribute to the macroscopic movement of the body [3,30].Thus, the total force acting on a dielectric body should be composed of only the time-derivative term in Eq. ( 17), i.e., Abraham's force.The movement of a torsion pendulum induced by the simultaneous application of low frequency time-dependent electric fields and static magnetic fields was measured in Ref. [97], and agreed with Abraham's force.The experiment and results reported in Ref. [98] are similar, but more detailed as it also covered the case of electromagnetic forces generated by static electric fields together with time-varying magnetic fields.At last, in Ref. [99], the pressure variation of a confined gas due to the presence of electromagnetic fields was measured to be compatible with Abraham's force.All these results are in agreement with MA formulation.

VI. DISCUSSION
Although, as mentioned earlier, a detailed discussion on the different electromagnetic formulations existing in the literature is out of our scope, we recall that every experiment of the long list cited in the last section can be quantitatively interpreted by using at least one of them -however, we emphasize that only MA, Einstein-Laub (EL) and Chu formulations can be consistent (at least partially) with all of the experiments.Nevertheless, there are strong theoretical issues with EL and Chu formulations which will be discussed now.As their force density is the same (see Appendix A), we can focus only on the former.
The EL force density can be written as [3] Comparing the above equation to our proposed force density, Eq. ( 17), we see their difference lies in the second term, the magnetostriction effect -which is a pure magnetic effect.Such difference is expected because the magnetization mechanism in the EL formulation is due to magnetic monopoles [24], while in MA formulation we adopt the classical Ampèrian current loops model.The magnetic force under the monopole model is [100] F m = (m • ∇)B, while for the current loop model we have Under the point dipole approximation we adopted, the spatial derivatives of m are neglected, so that F l = (m • ∇)B + m × (∇ × B) (exactly as in Eq. (C12)).Therefore, the models yield the same force only if ∇ × B = 0.More importantly, the monopole model has been convincingly observed to disagree with experimental results, while the loop model has shown very good agreement [26,28,29].Experimentally, there are unfortunately very few experiments probing the magnetostriction effect [101][102][103][104][105][106][107], and, to our knowledge, no experiment at all under optical excitation -mainly because thermal effects would in this case be dominant.Additionally, in these works the classical treatment of magnetostriction is phenomenological [108,109], given in terms of magneto-elastic coefficients, which is a significantly different approach from EL and MA descriptions.Thus, currently EL and MA force densities can not be directly discriminated on experimental basis, and do yield the same results for non-magnetic media.Nonetheless, the major concern about the EL formulation is based on the theoretical and experimental results just discussed, which we feel are compelling enough.Another important feature of the EL formulation is that it satisfies the so-called duality transformations of electromagnetic fields [26,110].Again, this is expected since in this formulation the polarization and magnetization mechanisms are treated in direct analogy.This is not the case with MA formulation, where these two mechanisms are distinct -therefore, not satisfying these transformations is actually expected for the MA force density.We notice, however, that a dual-symmetric classical electromagnetic theory can be built without invoking magnetic charges [111].
The force density presented in Eq. ( 17) must, of course, be theoretically compatible with the covariance requirements from special relativity.It was shown that to fulfil this condition there must be a coupled state of field and matter propagating through the material [83,112], which is described by the so-called mass-polariton (MP) quasi-particle.In this theory, the medium contribution -whose importance had been already noticed in earlier works in connection to the photon mass drag effect [7,[113][114][115] -is described as a mass density wave propagating along with the electromagnetic wave due to the atoms being driven forward by the optical force.This theory was shown to be covariant [112,116], to conserve angular momentum (from orbital and spin origins) [117], to be consistent with both classical field and photon descriptions [83,117] and, at last, to be in agreement with the celebrated recognition of Abraham's and Minkowski's momenta as kinetic and canonical momenta, respectively [80,118].The MP dynamics and its momentum transfer employ the Abraham force density, which consists of the three last terms of Eq. ( 17).In the present work, the force density additionally has the striction effects, which correspond to pure stresses and do not affect the overall momentum transfer inside the material [3,94].Therefore, we expect Eq. ( 17) is also compatible with the MP dynamics, which can potentially provide a complete microscopical description of energy and momentum transfer in linear, lossless and non-dispersive dielectrics, but a detailed discussion will be addressed in a future work.As the fundamental small mass transfer associated to the mass-polariton has not yet been observed, the comparison with the experimental investigations from last section are not affected.
Regarding the total forces acting on the medium, it is known that different formulations can consistently provide the same results [119][120][121] -in fact, it has been argued that, with the proper choice of material contribution, every existing formulation leads to the same results and thus choosing from one of them is just a matter of personal convenience [4].The possibility of experimental discrimination between distinct force densities must, therefore, somehow identify their particular spatio-temporal dependence within the material.This is, of course, a challenging experimental work, but it was reported recently in all-optical pump and probe photo-induced lensing experiments [61,93] and the results showed an excellent agreement with MA force density.In fact, these two pioneering works are very valuable as they should change the aforementioned theoretical beliefs that force densities are not unique and that only total forces can be observed.
The new radiation pressure equation for p polarization, Eq. ( 23), was derived for non-magnetic media, where f MA can be written as f MA = (P • ∇)E + ∂ t P × B (Eq. (C13) with M = 0).Although not in the explicit form of Eq. ( 23), the radiation pressure arising from this force density has already been addressed in the literature [120], providing consistent theoretical results regarding total momentum conservation.Additionally, we can, for example, take the thought-experiment considered in Ref. [122], where a p-polarized beam propagating in air enters and then exits a prism's wedges at Brewster angles θ B and θ ′ B , respectively.In this case, r p is zero at both interfaces, yielding (apart from a common multiplicative constant) the pressures at entrance and exit interfaces (respectively) , where θ B = tan −1 (n 2 /n 1 ) and θ B +θ ′ B = 90 • .It can be numerically verified that the sum of these pressures is zero, as required for conservation of total momentum in this situation -in fact, this result still holds regardless of the dielectric medium surrounding the prism.For general incidence angles, we must resort to computational techniques as the internal beam reflections will naturally make the analysis much more complicated.
Though our presented derivation of the electromagnetic force density is currently limited to the simplest type of dielectric materials, we notice that there are many possibilities that can be explored starting from it.One can, for instance, try to extend the theory to more complex materials, where effects such as dispersion, absorption, anisotropy and nonlinearities can take place.Consideration of non-conservative optical forces can also be of interest [123].The analysis of angular momentum distributions inside materials is also relevant, especially because light can have both spin and orbital angular momentum [124] -this has been simulated, for example, in Ref. [125].Additionally, one can search for a physically more fundamental formalism by fully working in the Quantum Mechanics regime, where effects such as field fluctuations and vacuum pressure can occur -there are already some theoretical works in this regard, contemplating, for example, the radiation pressure due to the quantum-mechanical Lorentz force [126], QED corrections to the Abraham force [127] and Casimir-like effects [128].

VII. CONCLUSIONS
In this work we presented a new equation for the electromagnetic force densities inside linear, isotropic, nondispersive and lossless dielectric material, called Microscopic Ampère formulation.This result is derived from the well-established dipolar approximation for electromagnetic sources.Among the other formulations existing in the literature, our proposal is the only one capable of explaining the vast majority of experimental works reported to date and to be simultaneously in accordance with the absence of magnetic monopoles in nature.Additionally, the proposed force density needs the inclusion of the hidden momentum contribution, whose origin was briefly discussed and clarified in the context of the Abraham-Minkowski controversy.A new expression for the radiation pressure in non-magnetic dielectrics under oblique incidences for p-polarized beams was also derived.It is consistent with momentum conservation, but more investigations are necessary -especially experimental ones.
Even though currently limited to classical Physics and to the simplest type of dielectric materials, the electromagnetic force density and pressure presented here cover a lot of practical applications and provide an important step towards obtaining the definitive knowledge of the behaviour of light inside matter.

ACKNOWLEDGMENTS
The research leading to these results received funding from CNPq (304738/2019-0), CAPES (Finance Code 001), Fundação Araucária, and FINEP.M.P. acknowledges funding from the Academy of Finland under Contract No. 349971.
Appendix A: Force densities according to different formulations For a modern derivation of the electromagnetic force densities listed below, we refer to the works [3,5,7,25,31].

Formulation
Force density a Considering the total stress-energy tensor of field plus matter.
b Valid for static or quasi-static fields.The force density from Eq. ( 17) can be derived for homogeneous media (where ∇ε = ∇µ = 0) from the Lagrangian approach as well.The potential energy of induced electric and magnetic dipoles is known to be [56] In the presence of an external electromagnetic force, the dipoles are put into motion, acquiring a velocity v given in the laboratory frame.Due to this motion, to first order in |v|/c the dipole moments are given in the laboratory as p ′ = p+v×m/c 2 and m ′ = m−v×p [46].Therefore, in the laboratory frame the non-relativistic Lagrangian reads where the fields E and B are assumed to be also given in the laboratory frame.The canonical momentum is then from where we identify the Röntgen and hidden momentum contributions (the two last terms, respectively), discussed in Sec.II.The force m v on the dipole is then obtained from the Euler-Lagrange equation, which yields Isolating m v and neglecting the terms linear in v, we obtain which, once divided by the volume and rearranged, leads to Eq. ( 17) for homogeneous media.Notice that a similar Lagrangian approach was reported in Ref. [129] for static fields and permanent dipoles -therefore, the derivation shown here is an extension to induced dipoles and timedependent fields.This approach is here convenient as it elucidates very clearly the origin of the Röntgen and hidden momentum contributions as transformations between inertial reference frames of the moving dipole moments -however, it does not address the possibility of a non-homogeneous medium, which is crucial to explain the radiation pressure experiments, as seen in Sec.V A.

Appendix C: Force density calculation
The following derivation can also be found in the Supplementary Material of Ref. [61].Here, we add it with some more details for completeness and convenience.
We start from the Lorentz force equation with the microscopic sources, namely As stated in Sec.IV, the integral is performed over a small volume of the dielectric, δV , which is microscopically large but still much smaller than the dielectric's macroscopic volume.We must now calculate Eq. (C1), which gives the total force on the dielectric.Each term of the integral will be treated separately.The first term is Writing explicitly the components of the term inside the brackets, we have Integrating each term by parts we obtain The last term of the above equation is related to the spatial extent of the dipole and under the point dipole approximation can be neglected as p does not depend on r -so, The second term of Eq. (C1) is integrated trivially: The last term of Eq. (C1) is where the auxiliary vector C has been implicitly defined for simplicity and î, ĵ, k denote the cartesian unit vectors.
Performing the cross product, we have The x component of F 3 is then where we have integrated by parts again and applied the point dipole approximation.The y and z components are obtained analogously, yielding To obtain a more compact result, notice that and, consequently, which completes our calculation.Summing the three contributions to the force and dividing by the dielectric's volume, we obtain the positionand time-dependent force density in the so called Microscopic Ampère (MA) formulation where the fields are the macroscopic ones evaluated at the location of the dipoles.Eq. (C13) was first given, to our knowledge, in Ref.
[27] -but it was not explored in the context of the Abraham-Minkowski controversy.It is compatible with a pure, point dielectric dipole located at the origin and at rest in its own frame.Explicitly, the time derivative term is as the dipole is assumed to be at rest.Now, we want to rewrite Eq. (C13) in a form which it is most conveniently interpreted in the context of the Abraham-Minkowski controversy.This task will require a lot of vector algebra.We start by using the vector property The radiation pressure for oblique incidence in nonmagnetic dielectrics is given in literature by Eq. ( 24).This equation is equivalent to [71] where α is the angle between the electric field and the plane of incidence.Thus, for α = 0 we have a p polarized beam, while for α = π/2 we have a s polarized beam.In the former case, we have then The use of the last term in the above equation as the average of the squared field component normal to the interface is not physically justified.Therefore, we expect that Eq. ( 23) is the correct one for radiation pressure due to p-polarized beams.For s-polarized beams, Eqs.(D1) and ( 24) yield the same results.

FIG. 1 .
FIG.1.Convention for permittivities used in radiation pressure equation, Eq. (21), according to beam incidence direction.In (a), the beam propagates in z direction, while in (b) the propagation direction is reversed.

TABLE I .
Force density according to different electromagnetic formulations.
[130]2 θ i + cos 2 θ t )T p .(D2)We can rewrite this equation by using the relationT = (n 2 cos θ t /n 1 cos θ i )|t 2 |, which is valid for both polarizations[130], obtaining By applying Snell's law and the relation n21 E z,i = n 2 2 E z,t to the last equation, it is possible to write it explicitly in terms of the field components as y + E z,i E z,t .(D5)