胶东地区金矿床成矿时代及其成矿地球动力学背景-兼论壳幔相互作用与成岩成矿

Jiaodong Peninsula is the largest repository of gold in China. Varieties of studies have been involved in the mechanism of metallogenesis. This thesis is a part of the project "Study of basic geology related to the prespecting of the supra-large deposits" which supported by National Climbing Program of China to Prof. Zhou. One of the key scientific problems is to study the age and metallogenic dynamics of ore deposit and to understand how interaction between mantle and crust constrains on metallogenesis and lithogenesis. As Jiaodong Peninsula to be study area, the Rb-Sr, Sm-Nd and Pb isotopic systematics of pyrite and altered rocks are measured to define the age and origin of gold. The elemental and Sr-Nd-Pb isotopic compositions of dikes and granites was studied to implicate the source and lithogenesis of the dike and granite and removal of lithosphere and the interaction between mantle and crust in the Jiaodong Peninsula. Considering the tectonic of Jiaodong Peninsula, basic on the time and space, this thesis gives a metallogenic dynamics of gold mineralization and discusses the constraints of the interaction between mantle and crust on the metallogenesis and lithogenesis. This thesis reports the first direct Rb-Sr dating of pyrites and ores using sub-sampling from lode gold deposit in Linglong, Jiaodong Peninsula and the results demonstrate this as a useful geochronological technique for gold mineralization with poor age constraint. The Rb-Sr data of pyrites yields an isochron age of (121.6-122.7) Ma, whereas, those of ore and ore-pyrite spread in two ranges from 120.0 to 121.8 Ma and 110.0-111.7 Ma. Studies of characteristic of gold deposit, microscopy of pyrite and quartz indicate that the apparent ages of ore and ore-pyrite are not isochron ages, it was only mixed by two end members, i.e., the primitive hydrothermal fluids and wall rocks. However, the isochron age of pyrite samples constrains the age of gold mineralization, i.e., early Cretaceous, which is in good consistence with the published U-Pb ages of zircon by using the SHRIMP technique. The whole rock Rb-Sr isochron age of altered rocks indicates that the age of gold mineralizing in the Xincheng gold deposit is 116.6 ± 5.3 Ma. The Sr, Nd and Pb isotopic compositions of pyrite and altered rocks indicate that the gold and relevant elements were derived from multi-sources, i.e. dikes derived from enriched lithospheric mantle and granites, granodiorites and metamorphic rocks outcropped on the crust. It also shows that the hydrothermal fluids derived from mantle magma degassing had play an important role in the gold mineralizing. The major and trace elements, Sr-Nd-Pb isotopic data of granites and granodiorites suggest that the Linglong Granite and Kunyushan Granite were derived from partial melting of basement rocks in the Jiaodong Peninsula at post-collision of North China Craton with South China Craton. Guojialing Granodiorite was considered to be derived from a mixture source, that is, mixed by magmas derived from an enriched lithospheric mantle and crust during the delamination of lithosphere induced by the subduction of Izanagi Plate and the movement of Tancheng-Lujiang Fault. There are kinds of dikes occurred in the Jiaodong Peninsula, which are accompanying with gold mineralization in time and space. The dikes include gabrro, diabase, pyroxene diorite, gabrrophyre, granite-porphyry, and aplite. The whole rock K-Ar ages give two age intervals: 120-124 Ma for the dikes that erupted at the gold mineralizing stage, and <120 Ma of the dikes that intruded after gold mineralizing. According to the age and the relationship between the dikes and gold mineralizing, the dikes could be divided into two groups: Group I (t = 120-124 Ma) and Group II (t < 120Ma). Group I dikes show the high Mg and K, low Ti contents, negative Nb anomalies and positive Eu anomalies, high ~(87)Sr/~(86)Sr and negative εNd(t) values and an enrichment in light rare earth elements, large ion lithosphile elements and a depletion in high field strength elements. Thus the elemental and isotopic characteristics of the Group I dikes indicate that they were derived from an enriched lithospheric mantle perhaps formed by metasomatism of the melt derived from the recycled crustal materials during the deep subduction of continent. In contrast, the Group II dikes have high Ti, Mg and K contents, no negative Nb anomalies, high ~(87)Sr/~(86)Sr and positive or little negative εNd(t) values, which indicate the derivation from a source like OIB-source. The geochemical features also give the tectonic constraints of dikes, which show that Group I dikes were formed at continental arc setting, whereas Group II dikes were formed within plate background. Considering the tectonic setting of Jiaodong Peninsula during the period of gold mineralizing, the metallogenic dynamics was related to the subduction of Izanagi Plate, movement of Tancheng-Lujiang Fault and removal of lithopheric mantle during Late Mesozoic Era.

Characteristics and Sub-Characteristics of Two-Dimensional Parabolized Stability Equations

特征分析表明：对原始扰动量的抛物化稳定性方程组(PSE),它在亚超音速区分别具有椭圆和抛物特性,给出PSE特征对马赫数的依赖关系,阐明PSE仅把信息对流-扩散传播特性抛物化,而保留了信息对流-扰动传播特性,因此PSE应称为扩散抛物化稳定性方程(DPSE)。

Latin hypercube sampling used in the calculation of the fracture probability

The paper proposes Latin hypercube sampling combined with the stratified sampling of variance reduction technique to calculate accurate fracture probability. In the compound sampling, the number of simulations is relatively small and the calculation error is satisfactory.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

The Intelligent Measuring Sub-System in the Computer Integrated and Flexible Laser Processing System

Based on the computer integrated and flexible laser processing system, develop the intelligent measuring sub-system. A novel model has been built to compensate the deviations of the main frame, a new-developed 3-D laser tracker system is applied to adjust the accuracy of the system. Analyzing the characteristic of all kinds of automobile dies, which is the main processing object of the laser processing system, classify the types of the surface and border needed to be measured and be processed. According to different types of surface and border, develop 2-D adaptive measuring method based on B?zier curve and 3-D adaptive measuring method based on spline curve. During the data processing, a new 3-D probe compensation method has been described in details. Some measuring experiments and laser processing experiments are carried out to testify the methods. All the methods have been applied in the computer integrated and flexible laser processing system invented by the Institute of Mechanics, CAS.

A surface energy model and application to mechanical behavior analysis of single crystals at sub-micron scale

Size effects of mechanical behaviors of materials are referred to the variation of the mechanical behavior due to the sample sizes changing from macroscale to micro-/nanoscales. At the micro-/nanoscale, since sample has a relatively high specific surface area (SSA) (ratio of surface area to volume), the surface although it is often neglected at the macroscale, becomes prominent in governing the energy effect, although it is often neglected at the macroscale, becomes prominent in governing the mechanical behavior. In the present research, a continuum model considering the surface energy effect is developed through introducing the surface energy to total potential energy. Simultaneously, a corresponding finite element method is developed. The model is used to analyze the axial equilibrium strain problem for a Cu nanowire at the external loading-free state. As another application of the model, from dimensional analysis, the size effects of uniform compression tests on the microscale cylinder specimens for Ni and Au single crystals are analyzed and compared with experiments in literatures. (C) 2009 Elsevier B.V. All rights reserved.

Improved sampling techniques for the direct simulation Monte Carlo method

The direct simulation Monte Carlo (DSMC) method is a widely used approach for flow simulations having rarefied or nonequilibrium effects. It involves heavily to sample instantaneous values from prescribed distributions using random numbers. In this note, we briefly review the sampling techniques typically employed in the DSMC method and present two techniques to speedup related sampling processes. One technique is very efficient for sampling geometric locations of new particles and the other is useful for the Larsen-Borgnakke energy distribution.

Spectral analysis and control to obtain sub-5 fs pulses by femtosecond filamentation

A spectral-filter method is numerically demonstrated to obtain sub-5 fs pulses by using femtosecond filamentation in fused silica. Instead of employing spectral phase compensation, by properly employing a high-pass filter to select the broadened high-frequency spectra that are located almost in phase in the tailing edge of the self-compressed pulses owing to self-steepening, as short as single-cycle pulses can be obtained. For instance, for an input pulse with a duration of 50 fs and energy 2.2 mu J, the minimum pulse duration can reach to similar to 4 fs (about 1.5 cycles) by applying a proper spectral filter. (C) 2008 Optical Society of America

Sub-half-wavelength atom localization via two standing-wave fields

We propose a scheme for sub-half-wavelength atom localization in a four-level ladder-type atomic system, which is coupled by two classical standing-wave fields. We find that one of the standing-wave fields can help in enhancing the localization precision, and the other is of crucial importance in increasing the detecting probability and leading sub-half-wavelength localization.

Single sub-100 attosecond pulse generation in a two-colour time-gating laser field

We theoretically propose a method of generating a single sub-100 attosecond (as) pulse with a two-colour time-gating laser field. The field is synthesized by an 8 fs/800 nm (three optical cycles) pulse and a 24 fs/2400 nm (three optical cycles) pulse with an optimal time delay between them. In our simulation, we obtain a supercontinuum with an extremely broad spectrum of 150 eV and generate an isolated attosecond pulse with 96 as pulse duration without any dispersion compensation.

Generation of sub-5 fs pulses through cascade of filamentation and hollow fiber compression

A pulse-compression scheme based on cascade of filamentation and hollow fiber has been demonstrated, Pulses with duration of sub-5 fs and energy of 0.2 mJ near 800 nm have been generated by compressing the similar to 40 fs pulses from a commercial laser system. This method is promising to generate near monocycle high energy pulses. [GRAPHICS] Measured autocorrelation curve of the final compressed pulses with duration of sub-5 fs (black solid) and the simulated autocorrelation curve of 4.6 fs pulse near 800 rim (red dash) (C) 2008 by Astro Ltd. Published exclusively by WILEY-VCH Verlag GmbH & Co. KGaA

(Ce,Cu)∶LiNbO3晶体非挥发全息记录理论研究与优化

根据双中心带输运模型,对(Ce,Cu)∶LiNbO3晶体双中心非挥发全息记录进行了理论研究与优化。推导了(Ce,Cu)∶LiNbO3晶体的微观参量,采用数值方法通过严格求解模拟双中心带输运方程来模拟全息记录过程。分析了记录过程中,记录与敏化光强、Ce和Cu掺杂浓度以及晶体微观参量对(Ce,Cu)∶LiNbO3晶体双中心全息记录的影响。发现(Ce,Cu)∶LiNbO3晶体非挥发全息记录中实现高衍射效率与固定效率的主导因素是深中心Cu,在记录过程中,深中心Cu建立起了很强的空间电荷场。数值模拟的结果经过实验验