Effects of Thermocapillary Convection on the Growth Process in Industrial 8-inch Silicon Czochralski Growth

Czochralski (CZ) crystal growth process is a widely used technique in manufacturing of silicon crystals and other semiconductor materials. The ultimate goal of the IC industry is to have the highest quality substrates, which are free of point defect, impurities and micro defect clusters. The scale up of silicon wafer size from 200 mm to 300 mm requires large crucible size and more heat power. Transport phenomena in crystal growth processes are quite complex due to melt and gas flows that may be oscillatory and/or turbulent, coupled convection and radiation, impurities and dopant distributions, unsteady kinetics of the growth process, melt crystal interface dynamics, free surface and meniscus, stoichiometry in the case of compound materials. A global model has been developed to simulate the temperature distribution and melt flow in an 8-inch system. The present program features the fluid convection, magnetohydrodynamics, and radiation models. A multi-zone method is used to divide the Cz system into different zones, e.g., the melt, the crystal and the hot zone. For calculation of temperature distribution, the whole system inside the stainless chamber is considered. For the convective flow, only the melt is considered. The widely used zonal method divides the surface of the radiation enclosure into a number of zones, which has a uniform distribution of temperature, radiative properties and composition. The integro-differential equations for the radiative heat transfer are solved using the matrix inversion technique. The zonal method for radiative heat transfer is used in the growth chamber, which is confined by crystal surface, melt surface, heat shield, and pull chamber. Free surface and crystal/melt interface are tracked using adaptive grid generation. The competition between the thermocapillary convection induced by non-uniform temperature distributions on the free surface and the forced convection by the rotation of the crystal determines the interface shape, dopant distribution, and striation pattern. The temperature gradients on the free surface are influenced by the effects of the thermocapillary force on the free surface and the rotation of the crystal and the crucible.

Perturbational finite volume scheme for the one-dimensional Navier-Stokes equations

Starting from the second-order finite volume scheme,though numerical value perturbation of the cell facial fluxes, the perturbational finite volume (PFV) scheme of the Navier-Stokes (NS) equations for compressible flow is developed in this paper. The central PFV scheme is used to compute the one-dimensional NS equations with shock wave.Numerical results show that the PFV scheme can obtain essentially non-oscillatory solution.

Statistics of Turbulence in Benthal Boundary Layers

This paper deals with turbulence behavior inbenthalboundarylayers by means of large eddy simulation (LES). The flow is modeled by moving an infinite plate in an otherwise quiescent water with an oscillatory and a steady velocity components. The oscillatory one aims to simulate wave effect on the flow. A number of large-scale turbulence databases have been established, based on which we have obtained turbulencestatisticsof the boundarylayers, such as Reynolds stress, turbulence intensity, skewness and flatness ofturbulence, and temporal and spatial scales of turbulent bursts, etc. Particular attention is paid to the dependences of those statistics on two nondimensional parameters, namely the Reynolds number and the current-wave velocity ratio defined as the steady current velocity over the oscillatory velocity amplitude. It is found that the Reynolds stress and turbulence intensity profile differently from phase to phase, and exhibit two types of distributions in an oscillatory cycle. One is monotonic occurring during the time when current and wave-induced components are in the same direction, and the other inflectional occurring during the time when current and wave-induced components are in opposite directions. Current component makes an asymmetrical time series of Reynolds stress, as well as turbulence intensity, although the mean velocity series is symmetrical as a sine/cosine function. The skewness and flatness variations suggest that the turbulence distribution is not a normal function but approaches to a normal one with the increasing of Reynolds number and the current-wave velocity ratio as well. As for turbulent bursting, the dimensionless period and the mean area of all bursts per unit bed area tend to increase with Reynolds number and current-wave velocity ratio, rather than being constant as in steady channel flows.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

Torsion fracture behavior of drawn pearlitic steel wires with different heat treatments

In this paper, torsion fracture behavior of drawn pearlitic steel wires with different heat treatments was investigated. Samples with different heat treatments was investigated. Samples with different heat treatment conditions were subjected to torsion and tensile tests. The shear strain along the torsion sample after fracture was measured. Fracture surface of wires was examined by Scanning Electron Microscopy. In addition, the method of Differential Scanning Calorimetry was used to characterize the thermodynamic process in the heat treatment. A numerical simulation via finite element method on temperature field evolution for the wire during heat treatment process was performed. The results show that both strain aging and recovery process occur in the material within the temperature range between room temperature and 435 degrees C. It was shown that the ductility measured by the number of twists drops at short heating times and recovers after further heating in the lead bath of 435 degrees C. On the other hand, the strenght of the wire increases at short heating times and decreases after further heating. The microstructure inhomogeneity due to short period of heat treatment, coupled with the gradient characteristics of shear deformation during torsion results in localized shear deformation of the wire. In this situation, shear cracks nucleate between lamella and the wire breaks with low number of twists.

Subsonic interface crack with crack face contact

A steady-state subsonic interface crack propagating between an elastic solid and a rigid substrate with crack face contact is studied. Two cases with respective to the contact length are considered, i.e., semi-infinite and finite crack face contact. Different from a stationary or an open subsonic interface crack, stress singularity at the crack tip in the present paper is found to be non-oscillatory. Furthermore, in the semi-infinite contact case, the singularity of the stress field near the crack tip is less than 1/2. In the finite contact case, no singularity exists near the crack tip, but less than 1/2 singularity does at the end of the contact zone. In both cases, the singularity depends on the linear contact coefficient and the crack speed. Asymptotic solutions near the crack tip are given and analyzed. In order to satisfy the contact conditions, reasonable region of the linear contact coefficient is found. In addition, the solution predicts a non-zero-energy dissipation rate due to crack face contact.

Floating zone convection

The microgravity research, as a branch of the advanced sciences and a spe- cialized field of high technology, has been made in China since the late 1980's. The research group investigating microgravity fluid physics consisted of our col- leagues and the authors in the Institute of Mechanics of the Chinese Academy of Sciences (CAS), and we pay special attention to the floating zone convection as our first research priority. Now, the research group has expanded and is a part of the National Microgravity Laboratory of the CAS, and the research fields have been extended to include more subjects related to microgravity science. Howev- er, the floating zone convection is still an important topic that greatly holds our research interests.

On the contact width in generalized plain strain JKR model for isotropic solids

Adhesive contact model between an elastic cylinder and an elastic half space is studied in the present paper, in which an external pulling force is acted on the above cylinder with an arbitrary direction and the contact width is assumed to be asymmetric with respect to the structure. Solutions to the asymmetric model are obtained and the effect of the asymmetric contact width on the whole pulling process is mainly discussed. It is found that the smaller the absolute value of Dundurs' parameter beta or the larger the pulling angle theta, the more reasonable the symmetric model would be to approximate the asymmetric one.

Two-dimensional Vortex-induced Vibrations of Submerged Floating Tunnel Tethers

A vortex-induced vibration (VIV) model is presented for predicting the nonlinear dynamic response of submerged floating tunnel (SFT) tethers which are subjected to wave, current and tunnel oscillatory displacements at their upper end in horizontal and vertical directions. A nonlinear fluid force formula is introduced in this model, and the effect of the nonlinearity of tether is investigated. First, the tunnel is stationary and the tether vibrates due to the vortices shedding. The calculated results show that the cross-flow amplitude of VIV decreases compared with the linear model. However the in-line amplitude of VIV increases. Next, the periodical oscillation of tunnel is considered. The oscillation caused by wave forces plays the roles of parametric exciter and forcing exciter to the VIV of tether. The time history of displacement of the tether mid-span is obtained by the proposed model. It is shown that the in-line amplitude increases obviously and the corresponding frequency is changed. The cross-flow amplitude exhibits a periodic behavior.

Two-colour laser-induced electron dynamics in two-level quantum system

To attempt to control the quantum state of a physical system with a femtosecond two-colour laser field, a model for the two-level system is analysed as a first step. We investigate the coherent control of the two-colour laser pulses propagating in a two-level medium. Based on calculating the influence of the laser field with various laser parameters on the electron dynamics, it is found the electronic state can be changed up and down by choosing the appropriate laser pulses and the coherent control of the two-colour laser pulses can substantially modify the behaviour of the electronic dynamics: a quicker change of two states can be produced even for small pulse duration. Moreover, the oscillatory structures around the resonant frequency and the propagation features of the laser pulses depend sensitively on the relative phase of the two-colour laser pulses. Finally, the influence of a finite lifetime of the upper level is discussed in brief.

基于U型塑料光纤和邻啡咯啉钌的溶解氧传感器

本文提出了一种基于U型塑料光纤和邻啡咯啉钌的溶解氧传感器。光纤传感头放在样品池中，测量了不同氮氧体积比情况下荧光强度和寿命。研究了邻啡咯啉钌的浓度、传感头的退火时间和U型光纤的弯曲半径对系统灵敏度的影响。我们提出了双层模型来解释实验结果不符合S-V方程的实验现象，并用双层模型分析了弯曲半径对系统灵敏度的影响。通过对溶解氧浓度和荧光寿命亚线性曲线的拟合，我们得出敏感层的厚度和S-V系数分别为0.59和0.61。

Er3+ and Yb3+ codoped phosphate laser glass for high power flashlamp pumping

A novel Vb(3+)-Er-(3+) codoped phosphate glass for high power flashlamp pumping and high repetition rate laser at 1.54 mu m, designated EAT5-2, is developed. The weight-loss rate of is 1.3 x 10(-5) gcm(-2) h(-1) in boiling water, which is comparable to Kigre's QX-Er glass. Some spectroscopic parameters are analysed by Judd-Ofelt theory and McCumber theory The emission cross section is calculated to be 0.73 x 10(-20) cm(2). The thermo-mechanical properties of EAT5-2 are modified after an ion-exchange chemical strengthening process in a KNO3/NaNO3 molten salt bath. The thresholds for optical damage from the flashlamp pumping are tested on glass rods. A repetition rate of 15 Hz is achieved for chemically strengthened glass. The laser experimental results at. 1.54 mu m from flashlamp pumping are also reported.

多层介质膜脉冲宽度压缩光栅与超短脉冲作用时的性能分析

基于啁啾脉冲放大技术的超短脉冲激光系统是提供超快、超强激光的重要途径，具有良好输出波形和高损伤阈值的多层介质膜脉冲宽度压缩光栅是获得高峰值功率脉冲激光的关键。基于傅里叶谱变换方法和严格模式理论，分析了多层介质膜光栅（MDG）在超短脉冲作用下的光学特性。结果表明，当MDG的反射带宽小于具有高斯分布的入射脉冲的频谱宽度时，-1级反射脉冲呈非对称高斯分布，其前沿出现振荡，并且-1级反射脉冲能量开始剧烈下降，讨论了MDG结构参数对其反射带宽的影响。分析了MDG与超短脉冲作用时的近场光分布，对提高其抗激光损伤特性具

Amplitude/frequency of spontaneous mEPSC correlates to the degree of long-term depression in the CA1 region of the hippocampal slice

Prior synaptic or cellular activity influences degree or threshold for subsequent induction of synaptic plasticity, a process known as metaplasticity. Thus, the continual synaptic activity, spontaneous miniature excitatory synaptic current (mEPSC) may correlate to the induction of long-teen depression (LTD). Here, we recorded whole-cell EPSC and mEPSC alternately in the Schaffer-CA1 synapses in brain slice of young rats, and found that this recording configuration affected neither EPSC nor mEPSC. Low frequency stimulation (LFS) induced variable magnitudes of LTD. Remarkably, larger magnitudes of LTD were significantly correlated to smaller amplitude/lower frequency of the basal mEPSC. Furthermore, under the conditions reduced amplitude/frequency of the basal mEPSC by exposure to behavioral stress immediately before slice preparation or low concentration of calcium in bath solution, the magnitudes of LTD were still inversely correlated to mEPSC amplitude/frequency. These new findings suggest that spontaneous mEPSC may reflect functional and/or structural aspects of the synapses, the synaptic history ongoing metaplasticity. (C) 2005 Elsevier B.V. All rights reserved.

Long-term depression in rat CA1-subicular synapses depends on the G-protein coupled mACh receptors

The subiculum, which is the primary target of CA1 pyramidal neurons and sending efferent fibres to many brain regions, serves as a hippocampal interface in the neural information processes between hippocampal formation and neocortex. Long-term depression (LTD) is extensively studied in the hippocampus, but not at the CA1-subicular synaptic transmission. Using whole-cell EPSC recordings in the brain slices of young rats, we demonstrated that the pairing protocols of low frequency stimulation (LFS) at 3 Hz and postsynaptic depolarization of -50 mVelicited a reliable LTD in the subiculum. The LTD did not cause the changes of the paired-pulse ratio of EPSC. Furthermore, it did not depend on either NMDA receptors or voltage-gated calcium channels (VGCCs). Bath application of the G-protein coupled muscarinic acetylcholine receptors (mAChRs) antagonists, atropine or scopolamine, blocked the LTD, suggesting that mAChRs are involved in the LTD. It was also completely blocked by either the Ca2+ chelator BAPTA or the G-protein inhibitor GDP-beta-S in the intracellular solution. This type of LTD in the subiculum may play a particular role in the neural information processing between the hippocampus and neocortex. (c) 2005 Elsevier Ireland Ltd and the Japan Neuroscience Society. All rights reserved.