2B_r^{+} \over \mu_0
% \Sigma} {\partial B_z \over \partial r} 
\right] \xi_z = - {2B_r^+ \over \mu_0 \Sigma} \left[ \left( {\partial
\Phi'_{{\rm M},d} \over \partial r} \right)_{z=0} - {2 \mu_d \delta
\over r^3} \right].
\label{mod1}$$

The potential $\Phi'_{{\rm M},d}$ can be expressed in term of $\xi_z$ and its derivative with respect to $r$ using eq. (\[mpot\]) and (\[bpzp\]). Therefore equation (\[mod1\]) gives a linear equation for determining the vertical displacement $\xi_z$ induced by the tilted external dipole. The term proportional to $\mu_d$ in this equation acts like a forcing term for the vertical displacement and it causes the disc to become warped and may also excite bending waves.

As noted above, since there is no Lorentz force acting in the disc interior, $\Omega$ in equation (\[mod1\]) is the Keplerian angular velocity.

We remark that equation (\[mod1\]) is the same as equation (14) of APT with $B_z=0$, but with an additional nonhomogeneous term proportional to $\mu_d$. As noted above, this is expected as here we calculate the response which is forced by the nonaligned dipole. We remark that this nonhomogeneous term is simply the pressure force due to the misaligned component of the dipole vertically integrated through the disc:

$$P'_{m,ext} = - \frac{2 B^+_r B'^+_{r,ext}}{\mu_0} = \frac{2
B^+_r}{\mu_0} \frac{2 \mu_d \delta}{r^3},$$

where we have used equation (\[bext\]) (see also eq. \[24\] of Aly [@Aly]).

WKB waves {#sec:dispersion}
---------

When $\mu_d$ is set to zero in equation (\[mod1\]), this equation has solutions corresponding to free bending waves. In the local limit, these can be found by assuming that $\xi_z \propto \exp(ikr),$ where $
2B_r^{+ } \over \mu_0 
% \Sigma } { \partial B_z \over \partial r } 
 \right ] \xi_z = - { 2B_r^+ \over \mu_0 \Sigma } \left [ \left ({ \partial 
 \Phi'_{{\rm M},d } \over \partial r } \right)_{z=0 } - { 2 \mu_d \delta 
 \over r^3 } \right ]. 
 \label{mod1}$$ 

 The potential $ \Phi'_{{\rm M},d}$ can be expressed in condition of $ \xi_z$ and its derived function with respect to $ r$ using eq.   (\[mpot\ ]) and   (\[bpzp\ ]). Therefore equality   (\[mod1\ ]) gives a linear equation for decide the vertical displacement $ \xi_z$ induce by the tilted external dipole. The term proportional to $ \mu_d$ in this equality acts like a forcing term for the vertical translation and it causes the disc to become warped and may also excite bending wave. 

 As noted above, since there is no Lorentz force act in the disc interior, $ \Omega$ in equation   (\[mod1\ ]) is the Keplerian angular velocity. 

 We note that equation   (\[mod1\ ]) is the same as equation   (14) of APT with $ B_z=0 $, but with an additional nonhomogeneous term proportional to $ \mu_d$. As note above, this is expected as here we calculate the response which is forced by the nonaligned dipole. We remark that this nonhomogeneous term is simply the pressure force due to the misaligned component of the dipole vertically integrated through the disc: 

 $ $ P'_{m, ext } = - \frac{2 B^+_r B'^+_{r, ext}}{\mu_0 } = \frac{2 
 B^+_r}{\mu_0 } \frac{2 \mu_d \delta}{r^3},$$ 

 where we have used equation   (\[bext\ ]) (see also eq.   \[24\ ] of Aly [ @Aly ]). 

 WKB waves { # sec: distribution } 
 --------- 

 When $ \mu_d$ is sic to zero in equation   (\[mod1\ ]), this equality have solution corresponding to free bending waves. In the local limit, these can be found by assuming that $ \xi_z \propto \exp(ikr),$ where $
2B_r^{+} \lver \mu_0
% \Sigma} {\partial B_z \over \partial r} 
\right] \ei_z = - {2B_d^+ \over \mj_0 \Sigma} \left[ \left( {\partial
\Phi'_{{\cm M},e} \ovee \partial r} \right)_{z=0} - {2 \mj_d \delta
\oner r^3} \rigyt].
\lauel{mod1}$$

The potentmzl $\Phi'_{{\rm M},d}$ czk be zx'ressed in term of $\xi_z$ ang its derivatiee wnth respect to $r$ using eq. (\[mpot\]) and (\[bpz[\]). Theregoge equation (\[mod1\]) givts w lihvav equation for determining the vsrtical displacement $\xi_z$ induced by the tilted extfrnal dipole. The tfrm proportuonaj to $\mu_d$ in tfis equation acts like a forcing term for the verticau dis'lacement abd it wauses the visc tj become warivd and kay alsp excite bendikg wates.

Aw noted above, since tiere is no Lorentz fjrce actitg in the disc inteeiir, $\Omgga$ it eqjqtiun (\[mkd1\]) ia the Nepmerian anghlar velociry.

We remark that eqiaeppn (\[mod1\]) is the same ws equation (14) of APT with $B_z=0$, but with an adgitjonal nonhomogeneous teem proportional to $\mu_f$. As noteq above, this is expected as here we calculate the respkvse wmich us forced by the nonaligned dipole. We remark thwf uhix nonhomogeneoms term is simply yhf ltessure force aue to ths misaligned compojent of the eipole vewticslly integrated through the disc:

$$P'_{m,ext} = - \drac{2 B^+_r B'^+_{r,ext}}{\mu_0} = \fxac{2
B^+_r}{\mu_0} \frac{2 \mu_d \delts}{r^3},$$

where we have used eqbation (\[gext\]) (see aldo eq. \[24\] of Zuy [@Aly]).

WKB waves {#rec:civpersion}
---------

When $\mu_d$ is set to zero in xquatnon (\[mod1\]), tfis gquatiog has soluhions corresponding to fgee bgnding waves. In hhe local limit, these can be fonid by assuminb dhan $\xi_z \pro'to \exi(ikr),$ where $
2B_r^{+} \over \mu_0 % \Sigma} {\partial B_z r} \xi_z = {2B_r^+ \over \mu_0 M},d} \partial r} \right)_{z=0} {2 \mu_d \delta r^3} \right]. \label{mod1}$$ The potential $\Phi'_{{\rm can be expressed in term of $\xi_z$ and its derivative with respect to using eq. (\[mpot\]) and (\[bpzp\]). Therefore equation (\[mod1\]) gives a linear equation for the displacement induced the tilted external dipole. The term proportional to $\mu_d$ in this equation acts like a forcing for the vertical displacement and it causes the to become warped and also excite bending waves. As above, there is Lorentz acting the disc interior, in equation (\[mod1\]) is the Keplerian angular velocity. We remark that equation (\[mod1\]) is the same as (14) of $B_z=0$, but an nonhomogeneous proportional to $\mu_d$. above, this is expected as here response which is forced by the nonaligned dipole. remark that nonhomogeneous term is simply the pressure due to the misaligned component of the dipole integrated through the disc: $$P'_{m,ext} = - \frac{2 B^+_r B'^+_{r,ext}}{\mu_0} = \frac{2 B^+_r}{\mu_0} \frac{2 \mu_d we have used equation (see also eq. of [@Aly]). waves --------- When is set to zero in equation (\[mod1\]), this equation has solutions to free bending waves. In the local limit, these can by that $\xi_z \propto where $
2B_r^{+} \over \mu_0
% \Sigma} {\partial B_z \oveR \partial r} 
\rIght] \xI_z = - {2B_R^+ \ovEr \Mu_0 \SiGma} \lEft[ \left( {\partial
\pHi'_{{\rm m},d} \over \partial r} \right)_{z=0} - {2 \mu_D \deltA
\oVEr r^3} \rIGhT].
\labeL{mod1}$$

The POtENTiaL $\PHi'_{{\Rm M},D}$ cAN bE exprEssEd in terM of $\xi_z$ and iTs dErIvative with rESpEct to $r$ usinG eq. (\[Mpot\]) and (\[bpzp\]). THerEfore eQuAtiON (\[mod1\]) gIveS a linEar equATion foR determinInG The verTIcal disPLAcEmenT $\xi_z$ induced by the tILtED external dipolE. The teRm PRoPORtiOnaL to $\mu_d$ in thIs EquatIOn acts lIKe A FORciNG term for the veRtical displACemEnt and It CauSEs the dIsc to BeCOme Warped and maY alsO excite beNding wAVes.

As noTEd above, Since tHerE is No LoREnTz ForCe ACtiNG iN thE DisC interioR, $\OMeGa$ in eQuatION (\[MOd1\]) is The keplErian Angular velociTy.

WE remARk tHat eqUatioN (\[mod1\]) Is The saMe as eqUatioN (14) oF APT with $B_z=0$, but wiTh an AdditionaL noNhOmoGeNeous TErm proPorTioNal to $\mu_D$. As noteD AboVe, THIS iS expected as here we cAlCULaTe the resPonse wHIcH iS Forced by ThE noNaliGNEd dipOle. WE ReMark that This noNHoMoGeneous TeRm is siMpLy tHe pRessuRE forCe due tO the misaLigneD Component of the DIpole verticalLY iNTEgRAted ThrOugh the disc:

$$p'_{m,exT} = - \Frac{2 b^+_r B'^+_{r,EXt}}{\Mu_0} = \fRAc{2
B^+_r}{\mU_0} \frac{2 \Mu_D \DeLTa}{r^3},$$

where we have used eQuAtion (\[bExt\]) (seE also eq. \[24\] of Aly [@ALy]).

WKB waves {#SEC:DispersiOn}
---------

WhEN $\mU_D$ is set to zero in EquatIon (\[mod1\]), this EQuation hAs solUtions coRrespondiNG To free beNdiNg wAveS. In THE lOcal limit, thesE CAn be FoUnd by asSumIng that $\Xi_z \ProPto \Exp(IkR),$ where $
2B_r^{+} \over \mu_0
% \Si gma} {\par tialB_z  \o ve r \p arti al r} 
\right] \xi_ z = - {2B_r^+ \over \m u_0 \ Si g ma}\ le ft[ \ left( { \ pa r t ial 
\ Ph i'_ {{ \ rm  M},d } \ over \p artial r}\ri gh t)_{z=0} - { 2  \ mu_d \delt a
\ over r^3} \r igh t].
\l ab el{ m od1}$ $
The p otenti a l $\Ph i'_{{\rmM} , d}$ ca n  be exp r e ss ed i n term of $\xi_z$ an d  its derivativ e with  r e sp e c t t o $ r$ using e q.  (\[m p ot\]) a n d( \ [ bpz p \]). Therefor e equation( \[m od1\])  g ive s  a lin ear e qu a tio n for deter mini ng the ve rtical displac e ment $\ xi_z$ind uce d by th etil te d  ex t er nal dip ole. The  t er m pro port i o n a l to  $\ mu_d $ inthis equation  ac ts l i kea for cingterm  f or th e vert icaldi splacement andit c auses the  di sc  to  b ecome warped  an d m ay also  excite ben di n g wa ves.

As noted abo ve , si nce ther e is n o  L or e ntz forc eact ingi n  thedisc in terior,$\Omeg a $in  equati on  (\[mo d1 \])  is  theK eple rian a ngular v eloci t y.

We remarkt hat equation( \[ m o d1 \ ]) i s t he same asequa t ion(14) of  AP T  with  $B_z =0 $ ,b ut with an addition al  nonho mogen eous term pro portionalt o $\mu_d$.  Asn ot e d above, thisis ex pected ash ere we c alcul ate theresponsew h ich is f orc edbythe n on aligned dipol e .  Were mark th atthis no nho mog ene ous  t erm is si mply the  p re ss ur e f orced ue to th emis al ign ed co m ponent  of t he d ip ol e  ve rticall y  i n t egra te dthro ugh  t he di sc:$$P '_{m,ex t} = - \f rac { 2 B^ +_ rB'^+_{r ,ext}}{\mu_0}  =  \frac{2
B ^+ _r} {\mu_0 } \frac{2\mu_d \delta}{r^3},$$
w here we  ha ve us ed e quation ( \[b ext\])  (s e e also  eq. \ [24\]  o f A l y  [@Al y ] ). 

W KB  waves {#s e c :di spers io n}
- ------- -

When $\mu_d$ is set  to zero in e qua tion ( \[ mod 1 \] ) , t hi s  eq u a tion has soluti ons corres po n di ng to free ben di ng wave s. In t he lo c al limi t, thesecan be fo un d by a ssu ming that$\xi_z \ propto \e x p(ikr ) ,$  wher e $ 
2B_r^{+} \over_\mu_0
% \Sigma}_{\partial B_z \over \partial_r} 
\right]_\xi_z_= -_{2B_r^+_\over \mu_0 \Sigma}_\left[ \left( {\partial
\Phi'_{{\rm_M},d} \over \partial r}_\right)_{z=0} - {2_\mu_d_\delta
\over r^3} \right].
\label{mod1}$$

The potential $\Phi'_{{\rm M},d}$ can be expressed in term of $\xi_z$ and_its_derivative with_respect_to_$r$ using eq. (\[mpot\]) and (\[bpzp\]). Therefore_equation (\[mod1\]) gives a linear equation_for determining_the vertical displacement $\xi_z$ induced by the tilted_external_dipole. The term_proportional to $\mu_d$ in this equation acts like a_forcing term for the vertical displacement_and it causes_the_disc_to become warped and_may also excite bending waves.

As noted_above, since there is no Lorentz_force acting in the disc interior, $\Omega$_in equation (\[mod1\]) is the Keplerian angular_velocity.

We remark that equation (\[mod1\]) is_the same_as equation (14) of APT with_$B_z=0$, but with_an additional_nonhomogeneous term proportional_to $\mu_d$. As noted above, this_is expected as_here we calculate the response which_is_forced by the_nonaligned_dipole._We remark_that this nonhomogeneous_term_is simply_the_pressure force due to the misaligned_component_of the dipole vertically integrated through the_disc:

$$P'_{m,ext} = - \frac{2_B^+_r_B'^+_{r,ext}}{\mu_0} = \frac{2
B^+_r}{\mu_0} \frac{2_\mu_d \delta}{r^3},$$

where we have used_equation (\[bext\]) (see also eq. \[24\] of Aly_[@Aly]).

WKB waves_{#sec:dispersion}
---------

When $\mu_d$_is set to zero in equation (\[mod1\]), this equation has solutions corresponding_to free bending waves. In the_local limit, these can_be found_by_assuming that $\xi_z_\propto_\exp(ikr),$ where_$