Calculation of the current noise spectrum in mesoscopic transport: A quantum master equation approach

Based on our recent work on quantum transport [X. Q. Li , Phys. Rev. B 71, 205304 (2005)], we show how an efficient calculation can be performed for the current noise spectrum. Compared to the classical rate equation or the quantum trajectory method, the proposed approach is capable of tackling both the many-body Coulomb interaction and quantum coherence on an equal footing. The practical applications are illustrated by transport through quantum dots. We find that this alternative approach is in a certain sense simpler and more straightforward than the well-known Landauer-Buttiker scattering matrix theory.

Influence of lateral propagating modes on laser output characteristics in selectively oxidized vertical cavity surface-emitting lasers with double oxide layers

The influence of lateral propagating modes on the threshold current and the spontaneous emission factor in selectively oxidized vertical cavity surface-emitting lasers (VCSELs) is investigated based on the mode behaviors of lateral propagating modes and the rate equation model. The numerical results show that the lateral propagating modes may be trapped in the aperture region for the selectively oxidized VCSEL with two oxide layers, one above and one below the active region. The output characteristics of VCSELs can be affected due to the reabsorption of the quasitrapped lateral propagating modes. A lower threshold current can be expected for a VCSEL with double oxide layers than that with a single oxide layer. The numerical results of rate equations also show that a larger spontaneous emission factor can be obtained by fitting the output-input curves for the VCSEL with double oxide layers. (C) 1999 American Institute of Physics. [S0021-8979(99)07919-0].

Measurement and analysis of microwave frequency response of semiconductor optical amplifiers - art. no. 682406

The measurement and analysis of the microwave frequency response of semiconductor optical amplifiers (SOAs) are proposed in this paper. The response is measured using a vector network analyzer. Then with the direct-subtracting method, which is based on the definition of scattering parameters of optoelectronic devices, the responses of both the optical signal source and the photodetector are eliminated, and the response of only the SOA is extracted. Some characteristics of the responses can be observed: the responses are quasi-highpass; the gain increases with the bias current; and the response becomes more gradient while the bias current is increasing. The multisectional model of an SOA is then used to analyze the response theoretically. By deducing from the carrier rate equation of one section under the steady state and the small-signal state, the expression of the frequency response is obtained. Then by iterating the expression, the response of the whole SOA is simulated. The simulated results are in good agreement with the measured on the three main characteristics, which are also explained by the deduced results. This proves the validity of the theoretical analysis.

Temperature-dependent modulation characteristics for 1.3 mu m InAs/GaAs quantum dot lasers

Temperature-dependent modulation characteristics of 1.3 mu m InAs/GaAs quantum dot (QD) lasers under small signals have been carefully studied at various bias currents. Based on experimental observations, it is found that the modulation bandwidth significantly increases when excited state (ES) lasing emerges at high temperature. This is attributed to additional photons emitted by ES lasing which contribute to the modulation response. A rate equation model including two discrete electron energy levels and the level of wetting layer has been used to investigate the temperature-dependent dynamic behavior of the QD lasers. Numerical investigations confirm that the significant jump for the small signal modulation response is indeed caused by ES photons. Furthermore, we identify how the electron occupation probabilities of the two discrete energy levels can influence the photon density of different states and finally the modulation rate. Both experiments and numerical analysis show that the modulation bandwidth of QD lasers at high temperature can be increased by injecting more carriers into the ES that has larger electron state degeneracy and faster carrier's relaxation time than the ground state.

Self-Heating Effect on the Two-State Lasing Behaviors in 1.3-mu m InAs-GaAs Quantum-Dot Lasers

Experimental and theoretical study of the self-heating effect on the two-state lasing behaviors in 1.3-mu m self-assembled InAs-GaAs quantum dot (QD) lasers is presented. Lasing spectra under different injected currents, light-current (L-I) curves measured in continuous and pulsed regimes as well as a rate-equation model considering the current heating have been employed to analyze the ground-state (GS) and excited-state (ES) lasing processes. We show that the self-heating causes the quenching of the GS lasing and the ES lasing by the increased carrier escape rate and the reduced maximum modal gain of GS and ES.

Verification of EL2 electronic absorption effect on charge transfer in semi-insulating GaAs

The electronic absorption of EL2 centers has been clarified to be related to the electron acid hole photoionizations, and the transition from its ground state to metastable state, respectively. Under an illumination with a selected photon energy in the near infrared region, these three processes with different optical cross sections will show different kinetics against the illumination time. It has recently been shown that the photosensitivity (measured under 1.25 eV illumination) of the local vibrational mode absorption induced by some deep defect centers in SI-GaAs is a consequence of the electron and hole photoionizations of EL2. This paper directly measures the kinetics of the electronic transition associated with EL2 under 1.25 eV illumination, which implies the expected charge transfer among different charge states of the EL2 center. A calculation based on a simple rate equation model is in good agreement with the experimental results.

Interplay effects of temperature and injection power on photoluminescence of InAs/GaAs quantum dot with high and low areal density

In this report, we have investigated the temperature and injection power dependent photoluminescence in self-assembled InAs/GaAs quantum dots (QDs) systems with low and high areal density, respectively. It was found that, for the high-density samples, state filling effect and abnormal temperature dependence were interacting. In particular, the injection power-induced variations were most obvious at the temperature interval where carriers transfer from small quantum dots (SQDs) to large quantum dots (LQDs). Such interplay effects could be explained by carrier population of SQDs relative to LQDs, which could be fitted well using a thermal carrier rate equation model. On the other hand, for the low density sample, an abnormal broadening of full width at half maximum (FWHM) was observed at the 15-100 K interval. In addition, the FWHM also broadened with increasing injection power at the whole measured temperature interval. Such peculiarities of low density QDs could be attributed to the exciton dephasing processes, which is similar to the characteristic of a single quantum dot. The compared interplay effects of high-and low-density QDs reflect the difference between an interacting and isolated QDs system.

On the propagation rate of the chemical waves observed during the course of CO oxidation on a Ag/Pt(110) composite surface

We have analyzed the propagation rate of the chemical waves observed during the course of CO oxidation on a Ag/Pt(I 10) composite surface that were reported in our previous papers [Surf Interface Anal. 2001, 32, 179; J. Phys. Chem. B 2002, 106, 5645]. In all cases, the propagation rate v can be adequately fitted as v = v(0) + D-0/d, in which v(0) and D-0 are constants, and d is the distance between the reaction front of the chemical wave and the boundary from which the chemical wave originates. We propose that the surface species responsible for the formation of the chemical wave comes from two paths: the adsorption of molecules in the gas phase on the surface and the migration from the adjacent surface with different catalytic activity. v(0) corresponds to the contribution from the surface species due to the adsorption, and D-0/d to that of the surface species that migrates from the adjacent surface. The rate equation clearly suggests that the observed chemical wave results from the coupling between adjacent surfaces with different catalytic activities during the course of heterogeneous catalysis. These results, together with our previous reports, provide a good fundamental understanding of spillover, an important phenomenon in heterogeneous catalysis.

Simulated morphological landscape of polymer single crystals by phase-field model

The novel phase field model with the "polymer characteristic" was established based on a nonconserved spatiotemporal Ginzburg-Landau equation (TDGL model A). Especially, we relate the diffusion equation with the crystal growth faces of polymer single crystals. Namely, the diffusion equations are discretized according to the diffusion coefficient of every lattice site in various crystal growth faces and the shape of lattice is selected based on the real proportion of the unit cell dimensions.

Submesoscale activity over the shelf of the northern South China Sea in summer: simulation with an embedded model

We applied a primitive equation ocean model to simulate submesoscale activities and processes over the shelf of the northern South China Sea (NSCS) with a one-way nesting technology for downscaling. The temperature and density fields showed that submesoscale activities were ubiquitous in the NSCS shelf. The vertical velocity was considerably enhanced in submesoscale processes and could reach an average of 58 m per day in the subsurface. At this point, the mixed layer depth also was deepened along the front, and the surface kinetic energy also increased with the intense vertical movement induced by submesoscale activity. Thus, submesoscale stirring/mixing is important for tracers, such as temperature, salinity, nutrients, dissolved organic, and inorganic carbon. This result may have implication for climate and biogeochemical investigations.

边坡工程的极限分析上限法

Theory of limit analysis include upper bound theorem and lower bound theorem. To deal with slope stability analysis by limit analysis is to approximate the real solution from upper limit and lower limit. The most used method of limit analysis is upper bound theorem, therefore it is often applied to slope engineering in many cases. Although upper bound approach of limit analysis can keep away from vague constitutive relation and complex stress analyses, it also can obtain rigorous result. Assuming the critical surface is circular slip surface, two kinematically admissible velocity fields for perpendicular slice method and radial slice method can be established according to the limit analysis of upper bound theorem. By means of virtual work rate equation and strength reduction method, the upper-bound solution of limit analysis for homogeneous soil slope can be obtained. A log-spiral rotational failure mechanism for homogeneous slope is discussed from two different conditions which represent the position of shear crack passing the toe and below the toe. In the dissertition, the author also establishes a rotational failure mechanics with combination of different logarithmic spiral arcs. Furthermore, the calculation formula of upper bound solution for inhomogeneous soil slope stability problem can be deduced based on the upper bound approach of rigid elements. Through calculating the external work rate caused by soil nail, anti-slide pile, geotechnological grid and retaining wall, the upper bound solution of safety factor of soil nail structure slope, slip resistance of anti-slide pile, critical height of reinforced soil slope and active earth pressure of retaining wall can be obtained by upper bound limit analysis method. Taking accumulated body slope as subject investigated, with study on the limit analysis method to calculate slope safety factor, the kinematically admissible velocity fields of perpendicular slice method for slope with broken slip surface is proposed. Through calculating not only the energy dissipation rate produced in the broken slip surfaces and the vertical velocity discontinuity, but also the work rate produced by self-weight and external load, the upper bound solution of slope with broken slip surface is deduced. As a case study, the slope stability of the Sanmashan landslide in the area of the Three Gorges reservoir is analyzed. Based on the theory of limit analysis, the upper bound solution for rock slope with planar failure surface is obtained. By means of virtual work-rate equation, energy dissipation caused by dislocation of thin-layer and terrane can be calculated; furthermore, the formulas of safety factor for upper bound approach of limit analysis can be deduced. In the end, a new computational model of stability analysis for anchored rock slope is presented after taking into consideration the supporting effect of rock-bolts, the action of seismic force and fissure water pressure. By using the model, not only the external woke-rate done by self-weight, seismic force, fissure water pressure and anchorage force but also the internal energy dissipation produced in the slip surface and structural planes can be totally calculated. According to the condition of virtual work rate equation in limit state, the formula of safety factor for upper bound limit analysis can be deduced.

三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.

On post-instability processes in adiabatic shear in hot rolled steel

EXPERIMENTS carried out using a split Hopkinson torsional bar have shown that only one shear band develops in specimens of hot rolled steel which break during testing. We observed, however, that in specimens which were not deformed to failure, several fine shear bands appeared. We believe that these formed during the loading cycle before the appearance of the final shear band and were not due to the effect of unloading. So we developed a numerical model to study the evolution of shear banding from several finite amplitude disturbances (FADs) in both temperature and strain rate. This numerical model reveals the detailed processes by which the FADs evolve into a fully developed shear band and suggests that beyond instability, the so-called shear banding process consists of two stages: inhomogeneous shearing and true shear-banding. The latter is characterized by the collapse of the stress and an abrupt increase of the local shear strain rate.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

掺镱双包层光纤放大器的放大特性

从掺镱(Yb)光纤放大器的功率传输方程出发,利用有限差分法对小模场面积(SMA)和大模场面积(LMA)掺镱双包层光纤放大器的放大特性进行了分析比较。采用模场直径(MFD)6.5 μm和20 μm的双包层掺镱光纤作为放大器增益介质进行窄线宽连续信号的放大,在915 nm激光抽运下模拟计算了大、小模场面积输出功率随输入信号功率、抽运光功率和光纤长度的变化特性,特别是对于大模场面积光纤放大器,最优光纤长度的选择至关重要;讨论了模场直径不同时的最优抽运功率和光纤长度的选择,得出4 m光纤放大时的临界抽运功率为4