基于Abel逆变换的激光诱导等离子体辐射特性研究

激光诱导击穿光谱(LIBS)作为一种新型的物质成分测量方法已经在越来越多的领域得到广泛应用，但是与传统的分析方法相比，LIBS技术的分析性能还需进一步提高。LIBS技术的理论基础是激光诱导等离子体，从物理机理上研究等离子体特性，对LIBS系统实验参数的优化具有指导作用，也为提高LIBS技术的检测能力奠定理论基础。激光诱导等离子体是一个与空间相关的非稳态辐射源，空间分辨光谱测量是探究等离子体物理特性的重要手段之一。为研究激光诱导等离子体的辐射特性，采用1 064 nm的Nd∶YAG调Q固体激光器烧蚀合金钢样品产生等离子体，利用空间分辨装置测量二维空间的等离子体辐射光谱信号，通过分析可知实验采集的光谱信号是信号探测器测量路径上的积分光谱强度，由此计算得到的等离子体参数也是观测路径上的平均值。为了深入研究等离子体由内层到外层的辐射规律，首先测量得到等离子体路径积分光谱强度的横向空间分布，然后以等离子体为光学薄和圆柱对称的前提条件，采用三次样条函数算法对路径积分光谱强度进行Abel逆变换，反演得到等离子体由内层到外层谱线辐射率的径向空间分布。选取等离子体辐射光谱中的原子谱线Fe Ⅰ: 374.55 nm和Mn Ⅰ: 403.08 nm为研究对象，分析等离子体辐射光谱的空间分布特征，研究结果表明，等离子体辐射路径积分光谱强度的横向分布呈现出中心位置强度大边缘位置强度小的特征，这是由于等离子体膨胀扩张的结果引起的；通过Abel逆变换得到等离子体光谱辐射率的径向分布，结果表明等离子体从内层到外层谱线的辐射率经过了先增加后减小的变化规律，等离子体中心处出现辐射率的极小值，造成这种现象的主要原因是由于等离子体辐射源中心区域具有较低的电子密度；选取等离子体辐射光谱中Fe元素的11条原子谱线，采用Boltzman法分别由谱线相应的积分光谱强度和辐射率计算等离子体温度，得到等离子体温度的横向空间和径向空间的二维分布，两者具有类似的变化规律；由等离子体温度的横向空间分布可以看出，随着离样品表面距离的增加，等离子体温度呈现单调减小的趋势，等离子体中心到边缘区域等离子体温度逐渐降低，这是由等离子体膨胀扩张以及与环境气体相互作用共同的结果；由等离子体温度的径向空间分布可以看出等离子体由内层到外层等离子体温度逐渐降低，这是由于等离子体膨胀扩张冷却引起的。由此可见，采用Abel逆变换能够实现等离子体由内层到外层的辐射特性分析，为深入理解等离子体产生和演变的物理机理提供实验依据，从而为提高激光诱导击穿光谱技术的分析性能奠定理论基础。     

GaAlAs-GaAs半导体激光器增益参数的研究

通过对半导体激光器发射光谱的计算可获得激光器的增益谱.本文研究了不同电流注入下激光器的增益特性;激光器的峰值增益系数Gmax在阈值电流以下随注入电流的增大而提高,也随结温的升高而下降.     

杨-巴克斯特自旋1/2链模型的量子关联研究

量子系统各部分间的量子关联可以作为量子信息应用研究的基础资源. 而量子失协是度量量子关联大小的物理量. 由此研究杨-巴克斯特自旋1/2链模型的量子关联情况. 首先利用两个杨-巴克斯特方程的解得到相应的杨-巴克斯特自旋1/2链模型. 然后, 计算分析热平衡时杨-巴克斯特自旋1/2链模型的量子失协、几何量子失协和量子纠缠随着温度和外磁场的变化情况. 结果表明对于杨-巴克斯特自旋1/2链模型, 量子失协和几何量子失协能够比量子纠缠更好地度量量子关联.     

溶液中球形晶体溶解的分析

从物质溶液浓度变化角度考虑了球形晶体在溶液中溶解随表面张力的变化,利用渐近方法求出了在溶液中球形晶体溶解的浓度和界面的近似解析解,能够计算出溶解过程中球形晶体溶解的浓度、界面演化形态.研究了表面张力对于溶液中溶质浓度分布、球晶界面形态和溶解速度的影响.结果表明,表面张力促进了球形晶体在溶液中的溶解.随着表面张力参数增大,溶液中在界面前沿的溶质浓度升高,球形晶体的溶解速度增大;随着时间的增加,溶解速度逐渐变大,球形晶体半径逐渐变小,直至溶解结束.     

Kinetic equation for electrons of a semiconductor in a strong electric field: Green's resolvent

Using the Laplace transform of Boltzmann's equation, a system of linking equations is obtained for the series expansion of the mapping of the distribution function in powers of the field. Since the form of the equations is identical, instead of direct integration it is possible to use the method of Green's functions for the solution of the system. Using the method of perturbances Green's resolvent is constructed in the form of an infinite series. Using the example of the evaluation of the distribution function in a high-frequency field, the convergence of the proposed method is demonstrated.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, Vol. 16, No. 7, pp. 22–25, July, 1973.     

复合场下优化产生粒子对能量分布宽度的特性研究

采用计算量子场论的方法,对振荡场加稳恒场的组合外场下真空中正反粒子对的产生特性进行了研究.通过一系列的对比得到当振荡场的宽度减小时,一方面可增加正反粒子对的产生量,另一方面也可减小正反粒子对的能量分布宽度从而得到能量单一性更好的粒子对.同时,通过分析产生量、能量分布宽度与振荡场宽度的关系可得出,仅在一定范围内减小振荡场的宽度可使能量分布更加集中,则能量分布宽度趋于某个极限值.因此,要得到产生量多且能量分布集中的正反粒子对应选择合适的参数,这可为今后的实验设计提供数据参考.     

圆形高斯镜平凹腔腔内外激光场研究

运用边界元素法把圆形高斯镜平凹腔的自再现模的衍射积分方程转化为有限阶矩阵方程。计算了圆形高斯镜取不同反射率膜斑半径和中心振幅反射率情况下基模的场强分布、相位分布和本征值。研究表明:腔内光场分布半径随反射率膜斑半径的减小而增大,中心振幅反射率不影响光场分布半径;反射率膜斑半径影响远场分布,当其数值较小时,在光束的远场分布主峰周围产生弱衍射环;当高斯镜中心全反射时,远场分布随反射率膜斑半径的变化无明显变化;模式本征值随反射率膜斑半径的增大而增大。     

Simulation of the destruction of a woodland massif by the shock wave generated by the impact of a large cosmic body

The interaction of shock waves with woodland massifs with consideration given to vegetation destruction is examined. The results of numerical simulations of the propagation of a perturbation through a vegetation massif based on solving the Euler equations for two-dimensional and axisymmetric cases are presented. This approach was used by the present author in simulating the propagation of wave fronts through woodland massifs without regard for forest stand destruction. The effect of vegetation is reproduced by introducing bulk forces in the volume occupied by forest. An initial perturbation of explosion-in-flight type was specified in the form of an increased-pressure region located at various heights from the ground surface. Simple threshold models are used to examine the effect of the dynamic pressure on vegetation destruction. A refined model was developed to investigate the change in the configuration of woodland massifs with consideration given to vegetation fall and the effect of this process on the rate of perturbation dissipation. The system of equations was solved by the Godunov method implemented using parallel programming. Numerical experiments were conducted with the use of two types of models of forest stand destruction, with various threshold values of the destruction factor, various initial perturbation intensities, and various heights above the vegetation massif. The effects arising during the interaction of a high-intensity perturbation with a woodland massif were identified. An analysis of the results showed that the destruction of vegetation occurs during its interaction with the fronts of the incident and reflected shock wave and with the jet core of the perturbation. In some cases, the model taking into account the accumulation of fallen vegetation predicts situations where destruction will occur only in the upper part of the forest stand, without touching the canopy. This effect can explain the existence of such regions in the area of the Tunguska event. It was also demonstrated that estimates of the sizes of the destruction zones obtained with the use of the models taking into account the influence of vegetation on the perturbation and destructibility differ substantially from those obtained within the framework of the models without regard for vegetation resistance. The highest range of perturbation dissipation was predicted by the models simulating the formation of obstructions.     

Dynamics of solidification in binary melts during rapid cooling I. The basic equations

In this paper we deal with the analysis of the general complex model which describes the solidification of the binary melt. Within the framework of this model we consider the mass and energy transport in the system and the kinetics of the phase transformation. The paper is divided into two parts. In the first part the set of equations is derived, which describe the evolution of the system within the framework of the theory of stochastic processes. The emphasis is laid on the rapid changes of temperature in the system. In the second part the analysis of the mentioned set of equations is presented. It results in the delimitation of the conditions under which the adiabatic approximation and the kinetic phase diagrams at the phase interface can be applied.     

废旧纺织品近红外光谱定量分析的新模型

根据废旧纺织品所含成分对它们做分类回收和处理可节省大量纺织原材料。目前，在废旧纺织品的回收过程中往往使用人工分拣方法。这种方法成本高且效率低。近红外光谱分析是21世纪发展最迅速的技术之一，可以在不破坏样本的情况下快速测定样本的成分及每种成分的含量。利用该技术对废旧纺织品进行分析，预先判断废旧纺织品所含的成分及各种成分的含量，可为废旧纺织品的大规模精细分类回收提供帮助。多模型方法通过将各子模型的预测值做加权平均得到最终的预测值，用该方法建立的近红外光谱分析模型一般具有较好的稳定性。以废旧纺织品样本的锦纶含量为例，先用多模型方法建立了锦纶含量的近红外光谱分析模型。方法如下：将反射率向量按照波长划分为15组。用每组数据建立一个近红外光谱分析子模型。对子模型的预测值做加权平均得出锦纶含量的最终预测值。然后在多模型方法基础上，根据锦纶含量预测值与实验值之间的近似线性关系，通过用变量代替常量并对变量做标准化处理，给出了一种便于优化的预测锦纶含量的近红外光谱分析新模型。优化后的每个子模型中的参数比优化前减少了6个，这样可防止模型过拟合。将上述两个模型与常见的用偏最小二乘法建立的模型进行了对比。交叉验证的结果表明：（优化后的）新模型的拟合优度的平均值为0.820 7，单纯使用多模型方法所建模型的拟合优度的平均值为0.769 1，用偏最小二乘法建立的模型的拟合优度的平均值为0.746 7。因此, 使用多模型方法建立的模型的预测效果好于用偏最小二乘法建立的模型的预测效果。新模型的预测效果明显好于其他两个模型的预测效果。该研究主要创新之处是新模型的建立和优化。文中建模方法有望用于废旧纺织品样本其他成分的含量预测。     

Determination of the mass spectrum of the bound state in the framework of the relativistic hamiltonian approach

We propose one a variant of calculation of the energy spectrum of bound state systems with relativistic corrections. In the framework of quantum field theory, an expression that takes into account relativistic corrections to the mass of the bound state with a known nonrelativistic pair interaction potential is proposed on the basis of calculating the asymptotic behavior of correlation functions of the corresponding field currents with the necessary quantum numbers. Excluding the time variables allows one to determine nonperturbative corrections to the interaction potential. The following results have been obtained in the framework of this approach. The nonperturbative corrections arising due to the relativistic nature of a system to the interaction Hamiltonian are determined. The dependence of the constituent mass of bound-state forming particles on the free state mass and on the orbital and radial quantum numbers is analytically derived. The energy level shift of muonic hydrogen taking into account relativistic corrections is calculated. The energy spectrum of a wide class of potentials that describe the Coulomb bound state is analytically derived with relativistic corrections. The mass spectrum of the glueballs and the constituent masses of the gluons are analytically calculated taking into account spin-orbit, spin—spin, and tensor interactions. Our numerical results have shown very good agreement with the lattice data. Taking into account nonperturbative and nonlocality characters of interactions, the mass spectrum of the mesons consisting of light-light and light-heavy quarks with orbital and radial excitations is determined. Our results show that good agreement with the experimental data for the slope and the intercept of the Regge trajectory can be obtained only taking into account the nonperturbative and the nonlocal characters of interactions. The dependences of the constituent masses of constituent particles on the masses of a free state are certain. When quarks are light, then the difference between current and valent masses of quarks is greater than valent masses of quarks, and when quarks are heavy, then the difference between these masses is insignificant. One of the alternative variants of taking nonlocality into account has been suggested for the definition of properties of hadrons at large distances. The dependence of the constituent masses of constituent particles on the radius of confinement is determined.     

Reconstruction of the characteristics of a light-scattering layer from its reflectance and transmittance: The neural network method

The problem of reconstructing the probability of photon survival, the optical layer thickness, and the parameter of asymmetry of the scattering indicatrix of an elementary volume from the dependence of the reflectance and the transmittance on the angle of incidence of radiation onto a layer of a medium with Fresnel reflection from the boundaries is considered. To solve this problem, a method based on the use of correlation neural networks is proposed. For training and control of the operation of neural networks, the reflectance and the transmittance of the layer calculated within the framework of the radiation transfer theory were used. For the simulation of the scattering characteristics of an elementary volume, the Mie theory was used. Estimates of errors in the reconstruction of the characteristics of the medium were performed.     

Entanglement of two atoms in two-mode Raman coupled model with intrinsic decoherence

Considering intrinsic decoherence, the two-atom two-mode Raman coupled model is investigated in this paper. Utilizing the constants of motion in this model, we obtain the analytic expressions of the density operator of the system for investigating the entanglement of two atoms. The speed of entanglement decay increases with the increasing of the coupling coefficient of one atom. The difference between the oscillation periods when the initial state parameter of atomic subsystem belongs to two intervals becomes smaller with the increasing of the coupling coefficient of one atom. The increasing of the initial photon number of the second field can hasten the vanishing of entanglement of atomic subsystem. The robustness of atomic entanglement against decoherence depends on the interval of the initial state parameter of atomic subsystem.     

Interpolation method for transient electromagnetic fields

Summary The method for the transient electromagnetic fields interpolation from the irregular to the uniform array is suggested. It is based on the approximation of the observed data by the field of the current plane disposed at some depth in the homogeneous conductive Earth. The interpolation problem is reduced to solving the Fredholm integral equation of the first kind for unknown distribution of the current in the plane. The discretization of the integral equation results in the systems of the linear equations which are solved by the method of the singular value decomposition. The distribution of the current in the plane is used not only for the interpolation of the observed field but also for the calculation of the linear transformations of the field components. The results of testing of the method on theoretical models are presented. To speed up publication, the proofs were not sent to the authors and were supervised by the Scientific Committee.     

Classical theory of homogeneous nucleation in superheated liquids and its experimental verification

This work summarizes the cycle of investigations devoted to determination of the character of superheated liquid boiling-up in a glass capillary at the boundary of the attainable superheating. The brief history of studies is followed by a brief introduction into the theory of homogeneous nucleation. The result of the targeted experiments determining stationarity of a random process of a supercritical embryo generation is considered. From the experiments it may be concluded that the process is unsteady. Based on a large selection of life spans of superheated liquids, the authors have made parametric and nonparametric evaluations of the functions of distribution and dependence of boiling-up frequency on time. The comparison of the obtained results with exponential distribution shows significant differences that also prove the nonstationarity of the studied random process. Special experiments and calculations for evaluation of homogeneity of the superheated liquid boilingup at the boundary of attainable superheating are considered. It is shown that in the experiments with glass capillary, the boiling-up occurs on the wall. As the most convincing evidence of heterogeneity of the superheated liquid boiling-up in a glass capillary the authors provide the results of high-speed video filming in a silylated and clean capillaries.     

Reconstruction of the form of complex motion of an object from the autodyne detection signal of a semiconductor laser

The possibility of reconstructing the shape of an interference signal on the basis of the autodyne detection signal of a semiconductor laser is demonstrated in cases of harmonic and nonharmonic vibrations of the object of study. The shape of the interference signal coincides with the shape of the signal of an interference system decoupled from the radiation source. It is shown that the form of nonharmonic vibrations can be reconstructed from the values of the autodyne signal function at four varying distances from the external vibrating reflector. Numerical simulation of the procedure of reconstruction of the shape of the interference signal was performed. The main limitations of the method under consideration are discussed.     

Numerical study on the splitting of a vapor bubble in the ultrasonic assisted EDM process with the curved tool and workpiece

Electrical discharge machining (EDM) is a powerful and modern method of machining. In the EDM process, a vapor bubble is generated between the tool and the workpiece in the dielectric liquid due to an electrical discharge. In this process dynamic behavior of the vapor bubble affects machining process. Vibration of the tool surface affects bubble behavior and consequently affects material removal rate (MRR). In this paper, dynamic behavior of the vapor bubble in an ultrasonic assisted EDM process after the appearance of the necking phenomenon is investigated. It is noteworthy that necking phenomenon occurs when the bubble takes the shape of an hour-glass. After the appearance of the necking phenomenon, the vapor bubble splits into two parts and two liquid jets are developed on the boundaries of the upper and lower parts of the vapor bubble. The liquid jet developed on the upper part of the bubble impinges to the tool and the liquid jet developed on the lower part of the bubble impinges to the workpiece. These liquid jets cause evacuation of debris from the gap between the tool and the workpiece and also cause erosion of the workpiece and the tool. Curved tool and workpiece affect the shape and the velocity of the liquid jets during splitting of the vapor bubble. In this paper dynamics of the vapor bubble after its splitting near the curved tool and workpiece is investigated in three cases. In the first case surfaces of the tool and the workpiece are flat, in the second case surfaces of the tool and the workpiece are convex and in the third case surfaces of the tool and workpiece are concave. Numerical results show that in the third case, the velocity of liquid jets which are developed on the boundaries of the upper and lower parts of the vapor bubble after its splitting have the highest magnitude and their shape are broader than the other cases.