高能量密度氘、氦的第一原理研究进展

极端条件下氘、氦的研究对于凝聚态物理、等离子体物理、天体物理以及惯性约束聚变的研究具有重要意义.文章综述了近年来温稠密区、热稠密区以及低温高压区的最新研究进展,归纳了尚待解决的科学问题,从而为进一步的研究提供参考.     

高能量密度氘、氦的第一原理研究进展

极端条件下氘、氦的研究对于凝聚态物理、等离子体物理、天体物理以及惯性约束聚变的研究具有重要意义.文章综述了近年来温稠密区、热稠密区以及低温高压区的最新研究进展,归纳了尚待解决的科学问题,从而为进一步的研究提供参考.     

Dynamical screening and ionic conductivity in water from ab initio simulations

We present a method to calculate ionic conductivities of complex fluids from ab initio simulations. This is achieved by combining density functional theory molecular dynamics simulations with polarization theory. Conductivities are then obtained via a Green-Kubo formula using time-dependent effective charges of electronically screened ions. The method is applied to two different phases of warm dense water. We observe large fluctuations in the effective charges; protons can transport effective charges greater than +e for ultrashort time scales. Furthermore, we compare our results with a simpler model of ionic conductivity in water that is based on diffusion coefficients. Our approach can be directly applied to study ionic conductivities of electronically insulating materials of arbitrary composition, e.g., complex molecular mixtures under such extreme conditions that occur deep inside giant planets.     

Anharmonic Raman spectra in high-pressure ice from ab initio simulations

We calculate from ab initio molecular dynamics the Raman scattering of high-pressure ice. To this effect we apply a new method based on the Berry phase theory of polarization. Our results are in agreement with recent and difficult experiments and are compatible with a picture in which ice VII is a proton-disordered system and in ice X the hydrogen bond is symmetric.     

Real-time ab initio simulations of excited carrier dynamics in carbon nanotubes

Combining time-dependent density functional calculations for electrons with molecular dynamics simulations for ions, we investigate the dynamics of excited carriers in a (3,3) carbon nanotube at different temperatures. Following an hnu=6.8 eV photoexcitation, the carrier decay is initially dominated by efficient coupling to electronic degrees of freedom. At room temperature, the excitation gap is reduced to nearly half its initial value after approximately 230 fs, where coupling to ionic motion starts dominating the decay. We show that the onset point and damping rate in the phonon regime change with initial ion velocities, a manifestation of temperature-dependent coupling between electronic and ionic degrees of freedom.     

Monte Carlo simulations of neon and argon using ab initio potentials

Gibbs ensemble Monte Carlo simulations of neon and argon have been performed with pair potentials taken from literature as well as with new ab initio potentials from just above the triple point to close to the critical point. The densities of the coexisting phases, their pair correlation functions, the vapour pressure and the enthalpy and entropy of vaporization have been calculated. The influence of the potential choice and of the addition of the Axilrod-Teller (AT) three-body potential on the above mentioned properties have been investigated. It turns out that an accurate ab initio two-body potential in connection with the AT potential yields very good results for thermodynamic properties of phase equilibria.     

Model of silica glass from combined classical and ab initio molecular-dynamics simulations

A Car-Parrinello (CP) molecular dynamics simulation of vitreous silica, combined with classical molecular dynamics, is presented. The equilibration of the liquid, quench and relaxation of the glass are performed classically using the van Beest, Kramer and van Santen (BKS) potential [12] and the resulting configuration (coordinates and velocities) is used as input for the CP simulation. A remarkable stability of the CP dynamics is observed justifying this procedure and validating the BKS potential. Several structural and electronic properties are calculated and compared with experiments. Received 23 March 1999 and Received in final form 15 July 1999     

Elastodynamic wave propagation in graded materials: simulations, experiments, phenomena, and applications

A two-dimensional numerical simulation model for the elastodynamic wave propagation in two linear elastic, isotropic, joint half-spaces is presented. The border between the two half-spaced is graded in a way, that the values of the elastic properties and the densities vary smoothly (sinusoidally) from the values of one continuum to the values of the other continuum within a transition zone of a defined thickness. It is demonstrated, that a graded layer leads to a frequency and wavelength dependent refraction and reflection behavior of elastodynamic waves. Numerical results show that wavelengths which are long compared with the transition layer thickness are dominantly reflected whereas short waves are dominantly transmitted, a phenomena which does not occur in the case of an infinitely thin transition layer. Furthermore the frequency dependent reflection and transmission behavior of elastodynamic waves is verified experimentally. There the interface between two vapor deposited films is graded due to intermetallic diffusion effects. These graded microstructures are analyzed with a short-pulse-laser-acoustic set-up. The corresponding frequencies of the elastodynamic waves which are filtered with these functionally graded microstructures are in the range of 0.5 THz.     

Experimental,semi-experimental and ab initio equilibrium structures

The determination of equilibrium structures of molecules by spectroscopic methods or by quantum mechanical calculations is reviewed. The following structures are described in detail: experimental equilibrium structures, empirical structures, semi-experimental structures and ab initio structures. The approximations made by the different methods are discussed and their accuracies are compared.     

Study of cerium phase transitions in shock wave experiments

Cerium has a complex phase diagram that is explained by the presence of structural phase transitions. Experiments to measure the sound velocities in cerium by two methods were carried out to determine the onset of cerium melting on the Hugoniot. In the pressure range 4–37 GPa, the sound velocity in cerium samples was measured by the counter release method using manganin-based piezoresistive gauges. In the pressure range 35–140 GPa, the sound velocity in cerium was measured by the overtaking release method using carbogal and tetrachloromethane indicator liquids. The samples were loaded with plane shock wave generators using powerful explosive charges. The onset of cerium melting on the Hugoniot at a pressure of about 13 GPa has been ascertained from the measured elastic longitudinal and bulk sound velocities.     

Cl-ClO催化循环圈的从头算研究

采用分子轨道从头算方法,在B3LYP/6-311 G(3df)和G2水平上研究了极地平流层臭氧损耗的一个基本过程.计算结果明显支持Cl-ClO催化循环圈机理,并且从能量角度解释了了臭氧破坏的基本原因.还对循环圈中各个反应的反应能,生成焓,相对吉布斯自由能做了计算,计算结果相互协调都说明了Cl-ClO催化循环圈破坏臭氧机理的正确性.     

Proceeding on the riddles of ketene pyrolysis by applying ab initio quantum chemical computational methods in a detailed kinetic modeling study

As ketene is a crucial intermediate for the high-temperature combustion of oxygenated hydrocarbons in general, an in-depth understanding of its chemistry is a fundamental requirement for the kinetic modeling of bio-based fuels. To gain a profound insight into the decomposition of ketene and subsequent reaction pathways high level ab initio methods were used. DSD-PBEP86/cc-pVTZ level of theory was applied for the geometries and frequencies, while single-point energies were determined at the CCSDT-1a level of theory extrapolated to the basis set limit. The reaction rate parameters for 38 reactions involved in the ketene chemistry including C1 to C4 species like acetylene, ethylene, propyne and allene were computed. For a total of 16 species, the thermochemistry were updated. The calculated rate parameters and the two new species cyclopropenone and 1,4-dioxo-1,3-butadiene were used to update the AramcoMech 3.0 base mechanism, which was then validated against speciation measurements during ketene pyrolysis. A reaction pathway analysis was performed to find the most prominent reactions at the investigated conditions and to discuss the simulation results. A significant improvement in the model’s prediction capability was found when applying the newly calculated reaction rate parameter and thermochemical data.     

Nuclear quantum effects in liquid water from path-integral simulations using an ab initio force-matching approach

We have applied path integral simulations, in combination with new ab initio based water potentials, to investigate nuclear quantum effects in liquid water. Because direct ab initio path integral simulations are computationally expensive, a flexible water model is parameterised by force-matching to density functional theory-based molecular dynamics simulations. Static and dynamic properties of liquid water at ambient conditions are presented and the role of nuclear quantum effects, exchange-correlation functionals and dispersion corrections are discussed in regards to reproducing the experimental properties of liquid water.     

Atomic structure and diffusivity in liquid Al80Ni20 by ab initio molecular dynamics simulations

The atomic structure and diffusivity in liquid Al₈₀Ni₂₀ are studied by ab initio molecular dynamics simulations. The local structures are analyzed by the pair correlation function, structure factor, coordinate number, Honneycutt–Anderson bond pair, and Voronoi tessellation methods. It is observed that the amount of icosahedral clusters increases, and the liquid becomes more ordered as the temperature decreases. The predicted self-diffusion coefficients of Al and Ni via the mean square displacements are very close to each other and agree well with the quasi-elastic neutron scattering measurements in the literature. The observation of equal self-diffusivity of Al and Ni is attributed to the formation of local solute-centered polyhedra, coupling the migration of Al and Ni. The Manning dynamic correlation factor is evaluated and found to be close to unity. The predicted interdiffusion coefficients using the Darken equation agree well with experimental data in the literature.     

Hydrogen bonding in liquid water: An ab initio QM/MM MD simulation study

Pattern and dynamics of hydrogen bonds in liquid water were investigated by a quantum mechanical/molecular mechanical molecular dynamics (QM/MM MD) simulation at Hartree–Fock (HF) level of theory. A large subregion of the whole system comprising two complete coordination shells was treated quantum mechanically in order to include all polarization and charge transfer effects and to obtain accurate data about structure and dynamics of the intermolecular bonds. The results of this investigation are in agreement with recent experimental findings and suggest that in liquid water every molecule forms in average 2.8, but almost as a rule less than four intermolecular hydrogen bonds.     

Artificial carrier heating due to the introduction of ab initio Coulomb scattering in Monte Carlo simulations

Accurately treating ionised impurity scattering in a way suitable to describe the influence of random dopant fluctuations on device characteristics is important in next generation MOSFETs. Statistical variations are unobservable using a continuous treatment of the doping, requiring a discrete representation of impurities. In particle-based simulations the P³M method, which resolves the Coulomb interaction into long- and short-range components, is in principle capable of describing Coulomb scattering through propagation in this accurately resolved potential. However, numerically the integration of the equations of motion is inaccurate and controlling the errors in practical simulations is vital. In this paper we investigate the effect of the choice of short-range correction strategy and integration time step on accuracy in a 3D self-consistent ensemble Monte Carlo simulations featuring random discrete dopants. We illustrate the importance of the ‘ab initio’ Coulomb scattering comparing the effect of a single trapped charge in drift-diffusion and Monte Carlo simulations.     

Intervalley scattering in GaAs: ab initio calculation of the effective parameters for Monte Carlo simulations

Interactions between excited electrons and short-wavelength (intervalley) phonons in GaAs are studied using density functional theory for the conduction bands, and density functional perturbation theory for phonon frequencies and matrix elements of the electron–phonon interaction. We have calculated the deformation potentials (DPs) and the average intervalley scattering time 〈τ〉. The integration of the scattering probabilities over all possible final states in the Brillouin zone has been performed without any ad hoc assumption about the behavior of the electron–phonon matrix elements nor the topology of the conduction band. For transitions from the L point to Γ valley (within the first conduction band), we find 〈τ〉_L to be 1.5 ps at 300 K, in good agreement with time-resolved photoluminescence experiment. We discuss the difference between our calculated DPs, and effective parameters used in Monte Carlo simulations of optical and transport properties of semiconductors. The latter are based on Conwell’s model, in which electron–phonon interaction is described by one single constant and a parabolic model is used for conduction bands. We deduce the effective DP from our 〈τ〉, and compare it to our calculated DPs. We conclude that only effective DPs obtained from a full calculation of 〈τ〉 are relevant parameters for Monte Carlo simulations. PACS 71.10-w; 72.10.Di; 71.55.Eq     

多孔脆性介质冲击波压缩破坏的细观机理和图像

本文采用一种具有良好定量性质的离散元模型研究了带孔洞的各向同性脆性介质在细观尺度上的压缩破坏特征. 通过对孤立孔洞、三种简单的孔洞排布方式和大量孔洞随机排布等几种情况的模拟, 认识到了剪切破坏和局域拉伸破坏是冲击波压缩下多孔介质的基本破坏模式; 孔洞之间的损伤贯通会促进孔洞在较低应力下发生塌缩, 但损伤区的应力松弛过程却会对一定范围内的介质起到损伤屏蔽作用; 不同区域中损伤促进和损伤屏蔽的综合效果是在多孔脆性介质中形成一种高损伤区与低损伤区间错排布的奇特损伤分布. 本文的研究结果为深入理解脆性材料冲击波压缩破坏的演化过程和机理提供了细观尺度上的初步物理图像.