An apparatus is provided for determining a target wavelength λ of a target photon beam. The apparatus includes a photon emitter, a pre-selection polarizer, a prism composed of a faraday medium, a post-selection polarizer, a detector and an analyzer. The photon emitter projects a monochromatic light beam at the target wavelength λ substantially parallel to a magnetic field having strength b. The target wavelength is offset from established wavelength λ′ as λ=λ′+Δλ by wavelength difference of Δλ<<λ. The faraday prism has Verdet value V. After passing through the pre-selection polarizer, the light beam passes through the prism and is incident to an interface surface at incidence angle θ0 to the normal of the surface and exits into a secondary medium as first and second circularly polarized light beams separated by target separation angle δ and having average refraction angle θ. The secondary medium has an index of refraction of n0. After passing the post-selection polarizer, the detector measures target pointer rotation angle Aw based on the target separation angle δ. The analyzer determines the target wavelength λ by calculating offset pointer rotation angle ΔAw=Aw−A′w from calibrated pointer rotation angle A′w based on established separation angle δ′ that corresponds to the established wavelength λ′, and by estimating the wavelength difference based on

Δ λ - 2 ɛπ n 0 Δ A w cos θ VB sin θ 0 ,
in which ε is an amplification factor. A method is provided incorporating operations described for the apparatus.

Patent
   8674688
Priority
Jun 06 2011
Filed
Jul 06 2011
Issued
Mar 18 2014
Expiry
Oct 15 2032
Extension
497 days
Assg.orig
Entity
Large
1
5
currently ok
4. A method for determining a target wavelength λ of a target photon beam, said method comprising:
emitting a monochromatic light beam at the target wavelength λ, substantially parallel to a magnetic field having strength b and offset from established wavelength λ′ as λ=λ′+Δλ by wavelength difference of Δλ<<λ;
pre-selection filtering of said monochromatic light beam as a pre-selection light beam;
refracting said pre-selection light beam through a prism composed of a faraday medium having Verdet value V such that said light beam is incident to an interface surface at incidence angle θ0 to a normal of said surface and exits into a secondary medium as first and second circularly polarized light beams separated by target separation angle δ and having average target refraction angle θ, said secondary medium having index of refraction of n0;
post-selection filtering of said polarized light beams as post-selection light beams;
measuring target pointer rotation angle Aw based on said target separation angle δ of said post-selection light beams; and
determining the target wavelength λ by calculating offset pointer rotation angle Δλ=Aw−A′w from calibrated pointer rotation angle A′w based on established separation angle δ′ that corresponds to said established wavelength λ′, and by estimating said wavelength difference based on
Δλ - 2 ɛπ n 0 Δ A w cos θ VB sin θ 0 ,
in which ε is an amplification factor, and θ′ is average established refraction angle.
1. An apparatus for determining a target wavelength λ of a target photon beam, said apparatus comprising:
a photon emitter for projecting a monochromatic light beam at the target wavelength λ substantially parallel to a magnetic field having strength b and offset from established wavelength λ′ as λ=λ′+Δλ by a wavelength difference of Δλ<<λ;
a pre-selection polarizer through which said light beam passes from said emitter as a pre-selection light beam;
a prism composed of a faraday medium having Verdet value V through which said pre-selection light beam passes from said pre-selection polarizer and is incident to an interface surface at incidence angle θ0 to a normal of said surface and exists into a secondary medium as first and second circularly polarized light beams separated by target separation angle δ and having average target refraction angle θ, said secondary medium having index of refraction of n0;
a post-selection polarizer through which said polarized light beams pass as post-selection light beams;
a detector for receiving said post-selection light beams and measuring target pointer rotation angle Aw based on said target separation angle δ; and
an analyzer for determining the target wavelength λ by calculating offset pointer rotation angle ΔAw=Aw−A′w from calibrated pointer rotation angle A′w based on established separation angle δ′ that corresponds to said established wavelength λ′, and by estimating said wavelength difference based on
Δλ - 2 ɛπ n 0 Δ A w cos θ VB sin θ 0 ,
in which ε is an amplification factor, and θ′ is average established refraction angle.
2. The apparatus according to claim 1, wherein said pointer rotation angle is
A w = ( θ + - θ - ) cos ɛ + [ ( θ + + θ - ) - 2 θ 0 ] sin ɛ 2 sin ɛ
in which θ+ and θ are respectively right- and left-polarized refraction angles with said average target refraction angle such that θ=½(θ+).
3. The apparatus according to claim 1, wherein said analyzer estimates divergence δ≈2εAw for determining the target wavelength.
5. The method according to claim 4, wherein said pointer rotation angle is
A w = ( θ + - θ - ) cos ɛ + [ ( θ + + θ - ) - 2 θ 0 ] sin ɛ 2 sin ɛ
in which θ+ and θ are respectively right- and left-polarized refraction angles with said average target refraction angle such that
θ = 1 2 ( θ + + θ - ) .
6. The method according to claim 4, further including estimating divergence δ≈2εAw to determine the target wavelength.

The invention is a Continuation-in-Part, claims priority to and incorporates by reference in its entirety U.S. patent application Ser. No. 13/134,486 filed Jun. 6, 2011 titled “Magnetic Field Detection Using Faraday Effect” and assigned Navy Case 99670.

The invention described was made in the performance of official duties by one or more employees of the Department of the Navy, and thus, the invention herein may be manufactured, used or licensed by or for the Government of the United States of America for governmental purposes without the payment of any royalties thereon or therefor.

The invention relates generally to detection of small wavelength differences. In particular, this invention relates to a quantum enhanced method for determining small differences in coherent photon wavelengths using weak value amplification of a Faraday Effect.

Detecting the presence of two nearly coincident wavelengths can become necessary for various operations, such as for detection of small Doppler shifts. However, conventional devices and processes exhibit deficiencies in cost and/or portability.

Conventional wavelength discrimination devices yield disadvantages addressed by various exemplary embodiments of the present invention. Various exemplary embodiments provide an apparatus for determining a target wavelength λ of a target photon beam wavelength λ using a prism composed of a Faraday medium having Verdet value V.

The apparatus includes a photon emitter, a pre-selection polarizer, a prism composed of a Faraday medium, a post-selection polarizer, a detector and an analyzer. The photon emitter projects a monochromatic light beam at the target wavelength λ substantially parallel to a magnetic field having strength B. The target wavelength is offset from established wavelength λ′ as λ=λ′+Δλ by wavelength difference of Δλ<<λ.

The light beam passes through the pre-selection polarizer and the prism. The beam is incident to an interface surface at incidence angle θ0 to the normal of the surface and is refracted into a secondary medium as first and second circularly polarized light beams separated by separation angle δ and having average refraction angle θ. The secondary medium has an index of refraction of n0.

The two circularly polarized light beams pass through the post-selection polarizer and reach the detector, which measures target pointer rotation angle Aw based on the separation angle δ. The analyzer determines the wavelength difference Δλ, first by calculating offset pointer rotation angle ΔAw=Aw−A′w from calibrated pointer rotation angle A′w based on the separation angle θ′ that corresponds to the established wavelength λ′, and second by estimating the wavelength difference based on

Δ λ - 2 ɛπ n 0 Δ A w cos θ VB sin θ 0 ,
in which ε is an amplification factor. Various exemplary embodiments also provide a method that incorporates operations described for the apparatus.

These and various other features and aspects of various exemplary embodiments will be readily understood with reference to the following detailed description taken in conjunction with the accompanying drawings, in which like or similar numbers are used throughout, and in which:

FIG. 1 is a diagram view of an optical diagram; and

FIG. 2 is a schematic view of a wavelength discriminator apparatus.

In the following detailed description of exemplary embodiments of the invention, reference is made to the accompanying drawings that form a part hereof, and in which is shown by way of illustration specific exemplary embodiments in which the invention may be practiced. These embodiments are described in sufficient detail to enable those skilled in the art to practice the invention. Other embodiments may be utilized, and logical, mechanical, and other changes may be made without departing from the spirit or scope of the present invention. The following detailed description is, therefore, not to be taken in a limiting sense, and the scope of the present invention is defined only by the appended claims.

In accordance with a presently preferred embodiment of the present invention, the components, process steps, and/or data structures may be implemented using various types of operating systems, computing platforms, computer programs, and/or general purpose machines. In addition, those of ordinary skill in the art will readily recognize that devices of a less general purpose nature, such as hardwired devices, or the like, may also be used without departing from the scope and spirit of the inventive concepts disclosed herewith. General purpose machines include devices that execute instruction code. A hardwired device may constitute an application specific integrated circuit (ASIC) or a floating point gate array (FPGA) or other related component.

This disclosure provides an overview of how a combined application of a recently discovered Faraday effect and weak value amplification can be used to measure Δλ=λ−λ′ (thereby detecting λ) when λ′ is known. The methodology has potential spectroscopic utility in such areas as measuring small Doppler shifts and detecting the presence of otherwise indistinguishable chemical or biological spectroscopic markers. Such measurements can be made using either classically intense laser light or single photon streams. For example, let λ and λ′ be two wavelengths such that λ=λ′+Δλ with Δλ being small.

A longitudinal magnetic field induces a circular differential refraction of a linearly polarized photon beam at the boundary between a Faraday medium and a medium with negligible Verdet constant as reported by A. Ghosh et al., “Observation of the Faraday effect . . . ” in Phys. Rev. A 76, 055402 (2007). See either http://www.rowland.harvard.edukjf/fischedimages/PRA76055402.pdf or else http://anciv.org/PS_cache/physics/pdf/0702/0702063v1.pdf for details.

This differential refraction is independent of the photon's pathlength through the Faraday medium and occurs within a few wavelengths of the boundary. The Verdet constant V represents an optical parameter that describes the strength of Faraday rotation from interaction between light and a magnetic field for a particular material, named for French physicist Émile Verdet.

FIG. 1 depicts an optical diagram 100 with accompanying coordinate reference frame 110 for a monochromatic linearly polarized coherent laser beam 120 of a wavelength λ. The beam 120 forms either a classically intense continuum or a stream of single photons that is incident from a Faraday medium 130 to a secondary medium 140 of negligible Verdet constant separated by an interface boundary 150. The Faraday and secondary media 130 and 140 have respective indices of refraction n± and n0. The subscripts plus (+) and minus (−) respectively correspond to right and left circular polarized radiation. The diagram 100 illustrates the paths taken by the light beam 120 at the interface between the Faraday medium 130 and the secondary medium 140 with negligible Verdet constant.

The beam 120 has an angle of incidence of θ0 from the normal to the interface 150. In the presence of a longitudinal magnetic field {right arrow over (B)} (having strength B), the beam 120 refracts at angles θ± from normal at the interface 150 into two circularly polarized beams having an angular divergence δ approximated as:

δ - λ sin θ 0 π n 0 cos θ VB , ( 1 )
where angle

θ = 1 2 ( θ + + θ - )
is the average of θ+ and θ, which are respectively the right- and left-circularized refraction angles, and V is the Verdet constant for the Faraday medium 130.

As shown, the incident beam's direction of propagation determines the y-axis of the reference frame 110. The x-axis is in the plane containing the beam 120 and the normal to the interface at the point of incidence. The origin of the reference frame 110 is defined by the perpendicular intersection of the x-axis with the y-axis at the interface 150. The usual z-axis (into the plane) with positive direction {circumflex over (z)}={circumflex over (x)}×ŷ completes the reference frame 110.

The longitudinal magnetic field B is assumed to be present and parallel to the positive y-axis. If the photon distribution of the incident beam 120 has a Gaussian distribution symmetric about the positive y-axis with mean value at x=0, then the refracted beams exhibit Gaussian distributions that are symmetric about their refracted paths which are along the vectors {right arrow over (u)}± in the x−y plane of the reference frame 110. More specifically, the refracted beams exhibit photon distribution mean values which are rotationally displaced around the z-axis through distinct angles θ±−θ0 from the positive y-axis in the direction of vectors {right arrow over (u)}± in the x-y plane, respectively.

This refraction process can be described from a quantum mechanical measurement perspective using the mean value of the intensity distribution profile produced by a detector as a measurement pointer. This description maintains validity for both a single photon stream and a classically intense beam.

In particular, an Hermitean operator  can be constructed and used to form a Hamiltonian operator Ĥ that describes a photon-interface interaction which produces the required geometry of the refraction process. Let |+custom character and |−custom character be the right and left circular polarization eigenstates, respectively, of the photon circular polarization operator {circumflex over (σ)} which obey the eigenvalue equation:
{circumflex over (σ)}|±custom character=±|±custom character  (2)
and have the orthogonality properties:
custom character±|±custom character=1
and
custom character±|∓custom character=0.  (3)

One can define the “which path” operator  as:
Â≡(θ+−θ0)|+custom charactercustom character+|+(θ−θ0)|−custom charactercustom character−|  (4)
and the associated interaction Hamiltonian H can be expressed as:
Ĥ=ÂĴzδ(t−t0).  (5)

Here the Dirac delta function δ(t−t0) encodes the fact that the refraction occurs within a few wavelengths of the interface 150 by modeling the refraction effectively as an impulsive interaction between a photon of the beam 120 and the interface 150 at time t0. The “which path” operator  accounts for the refractive angular displacements of the initial photon beam 120 at the interface 150. The operator Ĵz constitutes the measurement pointer's z-component of angular momentum, and couples the refractive angular displacements to the measurement pointer. One can note that:
Â,Ĵz┘=0,  (6)
and that |±custom character are eigenstates of  with respective eigenvalues (θ±−θ0).

FIG. 2 shows an elevation schematic 200 of an apparatus that employs this angular divergence. A laser 210 emits a photon beam 220 as analogous to the beam 120. The beam 220 passes through a pre-selection polarizer or polarizer 230 to reach a Faraday medium 240 (in the form of a prism), analogous to the medium 130. The refracted beam passes a post-selection polarizer or polarizer 250 to reach a detector 260 that measures the intensity distribution of the refracted polarization post-selected light beams.

For an initial photon polarization state |ψicustom character, i.e., the pre-selected state, and an initial (Gaussian) pointer state |φcustom character, the initial state of the combined pre-selected system and measurement pointer prior to the interaction at the interface 150 at time t0 constitutes the tensor product state |ψicustom charactercustom character. Note that the beam 220 has passed through the pre-selection filter 230 prior to its entry into the Faraday medium 240.

Immediately following the measurement's impulsive interaction, the combined system is in the state:

Ψ = - i H ^ t ψ i φ = - i A ^ J ^ z ψ i φ . ( 7 )
where use has been made of the fact that the integral of the delta function is:
∫δ(t−t0)dt=1.  (8)
Now let the initial polarization state be expressed as:
icustom character=a|+custom character+b|custom character,  (9)
in which a and b are complex numbers that satisfy the condition custom characterψiicustom character=1, and rewrite eqn. (7) as:

Ψ = - i A ^ J ^ z ( a + + b - ) φ . ( 10 )

Because of the orthogonality proportion of eqn. (9), the nth power of  assumes the form:
Ân=(θ+−θ0)n|+custom charactercustom character+|+(θ−θ0)n|−custom charactercustom character−|,n=0, 1, 2, . . .   (11)
Then the exponential term of the system state of eqn. (7) can be written as:

- i A ^ J ^ z = n = 0 [ - i A ^ J ^ z ] n n ! = n = 0 [ - i ( θ + - θ - ) J ^ z ] n n ! + + + n = 0 [ - i ( θ - - θ 0 ) J ^ z ] n n ! - - = - i ( θ + - θ 0 ) J ^ z + + + - i ( θ - - θ 0 ) J ^ z - - , ( 12 )
where i=√{square root over (−1)} is the imaginary unit and

= h 2 π
represents the reduced Planck constant. This result correlates refraction angle rotations with polarization.

The exponential operators constitute the rotation operators {circumflex over (R)}z:

- i ( θ ± - θ 0 ) J ^ z = R ^ z ( θ ± - θ 0 ) R ^ z ± . ( 13 )
These operators rotate the x- and y-axes through angles (θ±−θ0) around the z-axis of the reference frame 110. The rotation notation is consistent with the convention used by A. Messiah, Quantum Mechanics, v. 2, p. 1068 (1961).

This enables the system state in the {|{right arrow over (r)}custom character} representation to be rewritten as:
custom character{right arrow over (r)}|Ψcustom character=a|+custom charactercustom character{right arrow over (r)}|{circumflex over (R)}z+custom character+b|−custom charactercustom character{right arrow over (r)}|{circumflex over (R)}zcustom character.  (14)
The associated pointer state distribution in the {|{right arrow over (r)}custom character}-representation is then:
|custom character{right arrow over (r)}Ψcustom character|2=|a|2|custom character{right arrow over (r)}|{circumflex over (R)}z+custom characterφ|2+|b|2|custom character{right arrow over (r)}|{circumflex over (R)}z|custom characterφ|2,  (15)
and clearly corresponds to a sum of two Gaussian distributions |custom character{right arrow over (r)}|{circumflex over (R)}z±custom character|2 which are each symmetrically distributed about the vectors {right arrow over (u)}±, respectively.

A final photon polarization state |ψfcustom character that is post-selected can be expressed as:
fcustom character=c|+custom character+d|−custom character  (16)
in which c and d represent complex numbers that satisfy the condition custom characterψffcustom character=1. Note that the post-selection polarizer 250 receives the beam after refraction by the Faraday medium 240. From this, the resulting pointer state becomes:
custom charactercustom characterψfcustom character=ac*{circumflex over (R)}z+custom character+bd*{circumflex over (R)}zcustom character,  (17)
in which the asterisk denotes the complex conjugate, and its {|{right arrow over (r)}custom character}-representation distribution is:
|custom character{right arrow over (r)}|Φcustom character|2=|ac*|2|custom character{right arrow over (r)}|{circumflex over (R)}z+custom character|2+|bd*|2|custom character{right arrow over (r)}|{circumflex over (R)}zφcustom character|2+2Reac*bd*custom character{right arrow over (r)}|{circumflex over (R)}z+custom charactercustom character{right arrow over (r)}|{circumflex over (R)}2custom character*.  (18)

One may observe that although eqn. (18) constitutes a sum of two Gaussian distributions that are symmetrically distributed around vector {right arrow over (u)}±, unlike eqn. (15), this distribution also contains an interference term. Careful manipulation of this interference term can be described herein that produces the desired amplification effect.

In contrast to a strong measurement, a weak measurement of the “which path” operator  occurs when the uncertainty Δθ in the pointer's rotation angle is much greater than the separation between Â's eigenvalues, and when the interaction between a photon and the pointer is sufficiently weak so that the system remains essentially undisturbed by that interaction. In this case, the post-selected pointer state is represented as:

Φ = ψ f - i A ^ J ^ z ψ i φ ψ f ( 1 - i A ^ J ^ z ) ψ i φ ψ f ψ i - i A w J ^ z φ , ( 19 )
or else as:
custom character{right arrow over (r)}custom charactercustom characterψficustom charactercustom character{right arrow over (r)}|{circumflex over (R)}z(ReAw)|φcustom character,  (20)
where the quantity Aw is expressed as:

A w = ψ f A ^ ψ i ψ f ψ i , ( 21 )
and constitutes the weak value of operator Â. Note that rotation angle Aw is generally a complex value that can be directly calculated from the associated theory. One may also note that in response to |ψicustom character and |ψfcustom character being nearly orthogonal, the real value ReAw can lie far outside the spectrum of eigenvalues for Â.

The pointer state distribution for eqn. (20) is:
|custom character{right arrow over (r)}|Φcustom character|2custom characterψficustom character|2|custom character{right arrow over (r)}|{circumflex over (R)}z(ReÂw)|φcustom character|2,  (22)
and corresponds to a broad distribution that is symmetric around a vector in the x-y plane. That vector can be determined by a rotation of the x- and y-axes through an angle ReAw about the z-axis. In order that eqn. (20) be valid, both of the two following general weakness conditions for the uncertainty in the pointer rotation angle must be satisfied:

( a ) Δ θ » A w and ( b ) Δ θ » { min ( n = 2 , 3 , ) A w ( A n ) 1 n - 1 } - 1 , ( 23 )
as reported by I. M. Duck et al., “The sense in which ‘weak measurement’ of a spin-½ particle's spin component yields a value 100” in Phys. Rev. D 40, 2112-17 (1989). See http://prd.aps.org/pdf/PRD/v40/i6/p21121 for details.

Use of the above pre- and post-selected states |ψicustom character and |ψfcustom character—along with eqn. (11)—provides the following scalar expression for the weak value of the nth moment of “which path” operator Â:

( A n ) w = a c * ( θ + - θ 0 ) n + b d * ( θ - - θ 0 ) n a c * + b d * , ( 24 )
where c* and d* represent complex conjugates of c and d.

When n=1, then the first moment corresponds to the pointer's peak intensity. The first moment is:

A w = a c * ( θ + - θ 0 ) + b d * ( θ - - θ 0 ) a c * + b d * . ( 25 )

When the transmission axis of the pre-selection polarizer 230 is set so that:
a=sin φ,
and
b=cos φ,  (26)
and that of the post-selection polarizer 250 is set so that:
c=cos χ,
and
d=−sin χ,  (27)
then the quantity Aw becomes:

A w = Re A w = ( θ + - θ 0 ) sin ϕ cos χ - ( θ - - θ 0 ) cos ϕ sin χ sin ( ϕ - χ ) . ( 28 )

One can observe from this that the absolute value of the “which path” scalar |Aw| can be made arbitrarily large by choosing φ≈χ, i.e., separated by a small difference term ε. In particular, let χ=φ−ε and φ=π/4 (in which case the pre-selected state is linearly polarized in the x-direction). Consequently.

( a ) sin ϕ = cos ϕ = 2 / 2 , ( b ) cos χ = 2 2 ( cos ɛ + sin ɛ ) , ( c ) sin χ = 2 2 ( cos ɛ - sin ɛ ) , and ( d ) sin ( ϕ - χ ) = sin ɛ . ( 29 )

The previous relation from eqn. (28) for the amplified pointer rotation angle associated with the post-selected circularly polarized beams then becomes:

A w = ( θ + - θ 0 ) ( cos ɛ + sin ɛ ) - ( θ - - θ 0 ) ( cos ɛ - sin ɛ ) 2 sin ɛ , ( 30 )
which can be rewritten (by simplifying grouped terms) alternatively as:

A w = ( θ + - θ - ) cos ɛ + [ ( θ + + θ - ) - 2 θ 0 ] sin ɛ 2 sin ɛ . ( 31 )
This quantity is the pointer rotation angle, which can be conveniently related to the angular divergence δ and the difference term ε:

A w = δ 2 tan ɛ + ( θ - θ 0 ) . ( 32 )
Note that for small difference such that:
0<ε<<1,  (33)
then the rotation angle becomes arbitrarily large in magnitude and can be approximated as:

A w δ 2 ɛ , ( 34 )
and because of this, ε can be called the amplification factor

The weakness condition constraint follows when eqns. (24) and (25) can be used to obtain the associated weakness condition when incorporating 0<ε<<1 from eqn. (33) into inequalities eqn. (23), along with selections for a, b, c, d, φ, χ. These steps yield:

( a ) Δθ » δ 2 ɛ and ( b ) Δθ » 2 θ - θ 0 . ( 35 )
Here use is made of the fact that for 0<ε<<1 being sufficiently small, then:

min ( n = 2 , 3 , ) A w ( A n ) 1 n - 1 min ( n = 2 , 3 , ) ( θ + - θ 0 ) - ( θ - - θ 0 ) ( θ + - θ 0 ) n - ( θ - - θ 0 ) n 1 n - 1 = ( θ + - θ 0 ) - ( θ - - θ 0 ) ( θ + - θ 0 ) 2 - ( θ - - θ 0 ) 2 = ( θ + - θ 0 ) - ( θ - - θ 0 ) - 1 ( 2 θ - θ 0 ) - 1 . ( 36 )

Satisfaction of both conditions (a) and (b) of eqn. (35) requires that when 0<ε<<1 from eqn. (33), then the uncertainty Δθ greatly exceeds the absolute ratio value:

Δθ » δ 2 ɛ . ( 37 )
This condition can be satisfied by making the initial Gaussian pointer distribution width sufficiently large.

Thus, as per eqns. (22) and (34), the rotation of the initial photon distribution axis of symmetry provides an amplified measurement of the angular divergence δ via the weak value of the “which path” operator Â. For a known amplification ε and a measured mean value of the intensity distribution profile produced by the detector 260 corresponding to Aw, then angular divergence δ can be estimated from eqn. (34) as:
δ≈2εAw.  (38)

The real component of the complex operator Aw=ReAw corresponds to the angle between the direction of the resultant photon distribution peak and the positive y-axis is measured when there is sufficient knowledge of the value of the other parameters (e.g., θ, n0, V, etc.) appearing on the right hand side of this expression.

Consider the case where a target wavelength λ, can be expressed as the sum of a known wavelength λ′ and a difference wavelength Δλ:
λ=λ′+Δλ.  (39)
and similarly average refraction angle θ can be expressed as the sum of a known angle θ′ and a corresponding difference Δθ:
θ=θ′+Δθ,  (40)
so that approximation eqn. (1) can be written as:

δ - ( λ + Δλ ) sin θ 0 π n 0 cos ( θ + Δθ ) VB ( 41 )
or rewritten as the approximation:

δ - λ sin θ 0 π n 0 cos θ VB - Δλ sin θ 0 π n 0 cos θ VB δ + Δδ , ( 42 )
where, because Δθ is small, use has been made of the approximation
cos(θ′+Δθ)≈ cos θ′.  (43)

Using eqn. (42) in eqn. (34) renders:

A w δ 2 ɛ + Δδ 2 ɛ A w + Δ A w , ( 44 )
where A′w represents rotation angle corresponding to the mean value of the photon distribution profile associated with the measurement of the known wavelength λ′. Note that the rotation angle Aw=A′w when there is no wavelength difference or Δλ=0=ΔAw.

The apparatus represented by the diagram 200 in FIG. 2 can be used to detect λ=λ′+Δλ from eqn. (39) and estimate Δλ when λ′ is known. In order to accomplish this, the apparatus must first be calibrated so that its pointer value is A′w when the source is monochromatic with a known wavelength λ′. In particular, the photon distribution peak A′w is determined by enabling monochromatic light of wavelength λ′ to traverse the apparatus that comprises a Faraday medium with Verdet constant V.

In this example, the parameters θ0, n0, B are fixed and the polarizers are set per above values to provide an amplification factor ε. When light of wavelength λ (bichromatic or monochromatic) traverses the calibrated apparatus, then the pointer deviates from the calibrated pointer value A′w by the amount

Δδ 2 ɛ .

As numerical examples, consider an additional two instances in which the medium 140 with negligible Verdet constant is air. For the first example, let the Faraday medium 130 be terbium gallium garnet (formula Tb3Ga5O12) which has a Verdet constant V=−134 rad·Tesla−1·m−1 at known wavelength λ′=632.8 nm (red light) so that

Δ A w 21.33 ( Δλ · B ɛ ) rad · Tesla - 1 · m - 1 .

If for this first example, the amplification is ε=10−4, the wavelength difference is Δλ=1 pm and magnetic field strength is B=1 Tesla, then rotation angle difference is ΔAw≈21.33 μrad, provided that the associated weakness condition of minimum rotation angle difference

Δθ » 1.35 ɛ × 10 - 5 rad = 0.135 rad
is satisfied. If the detector 260 is 1 m (one meter) from the Faraday medium 240, then the pointer is translated by −21 μm in the detector plane. Thus, the 1 pm spectrum separation has been amplified by a factor of ˜106 at the detector 260.

For the second example, let the Faraday medium be MR3-2 Faraday rotator glass which has a Verdet constant V=−31.4 rad·Tesla−1·m−1 at known wavelength λ′=1064 nm (and 20° C. temperature) so that the rotation angle difference is

Δ A w 4.997 ( Δλ · B ɛ ) rad · Tesla - 1 · m - 1 .
If ε=10−3, Δλ=1 nm and B=1 Tesla, then ΔAw≈4.997 μrad provided that the associated weakness condition

Δθ » 5.317 ɛ × 10 - 6 rad = 5.317 m rad
is satisfied.

If the detector 260 is 1 m (one meter) distant from the Faraday medium, then the pointer translates by ˜5.0 μm in the detector plane. Thus, the 1 nm spectral separation has been amplified by a factor of ˜103 at the detector 260 for the medium 140 with negligible Verdet constant (and unitary refraction index) being air.

As described above, known beam wavelength λ′, magnetic field strength B and angle-of-incidence θ0 are established á priori. The refraction angle θ′ represents the average of the refraction angles for the circularly polarized beams determined for the known wavelength λ′ based on the indices of refraction n± of the Faraday medium 130.

For small differences such that Δλ≡λ−λ′<<λ between unknown and known wavelengths, the average angular refraction angle difference between the unknown (i.e., target) and known refractions Δθ≡θ−θ′ is small. This enables the average refraction angle to be reasonably approximated by the known value as θ≈θ′. The approximation cos θ≈ cos θ′ can be made, because of the relation:
cos(θ′+Δθ)=cos θ′ cos Δθ−sin θ′ sin Δθ≈ cos θ′−Δθ sin θ′≈ cos θ′.  (45)

The distribution peak rotation angle A′w corresponding to the known wavelength λ′ can be established from a calibration measurement. For a small amplification factor E such that 0<ε<<1, the rotation angle Aw corresponding to the unknown wavelength λ provides an estimate of the divergence δ from eqn. (38) as δ=2εAw. Similarly, the known rotator angle A′w the estimate for the corresponding divergence δ′.

In response to a small change in wavelength from the known λ′ to an unknown λvalue, the measured rotator angle becomes Aw=A′w+ΔAw with the difference ΔAw corresponding to the change in measured rotation angle due to the offset wavelength that provides a measure of change in dispersion angle from eqn. (42) as Δδ. Subtracting the calibrated values and rearranging terms enables the wavelength difference to be determined as:

Δλ = - π n 0 Δδcos θ VB sin θ 0 - 2 ɛπ n 0 Δ A w cos θ VB sin θ 0 . ( 46 )

While certain features of the embodiments of the invention have been illustrated as described herein, many modifications, substitutions, changes and equivalents will now occur to those skilled in the art. It is, therefore, to be understood that the appended claims are intended to cover all such modifications and changes as fall within the true spirit of the embodiments.

Spence, Scott E., Parks, Allen D.

Patent Priority Assignee Title
8970844, Aug 24 2011 The United States of America as represented by the Secretary of the Navy Signal amplifier by triple Mach-Zehnder interferometers
Patent Priority Assignee Title
5008622, Dec 21 1989 United States Department of Energy Superconductive imaging surface magnetometer
5109196, Mar 05 1991 Vanderbilt University Method and apparatus for magnetic identification and localization of flaws in conductors by canceling the field about the conductor with the field about a flawless conductor
5333142, Oct 26 1992 The United States of America as represented by the Secretary of the Navy Technique for intracavity sum frequency generation
6885192, Feb 06 2002 Regents of the University of California, The SQUID detected NMR and MRI at ultralow fields
7193413, Jun 23 2003 HITACHI HIGH-TECH CORPORATION Magnetic field measurement apparatus
///
Executed onAssignorAssigneeConveyanceFrameReelDoc
Jun 30 2011PARKS, ALLEN D UNITED STATES OF AMERICA, REPRESENTED BY SEC OF NAVYASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0266840247 pdf
Jun 30 2011SPENCE, SCOTT E UNITED STATES OF AMERICA, REPRESENTED BY SEC OF NAVYASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0266840247 pdf
Jul 06 2011The United States of America as represented by the Secretary of the Navy(assignment on the face of the patent)
Date Maintenance Fee Events
Jul 11 2017M1551: Payment of Maintenance Fee, 4th Year, Large Entity.
Nov 08 2021REM: Maintenance Fee Reminder Mailed.
Dec 03 2021M1552: Payment of Maintenance Fee, 8th Year, Large Entity.
Dec 03 2021M1555: 7.5 yr surcharge - late pmt w/in 6 mo, Large Entity.


Date Maintenance Schedule
Mar 18 20174 years fee payment window open
Sep 18 20176 months grace period start (w surcharge)
Mar 18 2018patent expiry (for year 4)
Mar 18 20202 years to revive unintentionally abandoned end. (for year 4)
Mar 18 20218 years fee payment window open
Sep 18 20216 months grace period start (w surcharge)
Mar 18 2022patent expiry (for year 8)
Mar 18 20242 years to revive unintentionally abandoned end. (for year 8)
Mar 18 202512 years fee payment window open
Sep 18 20256 months grace period start (w surcharge)
Mar 18 2026patent expiry (for year 12)
Mar 18 20282 years to revive unintentionally abandoned end. (for year 12)