等离子体湍流对电子的加速

本文讨论了等离子体湍流对电子加速的两种模型:（1）假定在空间中存在一个空间均匀的等离子体湍流区,当具有一定初始分布的电子束通过此湍流区时,研究湍流场对电子束的加速过程;（2）在某一封闭的区域中,存在着具有一定初始分布和空间均匀的等离子体,当某种类型的等离子体波突然传入此等离子体区,然后考察此区中电子的加速过程。在这两种模型中,可能存在着某种电子消失机制。假定湍谱是幂指数形式,我们给出了不同类型湍流扩散系数的普遍形式。利用较简单的数学方法,求解了包括消失过程的一维准线性动力学方程,对于给定的初始分布,得出了分布函数的解析解,并给出了平均能量时间关系的表达式。另外,对于特定的湍谱指数,解出了当平行电场和湍流同时存在时的分布函数。最后,对所得结果进行了数值分析和讨论。

近壁湍流脉动的概率分布函数

采用大涡模拟方法，模拟槽道湍流，获得了不同雷诺数情形下的槽道流大涡模拟数据库．在此基础上，获得了流向和垂向脉动速度的概率分布函数，并运用假设检验，分析了其与正态分布的定量差别．进一步计算了流向和垂向脉动速度的偏斜度、平坦度，讨论了二者在粘性子层、过渡区和对数律区的变化．同时，讨论了粘性子层、过渡区和对数律区流向和垂向脉动速度概率分布函数的特点及其与湍流猝发的高速流下扫和低速流喷发事件的关系．最后，分析了雷诺数对流向、垂向脉动速度分布的影响。

双晶晶界在高温下的热稳定性的分子动力学研究

晶界结构在高温下的热稳定性问题是一个长期争论而又未能解决的问题，其争论的焦点是：在远低于熔点的温度下，晶界结构是否发生了可观察到的无序化，即是否存在一个远低于熔点的结构转化温度。为了能澄清这一争论，本文系统地研究了晶界结构的热稳定性。为了消除相互作用势的影响和系统误差，本文首先采用Morse势和经验多体势分别对铝、铜单晶的熔化过程进行了分子动力学模拟。在平衡态下，通过计算表征结构无序化的静态结构因子、径向分布函数和单晶原子位形图，获得了铝、铜单晶的熔点，结果表明：多体势计算的铝和铜的单晶熔点更接近实验值。因此，采用经验多体势应用分子动力学方法分别模拟了铝、铜Σ3、Σ5、Σ9、Σ11、Σ19、Σ33六种对称倾侧双晶晶界晶界结构由有序向无序转化的过程，计算了平衡态下的表征结构无序化的静态结构因子、径向分布函数和晶界原子位形图并将多体势获得的铝、铜单晶熔点作为晶界结构转化温度的约化熔点，获得了铝、铜Σ3、Σ5、Σ9、Σ11、Σ19、Σ33六种对称倾侧双晶晶界结构的转化温度和熔点，结果表明：1．Σ5、Σ9、Σ11、Σ19、Σ33五种对称倾侧双晶晶界均在远低于单晶熔点温度时，晶界结构发生了可观察到的无序化，而且双晶晶界结构的转变温度相差不大，双晶晶界熔点也低于单晶熔点。2．Σ3晶界在温度远低于熔点时，其晶界结构没有发生可观察到的无序化；Σ3晶界的转化温度与单晶熔点接近。所以，可以认为Σ3晶界不存在转化温度。这是由于Σ3晶界为共格孪晶，具有较低的能量。综上所述，除Σ3共格孪晶外，在远低于熔点温度下，晶界结构发生了可观察到的无序化，即：存在一个远低于熔点的转化温度，此时其静态结构因子约为0.5左右；晶界结构的熔点均低于单晶熔点，此时其静态结构因子约为0.15左右。从全文模拟结果可以看出，静态结构因子、径向分布函数、晶界原子位形图三种方法在确定晶界的结构转化温度和熔点时，静态结构因子是最有效、最准确的定量方法。

Turbulent collision of inertial particles: point-particle based, hybrid simulations and beyond

Point-particle based direct numerical simulation (PPDNS) has been a productive research tool for studying both single-particle and particle-pair statistics of inertial particles suspended in a turbulent carrier flow. Here we focus on its use in addressing particle-pair statistics relevant to the quantification of turbulent collision rate of inertial particles. PPDNS is particularly useful as the interaction of particles with small-scale (dissipative) turbulent motion of the carrier flow is mostly relevant. Furthermore, since the particle size may be much smaller than the Kolmogorov length of the background fluid turbulence, a large number of particles are needed to accumulate meaningful pair statistics. Starting from the relative simple Lagrangian tracking of so-called ghost particles, PPDNS has significantly advanced our theoretical understanding of the kinematic formulation of the turbulent geometric collision kernel by providing essential data on dynamic collision kernel, radial relative velocity, and radial distribution function. A recent extension of PPDNS is a hybrid direct numerical simulation (HDNS) approach in which the effect of local hydrodynamic interactions of particles is considered, allowing quantitative assessment of the enhancement of collision efficiency by fluid turbulence. Limitations and open issues in PPDNS and HDNS are discussed. Finally, on-going studies of turbulent collision of inertial particles using large-eddy simulations and particle- resolved simulations are briefly discussed.

Large-eddy simulation of turbulent collision of heavy particles in isotropic turbulence

The small-scale motions relevant to the collision of heavy particles represent a general challenge to the conventional large-eddy simulation (LES) of turbulent particle-laden flows. As a first step toward addressing this challenge, we examine the capability of the LES method with an eddy viscosity subgrid scale (SGS) model to predict the collision-related statistics such as the particle radial distribution function at contact, the radial relative velocity at contact, and the collision rate for a wide range of particle Stokes numbers. Data from direct numerical simulation (DNS) are used as a benchmark to evaluate the LES using both a priori and a posteriori tests. It is shown that, without the SGS motions, LES cannot accurately predict the particle-pair statistics for heavy particles with small and intermediate Stokes numbers, and a large relative error in collision rate up to 60% may arise when the particle Stokes number is near St_K=0.5. The errors from the filtering operation and the SGS model are evaluated separately using the filtered-DNS (FDNS) and LES flow fields. The errors increase with the filter width and have nonmonotonic variations with the particle Stokes numbers. It is concluded that the error due to filtering dominates the overall error in LES for most particle Stokes numbers. It is found that the overall collision rate can be reasonably predicted by both FDNS and LES for St_K>3. Our analysis suggests that, for St_K<3, a particle SGS model must include the effects of SGS motions on the turbulent collision of heavy particles. The spectral analysis of the concentration fields of the particles with different Stokes numbers further demonstrates the important effects of the small-scale motions on the preferential concentration of the particles with small Stokes numbers.

Exact solutions to collisionless external gas flows

This paper presents exact density, velocity and temperature solutions for two problems of collisionless gas flows around a flat plate or a spherical object. At any point off the object, the local velocity distribution function consists of two pieces of Maxwellian distributions: one for the free stream which is characterized by free stream density, temperature and average velocity, n0, T0, U0; and the other is for the wall and it is characterized by density at wall and wall temperature, nw,Tw. Directly integrating the distribution functions leads to complex but exact flowfield solutions. To validate these solutions, we perform numerical simulations with the direct simulation Monte Carlo (DSMC) method. In general, the analytical and numerical results are virtually identical. The evaluation of these analytical solutions only requires less than one minute while the DSMC simulations require several days.

Nearly free molecular heat transfer from a sphere

Consider a sphere immersed in a rarefied monatomic gas with zero mean flow. The distribution function of the molecules at infinity is chosen to be a Maxwellian. The boundary condition at the body is diffuse reflection with perfect accommodation to the surface temperature. The microscopic flow of particles about the sphere is modeled kinetically by the Boltzmann equation with the Krook collision term. Appropriate normalizations in the near and far fields lead to a perturbation solution of the problem, expanded in terms of the ratio of body diameter to mean free path (inverse Knudsen number). The distribution function is found directly in each region, and intermediate matching is demonstrated. The heat transfer from the sphere is then calculated as an integral over this distribution function in the inner region. Final results indicate that the heat transfer may at first increase over its free flow value before falling to the continuum level.

适用于中等耦合等离子体的碰撞算子

完全电离等离子体中，当试探粒子分布函数fα是关于试探粒子速度vα的偶函数时，导出了一个新的动力学方程的碰撞算子．该碰撞算子同时包括了大角散射（库仑近碰撞）和小角散射（库仑远碰撞）的二体碰撞的贡献，因此，该碰撞算子同时适用于弱耦合（库仑对数ln∧≥10）和中等耦合（库仑对数2≤ln∧≤10）等离子体．而且经过修改的碰撞算子和Rosenbluth势有直接的联系，当试探粒子和场粒子满足条件mα＜mβ（如电子-离子碰撞或Lorentz气体模型）和｜vα｜〉｜vβ｜时，经约化的电子-离子碰撞算子同最初的Fokker

Direct visualization of the dynamics of membrane-anchor proteins in living cells

Dynamic properties of proteins have crucial roles in understanding protein function and molecular mechanism within cells. In this paper, we combined total internal reflection fluorescence microscopy with oblique illumination fluorescence microscopy to observe directly the movement and localization of membrane-anchored green fluorescence proteins in living cells. Total internal reflect illumination allowed the observation of proteins in the cell membrane of living cells since the penetrate depth could be adjusted to about 80 nm, and oblique illumination allowed the observation of proteins both in the cytoplasm and apical membrane, which made this combination a promising tool to investigate the dynamics of proteins through the whole cell. Not only individual protein molecule tracks have been analyzed quantitatively but also cumulative probability distribution function analysis of ensemble trajectories has been done to reveal the mobility of proteins. Finally, single particle tracking has acted as a compensation for single molecule tracking. All the results exhibited green fluorescence protein dynamics within cytoplasm, on the membrane and from cytoplasm to plasma membrane.

Vibrational pooling and constrained equilibration on surfaces

In this thesis, we provide a statistical theory for the vibrational pooling and fluorescence time dependence observed in infrared laser excitation of CO on an NaCl surface. The pooling is seen in experiment and in computer simulations. In the theory, we assume a rapid equilibration of the quanta in the substrate and minimize the free energy subject to the constraint at any time t of a fixed number of vibrational quanta N(t). At low incident intensity, the distribution is limited to one- quantum exchanges with the solid and so the Debye frequency of the solid plays a key role in limiting the range of this one-quantum domain. The resulting inverted vibrational equilibrium population depends only on fundamental parameters of the oscillator (ω_e and ω_eχ_e) and the surface (ω_D and T). Possible applications and relation to the Treanor gas phase treatment are discussed. Unlike the solid phase system, the gas phase system has no Debye-constraining maximum. We discuss the possible distributions for arbitrary N-conserving diatom-surface pairs, and include application to H:Si(111) as an example.

Computations are presented to describe and analyze the high levels of infrared laser-induced vibrational excitation of a monolayer of absorbed ¹³CO on a NaCl(100) surface. The calculations confirm that, for situations where the Debye frequency limited n domain restriction approximately holds, the vibrational state population deviates from a Boltzmann population linearly in n. Nonetheless, the full kinetic calculation is necessary to capture the result in detail.

We discuss the one-to-one relationship between N and γ and the examine the state space of the new distribution function for varied γ. We derive the Free Energy, F = NγkT − kTln(∑P_n), and effective chemical potential, μn ≈ γkT, for the vibrational pool. We also find the anti correlation of neighbor vibrations leads to an emergent correlation that appears to extend further than nearest neighbor.

Random propagation in complex systems : nonlinear matrix recursions and epidemic spread

This dissertation studies long-term behavior of random Riccati recursions and mathematical epidemic model. Riccati recursions are derived from Kalman filtering. The error covariance matrix of Kalman filtering satisfies Riccati recursions. Convergence condition of time-invariant Riccati recursions are well-studied by researchers. We focus on time-varying case, and assume that regressor matrix is random and identical and independently distributed according to given distribution whose probability distribution function is continuous, supported on whole space, and decaying faster than any polynomial. We study the geometric convergence of the probability distribution. We also study the global dynamics of the epidemic spread over complex networks for various models. For instance, in the discrete-time Markov chain model, each node is either healthy or infected at any given time. In this setting, the number of the state increases exponentially as the size of the network increases. The Markov chain has a unique stationary distribution where all the nodes are healthy with probability 1. Since the probability distribution of Markov chain defined on finite state converges to the stationary distribution, this Markov chain model concludes that epidemic disease dies out after long enough time. To analyze the Markov chain model, we study nonlinear epidemic model whose state at any given time is the vector obtained from the marginal probability of infection of each node in the network at that time. Convergence to the origin in the epidemic map implies the extinction of epidemics. The nonlinear model is upper-bounded by linearizing the model at the origin. As a result, the origin is the globally stable unique fixed point of the nonlinear model if the linear upper bound is stable. The nonlinear model has a second fixed point when the linear upper bound is unstable. We work on stability analysis of the second fixed point for both discrete-time and continuous-time models. Returning back to the Markov chain model, we claim that the stability of linear upper bound for nonlinear model is strongly related with the extinction time of the Markov chain. We show that stable linear upper bound is sufficient condition of fast extinction and the probability of survival is bounded by nonlinear epidemic map.

Test particle motion in a Lorentz gas

This is a two-part thesis concerning the motion of a test particle in a bath. In part one we use an expansion of the operator PLe^it(1-P)LLP to shape the Zwanzig equation into a generalized Fokker-Planck equation which involves a diffusion tensor depending on the test particle's momentum and the time.

In part two the resultant equation is studied in some detail for the case of test particle motion in a weakly coupled Lorentz Gas. The diffusion tensor for this system is considered. Some of its properties are calculated; it is computed explicitly for the case of a Gaussian potential of interaction.

The equation for the test particle distribution function can be put into the form of an inhomogeneous Schroedinger equation. The term corresponding to the potential energy in the Schroedinger equation is considered. Its structure is studied, and some of its simplest features are used to find the Green's function in the limiting situations of low density and long time.

Kinetic theory of normal quantum fluids

Close to equilibrium, a normal Bose or Fermi fluid can be described by an exact kinetic equation whose kernel is nonlocal in space and time. The general expression derived for the kernel is evaluated to second order in the interparticle potential. The result is a wavevector- and frequency-dependent generalization of the linear Uehling-Uhlenbeck kernel with the Born approximation cross section.

The theory is formulated in terms of second-quantized phase space operators whose equilibrium averages are the n-particle Wigner distribution functions. Convenient expressions for the commutators and anticommutators of the phase space operators are obtained. The two-particle equilibrium distribution function is analyzed in terms of momentum-dependent quantum generalizations of the classical pair distribution function h(k) and direct correlation function c(k). The kinetic equation is presented as the equation of motion of a two -particle correlation function, the phase space density-density anticommutator, and is derived by a formal closure of the quantum BBGKY hierarchy. An alternative derivation using a projection operator is also given. It is shown that the method used for approximating the kernel by a second order expansion preserves all the sum rules to the same order, and that the second-order kernel satisfies the appropriate positivity and symmetry conditions.

Influence of composition on the structure, electric and magnetic properties of Pd-Mn-P and Pd-Co-P amorphous alloys

The influence of composition on the structure and on the electric and magnetic properties of amorphous Pd-Mn-P and Pd-Co-P prepared by rapid quenching techniques were investigated in terms of (1) the 3d band filling of the first transition metal group, (2) the phosphorus concentration effect which acts as an electron donor and (3) the transition metal concentration.

The structure is essentially characterized by a set of polyhedra subunits essentially inverse to the packing of hard spheres in real space. Examination of computer generated distribution functions using Monte Carlo random statistical distribution of these polyhedra entities demonstrated tile reproducibility of the experimentally calculated atomic distribution function. As a result, several possible "structural parameters" are proposed such as: the number of nearest neighbors, the metal-to-metal distance, the degree of short-range order and the affinity between metal-metal and metal-metalloid. It is shown that the degree of disorder increases from Ni to Mn. Similar behavior is observed with increase in the phosphorus concentration.

The magnetic properties of Pd-Co-P alloys show that they are ferromagnetic with a Curie temperature between 272 and 399°K as the cobalt concentration increases from 15 to 50 at.%. Below 20 at.% Co the short-range exchange interactions which produce the ferromagnetism are unable to establish a long-range magnetic order and a peak in the magnetization shows up at the lowest temperature range . The electric resistivity measurements were performed from liquid helium temperatures up to the vicinity of the melting point (900°K). The thermomagnetic analysis was carried out under an applied field of 6.0 kOe. The electrical resistivity of Pd-Co-P shows the coexistence of a Kondo-like minimum with ferromagnetism. The minimum becomes less important as the transition metal concentration increases and the coefficients of ℓn T and T^2 become smaller and strongly temperature dependent. The negative magnetoresistivity is a strong indication of the existence of localized moment.

The temperature coefficient of resistivity which is positive for Pd- Fe-P, Pd-Ni-P, and Pd-Co-P becomes negative for Pd-Mn-P. It is possible to account for the negative temperature dependence by the localized spin fluctuation model and the high density of states at the Fermi energy which becomes maximum between Mn and Cr. The magnetization curves for Pd-Mn-P are typical of those resulting from the interplay of different exchange forces. The established relationship between susceptibility and resistivity confirms the localized spin fluctuation model. The magnetoresistivity of Pd-Mn-P could be interpreted in tenns of a short-range magnetic ordering that could arise from the Rudennan-Kittel type interactions.

Structure and properties of an amorphous ferromagnetic alloy

The structure and the electrical and magnetic properties of an amorphous alloy containing approximately 80 at .% iron, 13 at.% phos phorus and 7 at.% carbon (Fe_(80)Fe_(13)C_7) obtained by rapid quenching from the liquid state have been studied. Transmission electron diffraction data confirm the amorphous nature of this alloy. An analysis of the radial distribution function obtained from X-ray diffraction data indicates that the number of nearest neighbors is approximately seven, at a distance of 2.6A. The structure of the alloy can be related to that of silicate glasses and is based on a random arrangement of trigonal prisms of Fe_2P and Fe_3C types in which the iron atoms have an average ligancy of seven. Electrical resistance measurements show that the alloys are metallic. A minimum in the electrical resistivity vs. temperature curve is observed between 10° K to 50° K depending on the specimen, and the temperature at which the minimum occurs is related to the degree of local ordering. The Fe-P-C amorphous alloys are ferromagnetic. The Curie temperature measured by the induction method and by Mossbauer spectroscopy is 315° C. The field dependence of the magneto-resistance at temperatures from liquid helium to room temperature is similar to that found in crystalline iron. The ordinary Hall coefficient is approximately 10^(-11) volt-cm/amp-G. The spontaneous Hall coefficient is about 0.6 x 10^(-9) volt-cm/amp-G and is practically independent of temperature from liquid helium temperature up to 300° c.