1. Introduction
Resonant magneto-optics and the related field of atomic magnetometry have a history started in the late Fifties of the past century with research of Dehmelt [
1] and Bell and Bloom [
2,
3] prosecuting in the Sixties, particularly with the work Cohen-Tannoudji [
4,
5]. The following decades brought important progresses in the comprehension of optical pumping phenomena [
6,
7] and prepared the conditions for an important revival at the beginning of the current century, which was inspired by numerous attracting applications and facilitated by several elapsed technological progresses.
Most of the application fields envisaged for atomic magnetometers are identified on the basis of their excellent sensitivity, and, in the case of RF magnetometers, on their response to high-frequency fields. Among these application fields, emerge the detection of bio-magnetism, e.g., in the construction of magneto-encephalographs [
8,
9], magneto-cardiographs [
10,
11,
12], magneto-miographs [
13]. Another promising area is represented by nuclear magnetic resonance in ultra-low (e.g., at microtesla level) or even vanishing fields, where atomic magnetometers find use as non-inductive detectors [
14,
15,
16,
17,
18] also for imaging experiments [
19,
20,
21,
22]. RF magnetometers find other interesting applications in the detection and imaging of eddy currents [
23,
24], with implications in non-destructive test of materials [
25], security [
26], and biomedicine [
27]. Further perspectives originate from the possibility of producing miniaturized sensors [
28,
29,
30,
31] and arranging them in arrays [
32].
There exist a variety of atomic magnetometers, nowadays also produced as commercial devices. However, they all share several common features [
33]. Generally, the working principle consists in using (near) resonant light to orient atoms in a long lived state (or to induce alignment or higher order momenta); making atoms evolve in the magnetic field under measurement, which thus imprints its features in the subsequent state; and eventually interrogating the atomic sample by means of probe radiation that characterises its optical (absorptive/dispersive) behaviour.
The above mentioned three-step (pumping-evolution-interrogation) procedure [
34] can be performed in a sequence or simultaneously, numerous interaction geometries (e.g., relative orientation of field, pump beam and probe beam) can be considered, tailored time-dependent fields can be applied, various approaches (absorption, dispersion, polarimetry) can be used at the interrogation stage etc.: a
pletora of configurations can be proposed and this large variety comes with wide ranges of sensitivity, bandwidth, dynamical range, time and space resolution etc.
An interesting implementation consists in pumping atoms synchronously with their precession around a (nearly) static field. In this case the pumping radiation must be modulated (in terms of amplitude, polarisation or wavelength) and the probe radiation detects a time dependent optical behaviour. Such light-modulated setup was firstly proposed by Bell and Bloom [
2], who used amplitude-modulated light from a Caesium discharge lamp.
The ideal (optimal) configuration of a Bell-and-Bloom setup considers a static magnetic field oriented transversely with respect to the pump-radiation wave vector. This causes the resonantly induced macroscopic magnetisation to precess in a plane perpendicular to . This precession can be profitably analysed by probe radiation propagating in that plane. The Faraday rotation effect offers a favourable detection scheme based on polarimetric measurement: the probe radiation is linearly polarised and the polarisation plane is rotated by an angle that oscillates at the precession frequency. The amplitude of this oscillation is maximised when the pumping radiation is modulated at the Larmor frequency set by the static field, i.e., at an atomic magnetic resonance.
If the pump modulation frequency is fixed, small and slow variations of the static field bring the system in a near-resonant condition, with the typical effects of signal reduction and dephasing, the latter being more effectively and easily detected, thanks to its (nearly) linear dependence on the detuning.
This classical picture of Bell-and-Bloom behaviour applies when the field variation is small compared with the magnetic resonance line-width and slow with respect to the relaxation time of the atomic orientation. Most of literature reporting Bell-and-Bloom magnetometers interprets the detected signals in terms of such quasi-static and near-resonant condition [
35], furthermore considering optimal field orientation. In other terms, the field is usually assumed to be (nearly) perpendicular to the pump and probe propagation, (nearly) time-independent, and (nearly) matching the Larmor resonant condition.
This paper aims to provide a more general analysis of the behaviour of a Bell-and-Bloom magnetometer when it operates in the presence of a generic time dependent field that is parallel to the static (bias) one. The availability of a more accurate model is of relevance when the output of a magnetometer is used as an error-signal to implement feedback based field-stabilization systems [
36,
37,
38,
39,
40], as well as when designing feed-forward ones [
41] or when evaluating the effects of strong stray fields. The later is of particular relevance for systems operated in unshielded environment [
9,
42].
This paper is organised as follows. In
Section 2 we describe the mathematical approaches used to solve the Larmor equation under the generic conditions mentioned above. In
Section 3 we report the numerical findings of the model, while the main outcomes are discussed in
Section 4 and some conclusions are drawn in
Section 5.
2. Model
The starting point in modeling of the spin response is the Larmor equation for the magnetization vector
:
Here
represents a relaxation mechanism (assumed isotropic),
the action of the forcing term obtained through optical pumping (see [
43] for details),
is the gyromagnetic factor and
is the applied magnetic field.
In this paper we are interested in a configuration where the pumping light and the magnetic field are orthogonal. We take the
x axis along the pump (so
) and the
y axis along the magnetic field, which, in turn, is composed of a static
and a time-dependent
part. Introducing the frequencies
and
the relevant Larmor equation becomes
where
, and we are interested in
, which is the quantity related to the Faraday rotation of the probe laser collinear with the pump as, for example, in [
43,
44].
The solution of (
2) can be written as
where we introduced the Larmor angle associated with the time-dependent field
We assume that both the forcing and the time-dependent field are real-valued periodic functions:
and we can put
(otherwise redefine
) obtaining that also
is a periodic function which permits to write
where the
coefficients are complicated functions of
. It is worth noticing that, in case of pure sinusoidal time-dependent field, the
coefficients are related to the Bessel functions of first kind.
Now the steady-state (SS) solution, valid for
, of (
3) becomes
which is still a complicated expression that can be simplified noticing that in the usual experimental condition the pumping frequency is (nearly) resonant with the Larmor frequency set by the static field, that is
resulting in
where we introduced the detuning
, implicitly defined the
coefficients and the
complex function.
We would like to point out that the result (
8) has a more general validity: indeed it is valid also when the experimental apparatus track the component at fundamental frequency of the pumping modulation. After this demodulation the quantity monitored is the phase of
, which is still a periodic function
where, in general, the
coefficients can be evaluated only numerically as
A notable exception is the case of low intensity time-dependent field which is discussed in the following.
2.1. Low-Intensity Limit
When the time-dependent field is small with respect to the static one, which in our formalism means that the
coefficients are small quantities, the Larmor angle
is a small quantity too. It follows that
and
,
for
. Substituting in the first line of (
8) one obtains
Considering that the second term in parenthesis is
and that
for small
X, one obtains
and neglecting the first, uninteresting term, the final result is
which is identical to the expression reported in [
45] and can be recast as the output of a linear system
3. Results
The result reported in (
8) is quite general and applies to complicated time-dependence of
which would require to consider many
parameters, as well as to the case of a simple sinusoidal oscillation.
It is worth noticing that in the phase time signal (Equation (
9)) the coefficient
gives only a constant contribution which amounts to an experimental offset and can be safely neglected, resulting in an experimental advantage because one does not need to take into account the precise form of the pumping signal. This is not the case if one chooses to monitor the modulus of
where the
coefficient enters in a multiplicative way.
Because of the complicated dependence of the
coefficients on the
, it is difficult to draw conclusions valid for arbitrary values of the parameters. In such sense it is a nice result the analytical low-intensity limit of Equation (
12), which allows a discussion in term of linear time-invariant systems. However, in this paper we explore the complementary regime presenting results obtained as explained in the previous section.
To be more concrete we fix the time-dependence as a simple cosine, namely,
so that
are the Bessel functions of first kind and the coefficients
can be written as
and we expect that for not-so-small
also the higher-harmonics coefficients
will become important. In fact, the system is a parametric one and it will show a kind of “non-linear” behaviour. In other terms, even for a simple sinusoidal time-dependent field as “input”, higher order harmonics will appear in the monitored phase.
In the following pictures we present numerical results for the first Fourier coefficients of the phase for different values of the amplitude of the time-dependent field.
In the
Figure 1,
Figure 2,
Figure 3 and
Figure 4 the modulus of
(from left to right) is reported as a function of
for fixed values of
, respectively. Each curve represents a different value of the detuning
in units of
.
As can be seen in
Figure 1 where
the perturbative limit of Equation (
13) applies. In fact the
and
are two and three orders of magnitude smaller than
respectively. However, as can be seen in the
Figure 2,
Figure 3 and
Figure 4 and
become more and more comparable in magnitude with
increasing the amplitude
of the time-dependent field.
Let us consider the behaviour of
for different values of
. For small
the curve can be interpreted as the response of a low-pass filter, while increasing
it becomes resembling a band-pass filter centred in
. This point of view is reinforced in the low
regime by the result of Equation (
13). However, as can be seen in the pictures, this interpretation can be pushed also for higher values of
even if some structure starts appearing close to
.
The behaviour of is more structured. For is zero (we observed this for all the even coefficients ), then a “low-pass” profile evolves into a double peaked structure with the highest peak centred around and the other around . Increasing gives rise to a small structure close to .
Similar conclusions can be drawn for : here the number of visible peaks is three even for small . It seems reasonable to assume that this is the general behaviour of the coefficients , i.e., increasing increase both the number of “relevant” and for each the number of visible peaks centred at .
4. Discussion
The main scope of this paper is to investigate the dynamic response of a Bell-and-Bloom system to a time dependent field parallel to the bias one. This subject is often presented in terms of an approximation whose validity requires that the time dependent field is small, varies slowly with respect to relaxation time and that the pumping radiation is modulated in resonance with the precession frequency set by the dominant static field. Under these stringent but frequently fulfilled conditions, the systems permits a simplified treatment and behaves as a linear forced-damped oscillator with a low-pass (Butterworth filter) response [
8,
46].
We have previously analysed the behaviour of such a system , using a perturbative approach in the amplitude of the time dependent field [
45], in a more general configuration, as a generic orientation of the field was considered. Before discussing the results presented in Sec.
Section 3 it is worth recalling the main outcomes obtained with that approach.
4.1. Perturbative Results
An interesting feature, resulting at the first order perturbative approximation, emerges when the pump radiation is modulated at a frequency detuned from the exact resonance. Under this condition, the simplified low-pass response becomes inadequate and a band-pass behaviour is observed, with an enhanced response to AC components at frequencies close to the pump-frequency detuning
, as also observed in Ref. [
40].
In application, this feature may help tailor the system response to the detection of oscillating terms at frequencies around . For instance, a magnetometer with a resonance as narrow as 20Hz maybe used to detect with a good efficiency an AC field oscillating at 200Hz, provided that the pump light is modulated at a frequency 200Hz away from .
An interesting behaviour is obtained at an intermediate detuning, . This condition produces a nearly flat amplitude response up to a cut-off frequency set by . Such extended flat bandwidth is obtained at expenses of a slight reduction of the signal amplitude, if compared to the case. Concerning the phase of the magnetometric signal, a maximally extended constant-phase (non dispersive) response is instead achieved under the condition .
In this first-order approximation, the precessing spins do not respond to field variations along directions perpendicular to the static field, while such a response appears in the next order. Indeed, the perturbative treatment developed in [
45] shows that the second-order approximation predicts terms that
double and
mix the frequencies at which time-dependent fields are applied. It is worth noting that the frequency mixing among diverse components of the time-dependent field comes with several peculiar features and it does not occurs among all the couples of components, as detailed [
45].
Further increase of the time-dependent field makes also the second-order approximation inadequate and the numerical analysis developed in
Section 2 becomes necessary. For the sake of simplicity, in this work the time-dependent field is considered at arbitrarily high amplitude, but only along the bias field direction. Moreover the presented results are obtained with a single-frequency term, so that no frequency mixing occur and the non-linearities are pointed out in terms of a harmonic analysis.
4.2. Non-Perturbative Results
The numerical analysis performed at rather weak intensity is perfectly consistent with the first-order perturbative analysis recalled in
Section 4.1. In particular the leftmost plot in
Figure 1 shows that at resonance (
) the system has a low-pass response, which evolves into an extended-flat response when
and becomes a band-pass one for
. At large detunings (band-pass regime) the peak response is half the maximum observed at low frequency and vanishing detunings. A harmonic analysis shows that in the low-field regime (
) the anharmonicity of the response is negligible. A close inspection to the vertical scales in
Figure 1 points out that the second- and third-harmonic terms (
) are depressed by two and three order of magnitudes with respect to the fundamental one, respectively.
Figure 2 and
Figure 3 show the effects of a moderate increase in the amplitude of the time-dependent field. Under conditions in which
, the low-pass and band-pass behaviour does not change appreciably, while the anharmonicity becomes progressively more evident, e.g., for
the third harmonic peaks observed at large detunings are only one order of magnitude weaker than the fundamental one.
The results obtained at small and moderate amplitude of the time dependent field (
Figure 1,
Figure 2, and
Figure 3) show pretty similar features in the spectral response, and the main difference concerns only the anharmonicity level. The response observed for large
is characterised by a peak at
for all the harmonic terms. However, the second harmonic term (
) has also a secondary peak at
and an additional peak at
emerges for the third harmonic term,
.
The
Figure 4 shows new spectral features in the responses. Strong time-dependent fields cause extra peaks emerge in the plots. In particular, the leftmost plot (
) of
Figure 4 shows the appearance of an extra peak at
.