The calculation of photoelectron angular distributions with jet-like structure from scattering theory

Photoelectron angular distributions produced in above-threshold ionization (ATI) are analysed using a nonperturbative scattering theory. The numerical results are in good qualitative agreement with recent measurements. Our study shows that the origin of the jet-like structure arises from the inherent properties of the ATI process and not from the angular momentum of either the initial or the excited states of the atom.

Effect of phases on the dynamic evolution of a Lambda system with spontaneously generated coherence

The dynamic evolution of a A system coupled by two strong coherent fields is investigated by taking spontaneously generated coherence (SGC) into account. By numericaly simulation, it is shown that the relative phase of the two coherent fields affects significantly the time scale to the coherent population trapping state. In addition, an analytical expression to the evolution rate which is consistent with the numerical results is given. (c) 2005 Elsevier B.V. All rights reserved.

Ion cascade acceleration from the interaction of a relativistic femtosecond laser pulse with a narrow thin target

Particle-in-cell simulations are performed to study the acceleration of ions due to the interaction of a relativistic femtosecond laser pulse with a narrow thin target. The numerical results show that ions can be accelerated in a cascade by two electrostatic fields if the width of the target is smaller than the laser beam waist. The first field is formed in front of the target by the central part of the laser beam, which pushes the electron layer inward. The major part of the abaxial laser energy propagates along the edges to the rear side of the target and pulls out some hot electrons from the edges of the target, which form another electrostatic field at the rear side of the target. The ions from the front surface are accelerated stepwise by these two electrostatic fields to high energies at the rear side of the target. The simulations show that the largest ion energy gain for a narrow target is about four times higher than in the case of a wide target. (c) 2006 American Institute of Physics.

Geometric integration applied to moving mesh methods and degenerate Lagrangians

Moving mesh methods (also called r-adaptive methods) are space-adaptive strategies used for the numerical simulation of time-dependent partial differential equations. These methods keep the total number of mesh points fixed during the simulation, but redistribute them over time to follow the areas where a higher mesh point density is required. There are a very limited number of moving mesh methods designed for solving field-theoretic partial differential equations, and the numerical analysis of the resulting schemes is challenging. In this thesis we present two ways to construct r-adaptive variational and multisymplectic integrators for (1+1)-dimensional Lagrangian field theories. The first method uses a variational discretization of the physical equations and the mesh equations are then coupled in a way typical of the existing r-adaptive schemes. The second method treats the mesh points as pseudo-particles and incorporates their dynamics directly into the variational principle. A user-specified adaptation strategy is then enforced through Lagrange multipliers as a constraint on the dynamics of both the physical field and the mesh points. We discuss the advantages and limitations of our methods. The proposed methods are readily applicable to (weakly) non-degenerate field theories---numerical results for the Sine-Gordon equation are presented.

In an attempt to extend our approach to degenerate field theories, in the last part of this thesis we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for variational integration. Our main observation is that the evolution takes place on the primary constraint and the 'Hamiltonian' equations of motion can be formulated as an index 1 differential-algebraic system. We then proceed to construct variational Runge-Kutta methods and analyze their properties. The general properties of Runge-Kutta methods depend on the 'velocity' part of the Lagrangian. If the 'velocity' part is also linear in the position coordinate, then we show that non-partitioned variational Runge-Kutta methods are equivalent to integration of the corresponding first-order Euler-Lagrange equations, which have the form of a Poisson system with a constant structure matrix, and the classical properties of the Runge-Kutta method are retained. If the 'velocity' part is nonlinear in the position coordinate, we observe a reduction of the order of convergence, which is typical of numerical integration of DAEs. We also apply our methods to several models and present the results of our numerical experiments.

正弦变化的振幅型光瞳滤波器

光瞳滤波器作为实现超分辨的基本元件之一，它的设计以及制作都非常重要。设计了一种正弦变化的振幅型光瞳滤波器，可以通过调节少量的参量实现各种不同的超分辨模式。选择正弦函数的周期以及光瞳中心点透过率变化两个参量来实现对最终超分辨效果的调节。数值计算结果表明：在整个人射光瞳上的透过率瞳函数分布具有0．5～2个正弦振幅周期时，较为合适。小于0．5个周期将不会有任何超分辨效果，大于2个周期超分辨效果反而变差。数值计算中还注意到，当周期数为整数时，斯特雷尔比将保持0．25不变。对正弦变化振幅型光瞳滤波器的计算结果，显示

Nematic ordering pattern formation in the process of self-organization of microtubules in a gravitational field

Papaseit et al. (Proc. Nati. Acad. Sci. U.S.A. 97, 8364, 2000) showed the decisive role of gravity in the formation of patterns by assemblies of microtubules in vitro. By virtue of a functional scaling, the free energy for MT systems in a gravitational field was constructed. The influence of the gravitational field on MT's self-organization process, that can lead to the isotropic to nematic phase transition, is the focus of this paper. A coupling of a concentration gradient with orientational order characteristic of nernatic ordering pattern formation is the new feature emerging in the presence of gravity. The concentration range corresponding to a phase coexistence region increases with increasing g or NIT concentration. Gravity facilitates the isotropic to nernatic phase transition leading to a significantly broader transition region. The phase transition represents the interplay between the growth in the isotropic phase and the precipitation into the nematic phase. We also present and discuss the numerical results obtained for local NIT concentration change with the height of the vessel, order parameter and phase transition properties.

Multiscale model reduction methods for deterministic and stochastic partial differential equations

Partial differential equations (PDEs) with multiscale coefficients are very difficult to solve due to the wide range of scales in the solutions. In the thesis, we propose some efficient numerical methods for both deterministic and stochastic PDEs based on the model reduction technique.

For the deterministic PDEs, the main purpose of our method is to derive an effective equation for the multiscale problem. An essential ingredient is to decompose the harmonic coordinate into a smooth part and a highly oscillatory part of which the magnitude is small. Such a decomposition plays a key role in our construction of the effective equation. We show that the solution to the effective equation is smooth, and could be resolved on a regular coarse mesh grid. Furthermore, we provide error analysis and show that the solution to the effective equation plus a correction term is close to the original multiscale solution.

For the stochastic PDEs, we propose the model reduction based data-driven stochastic method and multilevel Monte Carlo method. In the multiquery, setting and on the assumption that the ratio of the smallest scale and largest scale is not too small, we propose the multiscale data-driven stochastic method. We construct a data-driven stochastic basis and solve the coupled deterministic PDEs to obtain the solutions. For the tougher problems, we propose the multiscale multilevel Monte Carlo method. We apply the multilevel scheme to the effective equations and assemble the stiffness matrices efficiently on each coarse mesh grid. In both methods, the $\KL$ expansion plays an important role in extracting the main parts of some stochastic quantities.

For both the deterministic and stochastic PDEs, numerical results are presented to demonstrate the accuracy and robustness of the methods. We also show the computational time cost reduction in the numerical examples.

Polaron effects on the optical absorptions in asymmetrical semi-parabolic quantum wells

The linear and nonlinear optical absorptions considering the weak-coupling electron-LO-phonon interaction in asymmetrical semiparabolic quantum wells are theoretically investigated. The numerical results for the typical GaAs/AlxGa1-xAs material show that the factors of Al content x, the relaxation time and the photon energy have great influence on the optical absorption coefficients. Moreover, the theoretical values of the optical absorptions are more than a factor of 2-3 higher than the one in the structure without considering the electron-LO-phonon interaction by calculating. (C) 2007 Elsevier B.V. All rights reserved.

Valuation of wind energy projects: A real options approach

27 p.

Nonlinear dynamics of negatively refracted light in a resonantly absorbing Bragg reflector

The evolution of nonlinear light fields traveling inside a resonantly absorbing Bragg reflector is studied by use of Maxwell-Bloch equations. Numerical results show that a pulse initially resembling a light bullet may effectively experience negative refraction and anomalous dispersion in the resonantly absorbing Bragg reflector. (c) 2007 Optical Society of America.

The dispersive properties of an excited-doublet four-level atomic system

We have investigated the dispersive properties of excited-doublet four-level atoms interacting with a weak probe field and an intense coupling laser field. We have derived an analytical expression of the dispersion relation for a general excited-doublet four-level atomic system subject to a one-photon detuning. The numerical results demonstrate that for a typical rubidium D1 line configuration, due to the unequal dipole moments for the transitions of each ground state to double excited states, generally there exists no exact dark state in the system. Close to the two-photon resonance, the probe light can be absorbed orgained and propagate in the so-called superluminal form. This system may be used as an optical switch.

体全息光栅透镜的设计和应用

设计了一种新型的体全息光栅透镜，在一块光学平板(体全息记录材料)内可以将输入光束产生横向传输并聚焦，或对输入光点产生横传的准直．它由一束平面波和一束球面波正交入射到光学平板上干涉形成的．研究了该体全息透镜的光栅间距变化情况，为设计和制备体全息光栅透镜及相关器件提供了理论依据．基于两光束耦合波理论，得到了该光栅透镜的耦合波方程，近似计算了该透镜的衍射效率及其达到高衍射效率时透镜的最佳尺寸．最后，讨论了该透镜在集成光学等领域中的应用．

强外加电场与大调制度下光折变动力学光栅形成研究

强外加电场与大调制度在光折变效应的研究中已经得到了广泛应用。采用PDECOL算法, 严格求解光折变带输运方程, 得到外加电场时不同调制度下光折变晶体中随时间变化的空间电荷场、载流子浓度, 并讨论了外加电场对它们的影响。通过将物质方程与耦合波方程联立数值求解, 可得到光折变光栅形成过程中两波耦合增益系数以及光束条纹相位的变化。模拟结果表明, 在强外加电场作用下, 两束记录光之间的光强与相位耦合都得到了增强, 而原有的解析式忽视了强外加电场与大调制度对空间电荷场相位耦合的影响, 此时不再适用。同时发现折射率光

高阶群速度色散引起的高斯超短脉冲宽度的展宽和形变

基于动量守恒和光参变过程中的三波耦合波方程, 和负单轴非线性光学晶体CsLiB6O10的色散方程, 研究了在光参变效应中超短激光脉冲由于群速度色散引起的展宽和形变。数值模拟显示, 在超短脉冲波形为双曲正割形和无啁啾调制时, 高阶群速度色散引起的超短脉冲为50 fs时, 晶体长度为10 mm, 紫外光213 nm作为基波入射时的脉冲展宽是波长为532 nm绿光在同等条件下的1.6倍。脉冲展宽程度与入射波长和晶体长度有关, 波长越短和晶体长度越长则脉冲展宽和波形变化越严重,高阶色散引起的超短高斯脉冲展宽, 将

(Ce,Cu)∶LiNbO3晶体非挥发全息记录理论研究与优化

根据双中心带输运模型,对(Ce,Cu)∶LiNbO3晶体双中心非挥发全息记录进行了理论研究与优化。推导了(Ce,Cu)∶LiNbO3晶体的微观参量,采用数值方法通过严格求解模拟双中心带输运方程来模拟全息记录过程。分析了记录过程中,记录与敏化光强、Ce和Cu掺杂浓度以及晶体微观参量对(Ce,Cu)∶LiNbO3晶体双中心全息记录的影响。发现(Ce,Cu)∶LiNbO3晶体非挥发全息记录中实现高衍射效率与固定效率的主导因素是深中心Cu,在记录过程中,深中心Cu建立起了很强的空间电荷场。数值模拟的结果经过实验验