Franck-Condon simulation of the photoelectron spectrum of SO2- including Duschinsky effects

A theoretical method to calculate multidimensional Franck-Condon factors including Duschinsky effects is described and used to simulate the photoelectron spectrum of the anion SO. Geometry optimizations and harmonic vibrational frequency calculations have been performed on the XA(1) state of SO2 and (XB1)-B-2 state of SO2. Franck-Condon analyses and spectral simulation were carried out on the first photoelectron band of SO2. The theoretical spectra obtained by employing CCSD(T)/6-31 I+G(2d,p) values are in excellent agreement with the experiment. In addition, the equilibrium geometric parameters, r(c)(OS) = 0.1508 +/- 0.0005 nm and theta(e)(O-S-0) = 113.5 +/- 0.5 degrees, of the (XB1)-B-2 state of SO2, are derived by employing an iterative Franck-Condon analysis procedure in the spectral simulation. (c) 2005 Elsevier B.V. All rights reserved.

Numerical Investigation of Richtmyer-Meshkov Instability Driven by Cylindrical Shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov(RM) instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied.

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

Molecular Dynamics Simulation of Peeling a DNA Molecule on Substrate

Molecular dynamics (MD) simulations are performed to study adhesion and peeling of a short fragment of single strand DNA (ssDNA) molecule from a graphite surface. The critical peel-off force is found to depend on both the peeling angle and the elasticity of ssDNA. For the short ssDNA strand under investigation, we show that the simulation results can be explained by a continuum model of an adhesive elastic band on substrate. The analysis suggests that it is often the peak value, rather than the mean value, of adhesion energy which determines the peeling of a nanoscale material.

Experimental investigation on the chaotic phenomena in the wake of a natural thermal convection flow

Chaotic phenomena in the wake of thermal convection flow fields above a heating flat plate were investigated experimentally. A newly developed electron beam fluorescence technique (EBF) was used to simultaneously measure density fluctuation at 7 points in a cross section above the plate. Correlation dimensions, intermittence coefficients, Fourier spectrum have been obtained for different Grashof numbers. Spatial distribution of correlation dimensions are presented. The experimental result shows that there is a certain relationship between the density fluctuation and the Gr number. And time-spacial characteristic of chaos evolution is also given.

Novel Investigation on Principle of Similarity for Microgravity Two-Phase Flow

根据两流体同心环状流线性稳定性分析的结果 ,对微重力气 /液两相流地面模拟实验所应遵循的相似准则进行了探讨 ,得到了一个新的重力无关性准则 ,即Bond数和环形区流体相的毛细数之比的绝对值不大于 1 .此外 ,微重力气 /液两相流模拟实验还必须满足两个条件 ,即流量比和气相表观Weber数应与所模拟的流动中对应数值相等 . In the present paper, the principle of similarity for two phase flows at microgravity is studied based on the results of the linear stability analysis of the two fluid concentric annular flow configuration. A new criterion of gravity independence, namely the absolute value of the ratio between the Bond number and the capillary number of the phase flowing in the annulus is no more than one, is achieved. It is also pointed out that the flowrate ratio and the gas superficial Weber number must have the same ...

Molecular Dynamics Simulation of Barnacle Cement

Barnacle cement is an underwater adhesive that is used for permanent settlement. Its main components are insoluble protein complexes that have not been fully studied. In present article, we chose two proteins of barnacle cement for study, 36-KD protein and Mrcp-100K protein. In order to investigate the characteristic of above two proteins, we introduced the method of molecular modeling. And the simulation package GROMACS was used to simulate the behavior of these proteins. In this article, before the simulations, we introduce some theories to predict the time scale for polymer relaxation. During the simulation, we mainly focus on two properties of these two proteins: structural stability and adhesive force to substrate. First, we simulate the structural stability of two proteins in water, and then the stability of 36-KD protein in seawater environment is investigated.We find that the stability varies in the different environments. Next, to study adhesive ability of two proteins, we simulate the process of peeling the two proteins from the substrate (graphite). Then, we analyze the main reasons of these results. We find that hydrogen bonds in proteins play an important role in the protein stability. In the process of the peeling, we use Lennard–Jones 12-6 potential to calculate the van der Waals interactions between proteins and substrate.

Application of Three-Dimensional Discrete Element Face-to-Face Contact Model with Fissure Water Pressure to Stability Analysis of Landslide in Panluo Iron Mine

Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical phi numerically calculated is less than the one calculated by use of the limit equilibrium method for the same C. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.

Numerical Investigation On Detonation Cell Evolution In A Channel With Area-Changing Cross Section

The two-dimensional cellular detonation propagating in a channel with area-changing cross section was numerically simulated with the dispersion-controlled dissipative scheme and a detailed chemical reaction model. Effects of the flow expansion and compression on the cellular detonation cell were investigated to illustrate the mechanism of the transverse wave development and the cellular detonation cell evolution. By examining gas composition variations behind the leading shock, the chemical reaction rate, the reaction zone length, and thermodynamic parameters, two kinds of the abnormal detonation waves were identified. To explore their development mechanism, chemical reactions, reflected shocks and rarefaction waves were discussed, which interact with each other and affect the cellular detonation in different ways.

Numerical Investigation On Evolution Of Cylindrical Cellular Detonation

Cylindrical cellular detonation is numerically investigated by solving two-dimensional reactive Euler equations with a finite volume method on a two-dimensional self-adaptive unstructured mesh. The one-step reversible chemical reaction model is applied to simplify the control parameters of chemical reaction. Numerical results demonstrate the evolution of cellular cell splitting of cylindrical cellular detonation explored in experimentas. Split of cellular structures shows different features in the near-field and far-field from the initiation zone. Variation of the local curvature is a key factor in the behavior of cell split of cylindrical cellular detonation in propagation. Numerical results show that split of cellular structures comes from the self-organization of transverse waves corresponding to the development of small disturbances along the detonation front related to detonation instability.

Space Experimental Investigation On Thermocapillary Migration Of Bubbles

Results from a space experiment on bubble thermocapillary migration conducted on board the Chinese 22nd recoverable satellite were presented. Considering the temperature field in the cell was disturbed by the accumulated bubbles, the temperature gradient was corrected firstly with the help of the temperature measurement data at six points and numerical simulation. Marangoni number (Ma) of single bubble migrating in the space experiment ranged from 98.04 to 9288, exceeding that in the previous experiment data. The experiment data including the track and the velocity of two bubble thermocapillary migration showed that a smaller bubble would move slower as it was passed by a larger one, and the smaller one would even rest in a short time when the size ratio was large enough.

Wind Tunnel Experimental Investigation Of Sand Velocity In Aeolian Sand Transport

Sand velocity in aeolian sand transport was measured using the laser Doppler technique of PDPA (Phase Doppler Particle Analyzer) in a wind tunnel. The sand velocity profile, probability distribution of particle velocity, particle velocity fluctuation and particle turbulence were analyzed in detail. The experimental results verified that the sand horizontal velocity profile can be expressed by a logarithmic function above 0.01 in, while a deviation occurs below 0.01 m. The mean vertical velocity of grains generally ranges from -0.2 m/s to 0.2 m/s, and is downward at the lower height, upward at the higher height. The probability distributions of the horizontal velocity of ascending and descending particles have a typical peak and are right-skewed at a height of 4 turn in the lower part of saltation layer. The vertical profile of the horizontal RMS velocity fluctuation of particles shows a single peak. The horizontal RMS velocity fluctuation of sand particles is generally larger than the vertical RMS velocity fluctuation. The RMS velocity fluctuations of grains in both horizontal and vertical directions increase with wind velocity. The particle turbulence intensity decreases with height. The present investigation is helpful in understanding the sand movement mechanism in windblown sand transport and also provides a reference for the study of blowing sand velocity. (C) 2007 Elsevier B.V All rights reserved.

Investigation on mechanical properties and size effect of nanocrystalline twin copper

The stress-strain relations of nanocrystalline twin copper with variously sized grains and twins are studied by using FEM simulations based on the conventional theory of mechanism-based strain gradient plasticity (CMSG). A model of twin lamellae strengthening zone is proposed and a cohesive interface model is used to simulate grain-boundary sliding and separation. Effects of material parameters on stress-strain curves of polycrystalline twin copper are studied in detail. Furthermore, the effects of both twin lamellar spacing and twin lamellar distribution on the stress-strain relations are investigated under tension loading. The numerical simulations show that both the strain gradient effect and the material hardening increase with decreasing the grain size and twin lamellar spacing. The distribution of twin lamellae has a significant influence on the overall mechanical properties, and the effect is reduced as both the grain size and twin lamellar spacing decrease. Finally, the FEM prediction results are compared with the experimental data.

Atomistic investigation of the effects of temperature and surface roughness on diffusion bonding between Cu and Al

Molecular dynamics (MD) simulations are carried out to analyze the diffusion bonding at Cu/Al interfaces. The results indicate that the thickness of the interfacial layer is temperature-dependent, with higher temperatures yielding larger thicknesses. At temperatures below 750 K, the interface thickness is found to increase in a stepwise manner as a function of time. At temperatures above 750 K, the thickness increases rapidly and smoothly. When surface roughness is present, the bonding process consists of three stages. In the first stage, surfaces deform under stress, resulting in increased contact areas. The second stage involves significant plastic deformation at the interface as temperature increases, resulting in the disappearance of interstices and full contact of the surface pair. The last stage entails the diffusion of atoms under constant temperature. The bonded specimens show tensile strengths reaching 88% of the ideal Cu/Al contact strength. (c) 2007 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Numerical simulation of coherent structure in two-dimensional compressible mixing layers

The coherent structure in two-dimensional mixing layers is simulated numerically with the compressible Navier-Stokes equations. The Navier-Stokes equations are discretized with high-order accurate upwind compact schemes. The process of development of flow structure is presented: loss of stability, development of Kelvin-Helmholtz instability, rolling up and pairing. The time and space development of the plane mixing layer and influence of the compressibility are investigated.