Quantum master equation scheme of time-dependent density functional theory to time-dependent transport in nanoelectronic devices

In this work a practical scheme is developed for the first-principles study of time-dependent quantum transport. The basic idea is to combine the transport master equation with the well-known time-dependent density functional theory. The key ingredients of this paper include (i) the partitioning-free initial condition and the consideration of the time-dependent bias voltages which base our treatment on the Runge-Gross existence theorem; (ii) the non-Markovian master equation for the reduced (many-body) central system (i.e., the device); and (iii) the construction of Kohn-Sham master equations for the reduced single-particle density matrix, where a number of auxiliary functions are introduced and their equations of motion (EOMs) are established based on the technique of spectral decomposition. As a result, starting with a well-defined initial state, the time-dependent transport current can be calculated simultaneously along with the propagation of the Kohn-Sham master equation and the EOMs of the auxiliary functions.

Dynamic behaviors of mode III interfacial crack under a constant loading rate

In an earlier study on intersonic crack propagation, Gao et al. (J. Mech. Phys. Solids 49: 2113-2132, 2001) described molecular dynamics simulations and continuum analysis of the dynamic behaviors of a mode II dominated crack moving along a weak plane under a constant loading rate. The crack was observed to initiate its motion at a critical time after the onset of loading, at which it is rapidly accelerated to the Rayleigh wave speed and propagates at this speed for a finite time interval until an intersonic daughter crack is nucleated at a peak stress at a finite distance ahead of the original crack tip. The present article aims to analyze this behavior for a mode III crack moving along a bi-material interface subject to a constant loading rate. We begin with a crack in an initially stress-free bi-material subject to a steadily increasing stress. The crack initiates its motion at a critical time governed by the Griffith criterion. After crack initiation, two scenarios of crack propagation are investigated: the first one is that the crack moves at a constant subsonic velocity; the second one is that the crack moves at the lower shear wave speed of the two materials. In the first scenario, the shear stress ahead of the crack tip is singular with exponent -1/2, as expected; in the second scenario, the stress singularity vanishes but a peak stress is found to emerge at a distance ahead of the moving crack tip. In the latter case, a daughter crack supersonic with respect to the softer medium can be expected to emerge ahead of the initial crack once the peak stress reaches the cohesive strength of the interface.

Chain propagation kinetics on melt grafting reaction

Graft copolymerization in the molten state is of fundamental importance as a probe of chemical modification and reactive compatibilization. However, few grafting kinetic studies on reactive extrusion have been carried out because of the inherent difficulties, as expected. In this work, we have studied chain propagation kinetics on melt grafting using pre-irradiated linear low density polyethylene (LLDPE) and three monomers, acrylic acid (AA), methacrylic acid (MAA), and methyl methacrylate (MMA), as the model system. We measured the apparent chain propagation rate coefficients of grafting (k(p,g)) and homopolymerization (k(p,h)) at an initial stage for the melt grafting by FT-IR spectroscopy and electron spin resonance spectroscopy. It was observed that the convective mixing affected the rate coefficients. The magnitude of k(p,h) and k(p,g) were in the same order, but k(p,h) was slightly larger than k(p,g) The k(p,g) of the three grafting systems increased in the order: LLDPE/MMA < LLDPE/MAA < LLDPE/AA. These results are explained in terms of phase separation, solubility, and inherent reactivity of the monomer.

Interparticle distance-temperature-strain rate equivalence for the brittle-tough transition in polymer blends

From the angle of energy transformation an equation was obtained for the brittle transition in polymer blends. The effects of interparticle distance, temperature and strain rate on the brittle-tough transition in polymer blends were characterized by this equation. The calculations show that, for this transition: (1) increasing temperature and decreasing interparticle distance are equivalent and the shift factor increases with increasing temperature; (2) decreasing strain rate and decreasing interparticle distance have equivalent effects on the transition; (3) the strain rate must be optimum in order to find the brittle-tough transition phenomena for a given temperature region. (C) 1998 Published by Elsevier Science Ltd. All rights reserved.

TEMPERATURE AND STRAIN-RATE DEPENDENCE OF FRACTURE-TOUGHNESS OF PHENOLPHTHALEIN POLYETHER KETONE

A strong strain-rate and temperature dependence was observed for the fracture toughness of phenolphthalein polyether ketone (PEK-C). Two separate crack-blunting mechanisms have been proposed to account for the fracture-toughness data. The first mechanism involves thermal blunting due to adiabatic heating at the crack tip for the high temperatures studied. In the high-temperature range, thermal blunting increases the fracture toughness corresponding to an effectively higher test temperature. However, in the low-temperature range, the adiabatic temperature rise is insufficient to cause softening and Jic increases with increasing temperature owing to viscoelastic losses associated with the p-relaxation there. The second mechanism involves plastic blunting due to shear yield/flow processes at the crack tip and this takes place at slow strain testing of the single-edge notched bending (SENB) samples. The temperature and strain-rate dependence of the plastic zone size may also be responsible for the temperature and strain-rate dependence of fracture toughness.

Comparison of photobioreactors and optimal light regime for rapid vegetative propagation of transgenic Undaria pinnatifida gametophytes

BACKGROUND: Previously, tachyplesin gene (tac) has been successfully transferred into Undaria pinnatifida gametophytes using the method of microprojectile bombardment transformation. The objectives of this study were to compare and evaluate the performance of bubble-column and airlift bioreactors to determine a preferred configuration of bioreactor for vegetative propagation of transgenic U. pinnatifida gametophytes, and to then investigate the influence of light on vegetative propagation of these gametophytes, including incident light intensity, photoperiod and light quality to resolve the problems of rapid vegetative propagation within the selected bioreactor. RESULTS: Experimental results showed that final dry cell density in the airlift bioreactor was 12.7% higher than that in the bubble-column bioreactor under the optimal aeration rate of 1.2 L air min(-1) L-1 culture. And a maximum final dry cell density of 2830 mg L-1 was obtained within the airlift bioreactor using blue light at 40 mu mol m(-2) s(-1) with a light/dark cycle of 14/10 (h). Polymerase chain reaction (PCR) analysis indicated that genes (bar and tac) were not lost during rapid vegetative propagation within the airlift bioreactor. CONCLUSION: The airlift bioreactor was shown to be much more suitable for rapid vegetative propagation of transgenic U. pinnatifida gametophytes than the bubble-column bioreactor in the laboratory. The use of blue light allows improvement of vegetative propagation of transgenic U. pinnatifida gametophytes. (C) 2009 Society of Chemical Industry

Microprojectile bombardment of Laminaria japonica gametophytes and rapid propagation of transgenic lines within a bubble-column bioreactor

A human acidic fibroblast growth factor gene, hafgf, was successfully transferred into Laminaria japonica (kelp) gametophytes via microprojectile bombardment using the biolistic PDS-1000/He gene gun. Following phosphinothricin screening, PCR detection and Southern blot analysis, transgenic L. japonica gametophytes were cultivated in an illuminated bubble-column bioreactor to optimize growth conditions. A maximal final dry cell density of 1,695 mg l(-1) was obtained in a batch culture having an initial dry cell density of 129.75 mg l(-1). This was achieved using an aeration rate of 1.08 l air min(-1) l(-1) culture in a medium containing 1.5 mM inorganic nitrate and 0.15 mM phosphate. In addition, the relationship between different nitrogen sources and growth of transgenic gametophytes indicated that both urea and sodium nitrate were effective nitrogen sources for cell growth, while ammonium ions inhibited growth of these gametophytes.

Effects of Irradiance and Temperature on the Photosynthesis and Vegetative Propagation of Caulerpa serrulata

The photosynthetic oxygen evolution of Caulerpa serrulata was determined with oxygen electrodes. The effects of light and temperature on the growth and regeneration of fragmented C. serrulata thalli were analyzed. The regenerating rate and establishment of different sizes and portions of C. serrulata were studied. The results showed that the light saturation point of C. serrulata was 200 mu mol photons/m(2) per s and the optimum growth temperature was 25-30 degrees C. Under these conditions, the maximum photosynthetic oxygen evolution rate was 15.1 +/- 0.29 mg O-2/mg Chl a/h, the growth rate and elongation rate reached the highest values, 4.67 +/- 0.09 mg FW/d and 0.78 +/- 0.01 mm/d, respectively. The fragmented C. serrulata thalli was regenerated at 20-35 degrees C and survived at 15 degrees C and 200 mu mol photons/m(2) per s. A different survival rate was detected according to fragment size. All of these results indicated that C. serrulata was a candidate to become an invasive species if introduced into a new place. Therefore, we should pay more attention to C. serrulata for its potential threat to marine ecosystem when it is sold for aquarium use.

Analysis on suspended sediment deposition rate for muddy coast of reclaimed land

A new expression for calculating suspended fine-sediment deposition rate is developed based on theoretic analysis and experiments. The resulting equation is applied to simulation of fine sediment deposition in the reclaimed land in the Hangzhou Bay, China. The hydrodynamic environment in this area is solved by use of a long wave model, which gives the 2D-velocity field and considers bathymetric changes due to fine sediment deposition. The expression is proved convenient to use in engineering practice, and the predicted deposition rate agrees with the annual data available from field measurements from the first year to the third year after the construction of the long groin as a reclaiming method.

Asymptotic periodicity of the Volterra equation with infinite delay

For some species, hereditary factors have great effects on their population evolution, which can be described by the well-known Volterra model. A model developed is investigated in this article, considering the seasonal variation of the environment, where the diffusive effect of the population is also considered. The main approaches employed here are the upper-lower solution method and the monotone iteration technique. The results show that whether the species dies out or not depends on the relations among the birth rate, the death rate, the competition rate, the diffusivity and the hereditary effects. The evolution of the population may show asymptotic periodicity, provided a certain condition is satisfied for the above factors. (c) 2006 Elsevier Ltd. All rights reserved.

Traveling wave front for the Fisher equation on an infinite band region

Instead of discussing the existence of a one-dimensional traveling wave front solution which connects two constant steady states, the present work deals with the case connecting a constant and a nonhomogeneous steady state on an infinite band region. The corresponding model is the well-known Fisher equation with variational coefficient and Dirichlet boundary condition. (c) 2006 Elsevier Ltd. All rights reserved.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.