TY - JOUR

T1 - Elasticity of magnesite and dolomite from a genetic algorithm for inverting Brillouin spectroscopy measurements

AU - Chen, Po Fei

AU - Chiao, Ling Yun

AU - Huang, Pao Hsien

AU - Yang, Yi Jong

AU - Liu, Lin Gun

N1 - Funding Information:
We thank C.R. Bina, C.C. Lin and two anonymous reviewers for improving the manuscript. The present work was supported by research grants from the NSC.

PY - 2006/4/14

Y1 - 2006/4/14

N2 - A hybrid numerical scheme that simultaneously retrieves single-crystal elastic constants, Cij, and wave-normal directions is applied to determine the elasticity of natural samples of both magnesite and dolomite, as measured by Brillouin spectroscopy at ambient conditions. The scheme incorporates the genetic algorithm (GA) as a global searching tool and a final local linearization to overcome the intrinsic difficulty of fitting non-linear Christoffel's equation. To compensate for the stochastic nature of GA, especially when the misfit function of the problem manifests highly rugged topography in the model space, the procedure was repeated for multi-runs. This is to assure that the best solution is captured by examining the statistics among the runs that yield the fittest solutions. The resultant elastic constants C11, C12, C13, C14, C33 and C44 are, respectively, 260.3(2.6), 82.9(4.5), 59.6(3.1), (-)20.1(1.3), 153.7(4.1) and 59.7(1.4) GPa for magnesite. For dolomite, C11, C12, C13, C14, C15, C33 and C44 are, respectively, 204.1(2.2), 68.5(3.4), 45.8(4.4), 20.6(1.3), 6.7(1.5), 97.4(5.3) and 39.1(1.5) GPa. Unfortunately, the results for dolomite are statistically non-unique in the numerical calculation. The above data were adopted because the cleavage plane of dolomite was used in the experiment. The results for both magnesite and dolomite are compatible with those determined earlier by ultrasonic methods. It is anticipated that the new method developed in the present study should be applicable to less symmetric crystals, since there are no particular assumptions on the crystal symmetry embedded within the scheme.

AB - A hybrid numerical scheme that simultaneously retrieves single-crystal elastic constants, Cij, and wave-normal directions is applied to determine the elasticity of natural samples of both magnesite and dolomite, as measured by Brillouin spectroscopy at ambient conditions. The scheme incorporates the genetic algorithm (GA) as a global searching tool and a final local linearization to overcome the intrinsic difficulty of fitting non-linear Christoffel's equation. To compensate for the stochastic nature of GA, especially when the misfit function of the problem manifests highly rugged topography in the model space, the procedure was repeated for multi-runs. This is to assure that the best solution is captured by examining the statistics among the runs that yield the fittest solutions. The resultant elastic constants C11, C12, C13, C14, C33 and C44 are, respectively, 260.3(2.6), 82.9(4.5), 59.6(3.1), (-)20.1(1.3), 153.7(4.1) and 59.7(1.4) GPa for magnesite. For dolomite, C11, C12, C13, C14, C15, C33 and C44 are, respectively, 204.1(2.2), 68.5(3.4), 45.8(4.4), 20.6(1.3), 6.7(1.5), 97.4(5.3) and 39.1(1.5) GPa. Unfortunately, the results for dolomite are statistically non-unique in the numerical calculation. The above data were adopted because the cleavage plane of dolomite was used in the experiment. The results for both magnesite and dolomite are compatible with those determined earlier by ultrasonic methods. It is anticipated that the new method developed in the present study should be applicable to less symmetric crystals, since there are no particular assumptions on the crystal symmetry embedded within the scheme.

KW - Dolomite

KW - Elasticity

KW - Genetic algorithm

KW - Magnesite

UR - http://www.scopus.com/inward/record.url?scp=33644915020&partnerID=8YFLogxK

U2 - 10.1016/j.pepi.2005.10.004

DO - 10.1016/j.pepi.2005.10.004

M3 - 期刊論文

AN - SCOPUS:33644915020

SN - 0031-9201

VL - 155

SP - 73

EP - 86

JO - Physics of the Earth and Planetary Interiors

JF - Physics of the Earth and Planetary Interiors

IS - 1-2

ER -