The present disclosure relates to optical circuits, and more particularly to optical flip flop circuits.
Optical flip-flops are the key elements of all-optical memory and buffering devices. The development of such devices remains in its infancy, even though flip-flops would be valuable in contemporary communication networks. The reason for this is that flip-flops are not yet fast, robust, and low-power enough to be used in the large numbers that these applications demand. The most common means of implementing optical flip-flops is the use of active semiconductor devices. Both semiconductor optical amplifiers and semiconductor lasers have been widely used for this purpose. An alternative would be to use the Kerr effect in a passive, bistable resonator. This second approach has two potential advantages. First, passive flip-flops do not require current injection and devices that require cascading of such elements could be implemented with lower power requirements. Second, the Kerr effect has an almost instantaneous material response. As a result, switching speeds are not limited by the nonlinear medium but by the photon lifetime of the cavity, which can, in principle, be engineered to be as fast as necessary. Implementation of Kerr flip-flops using microresonators has proven to be difficult. When the optical power in a resonator is large, material absorption significantly heats it. The resulting thermo-optic change in refractive index is much stronger and slower than the Kerr effect, and as a result, switching between the “on” and “off” states is limited to microsecond time scales. Techniques for getting around this limitation have included cooling the cavity to cryogenic temperatures and using pulsed input fields with temporal durations much smaller than the cavity's thermal response time.
Two-input bistability, for which the resonator is equally full of light in both stable states, has the potential to circumvent this limitation. In such a device the thermal effect would be equally strong in both ‘on’ and ‘off’ states, producing a background refractive index change that does not respond on a time scale comparable to that of the much faster Kerr nonlinearity.
The present invention can be embodied as a two-input bistable resonator that can be switched in a robust way by phase modulation of the input beams. An optical flip-flop of the present invention may be embodied as a two-input Kerr resonator. The disclosed switching mechanism is compatible with cross-phase modulation induced by set and reset pulses for realization of an ultrafast, passive, all-optical memory element.
For a fuller understanding of the nature and objects of the invention, reference should be made to the following detailed description taken in conjunction with the accompanying drawings, in which:
The present invention may be embodied as an optical flip-flop 10, shown schematically in
When both input fields (Ain(1) and Ain(2)) are equally intense and equally detuned from resonance, the asymmetry between SPM and XPM effects can lead to the existence of two stable states, as illustrated in
To quantify the performance of the two-input optical flip-flop of the present embodiment, the methodology disclosed in B. A. Daniel, D. N. Maywar, and G. P. Agrawal, “Dynamic mode theory of optical resonators undergoing refractive index changes,” J. Opt. Soc. Am. B 28, 2207-2215 (2011), incorporated herein by reference, may be used. Use of such methodology leads to the following two coupled nonlinear equations for the mode amplitudes a1 and a2:
where |ak(t)|2 represents the optical energy stored in mode k at time t, Δωk=ωk−ωr is the detuning from the cavity resonance ωr, and Ain(k) is the input field normalized so that |Ain(k)|2 is the optical power. Coupling of the input field into the resonator is governed by κ; it is given by κ=(2τph)−1/2 when cavity losses from absorption and scattering are relatively small. The SPM and XPM effects are included through the nonlinear parameters γkl given by:
where n2 is the Kerr coefficient, Vk is the effective mode volume, and nkl is a parameter which measures how well the resonator mode overlaps with the nonlinear medium. Typically nkl≈1 and Vk≈Vl. In what follows, the approximation γ12≈γ22≈γ11=γ is made and equations (1) and (2) are solved numerically. The power transmissivity of the two output ports is calculated using:
T
k
=|A
out
(k)
/A
in
(k)|2=|κak/Ain(k)|2. (4)
Before considering the phase-switching dynamics, the stable steady states of the optical flip-flop are identified by considering the continuous-wave (“CW”) case for which Ain(k)=√{square root over (Pk)} is a constant. The steady-state solutions are found by setting the time derivatives in equations (1) and (2) equal to zero. Introducing the mode energy Ek=|ak|2, the following set of two algebraic equations (k=1 or 2) are found:
E
k[(1/2τph)2+(AΔωk+γEk+2γE3-5)2]=κ2Pk. (5)
Considering a flip-flop which is biased such that the two input beams have the same power (P1=P2=P0) and the same detuning (Δω=Δω2=Δω0), the solutions of equation 5 can be divided into two categories: (1) symmetric (E1=E2) solutions, and (2) asymmetric (E1≠E2) solutions. The asymmetric solutions come in pairs since the roles of E1 and E2 can be reversed. Using known techniques,
Phase Modulation of Two-Input Bistable Resonators
Phase switching dictates that input phases of the two beams may change with time such that Ain(k)(t)=√{square root over (P0)}exp[iφk(t)]. A time-dependent phase is equivalent to imposing a frequency chirp which, in turn, modifies the biasing conditions in a transient fashion. Mathematically, an instantaneous change in the detuning of an input beam from resonance is given by:
This change in detuning modifies the available steady states towards which the system will evolve.
The phase of an input optical beam can be changed through a variety of techniques. For example, a phase modulator can be used for this purpose to switch a flip-flop electro-optically. An all-optical flip-flop can be switched using set and reset optical pulses at wavelengths that are different from those of the CW inputs used to bias the device. In devices and methods according to the present invention, the XPM phenomenon is used to modulate the phases of the two CW inputs. As an example, assuming a Gaussian shape for set and reset pulses of width T0, the XPM-induced phase shift can be written in the form:
φk(t)=φ0exp[(−(t−tp)2/T02], (7)
where φ0 is the maximum phase shift occurring at the location t=tp. The phase modulation induces a negative detuning of the CW bias beams for t<tp and a positive detuning for t>tp. The largest value of the detuning is approximately given by:
Δωkmax≈Δω0+0.86φ0/T0. (8)
If the ratio φ0=T0 is sufficiently large (˜1/τph), and if this detuning is maintained for a duration long enough (˜τph) that the resonator can respond, then the phase modulation will flip the device.
Numerical Simulation
An exemplary device of the present invention is shown to flipflop in a robust manner through numerical simulation by first solving equations (1) and (2) for phase-modulated input beams. The set and reset operations are performed by applying Gaussian phase modulations of the form in equation (7) to input beams 2 and 1, respectively. A pulse width of T0=2τph is assumed, for which simulations show that switching occurs over the range of maximum phase shifts 2.3<φ0<9, leaving a broad window for phase error.
In
The present invention may be embodied as a method 100 for phase-shift switching of a bistable optical resonator as a flip-flop, wherein the resonator has a first input and a second input. As such, the bistable resonator is switched between a first stable state and a second stable state. The resonator is biased 103 so as to have only two stable asymmetric states as described above. It should be noted that such biasing 103 may or may not make up a step of the present method 100 (e.g., the biasing 103 condition exists before implementation of the method 100). The method 100 comprises the step of modulating 106 the phase of the first input holding beam momentarily such that the resonator switches from the second stable state to the first stable state. The method 100 may further comprise the step of modulating 109 the phase of the second input holding beam momentarily such that the resonator switches from the first stable state to the second stable state.
The momentary modulation 106, 109 of the holding (input) beams can be considered as set/reset pulses of a flip-flop. The set and/or reset pulses of the method 100 may have a Gaussian shape.
The method 100 may further comprise the step of determining 112 whether the resonator is in the first stable state or the second stable state.
The present invention may be embodied as an flip-flop system 50 comprising a bistable optical resonator 52 having a first stable state and a second stable state. The resonator 52 has a first input 54 and a second input 56. The system 50 further comprises a optical source 58 (e.g., a laser source) for providing input holding beams to the first input 54 and the second input 56. The system 50 further comprises a modulator 60 configured to selectively modulate a phase of the first input 54 holding beam such that the resonator 52 is in the first stable state (either switching to the first stable state or remaining in the first stable state—depending on the initial state of the resonator 52). The modulator 60 may further be configured to selectively modulate a phase of the second input 56 holding beam such that the resonator 52 is in the second stable state (either switching to or remaining in the second stable state).
The resonator 52 may further comprise a first output 62, a second output 64, or both. The first output 62 and second output 64 may be used to determine the state of the resonator 52 (i.e., determine whether the resonator 52 is in the first stable state or the second stable state).
The model, further described below, makes no assumptions about the specific structure of a suitable optical cavity. The model describes exemplary embodiments of devices including, ring resonators, Fabry-Perot cavities illuminated by two input beams of the same frequency propagating at different angles (previously developed), and Fabry-Perot cavities with two resonator modes, having different resonance frequencies, propagating in the same direction. It is a dynamic model that can be used to describe the switching process, as well as the steady-state field behavior.
In this third configuration, a device (see, e.g.,
A. Dynamic-Mode Amplitude Equations
The model expands on the known dynamic-mode theory of resonators. The electric field in the cavity is written as a sum of the two resonant modes that are excited by the input fields
where ek(r) is the electric-field profile associated with the kth mode (k=1, 2) of the resonator. The mode amplitudes are normalized so that |ak(t)|2 is the electromagnetic energy stored in mode k at time t. The constant Nk is a normalization factor given by (10)
N
k=1/2∫∈0∈(r)|ek|2d3r, (10)
where the dielectric permittivity ∈(r) describes the structure of the cavity.
When the medium inside the cavity exhibits the Kerr effect, there is, in addition to the linear response described by the permittivity, a nonlinear response described by the third-order dipole-moment density
where xμαβγ(3) is the third-order susceptibility tensor. It has been shown that, for an itafly assumed solution of the form (9), Maxwell's equations imply the following equations for the mode amplitudes:
where τph is the photon lifetime of the cavity; κ is a coupling coefficient; and Ain(k) is the input field to mode k, which is normalized so that |Ain(k)|2 is its optical power. It is assumed for the sake of simplicity that τph and κ are the same for the two modes. Using equation (11) in equation (12) together with equation (9), we obtain the following two coupled nonlinear differential equations:
In deriving these equations, a number of terms have been neglected. Depending on which device configuration is being considered, there are different justifications for this neglect. For a configuration using two spectrally distinct modes, these additional terms are not frequency matched to the mode resonances, and hence they have a negligible influence on the mode amplitudes. In the case of modes that may have the same resonance frequency but are spatially distinct, such as counter-propagating modes in a ring resonator, these terms are vanishingly small as a result of the spatial phase structure of the modes.
The SPM and XPM terms appearing in equations (13) and (14) depend on a set of four nonlinear parameters given by (see also, equation (3))
In equation (15) Vk is an effective mode volume defined as
The parameter n2 that appears in equation (15) is the Kerr coefficient responsible for the intensity dependence of refractive index in the nonlinear medium. In general, n2 depends on the orientation of the electric field with respect to the crystallographic axes. In practice, n2 is chosen to be the value for some particular crystallographic direction. The third-order susceptibility xc(3) in this direction is related to n2 as
x
c
(3)=4/3∈0cn02n2, (17)
where n0 is the linear (low-intensity) refractive index of the medium. The parameter nkl that appears in equation (15) is a nonlinear overlap factor given by
Physically, the nonlinear overlap factors measure how effectively the modes interact through the third-order susceptibility. It is often a good approximation to take nkl≈1. It is also often a good approximation to take Vk≈Vcav, where Vcav is the volume of the cavity. If the two input fields additionally have nearly the same frequency, then γ11≈γ22≈γ12≈γ21≈γ, where
γ≈ω1cn2(n02Vcav) (19)
This approximate form of γ will be used throughout the rest of this section.
B. Steady-State Solutions
Equations (13) and (14) describe the behavior of the resonator for any input-field temporal profiles Ain(1)(t) and Ain(2)(t). For designing an optical flip-flop, we are interested in the stable, steady states of the resonator when two continuous-wave (CW) fields with constant intensities are launched into it. The input fields then take the form
where Bk are constants. For such input fields, the steady-state solutions of equations (13) and (14) take the form
Using equations (20) and (21) in equations (13) and (14), we obtain a pair of algebraic equations for the complex constants bk,
[−i(Δωk+γ|bk|2+2γ|b3-k|2)+1/2τph]bk=κBk, (22)
where Δωk=ω0k−ωk is the detuning of the kth input field from resonance. These equations result in the following pair of coupled equations for the mode energies Ek=|bk|2,
[(Δωk+γEk+2γE3-k)2+(1/2τph)2]Ek=|κ|2Pk, (23)
where Pk=|Bk|2 is the power of the kth input field. Once a solution is found for the mode energies by solving equation (23), equation (22) can be used to find the phases of bk.
Solving equation (23) does not, however, guarantee that the resulting solution represents a physically realizable state of the device. In order to be realizable, it is also necessary that the solution be stable. A stable solution is characterized by its being robust to small perturbations. If a small perturbation is applied to a stable state, it tends to die out and the system remains in that state. In contrast, if a small perturbation is applied to an unstable state, the system evolves away from that state and does not return to it. Stability of various solutions of equation (23) can be examined by performing a linear stability analysis of equations (13) and (14), as outlined in the following discussion.
Linear Stability Analysis
Consider the solutions of equation (22) for the mode amplitudes bk. These correspond to steady-state solutions of equations (13) and (14) of the form of equation (21). Now imagine that the fields are slightly perturbed, such as would occur regularly in a real device. After the perturbation the solutions of equations (13) and (14) can be written in the form
Where |ck|<<1. Using this form in equations (13) and (14), the fact that bk satisfy equation (22), and neglecting all terms higher than the first order in ck, we obtain a system of four differential equations that describe how the perturbation evolves in time. We can write them in matrix form as
where c=[c1 c2 c1*c2*]T is a column vector and the matrix is given by
Any solution of equation (25) can be written as a linear combination of its eigenmodes. The eigenmodes are constructed from the eigenvectors and corresponding eigenvalues of the matrix which satisfy
cm=λmcm. (27)
The eigenmode solutions are given by cmeiγλ
Im{λ
m}>0, (28)
which is the criterion used to determine the stability of steady-state solutions in this disclosure.
Flip-Flop Design Criteria
A. Biasing Conditions
The available stable states of the flip-flop depend on the power levels and detunings of the two input fields, as indicated by equation (23). The flip-flop could be designed such that these properties are different between the two fields, but this would be undesirable. As discussed above, thermal nonlinearities can be a problem for a flip-flop having different intracavity intensities in its two states. As a result, it is desirable to bias the device symmetrically so that the input fields have the same power and detuning from their respective resonances. In this situation, we can expect a pair of states to exist that exhibit the same intracavity energy, such as those depicted in
Δω1=Δω2Δω0, (29)
P
1
=P
2
=P
0, (30)
and then examining their stability using the linear stability analysis presented above.
It is known that the solutions of equation (23) under symmetric biasing conditions can be divided into two categories. The first category contains all symmetric solutions, characterized by the equality of their mode energies (E1=E2). The second category contains all asymmetric solutions for which E1≠E2. Because the biasing conditions are themselves symmetric, the asymmetric solutions come in pairs because the roles of E1 and E1 can always be reversed.
The solutions depicted in
B. Phase-Modulation Profile
We consider switching of the two-input Kerr resonator by pure phase modulations of the input fields. This phenomenon is modeled by solving equations (13) and (14) with input fields of the form
where φk (t) is a time-dependent phase imposed on the field by a control signal. This could be accomplished electrically by using an electro-optic phase modulator, or optically using XPM of an input field by a set or reset pulse. The result of this phase modulation is to temporarily change the detuning of the field's frequency from resonance as
If the phase modulation is slow enough for the resonator to respond, then its effect can be understood as temporarily modifying the biasing conditions and, hence, the available stable states toward which the system will evolve.
The phase switching of an optical flip-flop can be understood using
T
k
=|κa
k
/A
in
(k)|2. (33)
We assume that resonator loss is dominated by the two couplers so that |κ|2=1/2τph. If a phase modulation with a positive derivative is applied to input field 2 so that Δω2τph<−2, then
In practice, it may be necessary to turn off the signal that applies the modulation after a short time interval. As an example, we consider Gaussian phase shifts of the form
where φ0 is the maximum phase shift occurring at time tp, and T0 is a measure of the temporal duration of the phase modulation. Positive values of φ0 allow for switching to occur because the trailing edge of the modulation determines the final state of the device after the signal is gone. Thus, even though a positive value of φ0 increases the phase over the leading edge of the signal, which does not necessarily switch the device, the trailing edge creates a decreasing phase shift that can switch the device under the appropriate conditions. Maximum detuning can be derived using equations (32) and (34) and is found to be
Δωkmax≈Δω0+0.86φ0/T0. (35)
In order for the device to switch, this maximum detuning should be large enough to drive the device into the shaded switching region in
φ0>T0/τph. (36)
Equation (36) is not sufficient for a phase modulation to switch the device. Also, the modulation should occur over a long enough temporal duration that the resonator can respond. This leads to the following second criterion:
T
0>τph. (37)
The approximate criteria for the phase-modulation parameters in equations (36) and (37) are verified by rigorous numerical solutions of equations (13) and (14). The input fields are taken to be of the form of equation (31) with phase modulations of the form of equation (34). The flip-flop is biased using two CW fields with detunings Δω0τph=−2 and powers given by 4γ|κ|2τph3P0=2. As in
Initially, the flip-flop is in a state for which the transmission of input field 2 is high and the transmission of input field 1 is low. The impact of phase switching is then studied by applying three phase modulations of the form given in equation (34) at times tp=0, 50τph, and 100τph. The first modulation is applied to input field 2 (set operation) and the remaining two to field 1 (reset operations).
We consider first the role of the maximum phase shift φ0 and fix the duration of the phase modulation at T0=2τph.
As the maximum phase shift is further increased to π and 2π, as shown in
We consider next the role of the duration T0 of the phase modulation. For this purpose, we fix the maximum phase shift at a value of φ0=π.
Cavity Design Considerations
The analysis has so far assumed only one kind of material nonlinearity: the Kerr effect. There are, however, other effects associated with a given material system.
Semiconductors often exhibit nonlinear loss mechanisms such as two-photon absorption and subsequent free-carrier absorption, which are may prevent the device from functioning. If optical fields were employed with photon energies below the half-bandgap, this problem could be avoided so that the large Kerr nonlinearity of semiconductors could be leveraged for low-power operation. Other candidate material systems include silicon dioxide and silicon nitride. Neither of these exhibit significant nonlinear absorption near the 1.55 μm wavelength, and fabrication of resonators is technologically well developed for both of these materials.
There is one type of material nonlinearity that cannot be avoided. Even if the medium is nearly transparent, some optical power is always lost through material absorption. As discussed in above, this absorption heats the cavity and changes its refractive index via the thermo-optic effect. For a two-input flip-flop, this thermal index change does not affect the switching process because it has the same value in both of the stable states and responds on a time scale much longer than that required for switching. It does, however, present two practical problems. First, if the temperature change is large enough it can physically damage or destroy the resonator. Second, it presents a technical challenge for turning on the device. The reason for this can be understood as follows: in the cold-cavity state, before the input fields are turned on (i.e., when Ain(1)=Ain(2)=0) the kth mode's resonance frequency is ωk. In the operation mode, when both lasers are on and bistability has been achieved, the new resonance frequency is ωk′=ωk+ΔωTΔωKerr(k), where ΔωT is the resonance shift from the cavity's change in temperature, and ΔωKerr(k) is the resonance shift from the Kerr effect. When the lasers are first turned on, they might be so far from the cold-cavity resonances that neither input field is able to couple into the cavity and heat it up. As a result, simply turning the lasers on will not transition the device from the cold-cavity state into operation mode. In order to do this, it may be necessary to come up with a means of heating the cavity. This could be done using a thermo-electric temperature controller, or by sweeping the frequency of one of the input fields to “drag” the cold-cavity resonances close to the operation frequency. Either way it will be desirable that ΔωT, be small compared to ΔωKerr(k)). Its value will depend on the materials used to make the device as well as on the geometric structure of the resonator.
A. Simple Thermal Model
To analyze the influence of material properties on the temperature shifts that are responsible for ΔωT, we develop a thermal model of the device. Consider the influence of a weak material absorption on the cavity's photon lifetime. The overall photon lifetime is given by
where τphab, τphsc, and τphe represent contributions from material absorption, ph phph scattering losses, and coupling losses, respectively. Using the perturbative theory disclosed above, it can be shown that the contribution from material absorption is given by
where it has been assumed that the optical field is primarily confined to a single material with absorption coefficient αm and refractive index n0. By considering the rate at which electromagnetic energy is lost owing to this absorption, we find that the thermal energy ΔUT stored in the cavity changes with time as
where τT is the thermal lifetime.
The thermal energy is related to the thermo-optic change in refractive index as
where ∂n/∂T is the medium's thermo-optic coefficient, Cp is its specific heat capacity, and ρ is its density. Equations (40) and (41) indicate that ΔnT evolves with time as
Equation (42) describes how the optical field in the resonator influences the thermal index shift. The influence of the index shift ΔnT on the mode amplitudes is incorporated in the dynamic equations (13) and (14) by adding another term so that they become
The last term in this equation represents a thermally induced shift of the cavity's resonance frequency by an amount ΔωT=ωkΔnT/n0.
The thermal resonance shift of the device in operation mode can be calculated by considering the steady-state solution of equations (42) and (43), and it is found to be
where Ecav=E1+E2 is the optical energy stored in the cavity. Because the optical field in the cavity is dominated by one mode in each of the states (e.g., the kth mode), the optical energy can be approximated by that of this dominant mode (Ecav≈Ek). With this approximation, the Kerr-induced resonance shift of the dominant mode can be shown from equation (43) to be
where equation (19) has been used for γ. Equations (44) and (45) imply that the ratio of the thermal resonance shift to the Kerr-induced resonance shift is
Because γT depends primarily on material properties, it can be seen as a figure of merit for comparing candidate material systems out of which to construct the resonator. The only term in equation (46) that is not purely a material parameter is τT, which measures how quickly the cavity dissipates heat with the surrounding environment and has some dependence on its geometric structure. Active cooling of the cavity can reduce the effective value of τT and improve the figure of merit. Models that show the influence of material properties and cavity structure on τT have also been developed in for different types of cavities.
B. Influence of Resonator Structure
The resonator structure affects two parameters, the quality factor (Q) and cavity volume (Vcav), that influence device performance. Both of these should be selected so that the flip-flop can be made with a low enough bias power, low enough operating temperature, and fast enough speed.
The power of each of the input fields at the bias point that we have used in the preceding analysis is given by 2γτph2P0=2 when the resonator loss is dominated by coupling. Using this expression and equation (19) for γ leads to the following estimate of the required bias power:
This relation shows that the required power is proportional to the volume of the cavity. This makes sense because a smaller volume will require less input power to achieve the same intracavity intensity and therefore the same Kerr-induced resonance shift. Equation (47) also indicates that the needed bias power depends inversely on Q2. Physically, one factor of Q results from the fact that a larger quality factor implies a proportionately higher cavity enhancement of the optical power of an input field. The other factor of Q comes from the fact that a larger quality factor implies a proportionately smaller bandwidth of the cavity resonance. Bistable operation is achieved when the Kerr-induced shift of a resonance is comparable to its bandwidth. Because the Kerr-induced shift is proportional to the optical power inside the cavity, a smaller bandwidth results in a proportionately lower requirement for the bias power.
The temperature shift of the cavity when the flip-flop becomes operational is also important because too much heating can damage it. The temperature shift can be calculated by solving equation (40) for the steady-state thermal energy in the cavity and using the relation ΔT=ΔUT/ρCρVcav, where ΔT is the temperature shift. The result is
where the steady-state intracavity energy has been approximated by Ecav≈2/γτph and equation (19) has been used for γ. Equation (48) indicates that the operating temperature depends inversely on the quality factor, but that it does not depend on the cavity volume. The independence on cavity volume results from the fact that the required power scales with Vcav [equation (47)]. Thus, the intracavity intensity and, hence, the operating temperature do not depend on Vcav when the device is appropriately engineered. The inverse dependence of the operating temperature on Q occurs because a smaller Kerr-induced frequency shift is needed if the resonator has a higher quality factor. The Kerr-induced resonance shift and the thermal resonance shift are directly proportional [equation (46)], and the thermal shift is directly proportional to the change in temperature. A higher Q therefore implies a proportionately smaller temperature shift.
It might be concluded from equations (47) and (48) that it is desirable to design the cavity to have the largest Q possible. However, this is not necessarily the case because the quality factor also determines the device's switching speed. Numerical simulations in
Although the present invention has been described with respect to one or more particular embodiments, it will be understood that other embodiments of the present invention may be made without departing from the spirit and scope of the present invention. Hence, the present invention is deemed limited only by the appended claims and the reasonable interpretation thereof.