Ab Initio Calculations of the Electronic Structures of Copper Pyrites CuS2, CuSe2 and CuTe2

The electronic structures of CuS2, CuSe2 and CuTe2 with pyrite structures, within the framework the density-functional theory have been investigated. The calculated results explained the recent experimental results which show that there is no clear indication of strong electron correlations in the electronic properties of Cu pyrites, due to the dominant chalcogen p character rather than d characteristic of Cu at the Fermi level.

Elastic and Optoelectronic Properties of KCdF3:ab initio Calculations through LDA/GGA/TB-mBJ within FP-LAPW Method

Ab initio calculations are performed on the electronic, structural, elastic and optical properties of the cubic per- ovskite KCdF3. Tile Kohn Sham equations are solved by applying the full potential linearized augmented plane wave （FP-LAPW） method. The exchange correlation effects are included through the local density approximation （LDA ）, generalized gradient approximation （GGA） and modified Becke-Johnson （mBJ） exchange potential The calculated lattice constant is in good agreement with the experimental result. The elastic properties such as elastic constants, anisotropy factor, shear modulus, Young＇s modulus and Poisson＇s ratio are calculated. KCdF3 is ductile and elastically anisotropic. The calculations of the electronic band structure, density of states （DOS） and charge density show that this compound has an indirect energy band gap （M-F） with a mixed ionic and covalent bonding. The contribution of the different bands is analyzed from the total and partial density of states curves. Optical response of the dielectric functions, optical reflectivity, absorption coefficient, real part of optical conductivity, refractive index, extinction coefficient and electron energy loss, are presented for the energy range of O-40eV. The compound KCdF3 can be used for high-frequency optical and optoelectronic devices.

OPTIMAL SCHEME FOR SEQUENTIAL COMPUTATIONS OF F_m(z) INTEGRALS IN AB INITIO CALCULATIONS----COMBINATORY USE OF UPWARD AND DOWNWARD RECURSIONS

The quantitative rules of the transfer and variation of errors,when the Gaussian integral functions F.(z) are evaluated sequentially by recurring,have been expounded.The traditional viewpoint to negate the applicability and reliability of upward recursive formula in principle is amended.An optimal scheme of upward-and downward-joint recursions has been developed for the sequential F(z) computations.No additional accuracy is needed with the fundamental term of recursion because the absolute error of Fn(z) always decreases with the recursive approach.The scheme can be employed in modifying any of existent subprograms for Fn<z> computations.In the case of p-d-f-and g-type Gaussians,combining this method with Schaad's formulas can reduce,at least,the additive operations by a factor 40%;the multiplicative and exponential operations by a factor 60%.

Ab initio calculations of accurate dissociation energy and analytic potential energy function for the second excited state B^1∏ of ^7LiH

The reasonable dissociation limit of the second excited singlet state B1∏ of ^7LiH molecule is obtained. The accurate dissociation energy and equilibrium geometry of the B^∏ state are calculated using a symmetry-adaptedcluster configuration interaction method in full active space. The whole potential energy curve for the B1H state is obtained over the internuclear distance ranging from about 0.10 nm to 0,54 nm, and has a least-square fit to the analytic Murrell-Sorbie function form. The vertical excitation energy is calculated from the ground state to the B^1∏ state and compared with previous theoretical results. The equilibrium internuclear distance obtained by geometry optimization is found to be quite different from that obtained by single-point energy scanning under the same calculation condition. Based on the analytic potential energy function, the harmonic frequency value of the B^1∏ state is estimated. A comparison of the theoretical calculations of dissociation energies, equilibrium interatomic distances and the analytic potential energy function with those obtained by previous theoretical results clearly shows that the present work is more comprehensive and in better agreement with experiments than previous theories, thus it is an improvement on previous theories.

Ab initio calculations on the spectroscopic constants, vibrational levels and classical turning points for the 2^1Пu state of dimer ^7Li2

The accurate dissociation energy and harmonic frequency for the highly excited 2^1Пu state of dimer ^7Li2 have been calculated using a symmetry-adapted-cluster configuration-interaction method in complete active space. The calculated results are in excellent agreement with experimental measurements. The potential energy curves at numerous basis sets for this state are obtained over a wide internuclear separation range from about 2.4a0 to 37.0a0. And the conclusion is gained that the basis set 6-311＋＋G（d,p） is a most suitable one. The calculated spectroscopic constants De, Re, ωe, ωeχe, ae and Be at 6-311＋＋G（d,p） are 0.9670 eV, 0.3125 nm, 238.6 cm^-1, 1.3705 cm^-1, 0.0039 cm^-1 and 0.4921 cm^-1, respectively. The vibrational levels are calculated by solving the radial SchrSdinger equation of nuclear motion. A total of 53 vibrational levels are found and reported for the first time. The classical turning points have been computed. Comparing with the measurements, in which only the first nine vibrational levels have been obtained so far, the present calculations are very encouraging. A careful comparison of the present results of the parameters De and We with those obtained from previous theories clearly shows that the present calculations are much closer to the measurements than previous theoretical results, thus representing an improvement on the accuracy of the ab initio calculations of the potentials for this state.

Using Raman Spectroscopy and ab initio Calculations to Investigate Intermolecular Hydrogen Bonds in Binary Mixture (Tetrahydrofuran+Water)

We analyzed the properties and structures of the hydrogen-bonded complexes of tetrahydrofuran（THF） and water by means of experimental Raman spectra and ab initio calculations.The optimized geometries and vibrational frequencies of the neat THF molecule and its hydrogen-bonded complexes with water（THF/H2O） were calculated at the MP2/6-311＋G（d,p） level of theory.We found that the intermolecular hydrogen bonds which are formed from the binary mixtures of the neat THF and water with different molar ratios could explain the changes in wavenumber position and linewidth very well.The combination of ab initio calculations and experimental Raman spectral data provides an insight into the hydrogen bonds leading to the concentration dependent changes in the spectral features.

Ground-state mass of ^(22)Al and test of state-of-the-art ab initio calculations

The ground-state mass excess of the T_(z)=−2 drip-line nucleus ^(22)Al is measured for the first time as 18103(10)keV using the newly-developed Bρ-defined isochronous mass spectrometry method at the cooler storage ring in Lanzhou.The new mass excess value allowed us to determine the excitation energies of the two low-lying 1+states in ^(22)Al with significantly reduced uncertainties of 51 keV.When compared to the analogue states in its mirror nucleus ^(22)F,the mirror energy differences of the two 1^(+)states in the ^(22)Al-^(22)F mirror pair are determined to be−625(51)keV and−330(51)keV.The excitation energies and mirror energy differences are used to test the state-of-the-art ab initio valence-space in-medium similarity renormalization group calculations with four sets of interactions derived from the chiral effective field theory.The mechanism leading to the large mirror energy differences is investigated and attributed to the occupation of theπs_(1/2) orbital.

Identifying the crystal structure of T 1 precipitates in Al-Li-Cu alloys by ab initio calculations and HAADF-STEM imaging

Controversial experimental reports on the crystal structure of T 1 precipitates in Al-Li-Cu alloys have ex-isted for a long time,and all of them can be classified into five models.To clarify its ground-state atomic structure,herein,we have combined high-throughput first-principles calculations and CALPHAD,as well as aberration-corrected HAADF-STEM experiments.Employing the special quasi-random structure(SQS)and supercell approximation(SPA)methods to simulate the local disorder on Al-Cu sub-lattices,we find that none of the present models can satisfy the phase stability in Al-Li-Cu ternary system based on temperature-dependent convex hull analysis.Using the cluster expansion(CE)formulas,structural predic-tions derived from the five-frame models were performed.Subsequently,by introducing the vibrational contribution to the free energy at aging temperatures,we proposed a novel ground-state T 1 structure that maintains a coherent relationship with Al-matrix at the<112>Al orientation.The underlying phase transition between the variants of T 1 precipitates was further discussed.By means of ab initio molecular dynamics(AIMD)simulations,we resolved the controversy regarding the number of atomic layers con-stituting the T 1 phase and acknowledged the existence of Al-Li corrugated layers.The root cause of this structural distortion is triggered by atomic forces and bondings.Our work can have an positive impact on the novel fourth generation of Al-Cu-Li alloy designs by engineering the T 1 strengthening phase.

Ab initio calculations on accurate dissociation energy, equilibrium geometry, and analytic potential energy function for the 6~3Π state of ~7LiH molecule

The accurate dissociation energy and equilibrium geometry of the ball state of ^7LiH molecule is calculated using a symmetry-adapted-cluster configuration-interaction method in full active space. And the calculated results are 0.2580 eV and 0.1958 nm for the dissociation energy and equilibrium geometry, respectively. The whole potential energy curve for the b^3∏ state is also calculated over the internuclear separation range from about 0.10 to 0.54 nm. The results are fitted by the Murrell-Sorbie function. It is found that the Murrell-Sorbie function form, which is mainly used to fit the ground-state potential energy function, is well suitable for the excited triplet b^3∏ state. The vertical excitation energy from the ground state to the b^3∏ state is calculated to be 4.233 eV. Based on the analytic potential energy function, the harmonic frequency of 610.88 cm^-1 about this state is firstly estimated. Compared with other theoretical results, this work is the most complete effort to deal with the analytic potential energy function and the harmonic frequency of this state.

Rotational Isomeric State Model of Poly(vinylidene chloride) Based Upon ab Initio Electronic Structure Calculations

The rotational isomeric state（RIS） model was constructed for poly（vinylidene chloride）（PVDC） based on quantum chemistry calculations. The statistical weighted parameters were obtained from RIS representations and ab initio energies of conformers for model molecules 2,2,4,4-tetrachloropentane（TCP） and 2,2,4,4,6, 6-hexachloroheptane（HCH）. By employing the RIS method, the characteristic ratio C∞ was calculated for PVDC. The calculated characteristic ratio for PVDC is in good agreement with experiment result. Additionally, we studied the influence of the statistical weighted parameters on C∞ by calculating δC∞/δlnw. According to the values of δC∞/δlnw, the effects of second-order Cl-CH2 pentane type interaction and C1--C1 long range interaction on C∞ were found to be important. In contrast, first-order interaction is unimportant.

A charge-conserving approximation——method for ab initio calculations

A new method is presented for approximate ab initio calculations in quantum chemistry. It is called CCAM (charge conserving approximation method). The calculation method does not include the use of empirical parameters. We use Slater type orbitals as basis set, replacing STO's by STO-2G functions to evaluate three- and four-center integrals and making the STO-2G two- orbital charge distributions have the same total charge as STO. The results are presented for test calculations on five molecules. In view of these results, CCAM is better than ab initio calculations over STO-6G in the results on total energies, kinetic energies and occupied orbital energies. In atomic populations, dipole moments and unoccupied orbital energies, CCAM is also satisfactory. We estimate that CCAM would be as fast as ab initio calculations over STO-2G in evaluating molecular integrals.

Studies on nitramide and its methyl derivatives with ab initio calculations——Ⅳ.The harmonic force field and vibrational spectra of dimethylnitramine and its isotopic derivatives

The complete harmonic vibrational force field of dimethylnitramine has been calculated at the Hartree-Fock level using 4-21G basis set.The harmonic force field was then scaled with scale factors previously derived from N-methylnitramine,and the vibrational spectrum of dimethylnitramine was computed.This apriori prediction,made with no reference to observations on dimethylnitramine, agrees with the experimental IR spectrum in gas phase with a mean deviation of 8.4 cm^(-1).Some of the scale factors were reoptimized by fitting of the computed force field to experimental data.The new set of scale factors reduced the mean deviation to 4.5 cm^(-1),and was used to predict the vibrational spectrum of deuterated form of dimethylnitramine(-6D).Dipole moment derivatives were also cal- culated and used to predict infrared intensities which are comparable with experimental values.

Ab initio calculations on icosahedral molecules——Ⅰ.The implementation of symmetry reduction in calculation exemplified by B_(12)H_(12)^(2-)

The symmetry reduction of icosahedral group has been carried out by means of projective operator method.The irreducible bases and matrices were tabulated for generators C_2,C_3 and C_5 which would be sufficient for the implementation of ab initio or other computations in terms of symmetry-adapted linear combinations of atomic orbitals.As an initial application,the borane dianion B_(12)H_(12)~2-with I_h,symmetry has been calculated with a flexible s,p set,(9s5p/4s)being contracted into[3s2p/2s].The electron counting rule was then analyzed.

Investigation on ab initio calculations of the electron donor-acceptor complex H_3N·SiH_3Cl and the hypervalence of Si in it

The formation and the presence of hypervalent Si in the electron donor-acceptor com- plex H_3N·SiH_3 Cl have been investigated by ab initio calculation. The results show that there is a 0.707eV decrease of energy when the complex H_3N.SiH_3 Cl is formed from NH_3 and H_3SiCl, the interaction potential between the donor NH_3 and the acceptor H_3SiCl belongs to the Morse type, and the bond angle A(H-Si-Cl) versus bond length d(N-Si) presents a linear relation. The results also show that the interaction is mainly from giving the lone pair electrons in HOMO of NH_3 to LUMO of H_3SiCl, in which the 2P_z of N and the 3d_0 of Si play important role. Bond N-Si is a weak n-σ*type dative one.

Formation and dissociation of protonated cytosine–cytosine base pairs in i-motifs by ab initio quantum chemical calculations

Formation and dissociation mechanisms of C-C＋ base pairs in acidic and alkaline environments are investigated, employing ab initio quantum chemical calculations. Our calculations suggest that, in an acidic environment, a cytosine monomer is first protonated and then dimerized with an unprotonated cytosine monomer to form a C-C＋ base pair; in an alkaline environment, a protonated cytosine dimer is first unprotonated and then dissociated into two cytosine monomers. In addition, the force for detaching a C-C＋ base pair was found to be inversely proportional to the distance between the two cytosine monomers. These results provide a microscopic mechanism to qualitatively explain the experimentally observed reversible formation and dissociation of i-motifs.

Investigation of hydrogen bonding in neat dimethyl sulfoxide and binary mixture(dimethyl sulfoxide + water) by concentration-dependent Raman study and ab initio calculation

In this study, our vibrational spectroscopic analysis is made on hydrogen-bonding between dimethyl sulfoxide and water comprises both experimental Raman spectra and ab initio calculations on structures of various dimethyl sulfoxide/water clusters with increasing water content. The Raman peak position of the v（S=O） stretching mode of dimethyl sulfoxide serves as a probe for monitoring the degree of hydrogen-bonding between dimethyl sulfoxide and water. In addition, the two vibrational modes, namely, the CH3 symmetric stretching mode and the CH3 asymmetric stretching mode have been analysed under different concentrations. We relate the computational results to the experimental vibrational wavenumber trends that are observed in our concentration-dependent Raman study. The combination of experimental Raman data with ab initio calculation leads to a better knowledge of the nature of the hydrogen bonding and the structures of the hydrogen-bonded complexes studied.

Ab Initio Theoretical Investigation on the Geometrical and Electronic Structures of AlAu_n^(-/0)(n=2-4) Clusters

A systematic density functional theory and wave function theory investigation on the geometrical and electronic properties of A1Aun0/＂ （n = 2-4） clusters has been performed in this work. A1Aun- anions prove to possess ground states of the V-shaped C2v A1Au2, umbrella-shaped C3v A1Au3, and perfect tetrahedral Td A1Au4＂, while their neutrals favor the V-shaped CEv A1Au2, perfect planar triangular D3h A1Au3, and severely distorted Cs A1Au4, respectively. Aluminum aurides appear to be analogous to the corresponding aluminum hydrides, expect C~ A1Au4. Molecular orbitals （MOs） analyses also support this conclusion. Detailed orbital analyses indicate that Au 6s makes 94-96% and Au 5d makes 6-4% contribution to the Au-based orbitals in A1-Au bonds, which is smaller than the BAun0/- series, partially reflecting the relativistic effect of gold. The one-electron detachment energies of the anions and characteristic stretching vibrational frequencies of A1-Au bonds between 100-400 cm-1 have been calculated to facilitate future experimental characterization of these clusters.

Ab initio path-integral molecular dynamics and the quantum nature of hydrogen bonds

The hydrogen bond （HB） is an important type of intermolecular interaction, which is generally weak, ubiquitous, and essential to life on earth. The small mass of hydrogen means that many properties of HBs are quantum mechanical in nature. In recent years, because of the development of computer simulation methods and computational power, the influence of nuclear quantum effects （NQEs） on the structural and energetic properties of some hydrogen bonded systems has been intensively studied. Here, we present a review of these studies by focussing on the explanation of the principles underlying the simulation methods, i.e., the ab initio path-integral molecular dynamics. Its extension in combination with the thermodynamic integration method for the calculation of free energies will also be introduced. We use two examples to show how this influence of NQEs in realistic systems is simulated in practice.

Ab initio calculation of accurate dissociation energy, potential energy curve and dipole moment function for the A^1∑＋ state ^7LiH molecule

The reasonable dissociation limit of the A^1∑＋ state ^7LiH molecule is obtained. The accurate dissociation energy and the equilibrium geometry of this state are calculated using a symmetry-adapted-cluster configuration-interaction method in complete active space for the first time, The whole potential energy curve and the dipole moment function for the A^1∑＋ state are calculated over a wide internuclear separation range from about 0.1 to 1.4 nm. The calculated equilibrium geometry and dissociation energy of this potential energy curve are of Re=0.2487 nm and De=1.064eV, respectively. The unusual negative values of the anharmonicity constant and the vibration-rotational coupling constant are of ωeXe=-4.7158cm^-1 and αe=0.08649cm^-1, respectively. The vertical excitation energy from the ground to the A^1∑＋ state is calculated and the value is of 3.613eV at 0.15875nm （the equilibrium position of the ground state）. The highly anomalous shape of this potential energy curve, which is exceptionally flat over a wide radial range around the equilibrium position, is discussed in detail. The harmonic frequency value of 502.47cm^-1 about this state is approximately estimated. Careful comparison of the theoretical determinations with those obtained by previous theories about the A^1∑＋ state dissociation energy clearly shows that the present calculations are much closer to the experiments than previous theories, thus represents an improvement.