Thermoelastic Stresses in SiC Single Crystals Grown by the Physical Vapor Transport Method

A finite element-based thermoelastic anisotropic stress model for hexagonal silicon carbide polytype is developed for the calculation of thermal stresses in SiC crystals grown by the physical vapor transport method. The composite structure of the growing SiC crystal and graphite lid is considered in the model. The thermal expansion match between the crucible lid and SiC crystal is studied for the first time. The influence of thermal stress on the dislocation density and crystal quality is discussed.

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

Pressure Controlled Self-Assembly of High Quality Three-Dimensional Colloidal Photonic Crystals

A concise pressure controlled isothermal heating vertical deposition (PCIHVD) method is developed, which provides an optimal growing condition with better stability and reproducibility for fabricating photonic crystals (PCs) without the limitation of colloidal sphere materials and sizes. High quality PCs are fabricated with PCIHVD from polystyrene spheres with diameters ranging from 200 nm to 1 mu m. The deep photonic band gap and steep photonic band edge of the samples are most favorable for realizing ultrafast optical devices, photonic chips, and communications. This method makes a meaningful advance in the quality and diversity of PCs and greatly promotes their wide applications.

Silicon Surface Modification with a Mixed Silanes Layer to Immobilize Proteins for Biosensor with Imaging Ellipsometry

One kind of surface modification method on silicon wafer was presented in this paper. A mixed silanes layer was used to modify silicon surface and rendered the surface medium hydrophobic. The mixed silanes layer contained two kinds of compounds, aminopropyltriethoxysilane (APTES) and methyltriethoxysilane (NITES). A few of APTES molecules in the layer was used to immobilize covalently human immunoglobulin G (IgG) on the silicon surface. The human IgG molecules immobilized covalently on the modified surface could retain their structures well and bind more antibody molecules than that on silicon surface modified with only APTES. This kind of surface modification method effectively improved the sensitivity of the biosensor with imaging ellipsometry.

Modeling Ammonothermal Growth of GaN Single Crystals: The Role of Transport

Single crystal gallium nitride (GaN) is an important technological material used primarily for the manufacture of blue light lasers. An important area of contemporary research is developing a viable growth technique. The ammonothermal technique is an important candidate among many others with promise of commercially viable growth rates and material quality. The GaN growth rates are a complicated function of dissolution kinetics, transport by thermal convection and crystallization kinetics. A complete modeling effort for the growth would involve modeling each of these phenomena and also the coupling between these. As a first step, the crystallization and dissolution kinetics were idealized and the growth rates as determined purely by transport were investigated. The growth rates thus obtained were termed ‘transport determined growth rates’ and in principle are the maximum growth rates that can be obtained for a given configuration of the system. Using this concept, a parametric study was conducted primarily on the geometric and the thermal boundary conditions of the system to optimize the ‘transport determined growth rate’ and determine conditions when transport might be a bottleneck.

Immunosensor Interface Based on Physical and Chemical Immunogylobulin G Adsorption onto Mixed Self-Assembled Monolayers

An immunosensor interface based on mixed hydrophobic self-assembled monolayers (SAMs) of methyl and carboxylic acid terminated thiols with covalently attached human Immunoglobulin G (hIgG), is investigated. The densely packed and organised SAMs were characterised by contact angle measurements and cyclic voltammetry. The effect of the non-ionic surfactant, Tween 20, in preventing nonspecific adsorption is addressed by ellipsometry during physical and covalent hIgG immobilization on pure and mixed SAMs, respectively. It is clearly demonstrated that nonspecific adsorption due to hydrophobic interactions of hIgG on methyl ended groups is totally inhibited, whereas electrostatic/hydrogen bonding interactions with the exposed carboxylic groups prevail in the presence of surfactant. Results of ellipsometry and Atomic Force Microscopy, reveal that the surface concentration of covalently immobilized hIgG is determined by the ratio of COOH/CH3-terminated thiols in SAM forming solution. Moreover, the ellipsometric data demonstrates that the ratio of bound anti-hIgG/hIgG depends on the density of hIgG on the surface and that the highest ratio is close to three. We also report the selectivity and high sensitivity achieved by chronoamperometry in the detection of adsorbed hIgG and the reaction with its antibody.

From Hexagonally Arrayed Nanorods To Ordered Porous Film Through Controlling The Morphology Of Zno Crystals

The macrostructure can be changed by changing the morphology of its units. In this article, we use a colloidal template route, combined with hydrothermal growth method, to get the hexagonally arrayed ZnO nanorods on the polycrystalline ZnO substrate. More significantly, through controlling the morphology of ZnO crystals by adding structure-directing agent in the precursor solution, the highly ordered porous ZnO films were obtained instead of ZnO nanorods. This templated solvent-thermal method has great potential in micro/nano-fabrication. (C) 2008 Elsevier B.V. All rights reserved.

Numerical Study on Flow Field and Temperature Distribution in Growth Process of 200 mm Czochralski Silicon Crystals

The melt flow and temperature distribution in a 200 mm silicon Czochralski furnace with a cusp magnetic field was modeled and simulated by using a finite-volume based FLUTRAPP ( Fluid Flow and Transport Phenomena Program) code. The melt flow in the crucible was focused, which is a result of the competition of buoyancy, the centrifugal forces caused by the rotations of the crucible and crystal, the thermocapillary force on the free surfaces and the Lorentz force induced by the cusp magnetic field. The zonal method for radiative heat transfer was used in the growth chamber, which was confined by the crystal surface, melt surface, crucible, heat shield, and pull chamber. It was found that the cusp magnetic field could strength the dominant counter-rotating swirling flow cell in the crucible and reduce the flow oscillation and the pulling-rate fluctuation. The fluctuation of dopant and oxygen concentration in the growing crystal could thus be smoothed.

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 mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.

On the unstable stacking criterion for ideal and cracked copper-crystals

The unstable stacking criteria for an ideal copper crystal under homogeneous shearing and for a cracked copper crystal under pure mode II loading are analysed. For the ideal crystal under homogeneous shearing, the unstable stacking energy gamma(us) defined by Rice in 1992 results from shear with no relaxation in the direction normal to the slip plane. For the relaxed shear configuration, the critical condition for unstable stacking does not correspond to the relative displacement Delta = b(p)/2, where b(p) is the Burgers vector magnitude of the Shockley partial dislocation, but to the maximum shear stress. Based on this result, the unstable stacking energy Gamma(us) is defined for the relaxed lattice. For the cracked crystal under pure mode II loading, the dislocation configuration corresponding to Delta = b(p)/2 is a stable state and no instability occurs during the process of dislocation nucleation. The instability takes place at approximately Delta = 3b(p)/4. An unstable stacking energy Pi(us) is defined which corresponds to the unstable stacking state at which the dislocation emission takes place. A molecular dynamics method is applied to study this in an atomistic model and the results verify the analysis above.

The inhomogeneity of plastic-deformation in ductile single-crystals

A discrete slip model which characterizes the inhomogeneity of material properties in ductile single crystals is proposed in this paper. Based on this model rate-dependent finite element investigations are carried out which consider the finite deformation, finite rotation, latent hardening effect and elastic anisotropy. The calculation clearly exhibits the process from microscopic inhomogeneous and localized deformation to necking and the formation of LSBS and reveals several important features of shear localization. For example, the inhomogeneous deformation is influenced by the imperfections and initial non-uniformities of material properties. The inhomogeneous deformation may either induce necking which results in the lattice rotation and leads to geometrical softening, which in turn promotes the formation of CSBS, or induces heavily localized deformation. The microscopic localized deformation eventually develops into the LSBS and results in a failure. These results are in close agreement with experiment. Our calculations also find that the slip lines on the specimen's surface at necking become curved and also find that if the necking occurs before the formation of LSBS, this band must be misoriented from the operative slip systems. In this case, the formation of LSBS must involve non-crystallographic effects. These can also be indirectly confirmed by experiment. All these suggest that our present discrete slip model offers a correct description of the inhomogeneous deformation characterization in ductile crystals.

Growth of zero-mean-shear turbulent-mixed layer with suspended particles in a 2-layer fluid system

The growth behaviour of zero-mean-shear turbulent-mixed layer containing suspended solid particles has been studied experimentally and analysed theoretically in a two-layer fluid system. The potential model for estimating the turbulent entrainment rate of the mixed layer has also been suggested, including the results of the turbulent entrainment for pure two-layer fluid. The experimental results show that the entrainment behaviour of a mixed layer with the suspended particles is well described by the model. The relationship between the entrainment distance and the time, and the variation of the dimensionless entrainment rate E with the local Richardson number Ri1 for the suspended particles differ from that for the pure two-layer fluid by the factors-eta-1/5 and eta-1, respectively, where eta = 1 + sigma-0-DELTA-rho/DELTA-rho-0.

Crack extension and kinking in laminates and bicrystals

Singular fields at the tip of an interface crack in anisotropic solids are reviewed with emphasis on establishing a framework to quantify fracture resistance under mixed mode conditions. The concepts of mode mixity and surface toughness are unified by using generalized interface traction components. The similarity between the anisotropic theory and existing isotropic theory is shown. Explicit formulae are given for misoriented orthotropic bimaterials with potential applications envisioned including composite laminates and semiconductor crystals. Competition between crack extension along the interface and kinking into the substrate is investigated using a boundary layer formulation. Several case studies reveal the role of anisotropy. An explicit complex variable representation for orthotropic materials and a solution to a dislocation interacting with a crack are presented in two self-contained Appendices.