14.04.2013 Views

Analytical Expression for Hessian

Analytical Expression for Hessian

Analytical Expression for Hessian

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>Expression</strong> <strong>for</strong> A ij<br />

αβ :<br />

A ij<br />

αβ in expression 8 can be further simplified using<br />

∂χ mi<br />

2<br />

∂x j<br />

β<br />

= rmi α<br />

rmi ∂<br />

∂x j<br />

mi ∂ψ<br />

∂r<br />

β<br />

mi<br />

<br />

+ ∂ψmi<br />

∂rmi ∂<br />

∂x j<br />

mi rα r<br />

β<br />

mi<br />

<br />

,<br />

mi ψrr =<br />

rmi α<br />

+ ψmi<br />

rmi r rmi<br />

∂<br />

α<br />

∂rmi <br />

1<br />

rmi mi ∂r<br />

∂x j +<br />

β<br />

ψmi r<br />

rmi ∂rmi α<br />

∂x j ,<br />

β<br />

<br />

ψ<br />

=<br />

mi<br />

rr rmi α ψmi r<br />

−<br />

rmi rmi α<br />

(rmi ) 2<br />

<br />

∂rmi ∂x j +<br />

β<br />

ψmi r<br />

rmi ∂<br />

∂x j<br />

i<br />

xα − x<br />

β<br />

m α ,<br />

<br />

ψ<br />

=<br />

mi<br />

ji jm<br />

δ − δ + ψmi r<br />

rmi ji<br />

δ δαβ − δ jm <br />

δαβ ,<br />

where ψr ≡ ∂ψ<br />

∂r<br />

final <strong>for</strong>m as<br />

A ij<br />

αβ<br />

= ζ <br />

m=i<br />

χ mi<br />

1<br />

rr rmi α ψmi r<br />

−<br />

rmi rmi α<br />

(rmi ) 2<br />

<br />

mi rβ rmi etc. Thus Aij<br />

αβ in expression 8 can be re-written to obtain the<br />

<br />

ψmi rr<br />

(rmi ψmi r<br />

2 −<br />

) (rmi ) 3<br />

<br />

r mi<br />

α r mi ji jm<br />

β δ − δ <br />

+ ζδαβ χ<br />

m=i<br />

mi<br />

1<br />

A ij<br />

αβ reduces to the expression C4 of [1] <strong>for</strong> ζ = 1. Further, we can obtain the<br />

expression <strong>for</strong> the off-diagonal term from equation 10, by setting i = j as<br />

A ij<br />

<br />

ψ<br />

αβ = −ζχji 1<br />

ji<br />

<br />

r ji<br />

α rji<br />

<br />

β − ζδαβχ ji<br />

1<br />

rr<br />

(r ji )<br />

which can also be written as<br />

A ij<br />

<br />

ψ<br />

αβ = −ζχij 1<br />

ij<br />

rr<br />

(r ij )<br />

ψji r<br />

2 −<br />

(r ji ) 3<br />

ψij r<br />

2 −<br />

(r ij ) 3<br />

<br />

r ij<br />

α rij<br />

ψ<br />

β + δαβ<br />

ij<br />

r<br />

rij <br />

(9)<br />

ji ψr rji <br />

, (11)<br />

i = j, (12)<br />

which reduces to the expression C5 in [1] <strong>for</strong> ζ = 1. The expression <strong>for</strong> i = j<br />

can be obtained from equation 10<br />

<br />

<br />

ψmi r<br />

2 − r mi<br />

α rmi<br />

<br />

<br />

β + ζδαβ χ miψ<br />

1<br />

mi<br />

r<br />

, (13)<br />

rmi A ii<br />

<br />

αβ = ζ χ<br />

m=i<br />

mi<br />

1<br />

ψ mi<br />

rr<br />

(r mi )<br />

(r mi ) 3<br />

which is the generalisation of expression C6 in [1].<br />

3<br />

m=i<br />

mi ψr rmi ji jm<br />

δ − δ <br />

. (10)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!