Ab initio investigation of the structural and unusual electronic properties of α-CuSe （klockmannite）

In this article, a computational analysis has been performed on the structural properties and predominantly on the electronic properties of the α-CuSe （klockmannite） using density functional theory. The studies in this work show that the best structural results, in comparison to the experimental values, belong to the PBEsol-GGA and WC-GGA functionals. However, the best results for the bulk modulus and density of states （DOSs） are related to the local density approximation （LDA） functional. Through utilized approaches, the LDA is chosen to investigate the electronic structure. The results of the electronic properties and geometric optimization of α-CuSe respectively show that this compound is conductive and non-magnetic. The curvatures of the energy bands crossing the Fermi level explicitly reveal that major charge carriers in CuSe are holes, whose density is estimated to be 0.86×1022 hole/cm3. In particular, the Fermi surfaces in the first Brillouin zone demonstrate interplane conductivity between （001） planes. Moreover, the charge carriers among them are electrons and holes simultaneously. The conductivity in CuSe is mainly due to the hybridization between the d orbitals of Cu atoms and the p orbitals of Se atoms. The former orbitals have the dual nature of localization and itinerancy.     

Partial Wave Cross Sections for Collisions between Helium Atoms and Hydrogen Halide Molecules

The close-coupling method is utilized to calculate partial cross sections for elastic and inelastic scattering of He atoms with HX （X=F, Cl, Br） molecule based on the CCSD （T） potential energy surfaces obtained in the previous research. The calculation is performed at the incident energy of 200 me V. The rationality of our results has been confirmed by comparison with the available theoretical results. The tendency of the elastic and inelastic rotational excitation partial wave cross sections varying with the reduced mass of the three systems is obtained.     

Analytical Potential Energy Function for the Ground State X1∑+ of Lanthanum Monofluoride

The equilibrium geometry, harmonic frequency and bond dissociation energy of lanthanum monofluoride have been calculated using Density-Functional Theory (DFT), post-HF methods MP2 and CCSD(T) with the energyconsistent relativistic effective core potentials. The possible electronic state and reasonable dissociation limit of the ground state of LaF are determined based on atomic and molecular reaction statics. Potential energy curve scans for the ground state X 1∑+ have been performed at B3LYP and CCSD(T) levels, due to their better results of harmonic frequency and bond dissociation energy. We find that the potential energy calculated with CCSD(T) is about 0.6 eV larger than the bond dissociation energy, when the internuclear distance is as large as 0.8 nm. The problem that single-reference ab initio methods do not meet dissociation limit during calculations of lanthanide heavy-metal elements is analyzed. We propose the calculation scheme to derive the analytical Murrell-Sorbie potential energy function. Vibrotational spectroscopic constants Be, ωe, ωeχe, αe, βe, De and He obtained by the standard Dunham treatment coincide well with the results of rotational analyses on spectroscopic experiments.     

The geometry structures and electronic properties of Li_mB_n(m+n=12) clusters

The geometric structures, electronic properties, total and binding energies, harmonic frequencies, the highest occupied molecular orbital to the lowest unoccupied molecular orbital energy gaps, and the vertical ionization potential energies of small LimBn （m＋ n = 12） clusters were investigated by the density functional theory B3LYP with a 6-31 I＋G （2d, 2p） basis set. All the calculations were performed using the Gaussian09 program. For the study of the LimBn clusters, the global minimum of the B 12 cluster was chosen as the starting point and the boron atoms were gradually replaced by Li atoms. The results showed that as the number of Li atoms increased, the stability of the LimBn cluster decreased and the physical and chemical properties became more active. In addition, on average there was a large charge transfer from the Li atoms to the B atoms.     

Methyl orbital signatures in 2-amino-l-propanol

Electron density distributions of 2-aminoethanol （2AE） and 2-amino-l-propanol （2AP） are calculated in both the coordinate and the momentum spaces using the B3LYP/TZVP method. Using the dual space analysis, molecular orbital signatures of the methyl substituent in 2AP are identified with respect to 2AE. Relaxations of the geometry and the valence orbital in 2AP are found to be due to the insertion of the methyl group. Five orbitals, not four orbitals, are identified as the methyl signatures. They are orbital 5a in the core shell, orbitals 9a and 10a in the inner valence shell, and orbitals 15a and 16a in the outer valence. In the inner valence shell, the attachment of methyl to 2AE causes a splitting of its orbital 8a into orbitals 9a and 10a of 2AP, whereas in the outer valence shell, the methyl group results in the insertion of an additional orbital pair of 15a and 16a. The frontier molecular orbitals 21a, 20a, and 19a are found to have no significant role in the methylation of 2AE.     

Methyl orbital signatures in 2-amino-1-propanol

Electron density distributions of 2-aminoethanol (2AE) and 2-amino-1-propanol (2AP) are calculated in both the coordinate and the momentum spaces using the B3LYP/TZVP method. Using the dual space analysis, molecular orbital signatures of the methyl substituent in 2AP are identified with respect to 2AE. Relaxations of the geometry and the valence orbital in 2AP are found to be due to the insertion of the methyl group. Five orbitals, not four orbitals, are identified as the methyl signatures. They are orbital 5a in the core shell, orbitals 9a and 10a in the inner valence shell, and orbitals 15a and 16a in the outer valence. In the inner valence shell, the attachment of methyl to 2AE causes a splitting of its orbital 8a into orbitals 9a and 10a of 2AP, whereas in the outer valence shell, the methyl group results in the insertion of an additional orbital pair of 15a and 16a. The frontier molecular orbitals 21a, 20a, and 19a are found to have no significant role in the methylation of 2AE.     

Hydration effect on the electronic transport properties of oligomeric phenylene ethynylene molecular junctions

A first-principles computational method based on the hybrid density functional theory is developed to simulate the electronic transport properties of oligomeric phenylene ethynylene molecular junctions with H₂O molecules accumulated in the vicinity as recently reported by Na {\it et al.} [\wx{Nanotechnology}{18} 424001 (2007)]. The numerical results show that the hydrogen bonds between the oxygen atoms of the oligomeric phenylene ethynylene molecule and H₂O molecules result in the localisation of the molecular orbitals and lead to the lower transition peaks. The H₂O molecular chains accumulated in the vicinity of the molecular junction can not only change the electronic structure of the molecular junctions, but also open additional electronic transport pathways. The obvious influence of H₂O molecules on the electronic structure of the molecular junction and its electronic transport properties is thus demonstrated.     

Electron Momentum Distributions of the Highest Occupied Molecular Orbital of CF3Br： Delocalization of Halogen Lone-Pair Orbitals

The highest occupied molecular orbital (HOMO) of trifluorobromomethane (CF3Br) is studied by binary (e,2e)electron momentum spectroscopy. The experimental momentum profile of the HOMO is compared with the Hartree-Fock (HF) and density functional theory (DFT-B3LYP) calculations. The calculated results largely depend on the size of basis sets rather than theoretical methods. Both the HF and DFT calculations using the 6-311 G^** basis set give a good explanation to the experiment. Delocalization of halogen lone-pair orbitals in the series molecules CF3X (X=F, Cl, Br) has been investigated.     

Probing structural and electronic properties of divalent metal Mgn+1 and SrMgn (n=2–12) clusters and their anions

Divalent metal clusters have received great attention due to the interesting size-induced nonmetal-to-metal transition and fascinating properties dependent on cluster size,shape,and doping.In this work,the combination of the CALYPSO code and density functional theory(DFT)optimization is employed to explore the structural properties of neutral and anionic Mg_n+1 and SrMgn(n=2-12)clusters.The results exhibit that as the atomic number of Mg increases,Sr atoms are more likely to replace Mg atoms located in the skeleton convex cap.By analyzing the binding energy,second-order energy difference and the charge transfer,it can be found the SrMg9 cluster with tower framework presents outstanding stability in a studied size range.Further,bonding characteristic analysis reveals that the stability of SrMg9 can be improved due to the strong s-p interaction among the atomic orbitals of Sr and Mg atoms.     

Magnetization dynamics of mixed Co–Au chains on Cu(110) substrate:Combined ab initio and kinetic Monte Carlo study

We present an investigation of the one-dimensional ferromagnetism in Au–Co nanowires deposited on the Cu(110)surface. By using the density functional theory, the influence of the nonmagnetic copper substrate Cu(110) on the magnetic properties of the bimetallic Au–Co nanowires is studied. The results show the emergence of magnetic anisotropy in the supported Au–Co nanowires. The magnetic anisotropy energy has the same order of magnitude as the exchange interaction energy between Co atoms in the wire. Our electronic structure calculation reveals the emergence of new hybridized bands between Au and Co atoms and surface Cu atoms. The Curie temperature of the Au–Co wires is calculated by means of kinetic Monte Carlo simulation. The strong size effect of the Curie temperature is demonstrated.     

Adaptive grouped paired comparison: theory and selection of seeds

In the subjective sound quality evaluation based on paired comparison method, time cost grows with square of the number of samples, which makes the evaluation of large quantity of samples difficult. An adaptive grouped paired comparison （AGPC） design was proposed to greatly decrease time requirement and difficulties of subjective comparison test, and to avoid the reduction of the precision when improper reference samples were chosen in common grouped paired comparison by adapting the seeds. The mathematic theory and implementation method were given, including the experiential formulae. It was shown from computer simulations and the subjective sound quality evaluation that the results of the proposed AGPC method were in good agreement with the true merits of samples and those obtained from conventional paired comparison.     

First-principles simulation of Raman spectra and structural properties of quartz up to 5 GPa

The crystal structure and Raman spectra of quartz are calculated by using first-principles method in a pressure range from 0 to 5 GPa. The results show that the lattice constants(a, c, and V) decrease with increasing pressure and the a-axis is more compressible than the c axis. The Si–O bond distance decreases with increasing pressure, which is in contrast to experimental results reported by Hazen et al. [Hazen R M, Finger L W, Hemley R J and Mao H K 1989 Solid State Communications 725 507–511], and Glinnemann et al. [Glinnemann J, King H E Jr, Schulz H, Hahn T, La Placa S J and Dacol F 1992 Z. Kristallogr. 198 177–212]. The most striking changes are of inter-tetrahedral O–O distances and Si–O–Si angles. The volume of the Si O4-4tetrahedron decreased by 0.9%(from 0 to 5 GPa), which suggests that it is relatively rigid.Vibrational models of the quartz modes are identified by visualizing the associated atomic motions. Raman vibrations are mainly controlled by the deformation of the Si O4-4tetrahedron and the changes in the Si–O–Si bonds. Vibrational directions and intensities of atoms in all Raman modes just show little deviations when pressure increases from 0 to 5 GPa.The pressure derivatives(dνi/d P) of the 12 Raman frequencies are obtained at 0 GPa–5 GPa. The calculated results show that first-principles methods can well describe the high-pressure structural properties and Raman spectra of quartz. The combination of first-principles simulations of the Raman frequencies of minerals and Raman spectroscopy experiments is a useful tool for exploring the stress conditions within the Earth.     

Unified treatment of one-range addition theorems for integer and non-integer n-STO, -GTO and -generalized exponential type orbitals with hyperbolic cosine in position, momentum and four-dimensional spaces

Simpler formulas are derived for one-range addition theorems for the integer and noninteger n generalized exponential type orbitals, momentum space orbitals, and hyperspherical harmonics with hyperbolic cosine (GETO HC, GMSO HC, and GHSH HC) in position, momentum and four-dimensional spaces, respectively. The final results are expressed in terms of one-range addition theorems of complete orthonormal sets of ψα -exponential type orbitals, α - momentum space orbitals and z α -hyperspherical harmonics. We notice that the one-range addition theorems for integer and noninteger n-Slater type orbitals and Gaussian type orbitals in position, momentum and four dimensional spaces are special cases of GETO HC, GMSO HC, and GHSH HC. The theorems presented can be useful in the accurate study of the electronic structure of atomic and molecular systems.     

Accurate all-electron calculation on the vibrational and rotational spectra of ground states for O_2 and its ions

The potential energy curves, spectroscopic constants, and low-lying vibration–rotation levels of ground-state O_2 and its cation O_2~+ and anion O_2~- were calculated with the explicitly correlated multireference configuration interaction method.The zeroth-order reference wavefunction was treated with the complete active space multiconfigurational self-consistent field method, in which the active space was carefully selected, and an additional molecular orbital πu was added into the full valence active space.The electron correlation of the 1s core in the oxygen atom was considered in the computations.The Davidson correction on molecular energy was considered to account for higher electron excitation.The relativistic effects, including the scalar relativistic effect and spin–orbit coupling, were considered in the computation of potential energy curves.These physical effects on the spectroscopic constants were examined.The low-lying levels of vibration–rotation spectra of O2 and its ions were determined based on the computed potential energy curves.Comparisons with available experiments were made and excellent agreement was obtained for the vibrational and rotational parameters.The spectroscopic constants and vibration–rotation spectrum of O_2~-, which is sparse in experiments, were provided.Our study will shed some light on further theoretical and experimental studies on these simple but important molecular systems.     

Vibrational investigation of 1-cyclopentylpiperazine: A combined experimental and theoretical study

FT-IR and Raman spectra of 1-cyclopentylpiperazine(1cppp)have been experimentally examined in the region of 4000–200cm-1.The optimized geometric parameters,conformational equilibria,normal mode frequencies and corresponding vibrational assignments of 1cppp(C9H18N2)are theoretically examined by means of B3LYP hybrid density functional theory(DFT)method together with 6-31++G(d,p)basis set.On the basis of potential energy distribution(PED)reliable vibrational assignments have been made and the thermodynamics functions,highest occupied and lowest unoccupied molecular orbitals(HOMO and LUMO)of 1cppp have been predicted.Calculations are employed for four different conformations in C1 and Cs point groups of 1cppp in gas phase.Comparison between the experimental and theoretical results indicates that B3LYP method is able to provide satisfactory results for predicting vibrational frequencies and the structural parameters,vibrational frequencies and assignments.Furthermore,C1(equatorial-axial)point group has been found as the most stable conformer of 1cppp.     

Energetics and electronic structure of refractory elements in the dislocation of NiAl

Using DMol and the discrete variational method within the framework of the density functional theory,we study the alloying effects of Nb,Ti,and V in the [100](010) edge dislocation core of NiAl.We find that when Nb(Ti,V) is substituted for Al in the center-Al,the binding energy of the system reduces 3.00 eV(2.98 eV,2.66 eV).When Nb(Ti,V) is substituted for Ni in the center-Ni,the binding energy of the system reduces only 0.47 eV(0.16 eV,0.09 eV).This shows that Nb(Ti,V) exhibits a strong Al site preference,which agrees with the experimental and other theoretical results.The analyses of the charge distribution,the interatomic energy and the partial density of states show that some charge accumulations appear between the impurity atom and Ni atoms,and the strong bonding states are formed between impurity atom and neighbouring host atoms due mainly to the hybridization of 4d5s(3d4s) orbitals of impurity atoms and 3d4s4p orbitals of host Ni atoms.The impurity induces a strong pinning effect on the [100](010) edge dislocation motion in NiAl,which is related to the mechanical properties of the NiAl alloy.     

Assignment of terahertz vibrational modes of L-glutamine using density functional theory within generalized-gradient approximation

The crystal structure of L-glutamine is stabilized by a three-dimensional network of intermolecular hydrogen bonds.We utilize plane-wave density functional theory lattice-dynamics calculations within the generalized-gradient approximation(GGA), Perdew–Burke–Ernzerhof(PBE), PBE for solids(PBEsol), PBE with Wu–Cohen exchange(WC), and dispersion-corrected PBE, to investigate the effect of these intermolecular contacts on the absorption spectra of glutamine in the terahertz frequency range. Among these calculations, the solid-state simulated results obtained using the WC method exhibit a good agreement with the measured absorption spectra, and the absorption features are assigned with the help of WC. This indicates that the vibrational modes of glutamine were related to the combination of intramolecular and intermolecular motions, the intramolecular modes were dominated by rocking or torsion involving functional groups; the intermolecular modes mainly result from the translational motions of individual molecules, and the rocking of the hydrogenbonded functional groups.     

Electron momentum spectroscopy of NF3

The electronic structure of nitrogen trifluoride was investigated by combining the high-resolution electron momentum spectroscopy with the high-level calculations. The experimental binding energy spectra and the momentum distributions of each orbital were compared with the results of Hartree-Fock, density functional theory （DFT）, and symmetry-adapted- cluster configuration-interaction （SAC-CI） methods. SAC-CI and DFT-B3LYP with the aug-cc-pVTZ basis set can well reproduce the binding energy spectra and the observed momentum distributions of the valence orbitals except 1 a2 and 4e orbitals. It was found that the calculated momentum distributions using DFT-B3LYP are even better than those using the high-level SAC-CI method.     

Electron correlation in fast ion-impact single ionization of helium atoms

A four-body distorted-wave approximation is applied for theoretical analysis of the fully differential cross sections(FDCS)for proton-impact single ionization of helium atoms in their ground states.The nine-dimensional integrals for the partial amplitudes are analytically reduced to closed-form expressions or some one-dimensional integrals which can be easily calculated numerically.Calculations are performed in the scattering and perpendicular planes.The influence of the target static electron correlations on the process is investigated using a number of different bound-state wave functions for the ground state of the helium targets.An illustrative computation is performed for 75-ke V proton–helium collisions and the obtained results are compared with experimental data and other theoretical predictions.Although for small momentum transfers,the comparison shows a reasonable agreement with experiments in the scattering and perpendicular planes,some significant discrepancies are still present at large momentum transfers in these planes.However,our results are compatible and for some cases,better than those of the other sophisticated calculations.     

Two-dimensional arsenic monolayer sheet predicted from first-principles

Using first-principles calculations, we investigate the two-dimensional arsenic nanosheet isolated from bulk gray arsenic. Its dynamical stability is confirmed by phonon calculations and molecular dynamics analyzing. The arsenic sheet is an indirect band gap semiconductor with a band gap of 2.21 e V in the hybrid HSE06 functional calculations. The valence band maximum(VBM) and the conduction band minimum(CBM) are mainly occupied by the 4p orbitals of arsenic atoms,which is consistent with the partial charge densities of VBM and CBM. The charge density of the VBM G point has the character of a π bond, which originates from p orbitals. Furthermore, tensile and compressive strains are applied in the armchair and zigzag directions, related to the tensile deformations of zigzag and armchair nanotubes, respectively. We find that the ultimate strain in zigzag deformation is 0.13, smaller than 0.18 of armchair deformation. The limit compressive stresses of single-layer arsenic along armchair and zigzag directions are-4.83 GPa and-4.76 GPa with corresponding strains of-0.15 and-0.14, respectively.