Abstract:
A sensor for detecting both magnetic fields and electric fields can include at least one Sawyer-Tower (ST) circuit that can incorporate a multiferroic capacitor. An odd number of ST circuits coupled together in a ring configuration, so that for each ST circuit, the output of one ST circuit is an input to another of the ST circuits. The multiferroic capacitors can include a multiferroic layer that can be deposed on a substrate. For each multiferroic capacitor, the deposition process can cause an inherent amount of impurities in the multiferroic layer. The number of said odd number of ST circuits to be coupled together is chosen according to the amount of impurities, to “forgive” the impurities. The higher the level of impurities in the multiferroic layers, that more ST circuits that are required in the ring to achieve the same sensor sensitivity. BDFO can be chosen for the multiferroic material.

Description:
FEDERALLY-SPONSORED RESEARCH AND DEVELOPMENT 
     The United States Government has ownership rights in this invention. Licensing inquiries may be directed to Office of Research and Technical Applications, Space and Naval Warfare Systems Center, Pacific, Code 72120, San Diego, Calif., 92152; telephone (619)553-5118; email: ssc_pac_t2@navy.mil, referencing NC 102801. 
    
    
     FIELD OF THE INVENTION 
     The present invention pertains generally to sensors. More particularly, the present invention pertains to sensors for detecting both electric fields and magnetic fields. The present invention can be particularly, but not exclusively, useful as a sensor that can incorporate multiferroic materials in order to quantify magnetization as well as electric field by exploiting the nonlinear electrical behavior of the multiferroic element, using an applied electric field, instead of an applied (driving) electric current. 
     BACKGROUND OF THE INVENTION 
     Magnetometers, or sensors for detecting magnetic fields, are well known in the prior art. Such prior art magnetometers can consist of a small, magnetically susceptible core wound by two coils of wire. An alternating electrical current can be passed through one coil, which can induce an electrical current in the second coil, and this output current, mediated by the magnetically susceptible core, can be measured by a detector. In a magnetically neutral background, the input and output currents will match. However, when the core is exposed to a background field, it can be more easily saturated in alignment with that field and less easily saturated in opposition to it. Hence the alternating magnetic field, and the induced output current, will be out of phase with the input current. The extent to which this is the case will depend on the strength of the background magnetic field. Often, the current in the output coil can be integrated to yield an output analog voltage, which can be proportional to the magnetic field. But for these types of sensors, an applied current is required. 
     Multiferroics, or materials that simultaneously exhibit magnetic and ferroelectric orders, are also known in the prior art. These materials can often also be termed as magnetoelectrics, because the material magnetic and electric order parameters are coupled. Multiferroics can be technologically important, as they can have two or more switchable states, like a magnetization state that may be switched with an electric field, and a spontaneous electric polarization state that may be switched with a magnetic field. Such materials can play a vital role in the design of electric-field controlled ferromagnetic resonance devices, actuators, and variable transducers with magnetically-modulated piezoelectricity etc. Additionally, magnetoelectrics can also have tremendous potential for use in storage devices where writing and read-out can be carried out by both/either of electric or magnetic fields. In sum multiferroics can be important to any device where it is preferable to use an electric field, instead of electric current, to operate the device. 
     Interestingly, recent experimental and theoretical studies explicitly reveal novel behavior and exciting physics in multiferroic materials. Being of great interest, and being motivated by on-chip integration in microelectronic devices, nanostructured composites of ferroelectric and magnetic oxides deposition as a thin film on a substrate are being increasingly studied. Such initiatives are expected to lead to a better understanding of the basic nature of magneto-electric coupling so that the magneto-electric coupling can be used for specialized applications. 
     Research on the optimization of the performance of magnetometers that use ferromagnetic cores, and on nonlinear oscillators for electric field sensing based on ferroelectric capacitors continues to be ongoing in the prior art. However, up until now, minimization of the power demand for magnetometers (which is a key optimization feature) has been limited due to the intrinsic properties of the device, which must be “current driven” in order to ensure a proper magnetization of the core. Because the prior art magnetometers need an applied current to operate, the power budget for the prior art device cannot be reduced below a certain threshold, which can further place limits on both the size and the sensitivity of the device. 
     The availability of materials whose magnetization can be quantified not via an applied current but by using an electric field would represent a major breakthrough in the field given the inherent low power (“nearly zero power”) of this approach. This is the promise given by multiferroic materials that are therefore hysteretic in both the electric and the magnetic domain. An adequate knowledge of these materials and a suitable exploitation of their unique features could therefore lead to novel devices that can detect weak low frequency magnetic fields and that demand a negligible amount of power to operate. 
     In view of the above, it is an object of the present invention to provide a sensor that can incorporate multiferroic materials to detect electric fields or magnetic fields using the same underlying setup. Another object of the present invention is to provide a sensor that can incorporate multiferroic materials to measure magnetic fields using an applied electric field instead of an applied electric current. Still another object of the present invention is to provide a sensor that can incorporate multiferroic materials that have an extremely low power footprint to accomplish the measurement of electric fields and magnetic fields. Another object of the present invention to provide a sensor that can incorporate multiferroic materials, which can be easy to manufacture, and which can be used in a cost-efficient manner. Finally, the unique coupling configuration affords enhanced target signal resolution; in fact the resolution can be shown to improve with N the number of coupled circuit block. 
     SUMMARY OF THE INVENTION 
     A sensor for detecting magnetic fields and electric fields, and methods for detecting such fields, can include at least one Sawyer-Tower (ST) circuit that can incorporate a multiferroic capacitor. An odd number of ST circuits can be coupled together in a ring configuration, so that for each ST circuit, the output of one ST circuit can be an input to another of the ST circuits. The multiferroic capacitors can further include a multiferroic layer that can be deposed on a substrate. For the multiferroic capacitors, the deposition process can cause an inherent amount of impurities in the crystal structure of the multiferroic layer. These imperfections can lead to slight differences in the BDFO properties, which can be an undesirable result. However, the number of ST circuits to be coupled together can be chosen according to amount of impurities that can be caused by the deposition process, to “forgive” the impurities. The higher the level in impurities in the multiferroic layers, that more ST circuits are required in the ring to achieve the same ST circuit gain (sensor sensitivity) that is desired by the user. 
     The above structure and cooperation of structure can allow for a sensor that can detect magnetic fields using an applied electric field, instead of an electric current. One material that can be chosen for the multiferroic layer can be BDFO (a bulk material). Other materials that exhibit multiferroic properties could also be used. The multiferroic capacitor can further include a bias part and a sensing part. The sensing part can include a charge collector for sensing electric layer (e.g. PZT) and a magnetostrictive layer, e,g, Terfenol-D. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
       The novel features of the present invention will be best understood from the accompanying drawings, taken in conjunction with the accompanying description, in which similarly-referenced characters refer to similarly-referenced parts, and in which: 
         FIG. 1  is a graph of a theoretical modeling of the hysteresis loop of a representative multiferroic layer for the device of the present invention according to several embodiments; 
         FIG. 2  is circuit diagram of the sensor of the present invention; 
         FIG. 3  is a graph of voltage versus magnetic field amplitude, which can be used to estimate the sensitivity of the sensor of  FIG. 2 ; 
         FIG. 4  is a block diagram of a plurality of the sensors of  FIG. 2 , when in a ring configuration; 
         FIG. 5  is a graph of the voltage response of the ring configuration of  FIG. 4  for k=4; 
         FIG. 6  is the same graph as  FIG. 5 , but with k=40; 
         FIG. 7  is a top plan view of the sensor of  FIG. 2 , but when configured to simultaneously measure electrical field and magnetic fields; 
         FIG. 8  is a graph of the voltage response output signal of the ST circuit of the sensor of  FIG. 2  in absence of an external electric field; 
         FIG. 9  is the same graph as  FIG. 8 , but in the presence of an electric field; 
         FIG. 10  is a graph of voltage response the ring configuration of  FIG. 4  in the absence of an external electric field; 
         FIG. 11  is the same graph as  FIG. 10 , but in the presence of an electric field; and, 
         FIG. 12  is a block diagram of steps that can be taken to practice the methods of the present invention according to several embodiments. 
     
    
    
     DETAILED DESCRIPTION OF THE EMBODIMENTS 
     I. Introduction 
     In the prior art, a large class of materials can exhibit the ferroic orders of ferroelectricity, ferromagnetism or ferroelasticity. Materials in which at least two of the three ferroic orders exist can be called multiferroic (MF) materials. Of particular interest are magnetoelectric multiferroics, in which the ferroelectric and ferromagnetic orders not only coexist in the same material, but are also coupled to each other such that an E-field can induce a magnetization and a B-field can induce a polarization. The term magneto-electric effect, then, can be used to describe any form of cross-correlation between magnetic and dielectric properties. 
     Of the materials that exhibit magnetoelectricity, the perovskite phase of BiFeO 3  (BFO) can be unique. This is because BFO is an antiferromagnetic material, but the addition of Dy to the BFO can change the ordering from antiferromagnetic to ferromagnetic, and can result in a multiferroic material that can often exhibit magnetoelectric coupling at room temperature. In the recent prior art, a large spontaneous electric polarization (60 μCcm −2 ) in combination with a substantial thickness-dependent saturation magnetization (150 emucm −3 ), was observed above room temperature in thin films of BFO grown epitaxially on SrTiO 3  substrates. 
     To enhance the magnetoelectric coupling of BFO, several doping methods and materials can be used to modify BiFeO 3  by chemical substitutions at the Bi or Fe sites. In particular, Palkar and Prashanthi showed that a substitution of Dy at the Bi site in BiFeO 3  (BDFO) can induce ferromagnetism in BiFeO 3  without disturbing the ferroelectric behavior. The effect of an applied direct current (DC) B-field on the electric polarization was demonstrated through its effect on the ferroelectric hysteresis loop: both the slope and the saturation polarization increased with the B-field value; the converse effect was also demonstrated. At the same time, the inventors of the present invention according to several embodiments conceived of the idea of exploiting hysteretic behaviors of the multiferroic materials by using the multiferroic materials as ferroelectric capacitors for E-field sensing. To do this, a nonlinear ring oscillator underpinned by ferroelectric capacitors has been developed and characterized for its ability to measure weak electric fields. 
     II. Magnetic (B) Field Detection 
     The present invention according to several embodiments can focus on the change of the ferroelectric hysteresis with an applied B-field. The present invention can show that the magneto-electric effect of BDFO (BDFO is described as an example, but other multiferroic materials could be used) can be opportunistically exploited to detect and quantify an external B-field through the change that it produces in the electrical order. A theoretical phenomenological model of hysteretic behavior which can be based on the empirical data of BDFO material can be developed. The model can be used to quantify the response, to a B-field, of a single BDFO-based capacitive sensor. A single BDFO sensor, then, can underpins a more complex sensor comprising an odd number (N=3 elements can be used) of electrically (and unidirectional) coupled Sawyer-Tower (ST) circuits, underpinned by BDFO capacitors. The ring configuration can afford enhanced sensitivity (compared to a single ST circuit) to small changes in the ambient B-field (For purposes of this disclosure, an ST circuit can be defined as a circuit that can convert polarization to a voltage in order to quantify (measure) polarization). 
     For the present invention, and referring initially to  FIG. 1 , thin film BDFO devices can be subject to different values of DC B-fields (as perturbations to the ambient DC B-field) and the polarization (P) versus E-field (E) hysteresis loops can be obtained for each case. Graph  10  in  FIG. 1  illustrates prior art empirical data  12   a  for B=0 T, data  12   b  for B=0.01 T and data  12   c  for B=0.02 T. As shown in  FIG. 1 , good magnetoelectric coupling can be observed, i.e., the (ferroelectric) hysteresis loops characteristics change with the applied B-fields. 
     From the above, it can be seen that a “blueprint” can be developed for a detector that can exploit the dynamics of coupled nonlinear systems. The blueprint procedure can include the steps of: (1) Determining macroscopic dynamical equation for the ferroelectric component of the BDFO film, (2) Obtaining the material parameters in the (bistable) potential energy function by numerical fits to the experimental data; (3) Implementing a precisely crafted coupling of an odd (N&gt;1) number of devices using the obtained material parameters, a procedure that can be shown to yield enhanced sensitivity (over the N=1 case) to small B-field changes; and, (4) Simulating the device using a PSPICE circuit model. 
     To determine the model the polarization dynamics, standard Landau-Khalatnikov theory as known in the prior art could be used: 
                       τ   ⁢       ⅆ   P       ⅆ   t         =     -       ∂   U       ∂   P           ;       U   ⁡     (     P   ,   t     )       =         -     a   2       ⁢     P   2       +       b   4     ⁢     P   4       -     c   ⁢           ⁢     E   ⁡     (     t   -   ϕ     )       ⁢   P                 (   1   )               
Equation (1) can correspond to the particle-in-potential paradigm with U(P,t) being the (bistable) potential energy function, and the materials-based parameters (a, b) must be positive to ensure bistability. The parameters (a,b,c) can be phenomenological coefficients calculated via a fitting algorithm to ensure the good agreement between the experimentally obtained hysteresis and its theoretically model. “c” can be a materials-based parameter that quantifies the effect of the external electric field on the multiferroic polarization. “τ” can represent the device time-constant, and “φ” can be a phase-lag between the input and output (introduced by the prior art test equipment). In the absence of the driving term (i.e., E=0 which can cause “c” to be 0), the width of the hysteresis loop (between its intersections of 0) the stable minima of the potential function can occur at ±r i =±√{square root over (b/a)}, and the energy barrier height can be U 0 =a 2 /4b. The connection to the hysteresis loop is made by noting that coercivity E c =±√{square root over (4a 2 r i   2 /27)}, while the width of the hysteresis loop (between its intersections of the E axis) can be proportional to U 0 /2r. Thus, a simplex parameter identification procedure can be used to obtain a theoretical representation of the dynamics that underpin the experimentally obtained hysteresis loops by Palkar et al. in the prior art. The theoretical representations can be illustrated as hysteresis loops  14   a ,  14   b  and  14   c  in  FIG. 1 . Fig. shows the experimentally obtained hysteresis loops together with the numerical fit using Equation (1). The numerical fitting yielded φ≈π/2 for all three cases.
 
     From the above, it can be seen that the application of the Landau-Khalatnikov theory can yield the coefficients (a, b, c, τ) for the 3 values of the applied B-field (B=0 T, B=0.01 T and B=0.02 T). These parameters, listed in Table 1, change monotonically with the applied B-field and could be used to extend the model to extrapolate the effects of B-fields not specifically applied in the prior art. 
     
       
         
               
               
               
               
             
           
               
                 TABLE 1 
               
               
                   
               
               
                 Magnetic Field 
                 0T 
                 0.01T 
                 0.02T 
               
               
                   
               
             
             
               
                 a 
                 8.967 × 10 −7   
                 1.029 × 10 −2   
                 1.133 × 10 −2   
               
               
                 b 
                 3.936 × 10 −4   
                 1.421 × 10 −4   
                 1.260 × 10 −4   
               
               
                 c 
                 −2.250 × 10 −7    
                 −2.395 × 10 −7    
                 −2.459 × 10 −7    
               
               
                 T 
                 1.259 × 10 −2   
                 1.013 × 10 −2   
                 8.507 × 10 −3   
               
               
                   
               
             
          
         
       
     
     To implement the above model, polarization can be measured through a charge-to-voltage conversion obtained by using an above-defined ST circuit. Such a circuit  20  can be illustrated in  FIG. 2 . As shown in  FIG. 2 , ST circuit  20  can include a multiferroic capacitor  22 , represented by C MF . A MF capacitor can be a capacitor that includes a multiferroic layer. A multiferroic layer is a layer can be a material that is multiferroic, such as BDFO. Other multiferroic layers could be used, provided such materials exhibit the hysteresis behavior described above. In circuit  20 , multiferroic capacitor  22  can function as the active nonlinear element, C f  can represent the feedback capacitor  24  used to fix the working point of the circuit  20  and the resistor  26  R f  can be introduced to avoid drift in the circuit output. The frequency response of the ST circuit  20  can be given by 
                         V   out     ⁡     (   s   )           V   in     ⁡     (   s   )         =       -       C   MF       C   F         ⁢       sC   f       1   +     sC   f         ⁢     (       R   f       R   f       )               (   2   )               
Choosing R f &gt;&gt;(sC f ) −1 , it can be shown that
 
                       V   out     =       -       A   MF       C   f         ⁢   P       ,           (   3   )               
where A MF  can be the surface area of the electrode (the multiferroic layer) of the MF capacitor.
 
     With a sinusoidal voltage (amplitude 100 mV and frequency 10 Hz) applied at V in , the changes in the output voltage amplitude, for different values of the external B-field (B=0 T, B=0.01 T and B=0.02 T), can be shown in graph  30  in  FIG. 3 . The results displayed in graph  30  can illustrate a sensitivity of 4.2V/T for circuit  20  over the range considered (0≦B≦0.02 T). Sensitivity can improve to 4.5V/T, corresponding to a change of 33% in the output voltage amplitude for a change of 0.01 T in the external B-field, when smaller (around zero) changes of the target B-field are considered. 
     The implementation of a single ST circuit  20  shown in  FIG. 2  with the BDFO capacitor  22  can underpin a coupling scheme wherein N circuit “blocks” of circuit  20  can be coupled. Referring now to  FIG. 4  (In  FIG. 4 , N=3), in a ring configuration  40  of ST circuits  20  having cyclic boundary conditions can be shown. As used herein, the term “ring” or “ring configuration” shall mean that for each ST circuit  20  in ring configuration  40  the output of one ST circuit  20  in ring configuration  40  shall be the input of another ST circuit  20 . In  FIG. 4 , one can start with the dynamics of a single capacitor, for example multiferroic capacitor  22   a  in this coupled ring  40 , we readily see that the voltage in one element V in1  is connected, linearly, to the output voltage of the preceding element via the gain coefficient k, i.e. V in1 =kV out3 . Repeating this for each capacitor  22   a ,  22   b ,  22   c  of  FIG. 4 , the dynamics of the coupled ring  4  can be expressed as:
 
τ {dot over (P)}   i   =aP   i   =bP   i   3   +λP   i-1   :i= 1 . . . 3.  (4)
 
where, taking Equations (1) and (3) into account, the coupling strength can be expressed in terms of a (negative) gain k as
 
             λ   =     c   ⁢       A   MF       d   ⁢           ⁢     C   f         ⁢     k   .             
λ can therefore be a tunable parameter, where varying the c, A MF , (“d” is the thickness of the multiferroic layer) or C f  parameters can vary λ.
 
     It should be noted that as λ can change due to the change in the multiferroic material c parameter, which can further change as a consequence of the external magnetic field. This can be the case even if k is held constant. Equation (4) can be simulated using Simulink (other types of software known in the art for this purpose) by using the material a, b, c and τ parameters from Table 1 from the measurements taken on the MF samples. The Simulink results, therefore, can be highly representative of the real system behavior. The system coefficients can be A MF =100×100 μm 2  (i.e., the area of the multiferroic layer  22  when viewed in top plan), C f =10 nF, and d=250 nm, in order to be close to a realistic set of values for the multiferroic capacitor  22 . The application of an external target B-field corresponds to a change in the material parameters in accordance with the listings in Table 1. 
     Referring now to  FIGS. 5 and 6 , simulated voltage responses can be shown by graphs  50  and  60 , respectively. Graph  50  in  FIG. 5  is a simulation curves  52   a ,  52   b , and  52   c  of the voltage response for each circuit  20   a ,  20   b    20   c  in ring  40  in  FIG. 4  for B=0.02 T, for k=4.0, which further corresponds to λ≈λ c  (λ c  can correspond to the λ at which oscillatory behavior begins to occur). Graph  60  in  FIG. 6  can illustrate simulation curves  62   a ,  62   b  and  62   c  for the same circuits  20  in ring  40 , but for k=40, which can correspond to λ&gt;&gt;λ c . By cross-referencing  FIGS. 5 and 6 , it can be inferred that past a critical value λ c =a/2 of the control parameter λ, the solutions of Equations (4) can be oscillatory (although not sinusoidal, see  FIG. 5 ) unless λ&gt;&gt;λ. But when λ&gt;&gt;λ c , and as shown in  FIG. 6 , the solutions of Equation (4) can be both oscillatory and sinusoidal, with amplitude and period dependent on the parameters (a, b, λ). Thus, the oscillation characteristics might be useful to quantify the external magnetic field B, via the changes in their characteristics induced by changing the parameters listed in Table 1. 
     It should be noted that the “oscillations” in  FIGS. 5 and 6  can actually be switching events between the stable attractors of each element in Equation (4); this switching can be driven by the coupling term λ since, for λ≦λ c , the system is overdamped and no switching occurs. The oscillatory behavior also occurs only for odd N (N≧3). A system of the form in Equation (4) has been previously disclosed in U.S. Pat. No. 7,420,366, by Visarath In et al., for an invention entitled Coupled Nonlinear Sensor System. In et al. disclosed that the sequential (due to the unidirectional coupling) nature of the switching between stable attractors of each element P i  can allow for an analytic computation of the oscillation frequency. 
     From  FIGS. 5 and 6 , it can be seen that the amplitude of oscillations  52  can correspond roughly to one-half the separation of the saturation states. The oscillation period T i  for each element P i  can be obtained as
 
 T   i   =N √{square root over (2)}π[| f ( P   10 ) f ″( P   10 )| −1/2   +|f (− P   10 ) f ″(− P   10 )| −1/2 ],  (5)
 
where f(P 1 ) can be defined as aP 1 −bP 1   3 −λP 2m  (see Equation 4), P 2m =√{square root over ((a+λ)/b)}, P 10 =√{square root over (a/3b)} and the primes denote differentiation with respect to P 1 . The oscillations in the 3 elements P i  are separated in phase by 2π/N, and the period T can be shown to scale according to the relationship 1/√{square root over (|λ−λc|)} which can be characteristic of this class of bifurcation. The amplitude and frequency of the oscillations generated for three cases of applied DC B-field, where the case were obtained through numerical simulation of the coupled system with BDFO capacitors from Table 1 (a, b, c and τ) can be described in Table 2 below. For each value (0 T, 0.01 T, 0.02 T) of the “target” B-field in Table 2, the parameters in the system of Equation (4) can be changed in accordance with Table 1, with k=40. Table 2 can also show the values of λ corresponding to each value of the B-field.
 
     
       
         
               
               
               
               
               
             
               
               
               
               
               
             
           
               
                   
                 TABLE 2 
               
               
                   
                   
               
               
                   
                 Magnetic Field 
                 Coupling Factor 
                 Frequency 
                 Amplitude 
               
               
                   
                 B (T) 
                 λ (k = 40) 
                 (Hz) 
                 (mV pp ) 
               
               
                   
                   
               
             
             
               
                   
               
             
          
           
               
                   
                 0 
                 0.036 
                 0.4 
                 9 
               
               
                   
                 0.01 
                 0.038 
                 0.5 
                 16 
               
               
                   
                 0.02 
                 0.039 
                 0.6 
                 18 
               
               
                   
                   
               
             
          
         
       
     
     From the above, it can be seen that the change in B-field can be quantified through the change in the oscillation frequency and especially through the change in amplitude of curves  52 . In fact, the sensitivity of the output frequency and amplitude, respectively, with regard to an applied change in external B-field can be estimated as 
                       S   B   f     =       Δ   ⁢           ⁢   f       Δ   ⁢           ⁢   B         ,       S   B   A     =       Δ   ⁢           ⁢   A       Δ   ⁢           ⁢   B                 (   6   )               
Using the data above, and referring to the changes in the B-field around zero, a sensitivity of 7V/T with respect to output voltage can be obtained, which can correspond to a change in the output voltage of 77% for a change of 0.01 T in the B-field. This can confirm that the change in the oscillation amplitude might be the better indicator of a change in the applied magnetic signal (field). Referring briefly back to  FIG. 3 , and recalling that the sensitivity for a single sensor was in the 4.2V/T to 4.5V/T range, it can be inferred that the sensitivity obtained for the coupled system can be almost twice that provided by the single element. Moreover, the output amplitude and frequency can increase with the coupling gain k; therefore k can be a tunable parameter to make the system work in a desired range of frequency or amplitude.
 
     It should be reiterated that from the above, a “blueprint” for a small and cheap B-field sensor that exploits microscale and/or nanoscale materials (in this case the multiferroic BDFO) and recent ideas in the physics of coupled nonlinear systems. Using experimental (time-series) data that leads to the hysteresis behavior shown in the prior art, suitable system parameters (Table 1) can be identified by suitable fitting techniques. Once the system parameters have been determined, these parameters can be used ( FIG. 5  and Table 2) to establish a multiferroic capacitor  22 , which can be incorporated into a ST circuit to quantify (measure) a magnetic field without requiring the use of a driving current. Moreover, coupling the ST circuits in a ring  40  can yield better sensitivity over the single element configuration. 
     III. Electric (E) Field Detection 
     Based on previous work by the inventors, it can be seen that a proper arrangement of the device electrodes can allow one to obtain responses both for an external electric field and for the external magnetic field, at the same time and by using the same sensor. In fact, and referring now to  FIG. 7 , the multiferroic capacitor  22  can be divided into a bias portion  72  and a sensing portion  74 , which can be deposed onto substrate  76 . Sensing portion  74  can include a charge collector for sensing electric field, such as Lead zirconium titanate, or PZT, as well as a magnetostrictive layer such as an alloy of the formula Tb x Dy 1-x Fe 2  (x˜0.3), more commonly referred to as Terfenol-D, The addition of a charge collector can increase the sensitivity of the device to external electric fields, and can allow sensing of a target external electric field that will appear as an AC ripple over the output voltage of the ST circuit, while the magnetic field will be measured through the estimation of the peak value. In order to perform both the measures, possibly simultaneously, multiferroic devices have to be realized with the above-mentioned suitable layout, in which the top electrode is divided into a bias part and a sensing part.  FIG. 7  illustrates an example of such a layout. 
     With the above layout, the same device can be used to measure electric fields. If present, an AC electric field can modulate the output of the sensor, yielding a slow modulating wave, which can be superimposed on the magnetic field induced oscillations, and whose characteristics can depend on the external electric field (or on the external electric polarization). Measures of the magnetic field and of the electric field can be performed simultaneously, because the amplitude and the frequency of the carrier wave depend on the magnetic field, while the (modulating) wave ripple and its frequency depend on the electric field. 
       FIGS. 8 and 9  can illustrate these phenomena.  FIG. 8  depicts a representative output signal of the ST circuit  20  in absence of an external electric field, while  FIG. 9  depicts the output signal of ST circuit  20  in presence of an electric field. Simulations were performed for an external magnetic field of 0.01 T, imposing an external sinusoidal polarization of 5 mC/cm 2  amplitude and a frequency of 1 Hz. Simulations were performed to verify that the sensor output changes with an imposed target electric field. The Results are reported in Table 3, which can list value of the ripple amplitude (A r  in  FIG. 9 ) as a function of applied electrical polarization signal for a single (N=1) ST circuit  20 . 
     
       
         
               
               
               
             
               
               
               
             
           
               
                   
                 TABLE 3 
               
               
                   
                   
               
               
                   
                 Electric Polarization 
                   
               
               
                   
                 (mC/cm 2 ) 
                 Ripple Amplitude (mV) 
               
               
                   
                   
               
             
             
               
                   
               
             
          
           
               
                   
                 5 
                 29 
               
               
                   
                 10 
                 55 
               
               
                   
                 15 
                 77 
               
               
                   
                   
               
             
          
         
       
     
     Next, the behavior of ring  40  configuration of ST circuits  20  in the presence of an external electric field was simulated. Simulations were performed with Simulink software. The results are reported in  FIGS. 10 and 11 .  FIG. 10  illustrates the output signal of the coupled ring  40  in absence of an external electric field, while  FIG. 11  illustrate the output signal of ring  40  in presence of an electric field. Simulations were performed for an external magnetic field of 0.01 T, imposing an external sinusoidal polarization of 5 mC/cm 2  amplitude and a frequency of 1 Hz. 
     From  FIGS. 10 and 11 , it can again be observed that the external electric field can modulate the outputs of the multiferroic devices adding a ripple, whose amplitude and frequency depend on the amplitude and frequency of the target electric field. Table 4 reports the results obtained by imposing three different electric polarization values. As shown in Table 4, the ripple amplitude A r  can increase with the external electric field value. 
                                               TABLE 4                       Electric Polarization (mC/cm 2 )   Ripple Amplitude (mV)                                        5   22           10   42           15   60                        
Analogous to the case of a magnetic field sensor, when the target field is an electric polarization the sensitivity can be evaluated considering the following equation:
 
 S   P   Ar =100((Δ A   r   /A   rmax )/(Δ P/P   max ))  (7)
 
Where A r  is the ripple amplitude and P is the external polarization. Noting that the amplitude of the oscillations is larger in the coupled system than in the single system, we find a sensitivity of 95% in this case compared to 93% for the single circuit case Thus, the responsivity can be slightly larger in the coupled case.
 
     Referring now to  FIG. 12 , a block diagram that explains steps that can be taken to practice the methods according to several embodiments is shown and generally designated by reference character  80 . As shown method  80  can include the initial step  82  of determining the dynamical parameters for the ferroelectrical component multiferroic capacitor. Landau-Khalatnikov theory can be used to accomplish this step. Next, the methods can include the step of obtaining in material parameters of the multiferroic layer, as shown by step  84 . Once the parameters have been obtained, the methods can include the step of crafting the multiferroic capacitor for ST circuit  20 , as shown by block  86 . This step can be accomplished by epitaxially growing BDFO or similar type material on a SrTiO 3  substrate. The BDFO surface area A MF  for multiferroic capacitor  22  can be chosen according to the parameters of step  84 . The methods can further include the steps of incorporating the multiferroic capacitor  22  into an ST circuit  20  (steps  88 ) and coupling an odd number of crafted ST circuits  20  in a ring  40  configuration, as shown by step  90 . The odd number of crafted ST circuit can be chosen according to the level in impurities in the multiferroic layer of multiferroic capacitor  22  and the desired size and sensitivity of the device. 
     The use of the terms “a” and “an” and “the” and similar references in the context of describing the invention (especially in the context of the following claims) is to be construed to cover both the singular and the plural, unless otherwise indicated herein or clearly contradicted by context. The terms “comprising,” “having,” “including,” and “containing” are to be construed as open-ended terms (i.e., meaning “including, but not limited to,”) unless otherwise noted. Recitation of ranges of values herein are merely intended to serve as a shorthand method of referring individually to each separate value falling within the range, unless otherwise indicated herein, and each separate value is incorporated into the specification as if it were individually recited herein. All methods described herein can be performed in any suitable order unless otherwise indicated herein or otherwise clearly contradicted by context. The use of any and all examples, or exemplary language (e.g., “such as”) provided herein, is intended merely to better illuminate the invention and does not pose a limitation on the scope of the invention unless otherwise claimed. No language in the specification should be construed as indicating any non-claimed element as essential to the practice of the invention. 
     Preferred embodiments of this invention are described herein, including the best mode known to the inventors for carrying out the invention. Variations of those preferred embodiments may become apparent to those of ordinary skill in the art upon reading the foregoing description. The inventors expect skilled artisans to employ such variations as appropriate, and the inventors intend for the invention to be practiced otherwise than as specifically described herein. Accordingly, this invention includes all modifications and equivalents of the subject matter recited in the claims appended hereto as permitted by applicable law. Moreover, any combination of the above-described elements in all possible variations thereof is encompassed by the invention unless otherwise indicated herein or otherwise clearly contradicted by context.