Structure dependence of mode edges in photonic crystal waveguide with silicon on insulator

The mode edges of photonic crystal waveguide with triangular lattice based on a silicon-on-insulator slab are investigated by combination of the effective index method and two-dimensional plane wave expansion method. The variations of waveguide-mode edges with structure parameters of photonic crystal are deduced. When the ratio of the radius of air holes to the lattice constrant, r/Lambda, is fixed and the lattice constant of photonic crystal, Lambda, increases, the waveguide-mode edges shift to longer wavelengths. When Lambda is fixed and r/Lambda increases, the waveguide-mode edges shift to shorter wavelengths. Additionally, when r/Lambda and Lambda are both fixed, the radius of the two-row air holes adjacent to the waveguide increases, the waveguide-mode edges shift to shorter wavelengths.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.

Design of a compact polarizing beam splitter based on a photonic crystal ring resonator with a triangular lattice

We propose and simulate a new kind of compact polarizing beam splitter (PBS) based on a photonic crystal ring resonator (PCRR) with complete photonic bandgaps. The two polarized states are separated far enough by resonant and nonresonant coupling between the waveguide modes and the microring modes. Some defect holes are utilized to control the beam propagation. The simulated results obtained by the finite-difference time-domain method show that high transmission (over 95%) is obtained and the polarization separation is realized with a length as short as 3.1 mu m. The design of the proposed PBS can be flexible, thanks to the advantages of PCRRs.

Optical interleaver based on directional coupler in a 2D photonic crystal slab with triangular lattice of air holes

A small-size optical interleaver based on directional coupler in a 2D photonic crystal slab with triangular lattice of air holes is designed and theoretically simulated using plane wave expansion and finite-difference time-domain method. The interleaver is formed by two parallel and identical photonic crystal slab waveguides which are separated by three rows of air holes. The coupling region is designed below the light line to avoid vertical radiation. The simulated results show that the coupling coefficient is increased and the final length of the interleaver is decreased by enlarging the radius of the middle row of air holes. The transmission properties are analyzed after the interleaver's structure is optimized, and around 100 GHz channel spacing can be got when the length of the interleaver is chosen as 40.5 mu m. (c) 2006 Elsevier B.V. All rights reserved.

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

Rossby Vortex Simulation on a Paraboloidal Coordinate System Using the Lattice Boltzmann Method

In this paper, we apply our compressible lattice Boltzmann model to a rotating parabolic coordinate system to simulate Rossby vortices emerging in a layer of shallow water flowing zonally in a rotating paraboloidal vessel. By introducing a scaling factor, nonuniform curvilinear mesh can be mapped to a flat uniform mesh and then normal lattice Boltzmann method works. Since the mass per unit area on the two-dimensional (2D) surface varies with the thickness of the water layer, the 2D flow seems to be "compressible" and our compressible model is applied. Simulation solutions meet with the experimental observations qualitatively. Based on this research, quantitative solutions and many natural phenomena simulations in planetary atmospheres, oceans, and magnetized plasma, such as the famous Jovian Giant Red Spot, the Galactic Spiral-vortex, the Gulf Stream, and the Kuroshio Current, etc,, can be expected.

Simulation of the Load-Unload Response Ratio and Critical Sensitivity in the Lattice Solid Model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

Lattice instability at a fast moving crack tip

A molecular dynamics method is used to analyze the dynamic propagation of an atomistic crack tip. The simulation shows that the crack propagates at a relatively constant global velocity which is well below the Rayleigh wave velocity. However the local propagation velocity oscillates violently, and it is limited by the longitudinal wave velocity. The crack velocity oscillation is caused by a repeated process of crack tip blunting and sharpening. When the crack tip opening displacement exceeds a certain critical value, a lattice instability takes place and results in dislocation emissions from the crack tip. Based on this concept, a criterion for dislocation emission from a moving crack tip is proposed. The simulation also identifies the emitted dislocation as a source for microcrack nucleation. A simple method is used to examine this nucleation process. (C) 1996 American Institute of Physics.

Algorithm for simulating fracture processes in concrete by lattice modeling

To simulate fracture behaviors in concrete more realistically, a theoretical analysis on the potential question in the quasi-static method is presented, then a novel algorithm is proposed which takes into account the inertia effect due to unstable crack propagation and meanwhile requests much lower computational efforts than purely dynamic method. The inertia effect due to load increasing becomes less important and can be ignored with the loading rate decreasing, but the inertia effect due to unstable crack propagation remains considerable no matter how low the loading rate is. Therefore, results may become questionable if a fracture process including unstable cracking is simulated by the quasi-static procedure excluding completely inertia effects. However, it requires much higher computational effort to simulate experiments with not very high loading rates by the dynamic method. In this investigation which can be taken as a natural continuation, the potential question of quasi-static method is analyzed based on the dynamic equations of motion. One solution to this question is the new algorithm mentioned above. Numerical examples are provided by the generalized beam (GB) lattice model to show both fracture processes under different loading rates and capability of the new algorithm.

Splitting the fast and slow motions in molecular dynamics simulations based on the change of cold potential well bottom

Abstract. The atomic motion is coupled by the fast and slow components due to the high frequency vibration of atoms and the low frequency deformation of atomic lattice, respectively. A two-step approximate method was presented to determine the atomic slow motion. The first step is based on the change of the location of the cold potential well bottom and the second step is based on the average of the appropriate slow velocities of the surrounding atoms. The simple tensions of one-dimensional atoms and two-dimensional atoms were performed with the full molecular dynamics simulations. The conjugate gradient method was employed to determine the corresponding location of cold potential well bottom. Results show that our two-step approximate method is appropriate to determine the atomic slow motion under the low strain rate loading. This splitting method may be helpful to develop more efficient molecular modeling methods and simulations pertinent to realistic loading conditions of materials.

The impact of imperfect symmetry on band edge modes of a two-dimensional photonic crystal with square lattice

A theoretical analysis has been performed by means of the plane-wave expansion method to examine the dispersion properties of photons at high symmetry points of an InP based two-dimensional photonic crystal with square lattice. The Q factors are compared qualitatively. The mechanism of surface-emitting is due to the photon manipulation by periodic dielectric materials in terms of Bragg diffraction. A surface-emitting photonic crystal resonator is designed based on the phenomenon of slow light. Photonic crystal slabs with different unit cells are utilized in the simulation. The results indicate that the change of the air holes can affect the polarization property of the modes. So we can find a way to improve the polarization by reducing the symmetry of the structure.

A ministop band in a single-defect photonic crystal waveguide based on silicon on insulator

This paper reports that a two-dimensional single-defect photonic crystal waveguide in the F-K direction with triangular lattice on a silicon-on-insulator substrate is fabricated by the combination of electron beam lithography and inductively coupled plasma etching. A ministop band (MSB) is observed by the measurement of transmission characteristics. It results from the coupling between the two modes with the same symmetry, which is analysed from the stimulated band diagram by the effective index and the two-dimensional plane wave expansion methods. The parameter working on the MSB is the ratio of the radius of air holes to the lattice constant, r/a. It is obtained that the critical r/a value determining the occurrence or disappearance of MSB is 0.36. When r/a is larger than or equal to 0.36, the MSB occurs. However, when r/a is smaller than 0.36, the MSB disappears.

A lattice dynamical treatment for the total potential energy of single-walled carbon nanotubes and its applications: relaxed equilibrium structure, elastic properties, and vibrational modes of ultra-narrow tubes

In this paper, we propose a lattice dynamic treatment for the total potential energy of single-walled carbon nanotubes (SWCNTs) which is, apart from a parameter for the nonlinear effects, extracted from the vibrational energy of the planar graphene sheet. The energetics, elasticity and lattice dynamics are treated in terms of the same set of force constants, independently of the tube structures. Based upon this proposal, we have investigated systematically the relaxed lattice configuration for narrow SWCNTs, the strain energy, the Young's modulus and Poisson ratio, and the lattice vibrational properties with respect to the relaxed equilibrium tubule structure. Our calculated results for various physical quantities are nicely in consistency with existing experimental measurements. In particular, we verified that the relaxation effect makes the bond length longer and the frequencies of various optical vibrational modes softer. Our calculation provides evidence that the Young's modulus of an armchair tube exceeds that of the planar graphene sheet, and that the large diameter limits of the Young's modulus and Poisson ratio are in agreement with the experimental values of graphite; the calculated radial breathing modes for ultra-narrow tubes with diameters ranging between 2 and 5 angstrom coincide with the experimental results and the existing ab initio calculations with satisfaction. For narrow tubes with a diameter of 20 angstrom, the calculated frequencies of optical modes in the tubule's tangential plane, as well as those of radial breathing modes, are also in good agreement with the experimental measurements. In addition, our calculation shows that various physical quantities of relaxed SWCNTs can actually be expanded in terms of the chiral angle defined for the corresponding ideal SWCNTs.

Lattice polarity detection of InN by circular photogalvanic effect

We report an effective and nondestructive method based on circular photogalvanic effect (CPGE) to detect the lattice polarity of InN. Because of the lattice inversion between In- and N-polar InN, the energy band spin splitting is opposite for InN films with different polarities. Consequently under light irradiation with the same helicity, CPGE photocurrents in In- and N-polar layers will have opposite directions, thus the polarity can be detected. This method is demonstrated by our CPGE measurements in both n- and p-type InN films.

Design and Fabrication of a Photonic Crystal Channel Drop Filter Based on an Asymmetric Silicon-on-Insulator Slab

We report on the design and fabrication of a photonic crystal (PC) channel drop filter based on an asymmetric silicon-on-insulator (SOI) slab. The filter is composed of two symmetric stick-shape micro-cavities between two single-line-defect (W1) waveguides in a triangular lattice, and the phase matching condition for the filter to improve the drop efficiency is satisfied by modifying the positions and radii of the air holes around the micro-cavities. A sample is then fabricated by using electron beam lithography (EBL) and inductively coupled plasma (ICP) etching processes. The measured 0 factor of the filter is about 1140, and the drop efficiency is estimated to be 73% +/- 5% by fitting the transmission spectrum.