YV_(1-x)B_xO_(4-x):Eu~(3+)(或Dy~(3+),Bi~(3+))的合成及光谱性质的研究

本文报道了Eu~(3+)、D~(3+)在YV_(1-x)B_xO_(4-x)中的光谱性质。实验表明:B~(3+)取代V~(5+)时,可使基质和Dy~(3+)的发射强度及黄兰比下降,并使基质的吸收带短移。而B~(3+)取代Y~(3+)时,则可提高Dy~(3+)的黄兰比。同时我们还研究了Bi~(3+)和温度对YV_(0.69)B_(0.31)O_(3.69):Dy~(3+)中Dy~(3+)发射强度的影响。

CaF_2:Ce~(3+)晶体的光谱性质及γ线辐照效应

Ce~(3+)只有一个4f电子,基态能级为~2F_(5/2)和~2F(7/2),两能级差约为2200cm~(-1),Ce~(3+)不同于其它三价稀土离子,一般是d→f跃迁,特征发射为d→~2F_(5/2)和d→~2F_(7/2)跃迁引起的两宽带。CaF_2∶Ce~(3+)晶体中,当Ce~(3+)取代Ca~(2+)时,存在电荷补偿问题。Manthey讨论了CaF_2∶Ce~(3+)中间隙F~-充当电荷剂问题,Feofilov探讨了CaF_2∶Ce~(3+)的氧补偿问题,指出当有氧存在时,它

海洋聚磷菌Halomonas YSR-3的除磷特性研究

本论文报道从海洋中分离到的一株聚磷菌的分离、鉴定、在系统发育中的地位、除磷特性、菌体内多磷酸盐颗粒的研究、D-海因酶和核苷二磷酸激酶基因的克隆及序列分析，为海水系统的生物除磷提供部分基础资料。 从黄海海域分离到聚磷菌Halomonas sp. YSR-3，菌体呈杆状，大小为3.5 μm×1 μm，革兰氏阴性，好氧生长，能运动。透射电镜观察发现，菌体内有致密颗粒。经DAPI染色确定该致密颗粒是多磷酸盐，亦可称为异染粒、迂回体。16S rDNA鉴定结果表明，YSR-3与Halomonas属中的marine bacterium B5-7有较高的同源性，相似值99%。YSR-3的生理生化特性：对氯霉素和卡那霉素敏感；淀粉水解呈阳性；反硝化和几丁质降解呈阴性；能将葡萄糖作为唯一碳源和能源。 对YSR-3的培养条件进行优化。以海水2216培养基、24 ℃、180 rpm、pH 6.5的条件培养，更利于菌体生长和菌体内多磷酸盐的形成。 对YSR-3的除磷特性进行研究。无磷培养时，菌体不能生长；用磷酸钾盐作为磷源时，菌体生长较好，形成多磷酸盐的菌体比例较高；较适合YSR-3菌体生长和多磷酸盐形成的磷源是KH2PO4，较适磷浓度为1.5 mmol/L。pH的变化影响菌株的生长、多磷酸盐形成和除磷效果。pH值为5时，菌体的数量几乎不增加，体内多磷酸盐和培养基中磷含量变化不大；pH值为6、7和8时，菌体生长良好，95%以上的菌体内形成多磷酸盐，培养基中磷含量明显下降。YSR-3在不同培养基中除磷量和除磷率不同。在高磷培养基中除磷量为0.7 mmol/L(磷含量由1.84 mmol/L降到1.14 mmol/L)，除磷率为37.5%；在低磷培养基中除磷量为0.02 mmol/L(磷含量由0.028 mmol/L降到0.008 mmol/L)，除磷率为72.2%。 以海洋聚磷菌Halomonas sp. YSR-3的总DNA为模板，用PCR法扩增D-海因酶基因和核苷二磷酸激酶基因，将扩增片段克隆到pGM-T载体，转化E.coli TOP10菌株，经蓝白斑筛选、菌落PCR得到阳性克隆，测序后对序列进行Blast比对分析。得到的D-海因酶基因序列长度为1510 bp，与Pseudomonas entomophila L48的海因酶基因序列的相似性为77%。翻译后的序列与Pseudomonas fluorescens Pf-5，Marinomonas sp. MED121，Burkholderia vietnamiensis G4的海因酶蛋白序列相似性分别为75%，73%，70%。得到的核苷二磷酸激酶基因序列长度为420bp，翻译后的序列与Loktanella vestfoldensis SKA53，Jannaschia sp. CCS1，Roseobacter sp. CCS2的核苷二磷酸激酶蛋白序列相似性分别为89%，86%，85%。 聚磷菌能将外界环境中的磷吸收到体内，并以多磷酸盐的形式储存。多磷酸盐对于细胞的生存和生长有很重要的作用，但目前对于多磷酸盐的形成过程以及过程调控还不是很清楚。在今后可以通过构建高效表达的重组菌，提高与除磷相关的酶的纯度及活性。同时可以将相关酶的基因进行突变，对基因表达的调控以及酶的代谢以及功能结构等多方面进行基础研究，使聚磷菌在生物除磷中得到广泛应用。

景物的三维信息获取及其运动分析

采用模糊熵函数对图象象素分类作出整体最优分类评价,实现了区域分割.利用矩及其函数做为各区域的特征表达,构成以区域为基元的符号特征集并描述图象内容。根据立体图象对间的几何关系,解出各区域(基元)的相对三维坐标。与象索匹配相比较,它可以获得较高精度的三维信息和可描述的景物信息.通过获取不同时空的各区域(基元)三维信息,确定了它们的空间运动状态。联系这些状态,构造出景物中物体间的空间关系和近似模型,实现了对景物的3-D识别和描述。

波前重建法射线追踪及走时表压缩

As a fast and effective method for approximate calculation of seismic numerical simulation, ray tracing method, which has important theory and practical application value, in terms of seismic theory and seismic simulation, inversion, migration, imaging, simplified from seismic theory according to geometric seismic, means that the main energy of seismic wave field propagates along ray paths in condition of high-frequency asymptotic approximation. Calculation of ray paths and traveltimes is one of key steps in seismic simulation, inversion, migration, and imaging. Integrated triangular grids layout on wavefront with wavefront reconstruction ray tracing method, the thesis puts forward wavefront reconstruction ray tracing method based on triangular grids layout on wavefront, achieves accurate and fast calculation of ray paths and traveltimes. This method has stable and reasonable ray distribution, and overcomes problems caused by shadows in conventional ray tracing methods. The application of triangular grids layout on wavefront, keeps all the triangular grids stable, and makes the division of grids and interpolation of a new ray convenient. This technology reduces grids and memory, and then improves calculation efficiency. It enhances calculation accuracy by accurate and effective description and division on wavefront. Ray tracing traveltime table, which shares the character of 2-D or 3-D scatter data, has great amount of data points in process of seismic simulation, inversion, migration, and imaging. Therefore the traveltime table file will be frequently read, and the calculation efficiency is very low. Due to these reasons, reasonable traveltime table compression will be very necessary. This thesis proposes surface fitting and scattered data compression with B-spline function method, applies to 2-D and 3-D traveltime table compression. In order to compress 2-D (3-D) traveltime table, first we need construct a smallest rectangular (cuboidal) region with regular grids to cover all the traveltime data points, through the coordinate range of them in 2-D surface (3-D space). Then the value of finite regular grids, which are stored in memory, can be calculated using least square method. The traveltime table can be decompressed when necessary, according to liner interpolation method of 2-D (3-D) B-spline function. In the above calculation, the coefficient matrix is stored using sparse method and the liner system equations are solved using LU decomposition based on the multi-frontal method according to the sparse character of the least square method matrix. This method is practiced successfully in several models, and the cubic B-spline function can be the best basal function for surface fitting. It make the construction surface smooth, has stable and effective compression with high approximate accuracy using regular grids. In this way, through constructing reasonable regular grids to insure the calculation efficiency and accuracy of compression and surface fitting, we achieved the aim of traveltime table compression. This greatly improves calculation efficiency in process of seismic simulation, inversion, migration, and imaging.

特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

北祁连与块状硫化物矿床有关的早古生代火山岩的地球化学研究

Different conclusions from previous work are made from the geochemical study for the early Paleozoic volcanic rocks hosting massive sulfide deposits in the north Qilian Orogen. The main points are: (1)The geochemical characteristics of the basalts and rhyolites from the Baiyin deposit are not consistent with that of the volcanic rocks in the continental rift setting, but show the relationship with subduction. The basalts and rhyolites from the Baiyin deposit are probably individual tectonic slice piled by subduction, and there is no bimodal volcanic rock suite occurred in the Baiyin deposit. Zircon U-Pb dating constrains the magmatic emplacement of basalts and rhyolites at 475±10Ma and 453±12Ma, respectively. The basalts are characterized by enriched Th and Sr, and depleted Nb, Ta and Ti. They have relatively high Th/Nb ratios between 0.9 and 1.3. Their εNd(T) values vary from -1.2 to +3.4. The chemical and isotopic compositions display a typical subduction-related signature, and they suggest that an enriched component with the isotopic composition of EMII might have contributed to the generation of the Baiyin basalts. The basalts were likely formed in a mature island-arc or a volcanic arc built on comparatively young or thin continental crust in an active continental margin. The rhyoIites have low concentrations of LILE compared to the basalts. They do not seen to have a relationship with the basalts, because of their significantly higher εNd(T) values (+4.3～+7.7). The high and positive εNd(T) values also rule out their derivation from anatexis of the continental crust. A modeling study suggests that the source.of the Zhe-Huo and Xiaotieshan rhyolites is similar to boninite and IAT (island-arc tholeiite), and hence indicating an intra-oceanic arc environment. (2) The formation of the Shangliugou volcanic rocks from .Qilian area is also related to subduction. The basaltic andesite have low TiO_2(0.45～0.63%) and P_2O_5(0.04～0.09) content, and high Th/Nb ratios (0.3～0.6). They show flat REE patterns. Their εNd(T) values vary in a narrow range from +4.8 to +6.4. The chemical and isotopic compositions indicate that they are derived from a slightly depleted mantle source and are fromed in intra-island arc setting. The rhyolites show calc-alkaline trend. They show enriched LREE and fiat HREE patterns with obvious negative Eu anomaly. They have high Th/Ta ratios (5.0 ～ 11.7) and large negative εNd(T) values (-2.6 ～ -8.4). The rhyolites are formed in active continental margin and result from a mixed process of two endmembers, or crust assimilation. (3) The metal elements of the volcanic-hosted massive sulfide deposit have two sources, the copper and zinc are derived from rhyolitic magmas whereas the lead are probably related to old sediments overlying the rhyolites. (4) It is suggested here that the volcanic rocks hosting massive sulfide deposit in the north Qilian orogen, which are previously considered as a bimodal suite of Neo-proterozoic to middle Cambrian age in a continental rift, are virtually related to subduction magmatism in Ordovician age, and there might have no continental rift magmatism of Neo-proterozoic to middle Cambrian in the north Qilian.

东营凹陷北部陡坡带砂砾岩体性油气藏形成机理及分布规律研究

The north steep slope zone of Dongying Depression has great potential in oil resource and as the usage of 3-d seismic data in the last decade, the exploration of oil and gas has get into the stage of sandy glavel body lithological oil-gas pool exploration. In this thesis, writer take the north steep slope zone of Dongying Depression as target area and take Sha-III and Sha-[V Menber as purpose stratum, study on sequence stratigraphy, depositional system, reservoir description, emphasesing on analyzing of forming of lithological oil-gas pool, especially the dynamics principle of oil and gas preliminary movement from the source rock to the reservoir form lithological oil-gas pools. The aim of this work is to give some quantitatively explanation for the mechanism of lithological oil-gas pool forming, and set up the theory of pool form with characteristic terrestrial faulted basin. There are main conclusions and views as follow. 1. Applying with principle of sequence stratigrapgy, according to the depositional cycles of Dongying Depression, the sequence stratigraphical partition of Tertiary was finished, stressing on dismembering Sha-III and Sha-IV Menber as 5system tracts. 2. The structure of Dongying Depression especially of the north steep slope zone has accomplished, including the analyzing the structural cortroling to depositional condition of the north steep slope zone of Dongying Depression, discussed relationship between the structure of the north steep slope zone and the pool-forming. 3. The horizontal and vertical exchanges of ancient climates and ancient physiognomy of the all stratum units and studies on characteristic of depositional system distribution have been finished, found that there are five depositional systems in the north steep slope zone of Dongying Depression as fluvial, delta (tan-delta), sub-water fluvial fan lacustrine, gravitive flow, and seven formations of sandy glavel body, and forecasting of all kinds of sandy glavel body has been made. 4. Seismic stratigraphy and log stratigraphy have been made, described and forecasted all kinds of reservoir of objective stratum by means of physical geography method, setup a series means of sandy glavel body description suit to target area. 5. The pool-forming system has been studied, analyzing all the elements in petroleum sub-system of Sha-III and Sha-IV Menber of Dongying Depression with view of source controlling, estimated the petroleum system applying source rock potential index combining with distribution ofreservior. 6.Through studying types of pool, the controlling factors of pool-forming of sandy glavel body were discussed by deposition stages, formation types, structure ect. as a conclusion that the characteristics of pool forming in the north steep slope zone of Dongying Depression are, the controlling factor of the pools is mainly lithology, petrophysics of oil sands vary greatly, with a large heterogeneity, all kind of reservoir with different formation has different pool-forming conditions, and as a result, formed various pools of sandy glavel body along the steep slope with regular combination, distribution and constituted the multiple petroleum accumulative pattern. 7. It's the first time to cauculate and estimate the fluid pressure in source rock of Dongying Depression, set up the stratum fluid pressure in Dongying Depression, and firstly use equivalent charging pressure and reservoir forming index to quantitatively evaluate the pool-forming condition of lithological pool.8. Above all studies, follow up the scent of the exploration combined with practice a lot of explorative targets were found, and got geat economic and social benefit.

单（多）井抽灌对浅部地温场的影响研究

The technique of energy extraction using groundwater source heat pumps, as a sustainable way of low-grade thermal energy utilization, has widely been used since mid-1990's. Based on the basic theories of groundwater flow and heat transfer and by employing two analytic models, the relationship of the thermal breakthrough time for a production well with the effect factors involved is analyzed and the impact of heat transfer by means of conduction and convection, under different groundwater velocity conditions, on geo-temperature field is discussed.A mathematical model, coupling the equations for groundwater flow with those for heat transfer, was developed. The impact of energy mining using a single well system of supplying and returning water on geo-temperature field under different hydrogeological conditions, well structures, withdraw-and-reinjection rates, and natural groundwater flow velocities was quantitatively simulated using the finite difference simulator HST3D. Theoretical analyses of the simulated results were also made. The simulated results of the single well system indicate that neither the permeability nor the porosity of a homogeneous aquifer has significant effect on the temperature of the production segment provided that the production and injection capability of each well in the aquifers involved can meet the designed value. If there exists a lower permeable interlayer, compared with the main aquifer, between the production and injection segments, the temperature changes of the production segment will decrease. The thicker the interlayer and the lower the interlayer permeability, the longer the thermal breakthrough time of the production segment and the smaller the temperature changes of the production segment. According to the above modeling, it can also be found that with the increase of the aquifer thickness, the distance between the production and injection screens, and/or the regional groundwater flow velocity, and/or the decrease of the production-and-reinjection rate, the temperature changes of the production segment decline. For an aquifer of a constant thickness, continuously increase the screen lengths of production and injection segments may lead to the decrease of the distance between the production and injection screens, and the temperature changes of the production segment will increase, consequently.According to the simulation results of the single well system, the parameters, that can cause significant influence on heat transfer as well as geo-temperature field, were chosen for doublet system simulation. It is indicated that the temperature changes of the pumping well will decrease as the aquifer thickness, the distance between the well pair and/or the screen lengths of the doublet increase. In the case of a low permeable interlayer embedding in the main aquifer, if the screens of the pumping and the injection wells are installed respectively below and above the interlayer, the temperature changes of the pumping well will be smaller than that without the interlay. The lower the permeability of the interlayer, the smaller the temperature changes. The simulation results also indicate that the lower the pumping-and-reinjection rate, the greater the temperature changes of the pumping well. It can also be found that if the producer and the injector are chosen reasonably, the temperature changes of the pumping well will decline as the regional groundwater flow velocity increases. Compared with the case that the groundwater flow direction is perpendicular to the well pair, if the regional flow is directed from the pumping well to the injection well, the temperature changes of the pumping well is relatively smaller.Based on the above simulation study, a case history was conducted using the data from an operating system in Beijing. By means of the conceptual model and the mathematical model, a 3-D simulation model was developed and the hydrogeological parameters and the thermal properties were calibrated. The calibrated model was used to predict the evolution of the geo-temperature field for the next five years. The simulation results indicate that the calibrated model can represent the hydrogeological conditions and the nature of the aquifers. It can also be found that the temperature fronts in high permeable aquifers move very fast and the radiuses of temperature influence are large. Comparatively, the temperature changes in clay layers are smaller and there is an obvious lag of the temperature changes. According to the current energy mining load, the temperature of the pumping wells will increase by 0.7°C at the end of the next five years. The above case study may provide reliable base for the scientific management of the operating system studied.

岩浆－热液过渡阶段体系中不相容元素地球化学行为及其机制－以新疆阿尔泰3号伟晶岩脉研究为例

伟晶岩形成和演化中岩浆-热液过渡阶段物理化学性质、过铝质富挥发分岩浆体系岩浆液相不混溶、流体相出溶事件的相对时间尺度以及在岩浆-热液过渡阶段体系中稀土和其它微量元素地球化学行为(稀土“四重效应”机制和微量元素分异及控制因素等)是目前过铝质岩浆体系研究的焦点问题。对上述问题的研究和探讨，将有助于了解伟晶岩形成和演化的全过程，有助于了解挥发分(H_20、F、B、P)对过铝质岩浆体系物理性质以及对REE和其它微量元素地球化学行为影响，这对于理解过铝质岩浆体系成岩、成矿作用过程具有重大的理论和实际意义。本论文选择著名的新疆阿尔泰3号伟晶岩脉为研究对象，由于它具有完美的分带特征，使我们有可能通过各结构带矿物学、地球化学的研究揭示伟晶岩成因和演化过程以及讨论岩浆-热液过渡阶段体系中的稀土和微量元素地球化学行为。由于LCT型伟晶岩具有明显的分带性、矿物颗粒大小不均匀性等特点，因此利用全岩样品研究伟晶岩往往因取样没有代表性而被大多数岩石学家和地球化学家所否定。考虑到这点，本文系统采集阿尔泰3号伟晶岩脉各结构带磷灰石、电气石、碱性长石、石英矿物，部分结构带的石榴石、绿柱石、锂辉石矿物以及伟晶岩冷凝边和蚀变围岩的全岩样品，通过各结构带磷灰石矿物化学组成和~(87)Sr／~(86)Sr比值测定，电气石矿物的化学组成和硼同位素组成分析，碱性长石中磷含量的电子探针分析和碱性长石-石英矿物对的氧同位素体系研究以及蚀变围岩和冷凝边的化学组成分析，本文较为详尽系统地探讨以上所有的问题。通过本次多年的研究，我们获得了以下几点新的认识：f1) 富挥发分过铝质岩浆体系的样品存在稀土“四重效应”和Y-Ho、Zr-Hf、Nb-Ta、Sr-Eu元素对的显著分异，表明在过铝质岩浆演化过程中含氟流体相起着重要的作用；目前，稀土“四重效应”机制被认为是含水流体相与过铝质岩浆熔体相互作用的结果(Irber 1999；Bau 1997)，而不是独居石、磷钇矿、石榴石矿物等矿物的早期结晶引起的残余熔体稀土含量的异常变化。由于阿尔泰3号伟晶岩脉各带磷灰石以及与其共生的石榴石、绿柱石、碱性长石、锂辉石矿物均存存明显的稀土“四重效应”以及相同电价、相似离子半径的不相容元素间存在显著的分异，并结合最近赵振华等(1999)和Sba and Chappell(1999)报道S型花岗岩全岩和单矿物(磷灰石、独居石、长石、黄玉等)均存在稀土“四重效应”这一现象，本文研究提出，稀土“四重效应”是富挥发分过铝质岩浆体系的一个基本特征，其机制既不可能由含稀土的副矿物早期结晶引起残余熔体相中REE含量变化的结果，也不能定性地归因于流体相与熔体相相互作用过程中稀土元素在流体／熔体之间分异的结果，而很可能与伟晶岩岩浆形成之前某些过程密切相关，S型花岗岩岩浆在液相线以上存在硅酸盐熔体与高盐熔体(或卤水相)的不混溶液相分离有可能是过铝质岩浆体系产生稀土“四重效应”的主要原因。(2) 由于Sr~(2+)与Eu~(2+)具有相同的电价、相似的离子半径，根据经典的地球化学原理，它们具有相似的地球化学行为，在各种地质过程中其它们具有紧密的一致性。阿尔泰3号伟晶岩磷灰石中Sr／Eu比值具有二歧变化现象，其中I、Il、III和IV带磷灰石具有低的Sr／Eu比值(16:0～111.2)，而V、VI、VII带、核部带和IV带中的一个磷灰石则显示较大的Sr／Eu比值(主要在246．6～514:9范围，其中IV带的一个磷灰石该比值在10000以上)。对世界范围内的过铝质花岗岩统计表明，过铝质岩浆岩的全岩样品中的Sr／Eu也呈现二歧变化，由此本文提出Sr／Eu比值是指示过铝质岩浆体系中岩浆流体相出溶的有效地球化学参数，大于300的Sr／Eu比值是Eu~(2+)强烈分配进入流体相的地球化学标志。(3) 阿尔泰3号伟晶岩I、II、III带电气石矿物的硼同位素组成(δ~(11)B)分布在-41.11‰～-30.90‰之间，V、VI、VII带电气石矿物硼同位素组成(δ~(11)B)在-15.23‰～-9.20‰范围，IV带电气石硼同位素分布于-39.19～13.10‰范围，如此大的硼同位素分馏无法根据实验研究获得的熔体／流体相间的分馏系数进行合理解释。我们的研究表明岩浆成因电气石的B同位素组成与电气石化学组成具相关性，随电气石Y位上AI、Li阳离子数的增大，δ~(11)B值显示明显增大的趋势。阿尔泰3号伟晶岩最初形成的结构带中电气石的B同位素组成是目前所发现的最低值之一(-41.1l～39.01‰)，指示形成伟晶岩脉的初始岩浆可能由含非海相蒸发岩的泥质岩经变质深熔作用形成，或是伟晶岩岩浆在上侵定位过程中同化了含非海相蒸发岩地层的结果。本次研究利用电气石-白云母和电气石-四硼酸盐矿物间的同位索分馏以及熔体相-流体相间的同位索分馏，合理解释了阿尔泰3号伟晶岩脉电气石B同位素分馏过程，并首次获得伟晶岩形成、演化过程中B同位素组成演化的全程图谱。(4)伟晶岩岩浆形成和演化过程是否存在由过磷引起的液相不混溶现象未见确凿的岩石学和实验地球化学证据。阿尔泰3号伟晶岩脉糖粒状钠长石带(II带)明显由两种不同组构的岩性组成，而且该带是3号脉主要的Be矿化带，也是P、REE最重要的沉淀场所。阿尔泰3号伟晶岩脉梳状结构带和I带的碱性长石具较高P含量(平均值分别为0.42 wt％和0.22 wt％)，随着岩浆的演化进行，从III带到核部带，碱性长石中P含量逐渐增大(平均值由0.21 wt％增至0.50 wt％)，而II带碱性长石中异常低的P含量(平均值为0.08、wt％)是该带大量磷灰石矿物饱和结晶，导致残余熔体相中P极度亏损所造成的，而不是晚期出溶的岩浆热液交代形成的。根据Watson (1978)、Ryerson and Hess(1980)在液相线以上温度下实验研究获得的微量元素在不混溶的基性和酸性液相中分配特征，以及Webster et a1．(1997）对德国Ehrenfriedersdorf矿区伟晶岩石英中熔体包裹体的研究成果，我们推断在富磷岩浆体系中，在其演化的早期将出现由磷引起的液相不混溶作用，它对伟晶岩成岩和Be的成矿具有重大的影响。

Biomethane production from grass silage: laboratory assessment to maximise yields

On-farm biogas production is typically associated with forage maize as the biomass source. Digesters are designed and operated with the focus of optimising the conditions for this feedstock. Thus, such systems may not be ideally suited to the digestion of grass. Ireland has ca. 3.85 million ha of grassland. Annual excess grass, surplus to livestock requirements, could potentially fuel an anaerobic digestion industry. Biomethane associated with biomass from 1.1 % of grassland in Ireland, could potentially generate over 10 % renewable energy supply in transport. This study aims to identify and optimise technologies for the production of biomethane from grass silage. Mono-digestion of grass silage and co-digestion with slurry, as would occur on Irish farms, is investigated in laboratory trials. Grass silage was shown to have 7 times greater methane potential than dairy slurry on a fresh weight basis (107 m3 t-1 v 16 m3 t-1). However, comprehensive trace element profiles indicated that cobalt, iron and nickel are deficient in mono-digestion of grass silage at a high organic loading rate (OLR) of 4.0 kg VS m-3 d-1. The addition of a slurry co-substrate was beneficial due to its wealth of essential trace elements. To stimulate hydrolysis of high lignocellulose grass silage, particle size reduction (physical) and rumen fluid addition (biological) were investigated. In a continuous trial, digestion of grass silage of <1 cm particle size achieved a specific methane yield of 371 L CH4 kg-1 VS when coupled with rumen fluid addition. The concept of demand driven biogas was also examined in a two-phase digestion system (leaching with UASB). When demand for electricity is low it is recommended to disconnect the UASB from the system and recirculate rumen fluid to increase volatile fatty acid (VFA) and soluble chemical oxygen demand (SCOD) production whilst minimising volatile solids (VS) destruction. At times of high demand for electricity, connection of the UASB increases the destruction of volatiles and associated biogas production. The above experiments are intended to assess a range of biogas production options from grass silage with a specific focus on maximising methane yields and provide a guideline for feasible design and operation of on-farm digesters in Ireland.

Finite size effects in the presence of a chemical potential: A study in the classical nonlinear O(2) sigma model

In the presence of a chemical potential, the physics of level crossings leads to singularities at zero temperature, even when the spatial volume is finite. These singularities are smoothed out at a finite temperature but leave behind nontrivial finite size effects which must be understood in order to extract thermodynamic quantities using Monte Carlo methods, particularly close to critical points. We illustrate some of these issues using the classical nonlinear O(2) sigma model with a coupling β and chemical potential μ on a 2+1-dimensional Euclidean lattice. In the conventional formulation this model suffers from a sign problem at nonzero chemical potential and hence cannot be studied with the Wolff cluster algorithm. However, when formulated in terms of the worldline of particles, the sign problem is absent, and the model can be studied efficiently with the "worm algorithm." Using this method we study the finite size effects that arise due to the chemical potential and develop an effective quantum mechanical approach to capture the effects. As a side result we obtain energy levels of up to four particles as a function of the box size and uncover a part of the phase diagram in the (β,μ) plane. © 2010 The American Physical Society.

D-Amino acid peptide residualizing agents bearing N-hydroxysuccinimido- and maleimido-functional groups and their application for trastuzumab radioiodination.

INTRODUCTION: Proteins that undergo receptor-mediated endocytosis are subject to lysosomal degradation, requiring radioiodination methods that minimize loss of radioactivity from tumor cells after this process occurs. To accomplish this, we developed the residualizing radioiodination agent N(ϵ)-(3-[(*)I]iodobenzoyl)-Lys(5)-N(α)-maleimido-Gly(1)-D-GEEEK (Mal-D-GEEEK-[(*)I]IB), which enhanced tumor uptake but also increased kidney activity and necessitates generation of sulfhydryl moieties on the protein. The purpose of the current study was to synthesize and evaluate a new D-amino acid based agent that might avoid these potential problems. METHODS: N(α)-(3-iodobenzoyl)-(5-succinimidyloxycarbonyl)-D-EEEG (NHS-IB-D-EEEG), which contains 3 D-glutamates to provide negative charge and a N-hydroxysuccinimide function to permit conjugation to unmodified proteins, and the corresponding tin precursor were produced by solid phase peptide synthesis and subsequent conjugation with appropriate reagents. Radioiodination of the anti-HER2 antibody trastuzumab using NHS-IB-D-EEEG and Mal-D-GEEEK-IB was compared. Paired-label internalization assays on BT474 breast carcinoma cells and biodistribution studies in athymic mice bearing BT474M1 xenografts were performed to evaluate the two radioiodinated D-peptide trastuzumab conjugates. RESULTS: NHS-[(131)I]IB-D-EEEG was produced in 53.8%±13.4% and conjugated to trastuzumab in 39.5%±7.6% yield. Paired-label internalization assays with trastuzumab-NHS-[(131)I]IB-D-EEEG and trastuzumab-Mal-D-GEEEK-[(125)I]IB demonstrated similar intracellular trapping for both conjugates at 1h ((131)I, 84.4%±6.1%; (125)I, 88.6%±5.2%) through 24h ((131)I, 60.7%±6.8%; (125)I, 64.9%±6.9%). In the biodistribution experiment, tumor uptake peaked at 48 h (trastuzumab-NHS-[(131)I]IB-D-EEEG, 29.8%±3.6%ID/g; trastuzumab-Mal-D-GEEEK-[(125)I]IB, 45.3%±5.3%ID/g) and was significantly higher for (125)I at all time points. In general, normal tissue levels were lower for trastuzumab-NHS-[(131)I]IB-D-EEEG, with the differences being greatest in kidneys ((131)I, 2.2%±0.4%ID/g; (125)I, 16.9%±2.8%ID/g at 144 h). CONCLUSION: NHS-[(131)I]IB-D-EEEG warrants further evaluation as a residualizing radioiodination agent for labeling internalizing antibodies/fragments, particularly for applications where excessive renal accumulation could be problematic.

A multi-scale atomistic-continuum modelling of crack propagation in a two-dimensional macroscopic plate

A novel multi-scale seamless model of brittle-crack propagation is proposed and applied to the simulation of fracture growth in a two-dimensional Ag plate with macroscopic dimensions. The model represents the crack propagation at the macroscopic scale as the drift-diffusion motion of the crack tip alone. The diffusive motion is associated with the crack-tip coordinates in the position space, and reflects the oscillations observed in the crack velocity following its critical value. The model couples the crack dynamics at the macroscales and nanoscales via an intermediate mesoscale continuum. The finite-element method is employed to make the transition from the macroscale to the nanoscale by computing the continuum-based displacements of the atoms at the boundary of an atomic lattice embedded within the plate and surrounding the tip. Molecular dynamics (MD) simulation then drives the crack tip forward, producing the tip critical velocity and its diffusion constant. These are then used in the Ito stochastic calculus to make the reverse transition from the nanoscale back to the macroscale. The MD-level modelling is based on the use of a many-body potential. The model successfully reproduces the crack-velocity oscillations, roughening transitions of the crack surfaces, as well as the macroscopic crack trajectory. The implications for a 3-D modelling are discussed.