<?xml version='1.0'?>
<!DOCTYPE art SYSTEM 'http://www.biomedcentral.com/xml/article.dtd'>
<art>
	<ui>1556-276X-7-572</ui>
	<ji>1556-276X</ji>
	<fm>
		<dochead>Nano Express</dochead>
		<bibl>
			<title>
				<p>Energy transfer performance of mechanical nanoresonators coupled with electromagnetic fields</p>
			</title>
			<aug>
				<au id="A1" ca="yes"><snm>Javaheri</snm><fnm>Hooman</fnm><insr iid="I1"/><email>hooman@ccs.neu.edu</email></au>
				<au id="A2"><snm>Barbiellini</snm><fnm>Bernardo</fnm><insr iid="I2"/><email>bba@neu.edu</email></au>
				<au id="A3"><snm>Noubir</snm><fnm>Guevara</fnm><insr iid="I1"/><email>noubir@ccs.neu.edu</email></au>
			</aug>
			<insg>
				<ins id="I1"><p>College of Computer and Information Science, Northeastern University, Boston, MA, USA</p></ins>
				<ins id="I2"><p>Physics Department, Northeastern University, Boston, 02115, MA, USA</p></ins>
			</insg>
			<source>Nanoscale Research Letters</source>
			<section><title><p>Regular submissions</p></title></section><issn>1556-276X</issn>
			<pubdate>2012</pubdate>
			<volume>7</volume>
			<issue>1</issue>
			<fpage>572</fpage>
			<url>http://www.nanoscalereslett.com/content/7/1/572</url>
			<xrefbib><pubidlist><pubid idtype="doi">10.1186/1556-276X-7-572</pubid><pubid idtype="pmpid">23075029</pubid></pubidlist></xrefbib>
		</bibl>
		<history><rec><date><day>23</day><month>7</month><year>2012</year></date></rec><acc><date><day>6</day><month>9</month><year>2012</year></date></acc><pub><date><day>17</day><month>10</month><year>2012</year></date></pub></history>
		<cpyrt><year>2012</year><collab>Javaheri et al.; licensee Springer.</collab><note>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</note></cpyrt>
		<kwdg>
			<kwd>Nanoresonators</kwd>
			<kwd>Magnetic nanoparticles</kwd>
			<kwd>Magnetosomes</kwd>
			<kwd>Energy transfer</kwd>
			<kwd>Nanotube radio</kwd>
		</kwdg>
		<abs>
			<sec>
				<st>
					<p>Abstract</p>
				</st>
				<sec>
					<st>
						<p/>
					</st><p>We study the energy transfer performance in electrically and magnetically coupled mechanical nanoresonators. Using the resonant scattering theory, we show that magnetically coupled resonators can achieve the same energy transfer performance as for their electrically coupled counterparts or even outperform them within the scale of interest. Magnetic and electric coupling are compared in the <it>nanotube radio</it>, a realistic example of a nano-scale mechanical resonator. The energy transfer performance is also discussed for a newly proposed bio-nanoresonator composed of magnetosomes coated with a net of protein fibers.</p>
				</sec>
			</sec>
		</abs>
	</fm>
	<bdy>
		<sec>
			<st>
				<p>Background</p>
			</st><p>Mechanical nanoresonators exhibit resonance behavior involving the mechanical vibrations of the system elements. The natural frequencies of such resonances will, generally, be in the radio frequency range. Nano-scale mechanical resonators coupled with electromagnetic fields have been receiving significant attention recently <abbrgrp>
					<abbr bid="B1">1</abbr>
					<abbr bid="B2">2</abbr>
					<abbr bid="B3">3</abbr>
				</abbrgrp>. The ability to interact with electromagnetic fields allow such resonators to be essential parts of nano-scale systems. Imaging, sensing, and targeted actuation in nano-scale are among several emerging technologies that rely on efficient energy and information transfer.</p><p>In principle, nanoresonators may couple to electromagnetic fields by the charge distributions (electric coupling) or by the magnetic moment they carry (magnetic coupling). Traditionally, the energy transfer via electric coupling has received more attention since materials are mostly transparent to the magnetic field. Also, magnetic field intensity in electromagnetic radiations is significantly smaller than the electric field. Consequently, magnetic coupling of mechanical resonators with electromagnetic radiations becomes impractical unless the size of the system significantly decreases. A desirable magnetic coupling, however, can be achieved if the coupling occurs within the near-field range <abbrgrp>
					<abbr bid="B4">4</abbr>
				</abbrgrp>. Take the example of a mechanical nanoresonator operating in a biological environment. In this case, magnetic coupling holds important advantages over electric coupling. First, magnetically coupled systems can provide more <it>selective</it> and <it>localized</it> energy transfer that is due to the fact that magnetic fields, unlike electric fields, couple weakly with non-targeted surrounding media, which are often not magnetic <abbrgrp>
					<abbr bid="B5">5</abbr>
					<abbr bid="B6">6</abbr>
				</abbrgrp>. Therefore, magnetic signals suffer from considerably less attenuations while propagating in the surrounding biological media and can drive a targeted resonator inaccessible to electric signals with the same level of energy. In addition, magnetic dipoles are normally more stable than electric dipoles and do not require significant energy from outside to maintain their state.</p><p>This work revisits the interactions of radiofrequency electromagnetic fields with mechanical nanoresonators. In particular, we are interested in the quantitative assessment of the energy transfer in such nanoresonators. We use the same methodology presented by Hamam et al. <abbrgrp>
					<abbr bid="B7">7</abbr>
				</abbrgrp> and focus on low-dissipation conditions that permit resonance. The feasibility of achieving such conditions has been demonstrated in the literature <abbrgrp>
					<abbr bid="B1">1</abbr>
					<abbr bid="B8">8</abbr>
				</abbrgrp>. The outline of this paper is as follows. We first present a general model for mechanical nanoresonators including electric and magnetic coupling mechanisms and describe the dynamics of the model. Then, we compare the resonant energy transfer performance of the resonator for electric and magnetic coupling using resonant scattering theory. Finally, we sketch a roadmap for a new nanoresonator composed of a magnetite nanoparticle embedded in a net of protein fibers.</p>
		</sec>
		<sec>
			<st>
				<p>Methods</p>
			</st><p>In general, the mechanical structure of a nanoresonator consists of an elastic cantilever beam equipped with a specialized tip, which is responsible for electromagnetic interaction, vibrating in a low-viscosity fluid such as low-pressure air. The viscoelastic model of the nanoresonator includes the coefficient of mechanical elasticity, <it>k</it>, and the dissipation coefficient, <it>D</it>. For a cylindrical beam with a spherical tip, <it>k</it>&#8201;&#8773;&#8201; <it>E</it>
				<it>I</it>
				<sub>
					<it>c</it>
				</sub>/<it>L</it>
				<sup>3</sup>, where <it>E</it>, <it>I</it>
				<sub>
					<it>c</it>
				</sub>, and <it>L</it> are the Young&#8217;s modulus, second moment of cross-section, and the length of the beam, respectively. Moreover, as shown in <abbrgrp>
					<abbr bid="B9">9</abbr>
				</abbrgrp>, the combination of intrinsic (e.g., plastic deformation and surface effects) and extrinsic (e.g., viscous forces of the surrounding fluid) dissipation mechanisms determines the value of <it>D</it>. Because the size of the nanoresonator is much smaller than the wavelength of the external field, the energy transfer is in the form of interactions between the incoming field and dipole moment of the nanoresonator&#8217;s tip. As shown in Figure <figr fid="F1">1</figr>, we consider two nanoresonators that have identical mechanical structures, yet interact with electromagnetic fields via different coupling mechanisms: electric coupling (<inline-formula>
					<m:math name="1556-276X-7-572-i1" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi mathvariant="fraktur">E</m:mi>
</m:math>
				</inline-formula>) and <it>magnetic coupling</it>(<inline-formula>
					<m:math name="1556-276X-7-572-i2" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi mathvariant="fraktur">M</m:mi>
</m:math>
				</inline-formula>). In the case of electric coupling, an alternating current electric field, <it>E </it>=<it> E</it>cos(<it>&#969;t</it>), produces a force, <it>F</it>, on an electric charge distribution, <it>q</it>, placed at the tip of the nanoresonator and causes oscillatory deflections in the cantilever. A very similar model has been discussed in the <it>nanotube radio</it>
				<abbrgrp>
					<abbr bid="B1">1</abbr>
				</abbrgrp>. For the magnetic coupling, assume the tip of the resonator is made of a ferromagnetic material such as magnetite (Fe<sub>3</sub>O<sub>4</sub>) and has a magnetic moment of <it>&#956;</it>. An external magnetic field produces two different motions on such magnetic moment: a translational force produced by the field&#8217;s gradient and a rotational torque trying to align the magnetic moment with the external field. To produce a translational force that is greater than the thermal noise on a nanoresonator, a magnetic field with ultra-high gradient is required. Producing such magnetic fields drastically increases the complexity of the system; therefore, we choose to neglect the magnetic force and only consider the rotational torque exerted by the external field on the magnetic moment of the nanoresonator&#8217;s tip. An AC magnetic field, <it>B </it>=<it> B</it>cos(<it>&#969;t</it>), generates a magnetic torque, <inline-formula>
					<m:math name="1556-276X-7-572-i3" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mrow>
   <m:msub>
      <m:mi mathvariant="script">T</m:mi>
      <m:mi mathvariant="normal">m</m:mi>
   </m:msub>
   <m:mo>=</m:mo>
   <m:mi>&#956;</m:mi>
   <m:mo>&#215;</m:mo>
   <m:mi>B</m:mi>
</m:mrow>
</m:math>
				</inline-formula>, and rotates the tip leading to oscillatory beam deflections. A similar device has been built and used for <it>ultra-sensitive magnetic resonance force microscopy</it>
				<abbrgrp>
					<abbr bid="B2">2</abbr>
					<abbr bid="B8">8</abbr>
				</abbrgrp>.</p>
			<fig id="F1"><title><p>Figure 1</p></title><caption><p>An overview of nanoresonators with electric (right) and magnetic (left) coupling</p></caption><text>
   <p><b>An overview of nanoresonators with electric (right) and magnetic (left) coupling.</b> The viscoelastic properties of the resonators are identical.</p>
</text><graphic file="1556-276X-7-572-1"/></fig><p>The dynamics of the system can be expressed by a Langevin equation for the resonator tip. After linearization for small deflections, we have </p><p>
				<display-formula id="M1">
					<m:math name="1556-276X-7-572-i4" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi>m</m:mi>
<m:mover accent="true">
   <m:mrow>
      <m:mi>x</m:mi>
   </m:mrow>
   <m:mo>&#776;</m:mo>
</m:mover>
<m:mo>+</m:mo>
<m:mi>D</m:mi>
<m:mover accent="true">
   <m:mrow>
      <m:mi>x</m:mi>
   </m:mrow>
   <m:mo>&#775;</m:mo>
</m:mover>
<m:mo>+</m:mo>
<m:mtext mathvariant="italic">kx</m:mtext>
<m:mo>=</m:mo>
<m:mi>F</m:mi>
<m:mtext>cos</m:mtext>
<m:mo>(</m:mo>
<m:mi>&#969;</m:mi>
<m:mi>t</m:mi>
<m:mo>)</m:mo>
<m:mo>+</m:mo>
<m:mi>N</m:mi>
<m:mo>(</m:mo>
<m:mi>t</m:mi>
<m:mo>)</m:mo>
<m:mo>,</m:mo>
</m:math>
				</display-formula>
			</p><p>where <it>x</it> is the displacement at the tip of the beam; <it>m</it>, the effective mass of the system; and <it>D</it>, the dissipation coefficient. <it>F </it>=<it> qE</it> for electric coupling, while <inline-formula>
					<m:math name="1556-276X-7-572-i5" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mrow>
   <m:mi>F</m:mi>
   <m:mo>=</m:mo>
   <m:mi mathvariant="script">T</m:mi>
   <m:mo>/</m:mo>
   <m:mi>L</m:mi>
   <m:mo>=</m:mo>
   <m:mi>&#956;</m:mi>
   <m:mi>B</m:mi>
   <m:mo>/</m:mo>
   <m:mi>L</m:mi>
</m:mrow>
</m:math>
				</inline-formula> for magnetic coupling. The term <it>N</it>(<it>t</it>) is a stochastic force with the correlation of &#9001;<it>N</it>(<it>t</it>)<it>N</it>(<it>t</it> + &#916;<it>t</it>)&#9002;=2<it>D</it>
				<it>k</it>
				<sub>B</sub>
				<it>T&#948;</it>(&#916;<it>t</it>), where <it>k</it>
				<sub>B</sub> and <it>T</it> are the Boltzmann constant and temperature in Kelvin, respectively. In the systems we will consider, the amount of energy stored in the resonator is well above <it>k</it>
				<sub>B</sub>
				<it>T</it>. Therefore, we can omit the stochastic term, <it>N</it>(<it>t</it>), from (1). The system&#8217;s natural frequency and the quality factor are given by <inline-formula>
					<m:math name="1556-276X-7-572-i6" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mrow>
   <m:msub>
      <m:mi>&#969;</m:mi>
      <m:mn>0</m:mn>
   </m:msub>
   <m:mo>=</m:mo>
   <m:mrow>
      <m:msqrt>
         <m:mi>k</m:mi>
         <m:mo>/</m:mo>
         <m:mi>m</m:mi>
      </m:msqrt>
   </m:mrow>
</m:mrow>
</m:math>
				</inline-formula>and <inline-formula>
					<m:math name="1556-276X-7-572-i7" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mrow>
   <m:mi>Q</m:mi>
   <m:mo>=</m:mo>
   <m:mrow>
      <m:msqrt>
         <m:mtext mathvariant="italic">km</m:mtext>
         <m:mo>/</m:mo>
         <m:mi>D</m:mi>
      </m:msqrt>
   </m:mrow>
</m:mrow>
</m:math>
				</inline-formula>, respectively. The steady state solution of the system is as follows:</p><p>
				<display-formula id="M2">
					<m:math name="1556-276X-7-572-i8" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi>x</m:mi>
<m:mo>(</m:mo>
<m:mi>&#969;</m:mi>
<m:mi>t</m:mi>
<m:mo>)</m:mo>
<m:mo>=</m:mo>
<m:msub>
   <m:mrow>
      <m:mi>x</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">m</m:mi>
   </m:mrow>
</m:msub>
<m:mtext>cos</m:mtext>
<m:mo>(</m:mo>
<m:mi>&#969;</m:mi>
<m:mi>t</m:mi>
<m:mo>+</m:mo>
<m:mi>&#966;</m:mi>
<m:mo>)</m:mo>
<m:mo>,</m:mo>
</m:math>
				</display-formula>
			</p><p>where <it>x</it>
				<sub>m</sub> and <it>&#966;</it> are the maximum deflection of the tip and the phase shift, respectively, given by </p><p>
				<display-formula id="M3">
					<m:math name="1556-276X-7-572-i9" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msub>
   <m:mrow>
      <m:mi>x</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">m</m:mi>
   </m:mrow>
</m:msub>
<m:mo>(</m:mo>
<m:mi>&#969;</m:mi>
<m:mo>)</m:mo>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mi>F</m:mi>
      <m:mo>/</m:mo>
      <m:mi>m</m:mi>
   </m:mrow>
   <m:mrow>
      <m:msqrt>
         <m:mrow>
            <m:msup>
               <m:mrow>
                  <m:mo>(</m:mo>
                  <m:msup>
                     <m:mrow>
                        <m:mi>&#969;</m:mi>
                     </m:mrow>
                     <m:mrow>
                        <m:mn>2</m:mn>
                     </m:mrow>
                  </m:msup>
                  <m:mo>&#8722;</m:mo>
                  <m:msubsup>
                     <m:mrow>
                        <m:mi>&#969;</m:mi>
                     </m:mrow>
                     <m:mrow>
                        <m:mn>0</m:mn>
                     </m:mrow>
                     <m:mrow>
                        <m:mn>2</m:mn>
                     </m:mrow>
                  </m:msubsup>
                  <m:mo>)</m:mo>
               </m:mrow>
               <m:mrow>
                  <m:mn>2</m:mn>
               </m:mrow>
            </m:msup>
            <m:mo>+</m:mo>
            <m:msup>
               <m:mrow>
                  <m:mo>(</m:mo>
                  <m:mi>&#969;</m:mi>
                  <m:msub>
                     <m:mrow>
                        <m:mi>&#969;</m:mi>
                     </m:mrow>
                     <m:mrow>
                        <m:mn>0</m:mn>
                     </m:mrow>
                  </m:msub>
                  <m:mo>/</m:mo>
                  <m:mi>Q</m:mi>
                  <m:mo>)</m:mo>
               </m:mrow>
               <m:mrow>
                  <m:mn>2</m:mn>
               </m:mrow>
            </m:msup>
         </m:mrow>
      </m:msqrt>
   </m:mrow>
</m:mfrac>
</m:math>
				</display-formula>
			</p><p/>
			<p>
				<display-formula id="M4">
					<m:math name="1556-276X-7-572-i10" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi>&#966;</m:mi>
<m:mo>(</m:mo>
<m:mi>&#969;</m:mi>
<m:mo>)</m:mo>
<m:mo>=</m:mo>
<m:mo>arctan</m:mo>
<m:mfenced separators="" open="(" close=")">
   <m:mrow>
      <m:mfrac>
         <m:mrow>
            <m:mi>&#969;</m:mi>
            <m:msub>
               <m:mrow>
                  <m:mi>&#969;</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mn>0</m:mn>
               </m:mrow>
            </m:msub>
            <m:mo>/</m:mo>
            <m:mi>Q</m:mi>
         </m:mrow>
         <m:mrow>
            <m:msup>
               <m:mrow>
                  <m:mi>&#969;</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mn>2</m:mn>
               </m:mrow>
            </m:msup>
            <m:mo>&#8722;</m:mo>
            <m:msubsup>
               <m:mrow>
                  <m:mi>&#969;</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mn>0</m:mn>
               </m:mrow>
               <m:mrow>
                  <m:mn>2</m:mn>
               </m:mrow>
            </m:msubsup>
         </m:mrow>
      </m:mfrac>
   </m:mrow>
</m:mfenced>
<m:mi>.</m:mi>
</m:math>
				</display-formula>
			</p><p>A resonance can be achieved if<inline-formula>
					<m:math name="1556-276X-7-572-i11" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mrow>
   <m:mi>D</m:mi>
   <m:mo>&lt;</m:mo>
   <m:mrow>
      <m:msqrt>
         <m:mn>2</m:mn>
         <m:mtext mathvariant="italic">km</m:mtext>
      </m:msqrt>
   </m:mrow>
</m:mrow>
</m:math>
				</inline-formula>.</p><p>The dynamics of the system can also be expressed by the following Langevin equation for rotational oscillation <abbrgrp>
					<abbr bid="B10">10</abbr>
				</abbrgrp>: </p><p>
				<display-formula id="M5">
					<m:math name="1556-276X-7-572-i12" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi>I</m:mi>
<m:mover accent="true">
   <m:mrow>
      <m:mi>&#952;</m:mi>
   </m:mrow>
   <m:mo>&#776;</m:mo>
</m:mover>
<m:mo>+</m:mo>
<m:mi>C</m:mi>
<m:mover accent="true">
   <m:mrow>
      <m:mi>&#952;</m:mi>
   </m:mrow>
   <m:mo>&#775;</m:mo>
</m:mover>
<m:mo>+</m:mo>
<m:mi>&#954;</m:mi>
<m:mi>&#952;</m:mi>
<m:mo>=</m:mo>
<m:mi mathvariant="script">T</m:mi>
<m:mtext>cos</m:mtext>
<m:mo>(</m:mo>
<m:mi>&#969;</m:mi>
<m:mi>t</m:mi>
<m:mo>)</m:mo>
<m:mo>+</m:mo>
<m:mi>&#968;</m:mi>
<m:mo>(</m:mo>
<m:mi>&#969;</m:mi>
<m:mi>t</m:mi>
<m:mo>)</m:mo>
<m:mi>.</m:mi>
</m:math>
				</display-formula>
			</p><p>Here, <it>&#952; </it>=<it> x</it>/<it>L</it> is the angular displacement, <it>I</it>&#8201;&#8773;&#8201;<it> m</it>
				<it>L</it>
				<sup>2</sup> is the system&#8217;s second moment of inertia, <it>&#954;</it>&#8201;&#8773;&#8201;<it> k</it>
				<it>L</it>
				<sup>2</sup> is the rotational spring constant of the cantilever, and <it>&#968;</it> is the stochastic torque caused by the thermal noise. For the magnetic coupling, <inline-formula>
					<m:math name="1556-276X-7-572-i13" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mrow>
   <m:mi mathvariant="script">T</m:mi>
   <m:mo>=</m:mo>
   <m:mi>&#956;</m:mi>
   <m:mi>B</m:mi>
</m:mrow>
</m:math>
				</inline-formula>, while <inline-formula>
					<m:math name="1556-276X-7-572-i14" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mrow>
   <m:mi mathvariant="script">T</m:mi>
   <m:mo>=</m:mo>
   <m:mtext mathvariant="italic">qEL</m:mtext>
</m:mrow>
</m:math>
				</inline-formula> in the case of the electric coupling.</p>
		</sec>
		<sec>
			<st>
				<p>Discussion</p>
			</st>
			<sec>
				<st>
					<p>Energy analysis</p>
				</st><p>We now consider the total energy of the oscillator </p><p>
					<display-formula id="M6">
						<m:math name="1556-276X-7-572-i15" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi>U</m:mi>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mn>1</m:mn>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:mfrac>
<m:mi>m</m:mi>
<m:msup>
   <m:mrow>
      <m:mi>v</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msup>
<m:mo>+</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mn>1</m:mn>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:mfrac>
<m:mi>k</m:mi>
<m:msup>
   <m:mrow>
      <m:mi>x</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msup>
<m:mi>.</m:mi>
</m:math>
					</display-formula>
				</p><p>When <it>&#969;</it>=<it>&#969;</it>
					<sub>0</sub>, this quantity is time-independent and is given by </p><p>
					<display-formula id="M7">
						<m:math name="1556-276X-7-572-i16" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi>U</m:mi>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:msup>
         <m:mrow>
            <m:mi>F</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
      <m:msup>
         <m:mrow>
            <m:mi>Q</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
      <m:mi>k</m:mi>
   </m:mrow>
</m:mfrac>
<m:mi>.</m:mi>
</m:math>
					</display-formula>
				</p><p>One can think of <it>U</it> as the energy capacity of the resonator. An important observation is that <it>U</it> scales with <it>Q</it>
					<sup>2</sup>.</p><p>Next, we consider the energy absorbed by the nanoresonator during the relaxation time <it>&#964;</it>=<it>Q</it>/<it>&#969;</it>
					<sub>r</sub>. This quantity can be calculated by averaging the instantaneous power absorbed by the nanoresonator, <it>P</it>, over <it>&#964;</it>. For our system, <it>P</it> can be written as inner product of incident force and velocity of the resonator </p><p>
					<display-formula id="M8">
						<m:math name="1556-276X-7-572-i17" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi>P</m:mi>
<m:mo>=</m:mo>
<m:mi>F</m:mi>
<m:mo>&#183;</m:mo>
<m:mi>v</m:mi>
<m:mo>=</m:mo>
<m:mi>F</m:mi>
<m:mtext>cos</m:mtext>
<m:mo>(</m:mo>
<m:mi>&#969;</m:mi>
<m:mi>t</m:mi>
<m:mo>)</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mi mathvariant="normal">d</m:mi>
      <m:mi>x</m:mi>
      <m:mo>(</m:mo>
      <m:mi>&#969;</m:mi>
      <m:mi>t</m:mi>
      <m:mo>)</m:mo>
   </m:mrow>
   <m:mrow>
      <m:mtext mathvariant="italic">dt</m:mtext>
   </m:mrow>
</m:mfrac>
<m:mi>.</m:mi>
</m:math>
					</display-formula>
				</p><p>After some algebra, the average absorbed power, <inline-formula>
						<m:math name="1556-276X-7-572-i18" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mover accent="true">
   <m:mrow>
      <m:mi>P</m:mi>
   </m:mrow>
   <m:mo>&#772;</m:mo>
</m:mover>
</m:math>
					</inline-formula> is given by </p><p>
					<display-formula id="M9">
						<m:math name="1556-276X-7-572-i19" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mover accent="true">
   <m:mrow>
      <m:mi>P</m:mi>
   </m:mrow>
   <m:mo>&#772;</m:mo>
</m:mover>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mi>F</m:mi>
      <m:mi>&#969;</m:mi>
      <m:msub>
         <m:mrow>
            <m:mi>x</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mi mathvariant="normal">m</m:mi>
         </m:mrow>
      </m:msub>
      <m:mtext>sin</m:mtext>
      <m:mo>(</m:mo>
      <m:mi>&#966;</m:mi>
      <m:mo>)</m:mo>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:mfrac>
<m:mi>.</m:mi>
</m:math>
					</display-formula>
				</p><p>At the resonance frequency, <inline-formula>
						<m:math name="1556-276X-7-572-i20" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msub>
   <m:mrow>
      <m:mi>x</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">m</m:mi>
   </m:mrow>
</m:msub>
<m:mo>=</m:mo>
<m:mtext mathvariant="italic">FQ</m:mtext>
<m:mo>/</m:mo>
<m:mi>m</m:mi>
<m:msubsup>
   <m:mrow>
      <m:mi>&#969;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">r</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msubsup>
</m:math>
					</inline-formula> and <it>&#966; </it>=<it> &#960;</it>/2, which gives the average absorbed power of <inline-formula>
						<m:math name="1556-276X-7-572-i21" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mover accent="true">
   <m:mrow>
      <m:mi>P</m:mi>
   </m:mrow>
   <m:mo>&#772;</m:mo>
</m:mover>
<m:mo>=</m:mo>
<m:msup>
   <m:mrow>
      <m:mi>F</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msup>
<m:mi>Q</m:mi>
<m:mo>/</m:mo>
<m:mn>2</m:mn>
<m:mi>m</m:mi>
<m:msub>
   <m:mrow>
      <m:mi>&#969;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">r</m:mi>
   </m:mrow>
</m:msub>
</m:math>
					</inline-formula>. Thus, the energy deposited in the nanoresonator during the relaxation time <it>&#964;</it> is </p><p>
					<display-formula id="M10">
						<m:math name="1556-276X-7-572-i22" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi>&#916;</m:mi>
<m:msub>
   <m:mrow>
      <m:mi>U</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">r</m:mi>
   </m:mrow>
</m:msub>
<m:mo>=</m:mo>
<m:mover accent="true">
   <m:mrow>
      <m:mi>P</m:mi>
   </m:mrow>
   <m:mo>&#772;</m:mo>
</m:mover>
<m:mi>&#964;</m:mi>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:msup>
         <m:mrow>
            <m:mi>F</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
      <m:msup>
         <m:mrow>
            <m:mi>Q</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
      <m:mi>m</m:mi>
      <m:msubsup>
         <m:mrow>
            <m:mi>&#969;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mi mathvariant="normal">r</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msubsup>
   </m:mrow>
</m:mfrac>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:msup>
         <m:mrow>
            <m:mi>F</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
      <m:msup>
         <m:mrow>
            <m:mi>Q</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
      <m:mi>k</m:mi>
   </m:mrow>
</m:mfrac>
<m:mi>.</m:mi>
</m:math>
					</display-formula>
				</p><p>Note that the energy absorbed by the resonator over the relaxation time matches the resonator energy capacity. In general, the calculation of the force (or torque) exerted on the nanoresonator through electromagnetic coupling is not straightforward. As an alternative approach, one can use scattering theory <abbrgrp>
						<abbr bid="B7">7</abbr>
						<abbr bid="B11">11</abbr>
					</abbrgrp>, which allows to work with fluxes instead of forces, to estimate the energy deposited on the resonant system. The two approaches are equivalent since our theoretical model is solely based on dipole-dipole interactions. In the next section, we will use this more convenient method to study the resonant energy transfer.</p>
			</sec>
			<sec>
				<st>
					<p>Resonant scattering analysis</p>
				</st><p>The coupling between external fields and the nanoresonator consists of an absorption and a scattering process. According to the scattering theory, the power absorbed by the resonant system equals to <it>P</it>
					<sub>a</sub> = &#934;<it>&#963;</it>
					<sub>a</sub>, where &#934; is the incident electromagnetic power flux, and <it>&#963;</it>
					<sub>
						<it>a</it>
					</sub> is the absorption cross-section given by <abbrgrp>
						<abbr bid="B7">7</abbr>
						<abbr bid="B11">11</abbr>
					</abbrgrp>
				</p><p>
					<display-formula id="M11">
						<m:math name="1556-276X-7-572-i23" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msub>
   <m:mrow>
      <m:mi>&#963;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">a</m:mi>
   </m:mrow>
</m:msub>
<m:mo>(</m:mo>
<m:mi>&#969;</m:mi>
<m:mo>)</m:mo>
<m:mo>=</m:mo>
<m:mn>12</m:mn>
<m:mi>&#960;</m:mi>
<m:msup>
   <m:mrow>
      <m:mo>(</m:mo>
      <m:mfrac>
         <m:mrow>
            <m:mi>c</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mi>&#969;</m:mi>
         </m:mrow>
      </m:mfrac>
      <m:mo>)</m:mo>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msup>
<m:mfrac>
   <m:mrow>
      <m:msub>
         <m:mrow>
            <m:mi mathvariant="normal">&#915;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mi mathvariant="normal">a</m:mi>
         </m:mrow>
      </m:msub>
      <m:msub>
         <m:mrow>
            <m:mi mathvariant="normal">&#915;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mi mathvariant="normal">s</m:mi>
         </m:mrow>
      </m:msub>
   </m:mrow>
   <m:mrow>
      <m:msup>
         <m:mrow>
            <m:mo>(</m:mo>
            <m:mi>&#969;</m:mi>
            <m:mo>&#8722;</m:mo>
            <m:msub>
               <m:mrow>
                  <m:mi>&#969;</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mi mathvariant="normal">r</m:mi>
               </m:mrow>
            </m:msub>
            <m:mo>)</m:mo>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
      <m:mo>+</m:mo>
      <m:msup>
         <m:mrow>
            <m:mo>(</m:mo>
            <m:msub>
               <m:mrow>
                  <m:mi mathvariant="normal">&#915;</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mi mathvariant="normal">a</m:mi>
               </m:mrow>
            </m:msub>
            <m:mo>+</m:mo>
            <m:msub>
               <m:mrow>
                  <m:mi mathvariant="normal">&#915;</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mi mathvariant="normal">s</m:mi>
               </m:mrow>
            </m:msub>
            <m:mo>)</m:mo>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
      <m:mo>/</m:mo>
      <m:mn>4</m:mn>
   </m:mrow>
</m:mfrac>
<m:mi>.</m:mi>
</m:math>
					</display-formula>
				</p><p>Here, <it>c</it> is the speed of light; &#915;<sub>a</sub>, the absorption width; and &#915;<sub>s</sub>, the scattering width. The widths are the ratio of the power loss to the characteristic energy of the corresponding process. For process <it>i</it>, &#915;<sub>
						<it>i</it>
					</sub> = 1/<it>&#964;</it>
					<sub>
						<it>i</it>
					</sub> =<it> &#969;</it>/<it>Q</it>
					<sub>
						<it>i</it>
					</sub>, where <it>&#964;</it>
					<sub>
						<it>i</it>
					</sub> and <it>Q</it>
					<sub>
						<it>i</it>
					</sub> represent the relaxation time and the quality factor, respectively. The total energy absorbed by the resonant system during the resonant process is given by the following: </p><p>
					<display-formula id="M12">
						<m:math name="1556-276X-7-572-i24" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi>&#916;</m:mi>
<m:mi>U</m:mi>
<m:mo>(</m:mo>
<m:mi>&#969;</m:mi>
<m:mo>)</m:mo>
<m:mo>=</m:mo>
<m:msub>
   <m:mrow>
      <m:mi>P</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">a</m:mi>
   </m:mrow>
</m:msub>
<m:mo>(</m:mo>
<m:mi>&#969;</m:mi>
<m:mo>)</m:mo>
<m:mo>&#215;</m:mo>
<m:msub>
   <m:mrow>
      <m:mi>&#964;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">a</m:mi>
   </m:mrow>
</m:msub>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mi mathvariant="normal">&#934;</m:mi>
      <m:msub>
         <m:mrow>
            <m:mi>&#963;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mi mathvariant="normal">a</m:mi>
         </m:mrow>
      </m:msub>
      <m:mo>(</m:mo>
      <m:mi>&#969;</m:mi>
      <m:mo>)</m:mo>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">&#915;</m:mi>
   </m:mrow>
</m:mfrac>
<m:mo>,</m:mo>
</m:math>
					</display-formula>
				</p><p>where &#915; = &#915;<sub>a</sub> + &#915;<sub>s</sub> is the total width of the system. For nano-scale systems of interest, &#915;&#8201;&#8776;&#8201;&#915;<sub>a</sub> because &#915;<sub>s</sub>&#8201;&#8810;&#8201;&#915;<sub>a</sub>. The maximal energy transfer occurs at the resonant frequency and can be written as </p><p>
					<display-formula id="M13">
						<m:math name="1556-276X-7-572-i25" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msub>
   <m:mrow>
      <m:mi>&#916;</m:mi>
      <m:mi>U</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">r</m:mi>
   </m:mrow>
</m:msub>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mn>48</m:mn>
      <m:mi mathvariant="normal">&#934;</m:mi>
      <m:msup>
         <m:mrow>
            <m:mi>c</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:msubsup>
         <m:mrow>
            <m:mi>&#969;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mi mathvariant="normal">r</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msubsup>
   </m:mrow>
</m:mfrac>
<m:mfenced separators="" open="(" close=")">
   <m:mrow>
      <m:mfrac>
         <m:mrow>
            <m:msub>
               <m:mrow>
                  <m:mi mathvariant="normal">&#915;</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mi mathvariant="normal">s</m:mi>
               </m:mrow>
            </m:msub>
         </m:mrow>
         <m:mrow>
            <m:msubsup>
               <m:mrow>
                  <m:mi mathvariant="normal">&#915;</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mi mathvariant="normal">a</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mn>2</m:mn>
               </m:mrow>
            </m:msubsup>
         </m:mrow>
      </m:mfrac>
   </m:mrow>
</m:mfenced>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mn>48</m:mn>
      <m:mi mathvariant="normal">&#934;</m:mi>
      <m:msup>
         <m:mrow>
            <m:mi>c</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:msubsup>
         <m:mrow>
            <m:mi>&#969;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mi mathvariant="normal">r</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>4</m:mn>
         </m:mrow>
      </m:msubsup>
   </m:mrow>
</m:mfrac>
<m:mspace width="1em"/>
<m:msub>
   <m:mrow>
      <m:mi mathvariant="normal">&#915;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">s</m:mi>
   </m:mrow>
</m:msub>
<m:msubsup>
   <m:mrow>
      <m:mi>Q</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">a</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msubsup>
<m:mi>.</m:mi>
</m:math>
					</display-formula>
				</p><p>
					<it>Q</it>
					<sub>a</sub> is obtained from the steady state solution.</p><p>By definition, the width of the scattering process is equal to the inverse of the decay time of radiating dipole given by the following: </p><p>
					<display-formula id="M14">
						<m:math name="1556-276X-7-572-i26" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msub>
   <m:mrow>
      <m:mi mathvariant="normal">&#915;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">s</m:mi>
   </m:mrow>
</m:msub>
<m:mo>=</m:mo>
<m:mo>&#8722;</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mi mathvariant="normal">d</m:mi>
      <m:mtext>ln</m:mtext>
      <m:mi>U</m:mi>
      <m:mo>(</m:mo>
      <m:mi>t</m:mi>
      <m:mo>)</m:mo>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">d</m:mi>
      <m:mi>t</m:mi>
   </m:mrow>
</m:mfrac>
<m:mo>=</m:mo>
<m:mo>&#8722;</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mi mathvariant="normal">d</m:mi>
      <m:mi>U</m:mi>
      <m:mo>/</m:mo>
      <m:mi mathvariant="normal">d</m:mi>
      <m:mi>t</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi>U</m:mi>
   </m:mrow>
</m:mfrac>
</m:math>
					</display-formula>
				</p><p>where <it>U</it> is the energy of the resonator, and <it>P</it>
					<sub>r</sub> = d<it>U</it>/d<it>t</it> is the radiative power of the resonator&#8217;s dipole. For the electric dipole model (<inline-formula>
						<m:math name="1556-276X-7-572-i27" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi mathvariant="fraktur">E</m:mi>
</m:math>
					</inline-formula>), we have </p><p>
					<display-formula id="M15 ">
						<m:math name="1556-276X-7-572-i28" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msubsup>
   <m:mrow>
      <m:mi>P</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">r</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="fraktur">E</m:mi>
   </m:mrow>
</m:msubsup>
<m:mo>=</m:mo>
<m:mo>&#8722;</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mn>1</m:mn>
   </m:mrow>
   <m:mrow>
      <m:mn>4</m:mn>
      <m:mi>&#960;</m:mi>
      <m:msub>
         <m:mrow>
            <m:mi>&#949;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>0</m:mn>
         </m:mrow>
      </m:msub>
   </m:mrow>
</m:mfrac>
<m:mfrac>
   <m:mrow>
      <m:msubsup>
         <m:mrow>
            <m:mi>p</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>0</m:mn>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msubsup>
      <m:msup>
         <m:mrow>
            <m:mi>&#969;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>4</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:mn>3</m:mn>
      <m:msup>
         <m:mrow>
            <m:mi>c</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>3</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
</m:mfrac>
</m:math>
					</display-formula>
				</p><p/>
				<p>
					<display-formula id="M16">
						<m:math name="1556-276X-7-572-i29" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msup>
   <m:mrow>
      <m:mi>U</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="fraktur">E</m:mi>
   </m:mrow>
</m:msup>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mn>1</m:mn>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:mfrac>
<m:mi>k</m:mi>
<m:msubsup>
   <m:mrow>
      <m:mi>x</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">m</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msubsup>
<m:mo>,</m:mo>
</m:math>
					</display-formula>
				</p><p>where <it>p</it>
					<sub>0</sub> =<it> q</it>
					<it>x</it>
					<sub>m</sub> is the maximal amplitude for the electric moment of the resonator. Thus, one obtains the following scattering width for the system: </p><p>
					<display-formula id="M17">
						<m:math name="1556-276X-7-572-i30" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msubsup>
   <m:mrow>
      <m:mi mathvariant="normal">&#915;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">s</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="fraktur">E</m:mi>
   </m:mrow>
</m:msubsup>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:msup>
         <m:mrow>
            <m:mi>q</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:mn>4</m:mn>
      <m:mi>&#960;</m:mi>
      <m:msub>
         <m:mrow>
            <m:mi>&#949;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>0</m:mn>
         </m:mrow>
      </m:msub>
   </m:mrow>
</m:mfrac>
<m:mfrac>
   <m:mrow>
      <m:mn>2</m:mn>
      <m:msup>
         <m:mrow>
            <m:mi>&#969;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>4</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:mn>3</m:mn>
      <m:msup>
         <m:mrow>
            <m:mi>c</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>3</m:mn>
         </m:mrow>
      </m:msup>
      <m:mi>k</m:mi>
   </m:mrow>
</m:mfrac>
<m:mi>.</m:mi>
</m:math>
					</display-formula>
				</p><p>Replacing (17) in (13) results in </p><p>
					<display-formula id="M18">
						<m:math name="1556-276X-7-572-i31" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msubsup>
   <m:mrow>
      <m:mi>&#916;</m:mi>
      <m:mi>U</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">r</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="fraktur">E</m:mi>
   </m:mrow>
</m:msubsup>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:msup>
         <m:mrow>
            <m:mi>q</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:mn>4</m:mn>
      <m:mi>&#960;</m:mi>
      <m:msub>
         <m:mrow>
            <m:mi>&#949;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>0</m:mn>
         </m:mrow>
      </m:msub>
   </m:mrow>
</m:mfrac>
<m:mfrac>
   <m:mrow>
      <m:mn>32</m:mn>
      <m:mi mathvariant="normal">&#934;</m:mi>
      <m:msup>
         <m:mrow>
            <m:msub>
               <m:mrow>
                  <m:mi>Q</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mi mathvariant="normal">a</m:mi>
               </m:mrow>
            </m:msub>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:mtext mathvariant="italic">ck</m:mtext>
   </m:mrow>
</m:mfrac>
<m:mo>,</m:mo>
</m:math>
					</display-formula>
				</p><p>which confirms that the energy deposited scales as <it>Q</it>
					<sup>2</sup>/<it>k</it>.</p><p>Using classical electrodynamics <abbrgrp>
						<abbr bid="B12">12</abbr>
					</abbrgrp>, one can show that the the radiative power of an oscillating magnetic dipole of moment <it>&#956;</it>
					<sub>eff</sub> is given by replacing <it>p</it>
					<sub>0</sub> by <it>&#956;</it>
					<sub>eff</sub>/<it>c</it> in 15. In the case of the spherical MNP shown in Figure <figr fid="F1">1</figr>, if <it>&#952;</it>
					<sub>m </sub> is the maximum angular deflection, then the oscillating part of the magnetic dipole of moment is <inline-formula>
						<m:math name="1556-276X-7-572-i32" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msub>
   <m:mrow>
      <m:mi>&#956;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mtext>eff</m:mtext>
   </m:mrow>
</m:msub>
<m:mo>=</m:mo>
<m:mi>&#956;</m:mi>
<m:msub>
   <m:mrow>
      <m:mtext>sin</m:mtext>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msub>
<m:mo>(</m:mo>
<m:msub>
   <m:mrow>
      <m:mi>&#952;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">m</m:mi>
   </m:mrow>
</m:msub>
<m:mo>)</m:mo>
</m:math>
					</inline-formula>. Therefore, for the magnetic coupling model (<inline-formula>
						<m:math name="1556-276X-7-572-i33" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi mathvariant="fraktur">M</m:mi>
</m:math>
					</inline-formula>), one can follow the same derivation as in the electric coupling scheme (<inline-formula>
						<m:math name="1556-276X-7-572-i34" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi mathvariant="fraktur">E</m:mi>
</m:math>
					</inline-formula>) and obtain </p><p>
					<display-formula id="M19">
						<m:math name="1556-276X-7-572-i35" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msubsup>
   <m:mrow>
      <m:mi>P</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">r</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="fraktur">M</m:mi>
   </m:mrow>
</m:msubsup>
<m:mo>=</m:mo>
<m:mo>&#8722;</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mn>1</m:mn>
   </m:mrow>
   <m:mrow>
      <m:mn>4</m:mn>
      <m:mi>&#960;</m:mi>
      <m:msub>
         <m:mrow>
            <m:mi>&#949;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>0</m:mn>
         </m:mrow>
      </m:msub>
   </m:mrow>
</m:mfrac>
<m:mfrac>
   <m:mrow>
      <m:msup>
         <m:mrow>
            <m:mi>&#956;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
      <m:msup>
         <m:mrow>
            <m:mi>&#969;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>4</m:mn>
         </m:mrow>
      </m:msup>
      <m:msubsup>
         <m:mrow>
            <m:mi>&#952;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mi mathvariant="normal">m</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msubsup>
   </m:mrow>
   <m:mrow>
      <m:mn>3</m:mn>
      <m:msup>
         <m:mrow>
            <m:mi>c</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>5</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
</m:mfrac>
<m:mo>,</m:mo>
</m:math>
					</display-formula>
				</p><p/>
				<p>
					<display-formula id="M20">
						<m:math name="1556-276X-7-572-i36" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msup>
   <m:mrow>
      <m:mi>U</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="fraktur">M</m:mi>
   </m:mrow>
</m:msup>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mn>1</m:mn>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:mfrac>
<m:mi>&#954;</m:mi>
<m:msubsup>
   <m:mrow>
      <m:mi>&#952;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi>m</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msubsup>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:mn>1</m:mn>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:mfrac>
<m:mi>k</m:mi>
<m:msup>
   <m:mrow>
      <m:mi>L</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msup>
<m:msubsup>
   <m:mrow>
      <m:mi>&#952;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi>m</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msubsup>
<m:mi>.</m:mi>
</m:math>
					</display-formula>
				</p><p>Therefore, the corresponding radiation width and deposited energy are </p><p>
					<display-formula id="M21">
						<m:math name="1556-276X-7-572-i37" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msubsup>
   <m:mrow>
      <m:mi mathvariant="normal">&#915;</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">s</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="fraktur">M</m:mi>
   </m:mrow>
</m:msubsup>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:msup>
         <m:mrow>
            <m:mi>&#956;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:mn>4</m:mn>
      <m:mi>&#960;</m:mi>
      <m:msub>
         <m:mrow>
            <m:mi>&#949;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>0</m:mn>
         </m:mrow>
      </m:msub>
   </m:mrow>
</m:mfrac>
<m:mfrac>
   <m:mrow>
      <m:mn>2</m:mn>
      <m:msup>
         <m:mrow>
            <m:mi>&#969;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>4</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:mn>3</m:mn>
      <m:msup>
         <m:mrow>
            <m:mi>c</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>3</m:mn>
         </m:mrow>
      </m:msup>
      <m:mi>k</m:mi>
      <m:msup>
         <m:mrow>
            <m:mo>(</m:mo>
            <m:mtext mathvariant="italic">Lc</m:mtext>
            <m:mo>)</m:mo>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
</m:mfrac>
</m:math>
					</display-formula>
				</p><p/>
				<p>
					<display-formula id="M22">
						<m:math name="1556-276X-7-572-i38" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:msubsup>
   <m:mrow>
      <m:mi>&#916;</m:mi>
      <m:mi>U</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="normal">r</m:mi>
   </m:mrow>
   <m:mrow>
      <m:mi mathvariant="fraktur">M</m:mi>
   </m:mrow>
</m:msubsup>
<m:mo>=</m:mo>
<m:mfrac>
   <m:mrow>
      <m:msup>
         <m:mrow>
            <m:mi>&#956;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msup>
   </m:mrow>
   <m:mrow>
      <m:mn>4</m:mn>
      <m:mi>&#960;</m:mi>
      <m:msub>
         <m:mrow>
            <m:mi>&#949;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>0</m:mn>
         </m:mrow>
      </m:msub>
   </m:mrow>
</m:mfrac>
<m:msup>
   <m:mrow>
      <m:mfenced separators="" open="[" close="]">
         <m:mrow>
            <m:mfrac>
               <m:mrow>
                  <m:mn>1</m:mn>
               </m:mrow>
               <m:mrow>
                  <m:mtext mathvariant="italic">Lc</m:mtext>
               </m:mrow>
            </m:mfrac>
         </m:mrow>
      </m:mfenced>
   </m:mrow>
   <m:mrow>
      <m:mn>2</m:mn>
   </m:mrow>
</m:msup>
<m:mfrac>
   <m:mrow>
      <m:mn>32</m:mn>
      <m:mi mathvariant="normal">&#934;</m:mi>
      <m:msubsup>
         <m:mrow>
            <m:mi>Q</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mi mathvariant="normal">a</m:mi>
         </m:mrow>
         <m:mrow>
            <m:mn>2</m:mn>
         </m:mrow>
      </m:msubsup>
   </m:mrow>
   <m:mrow>
      <m:mtext mathvariant="italic">ck</m:mtext>
   </m:mrow>
</m:mfrac>
<m:mi>.</m:mi>
</m:math>
					</display-formula>
				</p><p>Note that <it>Q</it>
					<sub>a</sub> and <it>k</it> in Equations (18) and (22) only depend on the viscoelastic structure of the resonator and are independent from the coupling type (magnetic or electric). Given a similar viscoelastic structure, the energy absorption value for electric and magnetic coupling will be comparable if <it>&#956;</it>/<it>Lc</it>&#8201;&#8776;&#8201;<it>q</it>. For nano-scale systems of interest, the condition <it>Lc </it>&lt; 1,000 m<sup>2</sup>/s normally holds. By comparing 22 to 10, it is possible to derive an expression for the average magnetic force experienced by the nanoresonator over the resonance relaxation time, which is given by </p><p>
					<display-formula id="M23">
						<m:math name="1556-276X-7-572-i39" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mover accent="true">
   <m:mrow>
      <m:mi>F</m:mi>
   </m:mrow>
   <m:mo>&#772;</m:mo>
</m:mover>
<m:mo>=</m:mo>
<m:msqrt>
   <m:mrow>
      <m:mfrac>
         <m:mrow>
            <m:msup>
               <m:mrow>
                  <m:mi>&#956;</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mn>2</m:mn>
               </m:mrow>
            </m:msup>
         </m:mrow>
         <m:mrow>
            <m:mn>4</m:mn>
            <m:mi>&#960;</m:mi>
            <m:msub>
               <m:mrow>
                  <m:mi>&#949;</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mn>0</m:mn>
               </m:mrow>
            </m:msub>
         </m:mrow>
      </m:mfrac>
      <m:mfrac>
         <m:mrow>
            <m:mn>32</m:mn>
            <m:mi mathvariant="normal">&#934;</m:mi>
         </m:mrow>
         <m:mrow>
            <m:msup>
               <m:mrow>
                  <m:mi>c</m:mi>
               </m:mrow>
               <m:mrow>
                  <m:mn>3</m:mn>
               </m:mrow>
            </m:msup>
         </m:mrow>
      </m:mfrac>
   </m:mrow>
</m:msqrt>
<m:mi>.</m:mi>
</m:math>
					</display-formula>
				</p>
			</sec>
			<sec>
				<st>
					<p>Applications</p>
				</st><p>Having discussed the mechanical dynamics of the nanoresonator as well as the theoretical formulation for energy transfer performance of different coupling mechanism, we apply our analysis to a possible nanoresonator sketched in Figure <figr fid="F2">2</figr>, and we also discuss the feasibility of a bio-nanoresonator composed of protein coated Fe<sub>3</sub>O<sub>4</sub> nanoparticles.</p>
				<fig id="F2"><title><p>Figure 2</p></title><caption><p>Magnetosome arrangement in magnetotactic bacteria</p></caption><text>
   <p><b>Magnetosome arrangement in magnetotactic bacteria.</b> The magnified part shows how elastic protein fibers embed magnetite (<it>F</it><it>e</it><sub>3</sub><it>O</it><sub>4</sub>) crystals in the cytoskeleton. Interaction of the magnetic dipole of the crystal with external fields within its viscoelastic environment can be analyzed by our presented theoretical model as a torsional nanoresonator shown on the right hand side. Magnetic torque rotates the MNP around its center of mass. The rotational spring constant is given by <it>&#954; </it>=<it> kR</it>, where <it>k</it> is the aggregate rigidity of the connecting protein fibers, and <it>R</it> is the radius of the MNP. Since the Reynold number of the MNP is very small, the drag forces are given by Stokes&#8217; law. Therefore, the rotational damping coefficient is <it>C </it>= 6<it>&#960;&#8201;&#951;</it><it>R</it><sup>3</sup>, where <it>&#951; </it> is the viscosity of the surrounding fluid.</p>
</text><graphic file="1556-276X-7-572-2"/></fig><p>In our first example, we compare the energy transfer performance of the magnetic and electric coupling in the nanotube radio, a realistic example of a mechanical nanoresonator <abbrgrp>
						<abbr bid="B1">1</abbr>
					</abbrgrp>. We replace the electric dipole of the nanotube tip with a magnetic dipole in the form of spherical magnetite nanoparticle. According to the original study, a nanotube radio built from a cylindrical carbon nanotube of length <it>L</it>&#8201;&#8776;&#8201;1<it> &#956;</it>m holding a net charge of <it>q </it>= 200<it> e</it>
					<sup>&#8722;</sup> absorbs an amount of energy enough to detect radio signals from the electromagnetic radiation. To achieve the same amount of energy deposit, the magnetic moment of the replacement tip should be in the order of <it>&#956;</it>&#8201;&#8776;&#8201;<it>qLc </it>= 9.6&#8201;&#215;&#8201;10<sup>&#8722;15</sup> Am<sup>2</sup>, which can be obtained by placing a magnetite nanoparticle of radius <it>R</it> approximately 160 nm.</p><p>Another interesting application is the possibility of transmitting energy to magnetic nanoparticles in the biological setting. Biogenic magnetite nanoparticles called magnetosomes, first discovered in magnetotactic bacteria <abbrgrp>
						<abbr bid="B13">13</abbr>
					</abbrgrp>, are also found in the brain of many animals and are believed to participate in determining the orientation in several species such as migratory birds <abbrgrp>
						<abbr bid="B14">14</abbr>
					</abbrgrp>. Interestingly, magnetosomes consist of magnetite particles of radius 50 to 100 nm and are embedded in the cytoskeleton bound to a viscoelastic system formed by a net of protein fibers. Because magnetic nanoparticles of such size are single domain with high coercivity <abbrgrp>
						<abbr bid="B15">15</abbr>
						<abbr bid="B16">16</abbr>
					</abbrgrp>, the magnetosome can be represented as a torsional nanoresonator with magnetic coupling (see Figure <figr fid="F2">2</figr>). According to Winklhofer and Kirschvink <abbrgrp>
						<abbr bid="B17">17</abbr>
					</abbrgrp>, the rigidity of the cytoskeleton can be estimated by <it>&#954; </it>= 100<it>k</it>
					<sub>B</sub>
					<it>T</it>/ Rad per connecting filament. Thus, for magnetite particles with a density of <it>&#961; </it>= 5,200 Kg/m<sup>3</sup>, a radius of <it>R </it>= 100nm, and the number of the connecting filaments ranging from 1 to 1,000, the natural frequency of the oscillator will fall between 2 and 66 MHz. A resonance is in principle not possible if we adopt the standard viscosity of the cytoplasm <abbrgrp>
						<abbr bid="B18">18</abbr>
					</abbrgrp>. However, in a carefully engineered synthetic system, one could lower the drag forces in order to achieve resonance and higher quality factor up to <it>Q </it>= s100. Figure <figr fid="F3">3</figr> shows possible quality factor values for a nanoparticle of radius 100 nm, assuming that the viscosity and elastic constants could be controlled.</p>
				<fig id="F3"><title><p>Figure 3</p></title><caption><p>Quality factor of the resonance for reasonable range of values (color online)</p></caption><text>
   <p><b>Quality factor of the resonance for reasonable range of values (color online).</b> For the environment viscosity and the rotational spring constant of elastic environment (in terms of <it>k</it><sub>B</sub><it>T</it>). We assume the design includes a magnetite nanoparticle of radius 100 nm. Note that resonance is possible in the region above <it>Q </it>= 0.5line. It is shown how resonance of a given quality can be achieved in lower frequency by reducing the viscosity experienced by the resonator.</p>
</text><graphic file="1556-276X-7-572-3"/></fig><p>In order to achieve higher quality factor, the nanoresonator should experience smaller viscous resistance. For instance, one can reduce drag forces in the system by coating the magnetite nanoparticle with hydrophobic proteins or lipids. In this case, the hydrophobic coat acts as a lubricant <abbrgrp>
						<abbr bid="B19">19</abbr>
						<abbr bid="B20">20</abbr>
					</abbrgrp>. In a more sophisticated design, a multi-layer shell of hydrophobic proteins may be used to engulf the nanoresonator and repel water molecules <abbrgrp>
						<abbr bid="B21">21</abbr>
					</abbrgrp>. In order to aggressively reduce the rotational friction, the nanoresonator could be packed in an inorganic shell that completely excludes the system from the cytoplasm. The elastic protein fibers may be replaced by synthetic nanowires or nanotubes with carefully designed rigidity. For example, del Barco et al. <abbrgrp>
						<abbr bid="B22">22</abbr>
					</abbrgrp> have demonstrated the possibility to have free rotation of magnetic nanoparticle embedded in a solid matrix.</p><p>Assuming that high quality factor, <it>Q </it>= 100, can be achieved, one finds that an AC magnetic field of intensity <it>B </it>= 3.5mT, generating an electromagnetic flux of 10 W/m<sup>2</sup>, deposits a significant amount of energy &#916;<it>U </it>= 2,500<it>k</it>
					<sub>B</sub>
					<it>T</it> into the system over the resonance relaxation time. Since this field intensity is well below the coercivity field of the nanoparticle <abbrgrp>
						<abbr bid="B15">15</abbr>
						<abbr bid="B16">16</abbr>
					</abbrgrp>, we neglect the energy losses via magnetic reversal. If this energy was entirely manifested as heat, the temperature of the magnetosome would be increased by 0.5&#176;C during the relaxation time <it>&#964; </it>= 0.1<it>&#956;</it>s. As shown in Figure <figr fid="F3">3</figr>, <it>Q </it>= 10corresponds to <it>&#969;</it>
					<sub>0</sub> = 66MHz, in air (<it>&#951; </it>= 10<sup>&#8722;5</sup>Pa.s), while the same quality factor can be achieved at frequencies as low as about 1 MHz if the viscosity can be reduced by a factor of 100 compared to air. Magnetically coupled mechanical nanoresonators with high quality factor show good energy transfer performance while being tunable and may be useful in frequency selective heat production in the biological environment. The important contribution of mechanical motions in magnetic hypothermia has been experimentally shown in <abbrgrp>
						<abbr bid="B23">23</abbr>
					</abbrgrp>; however, these applications have not yet benefited from the resonant energy transfer since their quality factors are well below one.</p>
			</sec>
		</sec>
		<sec>
			<st>
				<p>Conclusion</p>
			</st><p>In conclusion, we have shown that carefully engineered magnetically coupled nanoresonators can match the energy transfer performance of its electrically coupled counterpart, while providing a more selective and robust interaction in biological environments. We have used a unifying framework of resonant energy transfer for electrically coupled and magnetically coupled mechanical nanoresonators and compared the performance for the two couplings. Our analysis suggests that if the interacting electric dipole of a small electrically coupled resonator is replaced by a magnetic dipole, a comparable amount of energy can still be deposited on the system. We have considered the example of nanotube radio, and we have shown that the strength of electromagnetic coupling remains the same using a magnetite nanoparticle of radius 160 nm instead of the charged tip. We have proposed a new resonator composed of magnetosomes embedded in a net of protein fibers and analyzed its energy transfer performance. We have discussed possible pathways to further improve the quality factor of the resonator. While this article focuses on quantitative aspect of energy transfer, our work also opens up new interesting questions on how to use efficient energy channels to transmit information to a nano-scale device or organism. Characterizing the transmission of information and the channel capacity <abbrgrp>
					<abbr bid="B24">24</abbr>
				</abbrgrp> will be discussed in future studies.</p>
		</sec>
		<sec>
			<st>
				<p>Competing interests</p>
			</st><p>The authors declare that they have no competing interests.</p>
		</sec>
		<sec>
			<st>
				<p>Authors&#8217; contributions</p>
			</st><p>HJ designed the theoretical model, carried out the analysis, and wrote the manuscript. BB participated in developing the model and helped in conceptualizing the ideas. GN defined the research theme and participated in writing the manuscript. All authors read and approved the final manuscript.</p>
		</sec>
	</bdy>
	<bm>
		<ack>
			<sec>
				<st>
					<p>Acknowledgements</p>
				</st><p>We would like to thank Professor Mark Dykman for his important suggestions and enlightening comments. We also thank the participants of <it>NSF Workshop on Biologically-enabled Wireless Networks</it> for stimulating discussions. This work is supported in part by the National Science Foundation under grant number NSF CNS-1051240 and the US Department of Energy, Office of Science, Basic Energy Sciences with contract numbers DE-FG02-07ER46352 and DE-FG02-08ER46540 (CMSN), and benefited from allocation of computer time at the NERSC and NU-ASCC computation centers.</p>
			</sec>
		</ack>
		<refgrp><bibl id="B1"><title><p>Nanotube radio</p></title><aug><au><snm>Jensen</snm><fnm>K</fnm></au><au><snm>Weldon</snm><fnm>J</fnm></au><au><snm>Garcia</snm><fnm>H</fnm></au><au><snm>Zettl</snm><fnm>A</fnm></au></aug><source>Nano lett</source><pubdate>2007</pubdate><volume>7</volume><issue>11</issue><fpage>3508</fpage><lpage>3511</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1021/nl0721113</pubid><pubid idtype="pmpid" link="fulltext">17973438</pubid></pubidlist></xrefbib></bibl><bibl id="B2"><title><p>Nanoscale magnetic resonance imaging</p></title><aug><au><snm>Degen</snm><fnm>CL</fnm></au><au><snm>Poggio</snm><fnm>M</fnm></au><au><snm>Mamin</snm><fnm>HJ</fnm></au><au><snm>Rettner</snm><fnm>CT</fnm></au><au><snm>Rugar</snm><fnm>D</fnm></au></aug><source>P Natl Acad Sci Usa</source><pubdate>2009</pubdate><volume>106</volume><issue>5</issue><fpage>1313</fpage><lpage>1317</lpage><xrefbib><pubid idtype="doi">10.1073/pnas.0812068106</pubid></xrefbib></bibl><bibl id="B3"><title><p>Spectrum of an oscillator with jumping frequency and the interference of partial susceptibilities</p></title><aug><au><snm>Dykman</snm><fnm>MI</fnm></au><au><snm>Khasin</snm><fnm>M</fnm></au><au><snm>Portman</snm><fnm>J</fnm></au><au><snm>Shaw</snm><fnm>SW</fnm></au></aug><source>Phys Rev Lett</source><pubdate>2010</pubdate><volume>105</volume><issue>23</issue><fpage>230601</fpage><xrefbib><pubid idtype="pmpid" link="fulltext">21231442</pubid></xrefbib></bibl><bibl id="B4"><title><p>Efficient wireless non-radiative mid-range energy transfer</p></title><aug><au><snm>Karalis</snm><fnm>A</fnm></au><au><snm>Joannopoulos</snm><fnm>J</fnm></au><au><snm>Soljacic</snm><fnm>M</fnm></au></aug><source>Ann Phys</source><pubdate>2008</pubdate><volume>323</volume><fpage>34</fpage><lpage>48</lpage><xrefbib><pubid idtype="doi">10.1016/j.aop.2007.04.017</pubid></xrefbib></bibl><bibl id="B5"><title><p>Magnetite in human tissues: a mechanism for the biological effects of weak ELF magnetic fields</p></title><aug><au><snm>Kirschvink</snm><fnm>JL</fnm></au><au><snm>Kobayashi-Kirschvink</snm><fnm>A</fnm></au><au><snm>Diaz-Ricci</snm><fnm>JC</fnm></au><au><snm>Kirschvink</snm><fnm>SJ</fnm></au></aug><source>Bioelectromagnetics</source><pubdate>1992,</pubdate><volume>Suppl 1</volume><fpage>101</fpage><lpage>13</lpage><xrefbib><pubid idtype="pmpid">1285705</pubid></xrefbib></bibl><bibl id="B6"><title><p>Optimal frequency for wireless power transmission into dispersive tissue</p></title><aug><au><snm>Poon</snm><fnm>A</fnm></au><au><snm>O&#8217;Driscoll</snm><fnm>S</fnm></au><au><snm>Meng</snm><fnm>T</fnm></au></aug><source>Antennas and Propagation, IEEE Trans</source><pubdate>2010</pubdate><volume>58</volume><issue>5</issue><fpage>1739</fpage><lpage>1750</lpage></bibl><bibl id="B7"><title><p>Coupled-mode theory for general free-space resonant scattering of waves</p></title><aug><au><snm>Hamam</snm><fnm>R</fnm></au><au><snm>Karalis</snm><fnm>A</fnm></au><au><snm>Joannopoulos</snm><fnm>J</fnm></au><au><snm>Solja&#269;i&#263;</snm><fnm>M</fnm></au></aug><source>Phys Rev A</source><pubdate>2007</pubdate><volume>75</volume><issue>5</issue><fpage>53801</fpage></bibl><bibl id="B8"><title><p>Magnetic dissipation and fluctuations in individual nanomagnets measured by ultrasensitive cantilever magnetometry</p></title><aug><au><snm>Stipe</snm><fnm>B</fnm></au><au><snm>Mamin</snm><fnm>H</fnm></au><au><snm>Stowe</snm><fnm>T</fnm></au><au><snm>Kenny</snm><fnm>T</fnm></au><au><snm>Rugar</snm><fnm>D</fnm></au></aug><source>Phys Rev Lett</source><pubdate>2001</pubdate><volume>86</volume><issue>13</issue><fpage>2874</fpage><lpage>2877</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1103/PhysRevLett.86.2874</pubid><pubid idtype="pmpid" link="fulltext">11290061</pubid></pubidlist></xrefbib></bibl><bibl id="B9"><title><p>A tunable carbon nanotube resonator</p></title><aug><au><snm>Sazonova</snm><fnm>V</fnm></au></aug><source>Dissertation</source><note>Cornell University, Physics Department; 2006</note></bibl><bibl id="B10"><title><p>Efficient magnetic torque transduction in biological environments using tunable nanomechanical resonators</p></title><aug><au><snm>Javaheri</snm><fnm>H</fnm></au><au><snm>Barbiellini</snm><fnm>B</fnm></au><au><snm>Noubir</snm><fnm>G</fnm></au></aug><source>Conf Proc IEEE Eng Med Biol Soc</source><pubdate>2011</pubdate><volume>2011</volume><fpage>1863</fpage><lpage>1866</lpage><xrefbib><pubid idtype="pmpid" link="fulltext">22254693</pubid></xrefbib></bibl><bibl id="B11"><aug><au><snm>Blatt</snm><fnm>JM</fnm></au><au><snm>Weisskopf</snm><fnm>VF</fnm></au></aug><source>Theoretical Nuclear Physics</source><publisher>New York: Wiley</publisher><pubdate>1952</pubdate></bibl><bibl id="B12"><aug><au><snm>Jackson</snm><fnm>J</fnm></au></aug><source>Classical Electrodynamics</source><publisher>New York: Wiley</publisher><pubdate>1967</pubdate></bibl><bibl id="B13"><title><p>Magnetotactic bacteria</p></title><aug><au><snm>Blakemore</snm><fnm>R</fnm></au></aug><source>Science</source><pubdate>1975</pubdate><volume>190</volume><issue>4212</issue><fpage>377</fpage><lpage>379</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1126/science.170679</pubid><pubid idtype="pmpid" link="fulltext">170679</pubid></pubidlist></xrefbib></bibl><bibl id="B14"><title><p>Biophysics of magnetic orientation: strengthening the interface between theory and experimental design</p></title><aug><au><snm>Kirschvink</snm><fnm>JL</fnm></au><au><snm>Winklhofer</snm><fnm>M</fnm></au><au><snm>Walker</snm><fnm>MM</fnm></au></aug><source>J R Soc Interface</source><pubdate>2010</pubdate><volume>7</volume><fpage>S179</fpage><lpage>S191</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1098/rsif.2009.0491.focus</pubid><pubid idtype="pmcid">2843999</pubid><pubid idtype="pmpid" link="fulltext">20071390</pubid></pubidlist></xrefbib></bibl><bibl id="B15"><title><p>Critical superparamagnetic/single-domain grain sizes in interacting magnetite particles: implications for magnetosome crystals</p></title><aug><au><snm>Muxworthy</snm><fnm>AR</fnm></au><au><snm>Williams</snm><fnm>W</fnm></au></aug><source>J R Soc Interface</source><pubdate>2009</pubdate><volume>6</volume><issue>41</issue><fpage>1207</fpage><lpage>1212</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1098/rsif.2008.0462</pubid><pubid idtype="pmcid">2817151</pubid><pubid idtype="pmpid" link="fulltext">19091684</pubid></pubidlist></xrefbib></bibl><bibl id="B16"><title><p>Magnetic particle hyperthermia: nanoparticle magnetism and materials development for cancer therapy</p></title><aug><au><snm>Hergt</snm><fnm>R</fnm></au><au><snm>Dutz</snm><fnm>S</fnm></au><au><snm>M&#252;ller</snm><fnm>R</fnm></au><au><snm>Zeisberger</snm><fnm>M</fnm></au></aug><source>J Phys: Condens Matter</source><pubdate>2006</pubdate><volume>18</volume><fpage>S2919</fpage><xrefbib><pubid idtype="doi">10.1088/0953-8984/18/38/S26</pubid></xrefbib></bibl><bibl id="B17"><title><p>A quantitative assessment of torque-transducer models for magnetoreception</p></title><aug><au><snm>Winklhofer</snm><fnm>M</fnm></au><au><snm>Kirschvink</snm><fnm>JL</fnm></au></aug><source>J R Soc Interface</source><pubdate>2010</pubdate><volume>7</volume><fpage>S273</fpage><lpage>S289</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1098/rsif.2009.0435.focus</pubid><pubid idtype="pmcid">2843997</pubid><pubid idtype="pmpid" link="fulltext">20086054</pubid></pubidlist></xrefbib></bibl><bibl id="B18"><title><p>Vibrational resonances in biological systems at microwave frequencies</p></title><aug><au><snm>Adair</snm><fnm>RK</fnm></au></aug><source>Biophys J</source><pubdate>2002</pubdate><volume>82</volume><issue>3</issue><fpage>1147</fpage><lpage>1152</lpage><xrefbib><pubidlist><pubid idtype="doi">10.1016/S0006-3495(02)75473-8</pubid><pubid idtype="pmcid">1301920</pubid><pubid idtype="pmpid" link="fulltext">11867434</pubid></pubidlist></xrefbib></bibl><bibl id="B19"><title><p>Protein-protein interactions: methods for detection and analysis</p></title><aug><au><snm>Phizicky</snm><fnm>EM</fnm></au><au><snm>Fields</snm><fnm>S</fnm></au></aug><source>Microbiol Rev</source><pubdate>1995</pubdate><volume>59</volume><fpage>94</fpage><lpage>123</lpage><xrefbib><pubidlist><pubid idtype="pmcid">239356</pubid><pubid idtype="pmpid" link="fulltext">7708014</pubid></pubidlist></xrefbib></bibl><bibl id="B20"><title><p>Protein-passivated Fe3O4 nanoparticles: low toxicity and rapid heating for thermal therapy</p></title><aug><au><snm>Samanta</snm><fnm>B</fnm></au><au><snm>Yan</snm><fnm>H</fnm></au><au><snm>Fischer</snm><fnm>NO</fnm></au><au><snm>Shi</snm><fnm>J</fnm></au><au><snm>Jerry</snm><fnm>DJ</fnm></au><au><snm>Rotello</snm><fnm>VM</fnm></au></aug><source>J Mater Chem</source><pubdate>2008</pubdate><volume>18</volume><issue>11</issue><fpage>1204</fpage><xrefbib><pubidlist><pubid idtype="doi">10.1039/b718745a</pubid><pubid idtype="pmcid">2593465</pubid><pubid idtype="pmpid">19122852</pubid></pubidlist></xrefbib></bibl><bibl id="B21"><title><p>Protein-based composite materials</p></title><aug><au><snm>Hu</snm><fnm>X</fnm></au><au><snm>Cebe</snm><fnm>P</fnm></au><au><snm>Weiss</snm><fnm>AS</fnm></au><au><snm>Omenetto</snm><fnm>F</fnm></au><au><snm>Kaplan</snm><fnm>DL</fnm></au></aug><source>Mater Today</source><pubdate>2012</pubdate><volume>15</volume><issue>5</issue><fpage>208</fpage><lpage>215</lpage><xrefbib><pubid idtype="doi">10.1016/S1369-7021(12)70091-3</pubid></xrefbib></bibl><bibl id="B22"><title><p>Free rotation of magnetic nanoparticles in a solid matrix</p></title><aug><au><snm>del Barco</snm><fnm>ED</fnm></au><au><snm>Asenjo</snm><fnm>J</fnm></au><au><snm>Zhang</snm><fnm>X</fnm></au><au><snm>Pieczynski</snm><fnm>R</fnm></au><au><snm>Julia</snm><fnm>A</fnm></au><au><snm>Tejada</snm><fnm>J</fnm></au><au><snm>Ziolo</snm><fnm>R</fnm></au><au><snm>Fiorani</snm><fnm>D</fnm></au><au><snm>Testa</snm><fnm>A</fnm></au></aug><source>Chem mater</source><pubdate>2001</pubdate><volume>13</volume><issue>5</issue><fpage>1487</fpage><lpage>1490</lpage><xrefbib><pubid idtype="doi">10.1021/cm001122t</pubid></xrefbib></bibl><bibl id="B23"><title><p>Heat production by bacterial magnetosomes exposed to an oscillating magnetic field</p></title><aug><au><snm>Alphand&#233;ry</snm><fnm>E</fnm></au><au><snm>Faure</snm><fnm>S</fnm></au><au><snm>Raison</snm><fnm>L</fnm></au><au><snm>Duguet</snm><fnm>E</fnm></au><au><snm>Howse</snm><fnm>PA</fnm></au><au><snm>Bazylinski</snm><fnm>DA</fnm></au></aug><source>J Phys Chem C</source><pubdate>2011</pubdate><volume>115</volume><fpage>18</fpage><lpage>22</lpage><xrefbib><pubid idtype="doi">10.1021/jp104580t</pubid></xrefbib></bibl><bibl id="B24"><title><p>Spin microscopy&#8217;s heritage, achievements, and prospects</p></title><aug><au><snm>Sidles</snm><fnm>JA</fnm></au></aug><source>P Natl Acad Sci Usa</source><pubdate>2009</pubdate><volume>106</volume><issue>8</issue><fpage>2477</fpage><lpage>2478</lpage><xrefbib><pubid idtype="doi">10.1073/pnas.0813322106</pubid></xrefbib></bibl></refgrp>
	</bm>
</art>