Search

Element-specific unoccupied electronic states of Cu(In, Ga)S2 were studied as a function of the In/Ga ratio by combining X-ray absorption spectroscopy with density functional theory calculations. The S absorption edge shifts with changing In/Ga ratio as expected from the variation of the band gap. In contrast, the cation edge positions are largely independent of composition despite the changing band gap. This unexpected behavior is well reproduced by our calculations and originates from the dependence of the electronic states on the local atomic environment. The changing band gap arises from a changing spatial average of these localized states with changing alloy composition. [ABSTRACT FROM AUTHOR]

We describe the calculation of Raman spectra for periodic systems via ab initio molecular dynamics (AIMD) utilizing the Gaussian and plane wave method in the program package CP2K. The electric -dipole-electric-dipole polarizability tensor has been implemented for an arbitrary shape of the simulation cell. In addition, a computationally efficient approach for its decomposition into local contributions is presented. As an example for the application of computational Raman spectroscopy to liquids, the Raman spectra of S-methyloxirane in the liquid phase have been calculated together with Raman spectra obtained from static calculations employing the double-harmonic approximation. The comparison to experimental data illustrates that a very good agreement between experiment and simulated spectra can be obtained employing AIMD, which takes into account anharmonicities and dynamical effects at ambient conditions. [ABSTRACT FROM AUTHOR]

Ferroelectric properties of tetragonal (BaZrO3)1/(PbZrO3)n and (SrZrO3)1/(PbZrO3)n superlattices are studied as the ratio of PbZrO3 was changed (n=1-3) in the superlattices using density functional theory (DFT) method as well as an electrostatic model. We show that PbZrO3-an antiferroelectric in bulk displays ferroelectric behavior in the superlattices combined with the paraelectric cubic (BaZrO3) or (SrZrO3). The electrostatic model and the DFT calculations show that the polarization and the tetragonality (c/a ratio) are reduced in the BaZrO3/PbZrO3 and SrZrO3/PbZrO3 superlattices with respect to the bulk tetragonal PbZrO3. Moreover, the polarization and the tetragonality of superlattices increase with increasing the fraction of PbZrO3 in superlattices. The estimated polarization of the (BaZrO3)1/(PbZrO3)1 superlattice is (36.21 μ/cm2), while it is (56.82 μ/cm2) for the (BaZrO3)1/(PbZrO3)3 superlattice. The estimated polarization of the (SrZrO3)1/(PbZrO3)1 superlattice is (55.28 μ/cm2), while it is (64.86 μ/cm2) for the (SrZrO3)1/(PbZrO3)3 superlattice. These ferroelectric superlattices have good lattice matching with many Heusler alloys and they could be very useful as ferroelectric substrates to these systems. [ABSTRACT FROM AUTHOR]

The electromechanical coupling at the silicon (100) and (111) surfaces was studied via density functional theory by calculating the response of the ionization potential and the electron affinity to different types of strain. We find a branched strain response of those two quantities with different coupling coefficients for negative and positive strain values. This can be attributed to the reduced crystal symmetry due to anisotropic strain, which partially lifts the degeneracy of the valence and conduction bands. Only the Si(111) electron affinity exhibits a monotonously linear strain response, as the conduction band valleys remain degenerate under strain. The strain response of the surface dipole is linear and seems to be dominated by volume changes. Our results may help to understand the mechanisms behind electromechanical coupling at an atomic level in greater detail and for different electronic and atomic structures. [ABSTRACT FROM AUTHOR]

We have studied the structural and electronic properties of YPd2Sn in the Heusler structure using a generalized gradient approximation of the density functional theory and the ab initio pseudopotential method. The electronic results indicate that the density of states at the Fermi level is primarily derived from Pd d states, which hybridize with Y d and Sn p states. Using our structural and electronic results, phonons and electron-phonon interactions have been studied by employing a linear response approach based on the density functional theory. Phonon anomalies have been observed for transverse acoustic branches along the [110] direction. This anomalous dispersion is merely a consequence of the strong coupling. By integrating the Eliashberg spectral function, the average electron-phonon coupling parameter is found to be λ=0.99. Using this value, the superconducting critical temperature is calculated to be 4.12 K, in good accordance with the recent experimental value of 4.7 K. [ABSTRACT FROM AUTHOR]

We report a first-principles prediction of the Raman shifts of parahydrogen (pH2) clusters of sizes N = 4-19 and 33, based on path integral ground-state simulations with an ab initio potential energy surface. The Raman shifts are calculated, using perturbation theory, as the average of the differencepotential energy surface between the potential energy surfaces for vibrationally excited and groundstate parahydrogen monomers. The radial distribution of the clusters is used as a weight function in this average. Very good overall agreement with experiment [G. Tejeda, J. M. Fernández, S. Montero, D. Blume, and J. P. Toennies, Phys. Rev. Lett. 92, 223401 (2004)] is achieved for p(H2)2-8,13,33. A number of different pair potentials are employed for the calculation of the radial distribution functions. We find that the Raman shifts are sensitive to slight variations in the radial distribution functions. [ABSTRACT FROM AUTHOR]

Time independent quantum mechanical (TIQM) scattering calculations have been carried out for the O(¹D) + CH4(X¹A1) → CH3(X²A2″) + OH(X²∏) atmospheric reaction, using an ab initio ground potential energy surface where the CH3 group is described as a pseudo-atom. Total and state-to-state reaction probabilities for a total angular momentum J = 0 have been determined for collision energies up to 0.5 eV. The vibrational and rotational state OH product distributions show no specific behavior. The rate coefficient has been calculated by means of the J-shifting approach in the 10-500 K temperature range and slightly depends on T at ordinary temperatures (as expected for a barrierless reaction). Quantum effects do not influence the vibrational populations and rate coefficient in an important way, and a rather good agreement has been found between the TIQM results and the quasiclassical trajectory and experimental ones. This reinforces somewhat the reliability of the pseudo-triatomic approach under the reaction conditions explored. [ABSTRACT FROM AUTHOR]

The (O3)2 dimer potential energy surface is thoroughly explored at the ab initio CCSD(T) computational level. Five minima are characterized with binding energies between 0.35 and 2.24 kcal/mol. The most stable may be characterized as slipped parallel, with the two O3 monomers situated in parallel planes. Partitioning of the interaction energy points to dispersion and exchange as the prime contributors to the stability, with varying contributions from electrostatic energy, which is repulsive in one case. Atoms in Molecules analysis of the wavefunction presents specific O·O bonding interactions, whose number is related to the overall stability of each dimer. All internal vibrational frequencies are shifted to the red by dimerization, particularly the antisymmetric stretching mode whose shift is as high as 111 cm-1. In addition to the five minima, 11 higher-order stationary points are identified. [ABSTRACT FROM AUTHOR]

Vibrationally resolved spectra provide a stringent test of the accuracy of theoretical calculations. We combine the thawed Gaussian approximation (TGA) with an on-the-fly ab initio (OTF-AI) scheme to calculate the vibrationally resolved emission spectra of oligothiophenes with up to five rings. The efficiency of the OTF-AI-TGA permits treating all vibrational degrees of freedom on an equal footing even in pentathiophene with 105 vibrational degrees of freedom, thus obviating the need for the global harmonic approximation, popular for large systems. Besides reproducing almost perfectly the experimental emission spectra, in order to provide a deeper insight into the associated physical and chemical processes, we also develop a novel systematic approach to assess the importance and coupling between individual vibrational degrees of freedom during the dynamics. This allows us to explain how the vibrational line shapes of the oligothiophenes change with increasing number of rings. Furthermore, we observe the dynamical interplay between the quinoid and aromatic characters of individual rings in the oligothiophene chain during the dynamics and confirm that the quinoid character prevails in the center of the chain. [ABSTRACT FROM AUTHOR]

Photodissociation pathways of nitromethane following π → π ∗ electronic excitation are reported. The potential energy surfaces for four lowest singlet states are explored, and structures of many intermediates, dissociation limits, transition states, and minimum energy conical intersections were determined using the automated searching algorism called the global reaction route mapping strategy. Geometries are finally optimized at CASSCF(14e,11o) level and energies are computed at CAS(14o,11e)PT2 level. The calculated preferable pathways and important products qualitatively explain experimental observations. The major photodissociation product CH3 and NO2 (2B2) is formed by direct dissociation from the S1 state. Important pathways involving S1 and S0 states for production of various dissociation products CH3NO + O (¹D), CH3O(X²E) + NO (X²Π), CH2NO + OH, and CH2O + HNO, as well as various isomerization pathways have been identified. Three roaming processes also have been identified: the O atom roaming in O dissociation from CH3NO2, the OH radical roaming in OH dissociation from CH2N(O)(OH), and the NO roaming in NO dissociation from CH3ONO. [ABSTRACT FROM AUTHOR]

Motivated by a recent experiment [H. Pan and K. Liu, J. Chem. Phys. 140, 191101 (2014)], we report a quasiclassical trajectory study of the O(3P) + CH4(vk = 0, 1) → OH + CH3 [k = 1 and 3] reactions on an ab initio potential energy surface. The computed angular distributions and cross sections correlated to the OH(v = 0, 1) + CH3(v = 0) coincident product states can be directly compared to experiment for O + CH4(v3 = 0, 1). Both theory and experiment show that the ground-state reaction is backward scattered, whereas the angular distributions shift toward sideways and forward directions upon antisymmetric stretching (v3) excitation of the reactant. Theory predicts similar behavior for the O + CH4(v1 = 1) reaction. The simulations show that stretching excitation enhances the reaction up to about 15 kcal/mol collision energy, whereas the O + CH4(vk = 1) reactions produce smaller cross sections for OH(v = 1) + CH3(v = 0) than those of O + CH4(v = 0) →OH(v = 0) + CH3 (v = 0). The former finding agrees with experiment and the latter awaits for confirmation. The computed cold OH rotational distributions of O + CH4(v = 0) are in good agreement with experiment. [ABSTRACT FROM AUTHOR]

A general methodology has been devised and implemented into the solid-state ab initio quantummechanical CRYSTAL program for studying the evolution under geophysical pressure of the elastic anisotropy of crystalline materials. This scheme, which fully exploits both translational and point symmetry of the crystal, is developed within the formal frame of one-electron Hamiltonians and atom-centered basis functions. Six silicate garnet end-members, among the most important rockforming minerals of the Earth's mantle, are considered, whose elastic anisotropy is fully characterized under high hydrostatic compressions, up to 60 GPa. The pressure dependence of azimuthal anisotropy and shear-wave birefringence of seismic wave velocities for these minerals are accurately simulated and compared with available single-crystal measurements. [ABSTRACT FROM AUTHOR]

A new global ab initio potential energy surface (called ZMB-a) for the 11A' state of the C(1D)+H2 reactive system has been constructed. This is based upon ab initio calculations using the internally contracted multireference configuration interaction approach with the aug-cc-pVQZ basis set, performed at about 6300 symmetry unique geometries. Accurate analytical fits are generated using many-body expansions with the permutationally invariant polynomials, except that the fit of the deep well region is taken from our previous fit. The ZMB-a surface is unique in the accurate description of the regions around conical intersections (CIs) and of van der Waals (vdW) interactions. The CIs between the 11A' and 21A' states cause two kinds of barriers on the ZMB-a surface: one is in the linear H-CH dissociation direction with a barrier height of 9.07 kcal/mol, which is much higher than those on the surfaces reported before; the other is in the C(1D) collinearly attacking H2 direction with a barrier height of 12.39 kcal/mol. The ZMB-a surface basically reproduces our ab initio calculations in the vdW interaction regions, and supports a linear C-HH vdW complex in the entrance channel, and two vdW complexes in the exit channel, at linear CH-H and HC-H geometries, respectively. [ABSTRACT FROM AUTHOR]

Excitation of HeNe by synchrotron light just below the frequency of the 1s→3p transition of isolated He has been recently shown to be followed by resonant interatomic Coulombic decay (ICD). The vibrationally resolved widths of the ICD states were extracted with high precision from the photoion spectra. In this paper, we report the results of ab initio calculations of these widths. We show that interaction between electronic states at about the equilibrium distance of HeNe makes dark states of He accessible for the photoexcitation and subsequent electronic decay.Moreover, the values of the calculated widths are shown to be strongly sensitive to the presence of the non-adiabatic coupling between the electronic states participating in the decay. Therefore, only by considering the complete manifold of interacting decaying electronic states a good agreement between the measured and computed ICD widths can be achieved. [ABSTRACT FROM AUTHOR]

An accurate full-dimensional global potential energy surface (PES) is developed for the title reaction. While the long-range interactions in the reactant asymptote are represented by an analytical expression, the interaction region of the PES is fit to more than 81 000 of ab initio points at the UCCSD(T)-F12b/AVTZ level using the permutation invariant polynomial neural network approach. Fully symmetric with respect to permutation of all four hydrogen atoms, the PES provides a faithful representation of the ab initio points, with a root mean square error of 1.8 meV or 15 cm-1. The reaction path for this exoergic reaction features an attractive and barrierless entrance channel, a submerged saddle point, a shallow H4O+ well, and a barrierless exit channel. The rate coefficients for the title reaction and kinetic isotope effect have been determined on this PES using quasi-classical trajectories, and they are in good agreement with available experimental data. It is further shown that the H2O+ rotational enhancement of reactivity observed experimentally can be traced to the submerged saddle point. Using our recently proposed Sudden Vector Projection model, we demonstrate that a rotational degree of freedom of the H2O+ reactant is strongly coupled with the reaction coordinate at this saddle point, thus unraveling the origin of the pronounced mode specificity in this reaction. [ABSTRACT FROM AUTHOR]

Starting from experimental findings and interface growth problems in Si/Ge superlattices, we have investigated through ab initio methods the concurrent and competitive behavior of strain and defects in the second-harmonic generation process. Interpreting the second-harmonic intensities as a function of the different nature and percentage of defects together with the strain induced at the interface between Si and Ge, we found a way to tune and enhance the second-harmonic generation response of these systems. [ABSTRACT FROM AUTHOR]

We present an analysis of the spin-rotation and absolute shielding constants of XF6 molecules (X = S, Se, Te, Mo, W) based on ab initio coupled cluster and four-component relativistic densityfunctional theory (DFT) calculations. The results show that the relativistic contributions to the spin- rotation and shielding constants are large both for the heavy elements as well as for the fluorine nuclei. In most cases, incorporating the computed relativistic corrections significantly improves the agreement between our results and the well-established experimental values for the isotropic spin-rotation constants and their anisotropic components. This suggests that also for the other molecules, for which accurate and reliable experimental data are not available, reliable values of spin-rotation and absolute shielding constants were determined combining ab initio and relativistic DFT calculations. For the heavy nuclei, the breakdown of the relationship between the spin- rotation constant and the paramagnetic contribution to the shielding constant, due to relativistic effects, causes a significant error in the total absolute shielding constants. [ABSTRACT FROM AUTHOR]

The treatment of dispersion interactions is ubiquitous but computationally demanding for seamless ab initio approaches. A highly popular and simple remedy consists in correcting for the missing interactions a posteriori by adding an attractive energy term summed over all atom pairs to standard density functional approximations. These corrections were originally based on atom pairwise parameters and, hence, had a strong touch of empiricism. To overcome such limitations, we recently proposed a robust system-dependent dispersion correction, dDsC, that is computed from the electron density and that provides a balanced description of both weak inter- and intramolecular interactions. From the theoretical point of view and for the sake of increasing reliability, we here verify if the self-consistent implementation of dDsC impacts ground-state properties such as interaction energies, electron density, dipole moments, geometries, and harmonic frequencies. In addition, we investigate the suitability of the a posteriori scheme for molecular dynamics simulations, for which the analysis of the energy conservation constitutes a challenging tests. Our study demonstrates that the post-SCF approach in an excellent approximation. [ABSTRACT FROM AUTHOR]

We present and test a new approximation for the exchange-correlation (xc) energy of Kohn-Sham density functional theory. It combines exact exchange with a compatible non-local correlation functional. The functional is by construction free of one-electron self-interaction, respects constraints derived from uniform coordinate scaling, and has the correct asymptotic behavior of the xc energy density. It contains one parameter that is not determined ab initio. We investigate whether it is possible to construct a functional that yields accurate binding energies and affords other advantages, specifically Kohn-Sham eigenvalues that reliably reflect ionization potentials. Tests for a set of atoms and small molecules show that within our local-hybrid form accurate binding energies can be achieved by proper optimization of the free parameter in our functional, along with an improvement in dissociation energy curves and in Kohn-Sham eigenvalues. However, the correspondence of the latter to experimental ionization potentials is not yet satisfactory, and if we choose to optimize their prediction, a rather different value of the functional's parameter is obtained. We put this finding in a larger context by discussing similar observations for other functionals and possible directions for further functional development that our findings suggest. [ABSTRACT FROM AUTHOR]

The effect of the spin-orbit interaction to photofragmentation is investigated in the mercury(II) bromide (HgBr2) molecule. Changes in the fragmentation between the two spin-orbit components of Hg 5d photoionization, as well as within the molecular-field-splitted levels of these components are observed. Dissociation subsequent to photoionization is studied with synchrotron radiation and photoelectron-photoion coincidence spectroscopy. The experimental results are accompanied by relativistic ab initio analysis of the photoelectron spectrum. [ABSTRACT FROM AUTHOR]

A systematic first-principles non-equilibrium Green's function study is conducted on the contact resistance between a series of metals (Au, Ag, Pt, Cu, Ni, and Pd) and graphene in the side contact geometry. Different factors such as the termination of the graphene edge, contact area, and point defect in contacted graphene are investigated. Notable differences are observed in structural configurations and electronic transport characteristics of these metal-graphene contacts, depending on the metal species and aforementioned influencing factors. It is found that the enhanced chemical reactivity of the graphene due to dangling bonds from either the unsaturated graphene edge or point defects strengthens the metal-graphene bonding, leading to a considerable contact resistance reduction for weakly interacting metals Au and Ag. For stronger interacting metals Pt and Cu, a slightly reduced contact resistance is found due to such influencing factors. However, the wetting metals Ni and Pd most strongly hybridize with graphene, exhibiting negligible dependence on the above influencing factors. This study provides guidance for the optimization of metal-graphene contacts at an atomic scale. [ABSTRACT FROM AUTHOR]

Based on the non-local van der Waals density functional (vdW-DF)+U scheme, we carry out the ab initio molecular dynamics (AIMD) study of the interaction dynamics for H2 impingement against the stoichiometric PuO2(111), the reduced PuO2(111), and the stoichiometric α-Pu2O3(111) surfaces. The hydrogen molecular physisorption states, which cannot be captured by pure DFT+U method, are obtained by employing the vdW-DF+U scheme. We show that except for the weak physisorption, PuO2(111) surfaces are so difficult of access that almost all of the H2 molecules will bounce back to the vacuum when their initial kinetic energies are not sufficient. Although the dissociative adsorption of H2 on PuO2(111) surfaces is found to be very exothermic, the collision-induced dissociation barriers of H2 are calculated to be as high as 3.2 eV and 2.0 eV for stoichiometric and reduced PuO2 surfaces, respectively. Unlike PuO2, our AIMD study directly reveals that the hydrogen molecules can penetrate into α-Pu2O3(111) surface and diffuse easily due to the 25% native O vacancies located along the <111> diagonals of α-Pu2O3 matrix. By examining the temperature effect and the internal vibrational excitations of H2, we provide a detailed insight into the interaction dynamics of H2 in α-Pu2O3. The optimum pathways for hydrogen penetration and diffusion, the corresponding energy barriers (1.0 eV and 0.53 eV, respectively) and rate constants are systematically calculated. Overall, our study fairly reveals the different interaction mechanisms between H2 and Pu-oxide surfaces, which have strong implications to the interpretation of experimental observations. [ABSTRACT FROM AUTHOR]

Powder and single crystal Raman spectra of the two most common phases of calcium carbonate are calculated with ab initio techniques (using a “hybrid” functional and a Gaussian-type basis set) and measured both at 80 K and room temperature. Frequencies of the Raman modes are in very good agreement between calculations and experiments: the mean absolute deviation at 80 K is 4 and 8 cm-1 for calcite and aragonite, respectively. As regards intensities, the agreement is in general good, although the computed values overestimate the measured ones in many cases. The combined analysis permits to identify almost all the fundamental experimental Raman peaks of the two compounds, with the exception of either modes with zero computed intensity or modes overlapping with more intense peaks. Additional peaks have been identified in both calcite and aragonite, which have been assigned to 18O satellite modes or overtones. The agreement between the computed and measured spectra is quite satisfactory; in particular, simulation permits to clearly distinguish between calcite and aragonite in the case of powder spectra, and among different polarization directions of each compound in the case of single crystal spectra. [ABSTRACT FROM AUTHOR]

The electronic absorption spectrum of liquid water was investigated by coupling a one-body energy decomposition scheme to configurations generated by classical and Born-Oppenheimer Molecular Dynamics (BOMD). A Frenkel exciton Hamiltonian formalism was adopted and the excitation energies in the liquid phase were calculated with the equation of motion coupled cluster with single and double excitations method. Molecular dynamics configurations were generated by different approaches. Classical MD were carried out with the TIP4P-Ew and AMOEBA force fields. The BLYP and BLYP-D3 exchange-correlation functionals were used in BOMD. Theoretical and experimental results for the electronic absorption spectrum of liquid water are in good agreement. Emphasis is placed on the relationship between the structure of liquid water predicted by the different models and the electronic absorption spectrum. The theoretical gas to liquid phase blue-shift of the peak positions of the electronic absorption spectrum is in good agreement with experiment. The overall shift is determined by a competition between the O-H stretching of the water monomer in liquid water that leads to a red-shift and polarization effects that induce a blue-shift. The results illustrate the importance of coupling many-body energy decomposition schemes to molecular dynamics configurations to carry out ab initio calculations of the electronic properties in liquid phase. [ABSTRACT FROM AUTHOR]

The cost of calculating nuclear hessians, either analytically or by finite difference methods, during the course of quantum chemical analyses can be prohibitive for systems containing hundreds of atoms. In many applications, though, only a few eigenvalues and eigenvectors, and not the full hessian, are required. For instance, the lowest one or two eigenvalues of the full hessian are sufficient to characterize a stationary point as a minimum or a transition state (TS), respectively. We describe here a method that can eliminate the need for hessian calculations for both the characterization of stationary points as well as searches for saddle points. A finite differences implementation of the Davidson method that uses only first derivatives of the energy to calculate the lowest eigenvalues and eigenvectors of the hessian is discussed. This method can be implemented in conjunction with geometry optimization methods such as partitioned-rational function optimization (P-RFO) to characterize stationary points on the potential energy surface. With equal ease, it can be combined with interpolation methods that determine TS guess structures, such as the freezing string method, to generate approximate hessian matrices in lieu of full hessians as input to P-RFO for TS optimization. This approach is shown to achieve significant cost savings relative to exact hessian calculation when applied to both stationary point characterization as well as TS optimization. The basic reason is that the present approach scales one power of system size lower since the rate of convergence is approximately independent of the size of the system. Therefore, the finite-difference Davidson method is a viable alternative to full hessian calculation for stationary point characterization and TS search particularly when analytical hessians are not available or require substantial computational effort. [ABSTRACT FROM AUTHOR]

Photoinduced charge-transfer processes play a key role for novel photovoltaic phenomena and devices. Thus, the development of ab initio methods that allow for an accurate and computationally inexpensive treatment of charge-transfer excitations is a topic that nowadays attracts a lot of scientific attention. In this paper we extend an approach recently introduced for the description of single and double excitations [M. Tassi, I. Theophilou, and S. Thanos, Int. J. Quantum Chem. 113, 690 (2013); J. Chem. Phys. 138, 124107 (2013)] to allow for the description of intermolecular chargetransfer excitations. We describe an excitation where an electron is transferred from a donor system to an acceptor one, keeping the excited state orthogonal to the ground state and avoiding variational collapse. These conditions are achieved by decomposing the space spanned by the Hartree-Fock (HF) ground state orbitals into four subspaces: The subspace spanned by the occupied orbitals that are localized in the region of the donor molecule, the corresponding for the acceptor ones and two more subspaces containing the virtual orbitals that are localized in the neighborhood of the donor and the acceptor, respectively. Next, we create a Slater determinant with a hole in the subspace of occupied orbitals of the donor and a particle in the virtual subspace of the acceptor. Subsequently we optimize both the hole and the particle by minimizing the HF energy functional in the corresponding subspaces. Finally, we test our approach by calculating the lowest charge-transfer excitation energies for a set of tetracyanoethylene-hydrocarbon complexes that have been used earlier as a test set for such kind of excitations. [ABSTRACT FROM AUTHOR]

A special case of predissociation, known as indirect or accidental predissociation observed in the Rydberg E1Π bound state of CO is discussed. We resort to ab initio potentials in order to determine the plausible mechanism for this predissociation. Values of the predissociation width for the valence k3Π state of CO, as obtained from Fermi's golden rule, are also reported. The predissociation width obtained for the mixed E1Π (v = 1, J = 7) state is 0.033 cm-1 compared to the experimental value of 0.034 cm-1. The mixed E − E' state with J = 28, v = 0 is found to be in near resonance condition with the k3Π (v = 4, J = 28) state, thus providing the means to indirect predissociation. [ABSTRACT FROM AUTHOR]

We have studied the ... 2A1–... 2B2 laser-induced fluorescence (LIF) spectrum of the jet-cooled F2BO radical for the first time. The transition consists of a strong 000 band at 446.5 nm and eight weak sequence bands to shorter wavelengths. Single vibronic level emission spectra obtained by laser excitation of individual levels of the ... state exhibit two electronic transitions: a very weak, sparse ... band system in the 450–500 nm region and a stronger, more extensive set of ... 2A1– ... 2B1 bands in the 580–650 nm region. We have also performed a series of high level ab initio calculations to predict the electronic energies, molecular structures, vibrational frequencies, and rotational and spin-rotation constants in the ... 2B2, Ã2B1 and ... 2A1 electronic states as an aid to the analysis of the experimental data. The theoretical results have been used as input for simulations of the rotationally resolved ... 2A1–... 2B2 000 LIF band and Franck-Condon profiles of the LIF and single vibronic level emission spectra. The agreement between the simulations obtained with purely ab initio parameters and the experimental spectra validates the geometries calculated for the ground and excited states and the conclusion that the radical has C2v symmetry in the ..., Ã, and ... states. The spectra provide considerable new information about the vibrational energy levels of the ... and à states, but very little for the ... state, due to the very restrictive Franck-Condon factors in the LIF spectra. [ABSTRACT FROM AUTHOR]

A full-dimensional, global ab initio potential energy surface (PES) for the Ar-HOCO system is presented. The PES consists of a previous intramolecular ab initio PES for HOCO [J. Li, C. Xie, J. Ma, Y. Wang, R. Dawes, D. Xie, J. M. Bowman, and H. Guo, J. Phys. Chem. A 116, 5057 (2012)], plus a new permutationally invariant interaction potential based on fitting 12 432 UCCSD(T)-F12a/aVDZ counterpoise-corrected energies. The latter has a total rms fitting error of about 25 cm-1 for fitted interaction energies up to roughly 12 000 cm-1. Two additional fits are presented. One is a novel very compact permutational invariant representation, which contains terms only involving the Ar-atom distances. The rms fitting error for this fit is 193 cm-1. The other fit is the widely used pairwise one. The pairwise fit to the entire data set has an rms fitting error of 427 cm-1. All of these potentials are used in preliminary classical trajectory calculations of energy transfer with a focus on comparisons with the results using the benchmark potential. [ABSTRACT FROM AUTHOR]

Classical density-functional theory is the most direct approach to equilibrium structures and free energies of inhomogeneous liquids, but requires the construction of an approximate free-energy functional for each liquid of interest. We present a general recipe for constructing functionals for small-molecular liquids based only on bulk experimental properties and ab initio calculations of a single solvent molecule. This recipe combines the exact free energy of the non-interacting system with fundamental measure theory for the repulsive contribution and a weighted density functional for the short-ranged attractive interactions. We add to these ingredients a weighted polarization functional for the long-range correlations in both the rotational and molecular-polarizability contributions to the dielectric response. We also perform molecular dynamics calculations for the free energy of cavity formation and the high-field dielectric response, and show that our free-energy functional adequately describes these properties (which are key for accurate solvation calculations) for all three solvents in our study: water, chloroform, and carbon tetrachloride. [ABSTRACT FROM AUTHOR]

We study charge transfer in bridged di- and triruthenium complexes from a theoretical and computational point of view. Ab initio computations are interpreted from the perspective of a simple empirical Hamiltonian, a chemically specific Mott-Hubbard model of the complexes' π electron systems. This Hamiltonian is coupled to classical harmonic oscillators mimicking a polarizable dielectric environment. The model can be solved without further approximations in a valence bond picture using the method of exact diagonalization and permits the computation of charge transfer reaction rates in the framework of Marcus' theory. In comparison to the exact solution, the Hartree-Fock mean field theory overestimates both the activation barrier and the magnitude of charge-transfer excitations significantly. For triruthenium complexes, we are able to directly access the interruthenium antiferromagnetic coupling strengths. [ABSTRACT FROM AUTHOR]

The density functional theory (DFT)+U method is an efficient and effective way to calculate the ground-state properties of strongly correlated transition metal compounds, with the effective U parameters typically determined empirically. Two ab initio methods have been developed to compute the U parameter based on either constrained DFT (CDFT) or unrestricted Hartree-Fock (UHF) theory. Previous studies have demonstrated the success of both methods in typical open-shell materials such as FeO and NiO. In this Communication we report numerical instability issues that arise for the CDFT method when applied to closed-shell transition metals, by using ZnO and Cu2O as examples. By contrast, the UHF method behaves much more robustly for both closed- and open-shell materials, making it more suitable for treating closed-shell transition metals, as well as main group elements. [ABSTRACT FROM AUTHOR]

The rotational spectra of six isotopologues of a complex of ethyne and silver acetylide, C2H2· · ·AgCCH, are measured by both chirped-pulse and Fabry-Perot cavity versions of Fouriertransform microwave spectroscopy. The complex is generated through laser ablation of a silver target in the presence of a gas sample containing 1% C2H2, 1% SF6, and 98% Ar undergoing supersonic expansion. Rotational, A0, B0, C0, and centrifugal distortion ΔJ and ΔJK constants are determined for all isotopologues of C2H2· · ·AgCCH studied. The geometry is planar, C2v and T-shaped in which the C2H2 sub-unit comprises the bar of the “T” and binds to the metal atom through its π electrons. In the r0 geometry, the distance of the Ag atom from the centre of the triple bond in C2H2 is 2.2104(10) Å. The r(HC≡CH) parameter representing the bond distance separating the two carbon atoms and the angle, ∠(CCH), each defined within the C2H2 sub-unit, are determined to be 1.2200(24) Å and 186.0(5)◦, respectively. This distortion of the linear geometry of C2H2 involves the hydrogen atoms moving away from the silver atom within the complex. The results thus reveal that the geometry of C2H2 changes measurably on coordination to AgCCH. A value of 59(4) Nm−1 is determined for the intermolecular force constant, kσ , confirming that the complex is significantly more strongly bound than hydrogen and halogen-bonded analogues. Ab initio calculations of the re geometry at the CCSD(T)(F12∗)/ACVTZ level of theory are consistent with the experimental results. The spectra of the 107Ag13C13CH and 109Ag13C13CH isotopologues of free silver acetylide are also measured for the first time allowing the geometry of the AgCCH monomer to be examined in greater detail than previously. [ABSTRACT FROM AUTHOR]

Semi-empirical Tight Binding (TB) is known to be a scalable and accurate atomistic representation for electron transport for realistically extended nano-scaled semiconductor devices that might contain millions of atoms. In this paper, an environment-aware and transferable TB model suitable for electronic structure and transport simulations in technologically relevant metals, metallic alloys, metal nanostructures, and metallic interface systems are described. Part I of this paper describes the development and validation of the new TB model. The new model incorporates intra-atomic diagonal and off-diagonal elements for implicit self-consistency and greater transferability across bonding environments. The dependence of the on-site energies on strain has been obtained by appealing to the Moments Theorem that links closed electron paths in the system to energy moments of angular momentum resolved local density of states obtained ab initio. The model matches self-consistent density functional theory electronic structure results for bulk face centered cubic metals with and without strain, metallic alloys, metallic interfaces, and metallic nanostructures with high accuracy and can be used in predictive electronic structure and transport problems in metallic systems at realistically extended length scales. [ABSTRACT FROM AUTHOR]

Several decades ago, Mo2B was assumed to have an Al2Cu-type structure with I4/mcm space group. Using ab initio phonon calculations, we identify the earlier proposed Al2Cu-type structure is dynamically unstable at ambient pressure. An energetically more favorable phase with the tetragonal I4/m structure was then predicted by employing frozen-phonon technique. The currently predicted I4/m phase is mechanically and dynamically stable and energetically more favorable than that of the earlier proposed Al2Cu-type structure. The electronic structures calculations indicate that Mo2B is a metal with several bands crossing the Fermi level. Our analysis indicates that the three-dimensional network of the covalent Mo-B bond is responsible for the ultra-incompressible property of Mo2B. [ABSTRACT FROM AUTHOR]

The first theoretical exploration of superhalogen properties of polynuclear structures based on pseudohalogen ligand is reported here via a case study on eight triply-bridged [Mg2(CN)5]− clusters. From our high-level ab initio results, all these clusters are superhalogens due to their high vertical electron detachment energies (VDE), of which the largest value is 8.67 eV at coupled-cluster single double triple (CCSD(T)) level. Although outer valence Green's function results are consistent with CCSD(T) in most cases, it overestimates the VDEs of three anions dramatically by more than 1 eV. Therefore, the combined usage of several theoretical methods is important for the accuracy of purely theoretical prediction of superhalogen properties of new structures. Spatial distribution of the extra electron of high-VDE anions here indicates two features: remarkable aggregation on bridging CN units and non-negligible distribution on every CN unit. These two features lower the potential and kinetic energies of the extra electron respectively and thus lead to high VDE. Besides superhalogen properties, the structures, relative stabilities and thermodynamic stabilities with respect to detachment of CN-1 were also investigated for these anions. The collection of these results indicates that polynuclear structures based on pseudohalogen ligand are promising candidates for new superhalogens with enhanced properties. [ABSTRACT FROM AUTHOR]

The melting curve of H2O from 49 to 72 GPa was determined by using a laser-heated diamond anvil cell. Double-sided CO2 laser heating technique was employed in order to heat the sample directly. Discontinuous changes of the heating efficiency attributed to the H2O melting were observed between 49 and 72 GPa. The obtained melting temperatures at 49 and 72 GPa are 1200 and 1410 K, respectively. We found that the slope of the melting curve significantly decreases with increasing pressure, only 5 K/GPa at 72 GPa while 44 K/GPa at 49 GPa. Our results suggest that the melting curve does not intersect with the isentropes of Uranus and Neptune, and hence, H2O should remain in the liquid state even at the pressure and temperature conditions found deep within Uranus and Neptune. [ABSTRACT FROM AUTHOR]

Using ab initio molecular dynamics, the atomic structure and transport properties of eutectic Ga-In and Ga-In-Sn are investigated. The Kubo-Greenwood (K-G) and the Ziman-Faber (Z-F) formulations and the Wiedemann-Franz (W-F) law are used for the electrical and electronic thermal conductivity. The species diffusivity and the viscosity are also predicted using the mean square displacement and the Stokes-Einstein (S-E) relation. Alloying Ga causes more disordered structure, i.e., broadening the atomic distance near the In and Sn atoms, which reduces the transport properties and the melting temperature. The K-G treatment shows excellent agreement with the experimental results while Z-F treatment formula slightly overestimates the electrical conductivity. The predicted thermal conductivity also shows good agreement with the experiments. The species diffusivity and the viscosity are slightly reduced by the alloying of Ga with In and Sn atoms. Good agreements are found with available experimental results and new predicted transport-property results are provided. [ABSTRACT FROM AUTHOR]

We have performed new direct-potential-fit (DPF) analyses of the rotationally resolved A1Πu(v' = 2, 3; J' = 1, 2) → X 1Σ+g (v" ϵ [0, 11]; J " ϵ [0, 3]) stimulated emission pumping spectra of Be2 [J. M. Merritt, V. E. Bondybey, and M. C. Heaven, Science 324, 1548 (2009)] using two quite different analytical potential energy functions that incorporate the correct theoretically known long-range behaviour in different ways. These functions are: the damped Morse/long-range potential [R. J. Le Roy, C. C. Haugen, J. Tao, and H. Li, Mol. Phys. 109, 435 (2011)], and the Chebyshev polynomial expansion potential [L. Busevica, I. Klincare, O. Nikolayeva, M. Tamanis, R. Ferber, V. V. Meshkov, E. A. Pazyuk, and A. V. Stolyarov, J. Chem. Phys. 134, 104307 (2011)]. In contrast with the expanded Morse oscillator potential determined in the original DPF analysis of Merritt et al. [Science 324, 1548 (2009)], both of these functions unambiguously support the existence of the v" = 11 last vibrational levels which is bound by only ∼0.5 cm-1, and they give equivalent, essentially exact predictions for this level when using the original data set which ended at v" = 10. These empirical potentials predict an equilibrium distance of re = 2.445(5) Å and a well depth of De = 934.9(0.4) cm-1, values which agree (within the uncertainties) with the best ab initio estimates of 2.444(10) Å and 935(10) cm-1, respectively [J. Koput, Phys. Chem. Chem. Phys. 13, 20311 (2011)]. [ABSTRACT FROM AUTHOR]

Based exclusively on high-level ab initio calculations, a new full-dimensional analytical potential energy surface (PES-2014) for the gas-phase reaction of hydrogen abstraction from methane by an oxygen atom is developed. The ab initio information employed in the fit includes properties (equilibrium geometries, relative energies, and vibrational frequencies) of the reactants, products, saddle point, points on the reaction path, and points on the reaction swath, taking especial caution respecting the location and characterization of the intermediate complexes in the entrance and exit channels. By comparing with the reference results we show that the resulting PES-2014 reproduces reasonably well the whole set of ab initio data used in the fitting, obtained at the CCSD(T) = FULL/aug-ccpVQZ// CCSD(T) = FC/cc-pVTZ single point level, which represents a severe test of the new surface. As a first application, on this analytical surface we perform an extensive dynamics study using quasi-classical trajectory calculations, comparing the results with recent experimental and theoretical data. The excitation function increases with energy (concave-up) reproducing experimental and theoretical information, although our values are somewhat larger. The OH rotovibrational distribution is cold in agreement with experiment. Finally, our results reproduce experimental backward scattering distribution, associated to a rebound mechanism. These results lend confidence to the accuracy of the new surface, which substantially improves the results obtained with our previous surface (PES-2000) for the same system. [ABSTRACT FROM AUTHOR]

Ab initio and density functional theory calculations are used to examine complexes formed between H2CO and H2CS with 1, 2, and 3 molecules of SO2. The nature of the interactions is probed by a variety of means, including electrostatic potentials, natural bond orbital, atoms in molecules, energy decomposition, and electron density redistribution maps. The dimers are relatively strongly bound, with interaction energies exceeding 5 kcal/mol. The structures are cyclic, containing both a O/S· · ·S chalcogen bond and a CH· · ·O H-bond. Addition of a second SO2 molecule leads to a variety of heterotrimer structures, most of which resemble the original dimer, where the second SO2 molecule engages in a chalcogen bond with the first SO2, and a C· · ·O attraction with the H2CX. Some cooperativity is apparent in the trimers and tetramers, with an attractive three-body interaction energy and shortened intermolecular distances. [ABSTRACT FROM AUTHOR]

Long known as a fully polarized band with a near vanishing depolarization ratio [ηs = 0.05, W. Holzer and R. Ouillon, Chem. Phys. Lett. 24, 589 (1974)], the 2ν5 Raman overtone of SF6 has so far been considered as of having a prohibitively weak anisotropic spectrum [D. P. Shelton and L. Ulivi, J. Chem. Phys. 89, 149 (1988)]. Here, we report the first anisotropic spectrum of this overtone, at room temperature and for 13 gas densities ranging between 2 and 27 amagat. This spectrum is 10 times broader and 50 times weaker than the isotropic counterpart of the overtone [D. Kremer, F. Rachet, and M. Chrysos, J. Chem. Phys. 138, 174308 (2013)] and its profile much more sensitive to pressure effects than the profile of the isotropic spectrum. From our measurements an accurate value for the anisotropy matrix-element |〈000020|Δα|000000〉| was derived and this value was found to be comparable to that of the mean-polarizability |〈000020| ...α |000000〉|. Among other conclusions our study offers compelling evidence that, in Raman spectroscopy, highly polarized bands or tiny depolarization ratios are not necessarily incompatible with large polarizability anisotropy transition matrix-elements. Our findings and the way to analyze them suggest that new strategies should be developed on the basis of the complementarity inherent in independent incoherent Raman experiments that run with two different incident-beam polarizations, and on concerted efforts to ab initio calculate accurate data for first and second polarizability derivatives. Values for these derivatives are still rarities in the literature of SF6. [ABSTRACT FROM AUTHOR]

We elucidate the role of room-temperature-induced instantaneous structural distortions in the Li K-edge X-ray absorption spectra (XAS) of crystalline LiF, Li2SO4, Li2O, Li3N, and Li2CO3 using high resolution X-ray Raman spectroscopy (XRS) measurements and first-principles density functional theory calculations within the eXcited electron and Core Hole approach. Based on thermodynamic sampling via ab initio molecular dynamics simulations, we find calculated XAS in much better agreement with experiment than those computed using the rigid crystal structure alone. We show that local instantaneous distortion of the atomic lattice perturbs the symmetry of the Li 1s core-excitedstate electronic structure, broadening spectral line-shapes and, in some cases, producing additional spectral features. The excellent agreement with high-resolution XRS measurements validates the accuracy of our first-principles approach to simulating XAS, and provides both accurate benchmarks for model compounds and a predictive theoretical capability for identification and characterization of multi-component systems, such as lithium-ion batteries, under working conditions. [ABSTRACT FROM AUTHOR]

We present a Nested Markov chain Monte Carlo (NMC) scheme for building equilibrium averages based on accurate potentials such as density functional theory. Metropolis sampling of a reference system, defined by an inexpensive but approximate potential, was used to substantially decorrelate configurations at which the potential of interest was evaluated, thereby dramatically reducing the number needed to build ensemble averages at a given level of precision. The efficiency of this procedure was maximized on-the-fly through variation of the reference system thermodynamic state (characterized here by its inverse temperature β0), which was otherwise unconstrained. Local density approximation results are presented for shocked states of argon at pressures from 4 to 60 GPa, where--depending on the quality of the reference system potential--acceptance probabilities were enhanced by factors of 1.2-28 relative to unoptimized NMC. The optimization procedure compensated strongly for reference potential shortcomings, as evidenced by significantly higher speedups when using a reference potential of lower quality. The efficiency of optimized NMC is shown to be competitive with that of standard ab initio molecular dynamics in the canonical ensemble. [ABSTRACT FROM AUTHOR]

Due to its importance for the function of organic optoelectronic devices, accurate simulations of the singlet exciton diffusion are crucial to predict the performance of new materials. We present a protocol which allows for the efficient directional analysis of exciton transport with high-level ab initio methods. It is based on an alternative to the frequently employed rate equation since the latter was found to be erroneous in some cases. The new approach can be used in combination with the master equation which is considerably faster than the corresponding Monte Carlo approach. The long-range character of the singlet exciton coupling is taken into account by an extrapolation scheme. The approach is applied to singlet exciton diffusion in those substances where these quantities are experimentally best established: naphthalene and anthracene. The high quality of the crystals, furthermore, diminish uncertainties arising from the geometrical structures used in the computations. For those systems, our new approach provides exciton diffusion lengths L for naphthalene and anthracene crystals which show an excellent agreement with their experimental counterparts. For anthracene, for example, the computed L value in a direction is computed to 58 nm while the experimental value is 60 ± 10 nm. [ABSTRACT FROM AUTHOR]

The adsorption of oxygen on bcc Fe–Cr(100) surfaces with two different alloy concentrations is studied using ab initio density functional calculations. Atomic-scale analysis of oxygen–surface interactions is indispensable for obtaining a comprehensive understanding of macroscopic surface oxidation processes. Up to two chromium atoms are inserted into the first two surface layers. Atomic geometries, energies and electronic properties are investigated. A hollow site is found to be the preferred adsorption site over bridge and on-top sites. Chromium atoms in the surface and subsurface layers are found to significantly affect the adsorption properties of neighbouring iron atoms. Seventy-one different adsorption geometries are studied, and the corresponding adsorption energies are calculated. Estimates for the main diffusion barriers from the hollow adsorption site are given. Whether the change in the oxygen affinity of iron atoms can be related to the chromium-induced charge transfer between the surface atoms is discussed. The possibility to utilize the presented theoretical results in related experimental research and in developing semiclassical potentials for simulating the oxidation of Fe–Cr alloys is addressed. [ABSTRACT FROM AUTHOR]

Electrostatics of depolarization field Ed in relation to the polarization is studied. In particular, the value of permittivity for Ed (εd) in prototypical situations of ferroelectrics, including Mehta formula, is examined by ab initio calculations. By using spontaneous polarization PS corresponding to accurate experiment ones, we show εd = 1, which suggests that the results of εd ≫ 1 indicate hidden mechanisms; εd = 1 suggests that the effect of Ed is significant to induce intriguing important phenomena overlooked by εd ≫ 1. A bridge between εd = 1 and εd ≫ 1, i.e. the consistency of εd = 1 with conventional results is presented. The exact electrostatic equality of head-to-head–tail-to-tail domains to free-standing ferroelectrics is deduced. Hence, most stoichiometric clean freestanding monodomain ferroelectrics and head-to-head–tail-to-tail domains are shown unstable regardless of size, unless partially metallic. This verifies the previous results in a transparent manner. This conclusion is shown consistent with a recent hyperferroelectric LiBeSb and "freestanding" monolayer ferroelectrics, of which origin is suggested to be adsorbates. In addition, this restriction is suggested to break in externally strained ultrathin ferroelectrics. The macroscopic formulas of Ed are found valid down to a several unit-cells, when electronic and atomic-scale surface effects are unimportant and accurate PS is used. [ABSTRACT FROM AUTHOR]

Using accurate dissipative DFT-NEGF atomistic-simulation techniques within the Wannier-Function formalism, we give a fresh look at the possibility of sub-10-nm scaling for high-performance complementary metal oxide semiconductor (CMOS) applications. We show that a combination of good electrostatic control together with high mobility is paramount to meet the stringent roadmap targets. Such requirements typically play against each other at sub-10-nm gate length for MOS transistors made of conventional semiconductor materials like Si, Ge, or III–V and dimensional scaling is expected to end ~12 nm gate-length (pitch of 40 nm). We demonstrate that using alternative 2D channel materials, such as the less-explored HfS2 or ZrS2, high-drive current down to ~6 nm is, however, achievable. We also propose a dynamically doped field-effect transistor concept, that scales better than its MOSFET counterpart. Used in combination with a high-mobility material such as HfS2, it allows for keeping the stringent high-performance CMOS on current and competitive energy-delay performance, when scaling down to virtually 0 nm gate length using a single-gate architecture and an ultra-compact design (pitch of 22 nm). The dynamically doped field-effect transistor further addresses the grand-challenge of doping in ultra-scaled devices and 2D materials in particular. [ABSTRACT FROM AUTHOR]

Higher vibrational levels of the C2 molecule than those observed so far were investigated for the X1Σg+, A1Πu, a3Πu, c3Σu+, and d3Πg states through the Phillips, Swan, and d3Πg-c3Σu+ band systems under a jet-cooled condition. The term values and the molecular constants for 21 new vibronic levels were determined from rotationally resolved excitation spectra. The determined term values and rotational constants were compared to those derived from high-level ab initio potential curves. Perturbations identified in low J levels of the d3Πg(v = 8) state are most likely to be caused by the 15Πg(v = 3) state. [ABSTRACT FROM AUTHOR]

Using state-of-the-art theoretical methods, we investigate the lowest electronic states of singlet and triplet spin multiplicities of HSNO. These computations are done using configuration interaction ab initio methods and the aug-cc-pV5Z basis set. One-dimensional cuts of the six-dimensional potential energy surfaces of these electronic states along the HS, SN stretches and HSN, SNO bending and torsion coordinates are calculated. Several avoided crossings and conical intersections are found. We computed also radiative lifetimes and spin-orbit couplings of these electronic states. Our work shows that the dynamics on these excited states is very complex, and suggest that multi-step mechanisms will populate the ground state via radiationless processes or lead to predissociation or intramolecular isomerization. For instance, these potentials are used to propose mechanisms for the IR, Vis, and UV light-induced cis-trans interconversions of HSNO and reactivity towards HS + NO and H + SNO products. Our findings are in good agreement with previous experimental studies on the photochemistry of HSNO. The atmospheric implication of HSNO is also discussed. [ABSTRACT FROM AUTHOR]