XMVB: a program for ab initio nonorthogonal valence bond computations.
Song, Lingchun; Mo, Yirong; Zhang, Qianer; Wu, Wei
2005-04-15
An ab initio nonorthogonal valence bond program, called XMVB, is described in this article. The XMVB package uses Heitler-London-Slater-Pauling (HLSP) functions as state functions, and calculations can be performed with either all independent state functions for a molecule or preferably a few selected important state functions. Both our proposed paired-permanent-determinant approach and conventional Slater determinant expansion algorithm are implemented for the evaluation of the Hamiltonian and overlap matrix elements among VB functions. XMVB contains the capabilities of valence bond self-consistent field (VBSCF), breathing orbital valence bond (BOVB), and valence bond configuration interaction (VBCI) computations. The VB orbitals, used to construct VB functions, can be defined flexibly in the calculations depending on particular applications and focused problems, and they may be strictly localized, delocalized, or bonded-distorted (semidelocalized). The parallel version of XMVB based on MPI (Message Passing Interface) is also available. PMID:15704237
Ab initio study of the bonding in diatomic nickel
Noell, J.O.; Newton, M.D.; Hay, P.J.; Martin, R.L.; Bobrowicz, F.W.
1980-09-01
Hartree--Fock, GVB, and configuration interaction calculations were performed for diatomic nickel using an ab initio effective core potential. A basis set specifically optimized for the /sup 3/D state of atomic nickel is found to be far superior to the more common basis obtained from the /sup 3/F atomic state. Correlation effects are found to be significant in determining the bond energy. In particular, the two electrons of the s--s bond must be appropriately correlated. In addition, correlation effects which one would interpret as being principally intra-atomic in character are found to have a marked effect on the molecular properties. The theoretically predicted bond dissociation energy (D/sub e/) of 43.4 kcal/mol is significantly lower than the experimental estimate of 55 +- 5 kcal/mol. However, molecular partition functions calculated using the present results indicate that the experimental value should be revised downward to a value of approx.46 +- 5 kcal/mol, in good agreement with our calculations. An interatomic distance of 4.27 bohr is computed and compared with experimental estimates. Spectroscopic parameters for dipole-allowed transitions from the ground state were determined from SCF and GVB calculations and discussed in relation to the experimentally observed visible and ultraviolet spectra attributed to Ni/sub 2/.
Ab initio SCF calculations on hydrogen bonded cresol isomers
NASA Astrophysics Data System (ADS)
Pohl, M.; Kleinermanns, K.
1988-12-01
Ab initio GAUSSIAN 80 calculations with two different basis sets (STO-3G and 4 31 G*) were performed on hydrogen bonded cresol isomers for comparison with experimental data from free jet fluorescence excitation spectroscopy. For m-cresol, the calculated barriers for hindered internal rotation of the OH-group and the CH3-group are in good agreement with experiment. The calculations show the trans-linear configuration of p-cresol· B-clusters ( B = H2O, CH3OH) to be more stable than the all-planar configuration. This agrees with CI calculations and microwave spectroscopic investigations of the water dimer. Calculations of both the intermolecular stretch and bend frequencies of p-cresol· B-clusters show little dependence on the all-planar or trans-linear configuration but a strong dependence on the choice of the basis set. With the minimal basis set STO-3G, the vibrational energies are generally too high. The agreement between the calculated vibrational frequencies from the 4 31 G* basis set and the experimental values is fair.
Patil, Amol Baliram; Bhanage, Bhalchandra Mahadeo
2016-06-21
The nature of bonding interactions between the cation and the anion of an ionic liquid is at the heart of understanding ionic liquid properties. A particularly interesting case is a special class of ionic liquids known as protic ionic liquids. The extent of proton transfer in protic ionic liquids has been observed to vary according to the interacting species. Back proton transfer renders protic ionic liquids volatile and to be considered as inferior ionic liquids. We try to address this issue by employing modern ab initio valence bond theory calculations. The results indicate that the bonding in the cation and the anion of a prototypical ionic liquid, ethylammonium nitrate, is fundamentally different. It is neither characteristic of covalent/polar covalent bonding nor ionic bonding but rather charge shift bonding as a resonance hybrid of two competing ionic molecular electronic structure configurations. An investigation of other analogous protic ionic liquids reveals that this charge shift bonding seems to be a typical characteristic of protic ionic liquids while the ionic solid analogue compound ammonium nitrate has less charge shift bonding character as compared to protic ionic liquids. Further the extent of charge shift bonding character has been found to be congruent with the trends in many physicochemical properties such as melting point, conductivity, viscosity, and ionicity of the studied ionic liquids indicating that percentage charge shift character may serve as a key descriptor for large scale computational screening of ionic liquids with desired properties. PMID:27229870
Red-Shifting versus Blue-Shifting Hydrogen Bonds: Perspective from Ab Initio Valence Bond Theory.
Chang, Xin; Zhang, Yang; Weng, Xinzhen; Su, Peifeng; Wu, Wei; Mo, Yirong
2016-05-01
Both proper, red-shifting and improper, blue-shifting hydrogen bonds have been well-recognized with enormous experimental and computational studies. The current consensus is that there is no difference in nature between these two kinds of hydrogen bonds, where the electrostatic interaction dominates. Since most if not all the computational studies are based on molecular orbital theory, it would be interesting to gain insight into the hydrogen bonds with modern valence bond (VB) theory. In this work, we performed ab initio VBSCF computations on a series of hydrogen-bonding systems, where the sole hydrogen bond donor CF3H interacts with ten hydrogen bond acceptors Y (═NH2CH3, NH3, NH2Cl, OH(-), H2O, CH3OH, (CH3)2O, F(-), HF, or CH3F). This series includes four red-shifting and six blue-shifting hydrogen bonds. Consistent with existing findings in literature, VB-based energy decomposition analyses show that electrostatic interaction plays the dominating role and polarization plays the secondary role in all these hydrogen-bonding systems, and the charge transfer interaction, which denotes the hyperconjugation effect, contributes only slightly to the total interaction energy. As VB theory describes any real chemical bond in terms of pure covalent and ionic structures, our fragment interaction analysis reveals that with the approaching of a hydrogen bond acceptor Y, the covalent state of the F3C-H bond tends to blue-shift, due to the strong repulsion between the hydrogen atom and Y. In contrast, the ionic state F3C(-) H(+) leads to the red-shifting of the C-H vibrational frequency, owing to the attraction between the proton and Y. Thus, the relative weights of the covalent and ionic structures essentially determine the direction of frequency change. Indeed, we find the correlation between the structural weights and vibrational frequency changes. PMID:27074500
Hydrogen Bonds in Crystalline Imidazoles Studied by 15N NMR and ab initio MO Calculations
NASA Astrophysics Data System (ADS)
Ueda, Takahiro; Nagatomo, Shigenori; Masui, Hirotsugu; Nakamura, Nobuo; Hayashi, Shigenobu
1999-07-01
Intermolecular hydrogen bonds of the type N-H...N in crystals of imidazole and its 4-substituted and 4,5-disubstituted derivatives were studied by 15N CP/MAS NMR and an ab initio molecular orbital (MO) calculation. In the 15N CP/MAS NMR spectrum of each of the imidazole derivatives, two peaks due to the two different functional groups, >NH and =N-, were observed. The value of the 15N isotropic chemical shift for each nitrogen atom depends on both the length of the intermolecular hydrogen bond and the kind of the substituent or substituents. It was found that the difference between the experimen-tal chemical shifts of >NH and =N-varies predominantly with the hydrogen bond length but does not show any systematic dependence on the kind of substituent. The ab initio MO calculations suggest that the hydrogen bond formation influences the 15N isotropic chemical shift predominantly, and that the difference between the 15N isotropic chemical shift of >NH and =N-varies linearly with the hydrogen bond length.
An ab initio molecular dynamics study on hydrogen bonds between water molecules
NASA Astrophysics Data System (ADS)
Pan, Zhang; Chen, Jing; Lü, Gang; Geng, Yi-Zhao; Zhang, Hui; Ji, Qing
2012-04-01
The quantitative estimation of the total interaction energy of a molecular system containing hydrogen bonds (H bonds) depends largely on how to identify H bonding. The conventional geometric criteria of H bonding are simple and convenient in application, but a certain amount of non-H bonding cases are also identified as H bonding. In order to investigate the wrong identification, we carry out a systematic calculation on the interaction energy of two water molecules at various orientation angles and distances using ab initio molecular dynamics method with the dispersion correction for the Becke-Lee-Yang-Parr (BLYP) functionals. It is shown that, at many orientation angles and distances, the interaction energies of the two water molecules exceed the energy criterion of the H bond, but they are still identified as H-bonded by the conventional "distance-angle" criteria. It is found that in these non-H bonding cases the wrong identification is mainly caused by short-range interaction between the two neighbouring water molecules. We thus propose that, in addition to the conventional distance and angle criteria of H bonding, the distance dHṡṡṡH between the two neighbouring hydrogen atoms of the two water molecules should also be taken as a criterion, and the distance rOṡṡṡH between the hydrogen atom of the H-bond donor molecule and the oxygen atom of the acceptor molecule should be restricted by a lower limit. When dHṡṡṡH and rOṡṡṡH are small (e.g., dHṡṡṡH < 2.0 Å and rOṡṡṡH < 1.62 Å), the repulsion between the two neighbouring atoms increases the total energy of the two water molecules dramatically and apparently weakens the binding of the water dimer. A statistical analysis and comparison of the numbers of the H bonds identified by using different criteria have been conducted on a Car-Parrinello ab initio molecular dynamics simulation with dispersion correction for a system of 64 water molecules at near-ambient temperature. They
Probing Defects and Correlations in the Hydrogen-Bond Network of ab Initio Water.
Gasparotto, Piero; Hassanali, Ali A; Ceriotti, Michele
2016-04-12
The hydrogen-bond network of water is characterized by the presence of coordination defects relative to the ideal tetrahedral network of ice, whose fluctuations determine the static and time-dependent properties of the liquid. Because of topological constraints, such defects do not come alone but are highly correlated coming in a plethora of different pairs. Here we discuss in detail such correlations in the case of ab initio water models and show that they have interesting similarities to regular and defective solid phases of water. Although defect correlations involve deviations from idealized tetrahedrality, they can still be regarded as weaker hydrogen bonds that retain a high degree of directionality. We also investigate how the structure and population of coordination defects is affected by approximations to the interatomic potential, finding that, in most cases, the qualitative features of the hydrogen-bond network are remarkably robust. PMID:26881726
NASA Astrophysics Data System (ADS)
Maes, G.; Smets, J.; Adamowicz, L.; McCarthy, W.; Van Bael, M. K.; Houben, L.; Schoone, K.
1997-06-01
Correlations between selected ab initio predicted and experimentally observed properties of 1:1 H-bonded complexes of pyridines, pyrimidines, and imidazoles with water are investigated. Relationships are found between the experimental properties of proton affinity and water frequency shift, and the ab initio calculated bond distances, interaction energies and water frequency shifts. It is also found that well-defined relations can be established between calculated and observed properties for the pyridine complexes, but these cannot be reliably extended to the other N-base systems. The similarities demonstrate that the presently available ab initio methods are useful in predicting the experimental behaviour of H-bonded systems, but only for closely related molecules.
Klevets, Ivan; Bryk, Taras
2014-12-07
Electron-ion structure factors, calculated in ab initio molecular dynamics simulations, are reported for several binary liquids with different kinds of chemical bonding: metallic liquid alloy Bi–Pb, molten salt RbF, and liquid water. We derive analytical expressions for the long-wavelength asymptotes of the partial electron-ion structure factors of binary systems and show that the analytical results are in good agreement with the ab initio simulation data. The long-wavelength behaviour of the total charge structure factors for the three binary liquids is discussed.
Revealing halogen bonding interactions with anomeric systems: an ab initio quantum chemical studies.
Lo, Rabindranath; Ganguly, Bishwajit
2015-02-01
A computational study has been performed using MP2 and CCSD(T) methods on a series of O⋯X (X=Br, Cl and I) halogen bonds to evaluate the strength and characteristic of such highly directional noncovalent interactions. The study has been carried out on a series of dimeric complexes formed between interhalogen compounds (such as BrF, BrCl and BrI) and oxygen containing electron donor molecule. The existence and consequences of the anomeric effect of the electron donor molecule has also been investigated through an exploration of halogen bonding interactions in this halogen bonded complexes. The ab initio quantum chemical calculations have been employed to study both the nature and directionality of the halogen molecules toward the sp(3) oxygen atom in anomeric systems. The presence of anomeric nO→σ*CN interaction involves a dominant role for the availability of the axial and equatorial lone pairs of donor O atom to participate with interhalogen compounds in the halogen-bonded complexes. The energy difference between the axial and equatorial conformers with interhalogen compounds reaches up to 4.60 kJ/mol, which however depends upon the interacting halogen atoms and its attaching atoms. The energy decomposition analysis further suggests that the total halogen bond interaction energies are mainly contributed by the attractive electrostatic and dispersion components. The role of substituents attached with the halogen atoms has also been evaluated in this study. With the increase of halogen atom size and the positive nature of σ-hole, the halogen atom interacted more with the electron donor atom and the electrostatic contribution to the total interaction energy enhances appreciably. Further, noncovalent interaction (NCI) studies have been carried out to locate the noncovalent halogen bonding interactions in real space. PMID:25522359
An ab initio investigation of some hydrogen-bonded complexes of methanethiol
NASA Astrophysics Data System (ADS)
Bricknell, Bradley C.; Ford, Thomas A.
2010-07-01
The hydrogen-bonded complexes formed between methanethiol, as proton donor, and trimethyl amine, dimethyl ether, methyl fluoride, trimethyl phosphine, dimethyl sulphide and methyl chloride, as proton acceptors, have been studied by means of ab initio calculations at the second order level of Møller-Plesset perturbation theory and using the 6-311++G(d,p) basis set. The hydrogen bond energies were found to correlate with the gas phase basicities of the proton acceptors. The increases of the SH bond lengths, the shifts of the wavenumbers of the SH stretching modes and the increases of their infrared intensities were also found to be dependent on the interaction energies. The results for the CH 3SH·S(CH 3) 2 complex are anomalous, and this is attributed to deficiencies with the basis set used for sulphur. The data for the CH 3SH complex closely mirror the behaviour of the corresponding complexes of CH 3OH and (CH 3) 2NH, with the perturbations of the various properties varying consistently in the order CH 3OH > (CH 3) 2NH > CH 3SH.
Ab initio path-integral molecular dynamics and the quantum nature of hydrogen bonds
NASA Astrophysics Data System (ADS)
Yexin, Feng; Ji, Chen; Xin-Zheng, Li; Enge, Wang
2016-01-01
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. Project supported by the National Natural Science Foundation of China (Grant Nos. 11275008, 91021007, and 10974012) and the China Postdoctoral Science Foundation (Grant No. 2014M550005).
Matrix isolation infrared and ab initio study of the hydrogen bonding between formic acid and water
NASA Astrophysics Data System (ADS)
George, Lisa; Sander, Wolfram
2004-11-01
The infrared spectra of the formic acid-water complexes isolated in argon matrices are reported. Both supersonic jet expansion and a conventional effusive source followed by trapping in solid argon at 10 K are used to obtain the matrices. The experimental IR spectra are compared to the data obtained from high level ab initio (MP2) and DFT (B3LYP) calculations with 6-311++G(d,p) and aug-cc-pVTZ basis sets. The complex formation results in red shifts in the CO and OH stretching vibrations and a blue shift in the CO stretching vibration of formic acid. The OH stretching modes of water also exhibit pronounced red shifts. Both the MP2 and B3LYP calculations located three minima corresponding to cyclic HCOOH⋯H 2O complexes with two hydrogen bond interactions. The binding energies are -10.3, -5.1, and -3.5 kcal mol -1, respectively, for the three complexes at the MP2/ aug-cc-pVTZ level, corrected for the basis set superposition error (BSSE) using the Boys-Bernardi counterpoise scheme. Comparison of the calculated frequencies of the three complexes with the matrix IR spectrum reveals that the lowest energy complex is formed. In addition, a complex of formic acid with two water molecules is observed.
Bauzá, Antonio; Frontera, Antonio
2016-07-27
In this manuscript the ability of hydrogen and halogen bonding interactions, as well as metal coordination to enhance FF interactions involving fluorine substituted aromatic rings has been studied at the RI-MP2/def2-TZVPD level of theory. We have used 4-fluoropyridine, 4-fluorobenzonitrile, 3-(4-fluorophenyl)propiolonitrile and their respective meta derivatives as aromatic compounds. In addition, we have used HF and IF as hydrogen and halogen bond donors, respectively, and Ag(i) as the coordination metal. Furthermore, we have also used HF as an electron rich fluorine donor entity, thus establishing FF interactions with the above mentioned aromatic systems. Moreover, a CSD (Cambridge Structural Database) search has been carried out and some interesting examples have been found, highlighting the impact of FF interactions involving aromatic fluorine atoms in solid state chemistry. Finally, cooperativity effects between FF interactions and both hydrogen and halogen bonding interactions have been analyzed and compared. We have also used Bader's theory of "atoms in molecules" to further describe the cooperative effects. PMID:27401347
Wikfeldt, K. T.; Michaelides, A.
2014-01-28
Ab initio simulations that account for nuclear quantum effects have been used to examine the order-disorder transition in squaric acid, a prototypical H-bonded antiferroelectric crystal. Our simulations reproduce the >100 K difference in transition temperature observed upon deuteration as well as the strong geometrical isotope effect observed on intermolecular separations within the crystal. We find that collective transfer of protons along the H-bonding chains – facilitated by quantum mechanical tunneling – is critical to the order-disorder transition and the geometrical isotope effect. This sheds light on the origin of isotope effects and the importance of tunneling in squaric acid which likely extends to other H-bonded ferroelectrics.
NASA Astrophysics Data System (ADS)
Seller, Harrell; Shustorovich, Evgeny
1996-02-01
Binding energies for sulfur oxides, SO x, x = 1-3, have been determined for several coordination modes on silver, gold and palladium surfaces employing ab initio quantum chemical methods and the bond order conservation Morse potential (BOC-MP) method. SO 2 coordination was studied in the most detail. In general the agreement between the BOC-MP and ab initio binding energies is good for the (111) surfaces of silver and palladium with both methods predicting that, in the zero coverage limit, di-coordination via S,O and O,O will be more favorable energetically than mono-coordination via S. In the case of chemisorption on the Pd (110) surface the two methods agree well for the cases in which there are formulas for the BOC-MP binding energies. In going from the (111) surfaces to the (110) surfaces of silver and palladium the ab initio calculations predict that the preferred chemisorption site shifts from the bridge site to the hollow site. On the silver surfaces the net charge transferred to the adsorbate as judged from the Mulliken populations correlates roughly with the binding energy. No significant charge transfer was found on the palladium surfaces. Our SO 2 chemisorption calculations indicate that the work functions of the metal surfaces examined should increase upon mono-S adsorption, increase to a lesser extent upon di S,O adsorption and may even decrease upon di O,O adsorption. Ab initio calculations provide evidence of the existence of SO 2 surface dimers. The binding energy predicted by the BOC-MP model for SO 3 in the bridging site agrees well with the ab initio result for SO 3 di-coordinated in the long bridge of the Ag(110) surface. The methods yield similar predictions for the case of SO on silver. Our modeling provides a coherent picture consistent with many aspects of the experimental literature. We present some model predictions, particularly the di O,O coordination mode for SO 2, that require verification experimentally.
NASA Astrophysics Data System (ADS)
Pati, Ranjit; Karna, Shashi P.
2002-01-01
The dependence of electron transfer (ET) coupling element, VAB, on the length of rigid-rod-like systems consisting of bicyclo[1.1.1]pentane (BCP), cubane (CUB), and bicyclo[2.2.2]octane (BCO) monomers, has been investigated with the use of ab initio Hartree-Fock (HF) method employing Marcus-Hush two-state (TS) model. The value of VAB decreases exponentially with increase in the number of the cage units of the σ-bonded molecules. The calculated decay constant, β, shows good agreement with previously reported data. For molecular length⩾15 Å, the value of VAB becomes negligibly small, suggesting complete suppression of the through bond direct tunneling contribution to ET process.
Walker, Brent; Michaelides, Angelos
2010-11-01
The impact of quantum nuclear effects on hydrogen (H-) bond strength has been inferred in earlier work from bond lengths obtained from path integral molecular dynamics (PIMD) simulations. To obtain a direct quantitative assessment of such effects, we use constrained-centroid PIMD simulations to calculate the free energy changes upon breaking the H-bonds in dimers of HF and water. Comparing ab initio simulations performed using PIMD and classical nucleus molecular dynamics (MD), we find smaller dissociation free energies with the PIMD method. Specifically, at 50 K, the H-bond in (HF)(2) is about 30% weaker when quantum nuclear effects are included, while that in (H(2)O)(2) is about 15% weaker. In a complementary set of simulations, we compare unconstrained PIMD and classical nucleus MD simulations to assess the influence of quantum nuclei on the structures of these systems. We find increased heavy atom distances, indicating weakening of the H-bond consistent with that observed by direct calculation of the free energies of dissociation. PMID:21054031
Liu, Zhong-Li; Cai, Ling-Cang; Chen, Xiang-Rong; Wu, Qiang; Jing, Fu-Qian
2009-03-01
We report a detailed ab initio study for body-centered-cubic (bcc) Ta within the framework of the quasiharmonic approximation (QHA) to refine its thermal equation of state and thermodynamic properties. Based on the excellent agreement of our calculated phonon dispersion curve with experiment, the accurate thermal equations of state and thermodynamic properties are well reproduced. The thermal equation of state (EOS) and EOS parameters are considerably improved in our work compared with previous results by others. Furthermore, at high temperatures, the excellent agreement of our obtained thermal expansion and Hugoniot curves with experiments greatly verifies the validity of the quasiharmonic approximation at higher temperatures. It is known that pressure suppresses the vibrations of atoms from their equilibrium positions, i.e. the bondings among atoms are strengthened by pressure; for the same temperature, anharmonicity becomes less important at high pressure. Thus the highest valid temperature of the QHA can be reasonably extended to the larger range. PMID:21817394
Feng, Ya-Juan; Huang, Teng; Wang, Chao; Liu, Yi-Rong; Jiang, Shuai; Miao, Shou-Kui; Chen, Jiao; Huang, Wei
2016-07-14
Molecular level insight into the interaction between volatile organic compounds (VOCs) and aerosols is crucial for improvement of atmospheric chemistry models. In this paper, the interaction between adsorbed toluene, one of the most significant VOCs in the urban atmosphere, and the aqueous surface of aerosols was studied by means of combined molecular dynamics simulations and ab initio quantum chemistry calculations. It is revealed that toluene can be stably adsorbed on the surface of aqueous droplets via hydroxyl-π hydrogen bonding between the H atoms of the water molecules and the C atoms in the aromatic ring. Further, significant modifications on the electrostatic potential map and frontier molecular orbital are induced by the solvation effect of surface water molecules, which would affect the reactivity and pathway of the atmospheric photooxidation of toluene. This study demonstrates that the surface interactions should be taken into consideration in the atmospheric chemical models on oxidation of aromatics. PMID:27280740
Chemical Bonding in Si52- and NaSi5- via Photoelectron Spectroscopy and Ab Initio Calculations
Yu, Dmitry; Boldyrev, Zubarev; Boldyrev, Alexander I.; Li, Xi; Cui, Lifeng; Wang, Lai S.
2005-12-22
Photoelectron spectroscopy and ab initio calculations are used to investigate the electronic structure and chemical bonding of Si5 - and Si5 2- in NaSi5 -. Photoelectron spectra of Si5 - and NaSi5 - are obtained at several photon energies and are compared with theoretical calculations at four different levels of theory, TD-B3LYP, R(U)OVGF, UCCSD(T), and EOM-CCSD(T), all with 6-311+G(2df) basis sets. Excellent agreement is observed between experiment and theory, confirming the obtained ground-state structures for Si5 - and Si5 2-, which are both found to be trigonal bipyramid with D3h symmetry at several levels of theory. Chemical bonding in Si5, Si5 -, and Si5 2- is analyzed using NPA, molecular orbitals, ELF, and NICS indices. The bonding in Si5 2- is compared with that in the isoelectronic and isostructural B5H5 2- species, but they are found to differ due to the involvement of electron densities, which are supposed to be lone pairs in the skeletal bonding in Si5 2-.
Gorenstein, D G; Taira, K
1984-01-01
Ab initio molecular orbital calculations have been performed on the reaction profile for the addition/elimination reaction between ammonia and formic acid, proceeding via a tetrahedral intermediate: NH3 + HCO2H----H2NCH(OH)2----NH2CHO + H2O. Calculated transition state energies for the first addition step of the reaction revealed that a lone pair on the oxygen of the OH group, which is antiperiplanar to the attacking nitrogen, stabilized the transition state by 3.9 kcal/mol, thus supporting the hypothesis of stereoelectronic control for this reaction. In addition, a secondary, counterbalancing stereoelectronic effect stabilizes the second step, water elimination, transition state by 3.1 kcal/mol if the lone pair on the leaving water oxygen is not antiperiplanar to the C-N bond. The best conformation for the transition states was thus one with a lone pair antiperiplanar to the adjacent scissile bond and also one without a lone-pair orbital on the scissile bond oxygen or nitrogen antiperiplanar to the adjacent polar bond. The significance of these stereoelectronic effects for the mechanism of action of serine proteases is discussed. PMID:6394065
Song, Lingchun; Gao, Jiali
2009-01-01
A theoretical model is presented for deriving effective diabatic states based on ab initio self-consistent field valence bond (VBSCF) theory by reducing the multi-configurational VB Hamiltonian into an effective two-state model. We describe two computational approaches for the optimization of the effective diabatic configurations, resulting in two ways of interpreting such effective diabatic states. In the variational diabatic configuration (VDC) method, the energies of the individual diabatic states are variationally minimized. In the consistent diabatic configuration (CDC) method, both the configuration coefficients and orbital coefficients are simultaneously optimized to minimize the adiabatic ground-state energy in VBSCF calculations. In addition, we describe a mixed molecular orbital and valence bond (MOVB) approach to construct the CDC diabatic and adiabatic states for a chemical reaction, whereas the VDC-MOVB method has been described previously. Employing the symmetric SN2 reaction between NH3 and CH3NH3+ as a test system, we found that the results from ab initio VBSCF and from MOVB calculations are in good agreement, suggesting that the computationally efficient MOVB method is a reasonable model for VB simulations of condensed phase reactions. The results indicate that CDC and VDC diabatic states converge, respectively, to covalent and ionic states as the molecular geometries are distorted from the minimum of the respective diabatic state along the reaction coordinate. Furthermore, the resonance energy that stabilizes the energy of crossing between the two diabatic states, resulting in the transition state of the adiabatic ground-state reaction, has a strong dependence on the overlap integral between the two diabatic states and is a function of both the exchange integral and the total diabatic ground-state energy. PMID:18828577
Wan, Quan; Spanu, Leonardo; Galli, Giulia A; Gygi, François
2013-09-10
We report the first ab initio simulations of the Raman spectra of liquid water, obtained by combining first principles molecular dynamics and density functional perturbation theory. Our computed spectra are in good agreement with experiments, especially in the low frequency region. We also describe a systematic strategy to analyze the Raman intensities, which is of general applicability to molecular solids and liquids, and it is based on maximally localized Wannier functions and effective molecular polarizabilities. Our analysis revealed the presence of intermolecular charge fluctuations accompanying the hydrogen bond (HB) stretching modes at 270 cm(-1), in spite of the absence of any Raman activity in the isotropic spectrum. We also found that charge fluctuations partly contribute to the 200 cm(-1) peak in the anisotropic spectrum, thus providing insight into the controversial origin of such peak. Our results highlighted the importance of taking into account electronic effects in interpreting the Raman spectra of liquid water and the key role of charge fluctuations within the HB network; they also pointed at the inaccuracies of models using constant molecular polarizabilities to describe the Raman response of liquid water. PMID:26592405
Ab initio study of the one-dimensional H-bonded ferroelectric CsH2PO4
NASA Astrophysics Data System (ADS)
Lasave, J.; Abufager, P.; Koval, S.
2016-04-01
We studied the microscopic mechanism of the paraelectric-ferroelectric (PE-FE) phase transition of CsH2PO4 (CDP) by means of first-principles electronic structure calculations. The calculated structural parameters in the PE and FE phases as well as the total spontaneous polarization Ps obtained with the Berry phase formalism for CDP are in good agreement with experiments. The main contribution to Ps originates from a large y x component of the calculated Born effective-charge tensor for the disordered protons. Moreover, this component is ≈2.7 times larger than the z x component of the proton effective-charge tensor relevant to the polarization in the H-bonded FE KH2PO4 (KDP). This is the main feature that compensates the different number of protons per formula unit involved in the phase transitions for CDP and KDP to give close values for their measured Ps. Correlations among protons and heavy atoms along chains in the b direction lead to larger instabilities for the global and local FE distortions in CDP and its deuterated counterpart DCDP. We conclude that the tunneling particle is a dressed proton (deuteron) in agreement with experiments and with recent ab initio results for KDP.
An unconventional halogen bond with carbene as an electron donor: An ab initio study
NASA Astrophysics Data System (ADS)
Li, Qingzhong; Wang, Yilei; Liu, Zhenbo; Li, Wenzuo; Cheng, Jianbo; Gong, Baoan; Sun, Jiazhong
2009-02-01
An unconventional halogen bond has been proved to exist in H2C-BrH complex. The halogen bond energy of H2C-BrH complex is calculated at four levels of theory [MP2, MP4, CCSD, and CCSD(T)]. The result shows that the carbene is a better electron donor. The substitution effect is prominent in this interaction. For example, the interaction energy in H2C-BrCN complex is increased by more than 300% relative to H2C-BrH complex. The analyses of NBO, AIM, and energy components were used to unveil the nature of the interaction. The results show that this novel halogen bond has similar characteristics to hydrogen bonds.
Ab initio investigations of the electronic structure and chemical bonding of Li{sub 2}ZrN{sub 2}
Matar, S.F.; Poettgen, R.; Al Alam, A.F.; Ouaini, N.
2012-06-15
The electronic structure of the ternary nitride Li{sub 2}ZrN{sub 2} is examined from ab initio with DFT computations for an assessment of the properties of chemical bonding. The compound is found insulating with 1.8 eV band gap; it becomes metallic and less ionic upon removal of one equivalent of Li. The chemical interaction is found mainly between Zr and N on one hand and Li and N on the other hand. While all pair interactions are bonding, antibonding N-N interactions are found dominant at the top of the valence band of Li{sub 2}ZrN{sub 2} and they become less intense upon removal of Li. From energy differences the partial delithiation leading to Li{sub 2-x}ZrN{sub 2} (x={approx}1) is favored. - Graphical abstract: Trigonal structure of Li{sub 2}ZrN{sub 2} showing the Zr-N-Li layers along the c-axis. Highlights: Black-Right-Pointing-Pointer Li{sub 2}ZrN{sub 2} calculated insulating with a 1.8 eV gap in agreement with its light green color. Black-Right-Pointing-Pointer Lithium de-intercalation is energetically favored for one out of two Li equivalents. Black-Right-Pointing-Pointer Li plays little role in the change of the structure, ensured by Zr and N binding. Black-Right-Pointing-Pointer Similar changes in the electronic structure as for various intercalated phases of ZrN.
Kokalj, Anton
2015-01-01
The bonding of benzotriazole-an outstanding corrosion inhibitor for copper-on reduced and oxidized copper surfaces is discussed on the basis of density functional theory (DFT) calculations. Calculations reveal that benzotriazole is able to bond with oxide-free and oxidized copper surfaces and on both of them it bonds significantly stronger to coordinatively unsaturated Cu sites. This suggests that benzotriazole is able to passivate the reactive under-coordinated surface sites that are plausible microscopic sites for corrosion attack. Benzotriazole can adsorb in a variety of different forms, yet it forms a strong molecule-surface bond only in deprotonated form. The bonding is even stronger when the deprotonated form is incorporated into organometallic adcomplexes. This is consistent with existing experimental evidence that benzotriazole inhibits corrosion by forming protective organometallic complexes. It is further shown that adsorption of benzotriazole considerably reduces the metal work function, which is a consequence of a large permanent molecular dipole and a properly oriented adsorption structure. It is argued that such a pronounced effect on the work function might be relevant for corrosion inhibition, because it should diminish the anodic corrosion reaction, which is consistent with existing experimental evidence that benzotriazole, although a mixed type inhibitor, predominantly affects the anodic reaction. PMID:25955130
Intramolecular hydrogen bonds: ab initio Car Parrinello simulations of arylamide torsions
NASA Astrophysics Data System (ADS)
Doerksen, Robert J.; Chen, Bin; Klein, Michael L.
2003-10-01
Gas-phase, room temperature Car-Parrinello molecular dynamics simulations using the HCTH density functional are reported for the arylamides acetanilide ( 1) and ortho-methylthioacetanilide ( 2). The simulations show that in 1, rotation around the ring-amide bond is relatively unrestricted. By contrast, in 2 the methylthio side chain encourages the amide to be directed with N-H pointing toward S, not to flip by 360°, and furthermore to remain close to coplanar with the benzene ring. Because of an intramolecular N-H⋯S hydrogen bond, the N-H stretch frequency of 2 is red-shifted by ˜78 cm -1 compared to that of 1.
An eclipsed Csp3-CH3 bond? An ab initio investigation of an atypical rotation barrier
NASA Astrophysics Data System (ADS)
Vacek, George
In 1987, Seiler et al. (1987, Angew. Chem. Int. Edn Engl., 26, 1175) reported the observation of an eclipsed Csp3-CH3 bond in the low-temperature X-ray analysis of a crystalline trihydrate tricyclic orthoamide. This contrasts with the observed staggered conformation in an anhydrous crystal of the tricyclic orthoamide. Modest quantum chemical calculations estimated that a 5.5 kcal-mol-1 energy difference between the eclipsed form and the expected staggered form could be accounted for by unprecedented hydrogen bonding interaction strengths due to cooperative effects. In 1993, Novoa et al. (1993, Angew. Chem. Int. Edn Engl, 32, 588) showed that the earlier calculations were in error since they did not correct for basis set superposition error (BSSE). The BSSE corrected results of Novoa et al. could not account for the necessary 5.5 kcal mol-1 and this led them to question the original X-ray analysis of Seiler et al. Both sets of calculations suffered from the use of insufficiently accurate levels of theory and an inappropriate choice of model system. A current best estimate for the barrier to rotation is 1.4 kcal mol-1, far below that of the earlier studies. The stabilization of the eclipsed geometry due to hydrogen bonding after correction for BSSE is enough to match this small barrier. At today's best level of theory, the experimentally observed structure is 0.2 kcal mol-1 more stable than the staggered geometry.
Electronic structure and anisotropic chemical bonding in TiNF from ab initio study
Matar, Samir F.
2012-01-15
Accounting for disorder in anatase titanium nitride fluoride TiNF is done through atoms re-distributions based on geometry optimizations using ultra soft pseudo potentials within density functional theory DFT. The fully geometry relaxed structures are found to keep the body centering of anatase (I4{sub 1}/amd No. 141). The new structural setups are identified with space groups I-4m2 No. 119 and Imm2 No. 44 which obey the 'group to subgroup' relationships with respect to anatase. In the ground state Imm2 structure identified from energy differences, TiNF is found semi-conducting with similar density of states features to anatase TiO{sub 2} and a chemical bonding differentiated between covalent like Ti-N versus ionic like Ti-F. Inter-anion N-F bonding is also identified. - Graphical Abstract: The geometry optimized ground state anatase derived TiNF structure with arrangement of open faceted TiN3F3 distorted octahedra. The insert shows the arrangement of octahedra in anatase TiO{sub 2}. Highlights: Black-Right-Pointing-Pointer Original approach of TiNF structure for addressing the electronic band structure. Black-Right-Pointing-Pointer Based on anatase, two different ordering scheme models with geometry optimization. Black-Right-Pointing-Pointer New structures obeying the group{yields}subgroup relationships with Imm2 ground state from energy. Black-Right-Pointing-Pointer In the ground state TiNF is found semi-conducting with similar density of states to anatase TiO{sub 2}. Black-Right-Pointing-Pointer Chemical bonding differentiated between covalent like Ti-N and ionic Ti-F.
Tetrel bond cooperativity in open-chain (CH3CN)n and (CH3NC)n clusters (n = 2-7): An ab initio study
NASA Astrophysics Data System (ADS)
Esrafili, Mehdi D.; Mohammadirad, Nafiseh; Solimannejad, Mohammad
2015-05-01
In this work, ab initio calculations are performed to characterize cooperative effects between tetrel bond interactions in linear (CH3CN)2-7 and (CH3NC)2-7 clusters. It is found that increasing the size of the clusters brings about an important reduction in the tetrel bond distances, with the effect being most significant in the central units. To shed more light on the cooperative phenomenon, molecular electrostatic potential and electron density analyses are performed. According to interaction energy decomposition analysis, cooperative effects enhance the covalency of C⋯N and C⋯C interactions.
NASA Astrophysics Data System (ADS)
Oguri, Tomoya; Shimamura, Kohei; Shibuta, Yasushi; Shimojo, Fuyuki; Yamaguchi, Shu
2014-03-01
Dissociation of ethanol on a nickel cluster is investigated by ab initio molecular dynamics simulation to reveal the bond dissociation mechanism of carbon source molecules during carbon nanotube synthesis. C-C bonds in only CHxCO fragments are dissociated on the nickel cluster, whereas there is no preferential structure among the fragments for C-O bond dissociation. The dissociation preference is uncorrelated with the bond dissociation energy of corresponding bonds in freestanding molecules but is correlated with the energy difference between fragment molecules before and after dissociation on the nickel surface. Moreover, carbon-chain formation occurs after C-C bond dissociation in a continuous simulation. What determines the chirality of CNTs? What happens at the dissociation stage of carbon source molecules? Regarding the former question, many researchers have pointed out the good epitaxial relationship between a graphite network and a close-packed facet (i.e., fcc(1 1 1) or hcp(0 0 0 1)) of transition metals [17-19]. Therefore, the correlation between the chirality of CNTs and the angle of the step edge on metal (or metal carbide) surfaces has been closely investigated [20-22]. In association with this geometric matching, the epitaxial growth of graphene on Cu(1 1 1) and Ni(1 1 1) surfaces has recently been achieved via CCVD technique [23-25], which is a promising technique for the synthesis of large-area and monolayer graphene.Regarding the latter question, it is empirically known that the yield and quality of CNT products strongly depend on the choice of carbon source molecules and additives. For example, it is well known that the use of ethanol as carbon source molecules yields a large amount of SWNTs without amorphous carbons (called the alcohol CCVD (ACCVD) technique) compared with the CCVD process using hydrocarbons [4]. Moreover, the addition of a small amount of water dramatically enhances the activity and lifetime of the catalytic metal (called the
NASA Astrophysics Data System (ADS)
Petrongolo, Carlo; Ranghino, Graziella; Scordamaglia, Raimondo
1980-01-01
The potential energy surface of the β-lactam + OH - reaction, related to the mode of action of β-lactam antibiotics, was investigated using the ab initio Hartree—Fock method with the STO-3G basis set. Three possible reaction paths for the B A C2 breaking of the amidic CN bond were obtained and discussed. The minimum-energy reaction path is characterized by the following processes: (1) the formation of a tetrahedral intermediate, ≈ 121 kcal mol -1 more stable than the reagents; (2) a barrier, ≈ 15 kcal mol -1 above the intermediate, which is mainly due to the partial breaking of the amidic bond; (3) the complete breaking of the amidic bond concerted with a proton transfer till the formation of the final product, ≈ 34 kcal mol -1 more stable than the intermediate. The evolution of some molecular orbitals and of the electron population along the reaction path was also discussed.
Gopi, R; Ramanathan, N; Sundararajan, K
2014-07-24
The 1:1 hydrogen-bonded complex of fluoroform and hydrogen chloride was studied using matrix-isolation infrared spectroscopy and ab initio computations. Using B3LYP and MP2 levels of theory with 6-311++G(d,p) and aug-cc-pVDZ basis sets, the structures of the complexes and their energies were computed. For the 1:1 CHF3-HCl complexes, ab initio computations showed two minima, one cyclic and the other acyclic. The cyclic complex was found to have C-H · · · Cl and C-F · · · H interactions, where CHF3 and HCl sub-molecules act as proton donor and proton acceptor, respectively. The second minimum corresponded to an acyclic complex stabilized only by the C-F · · · H interaction, in which CHF3 is the proton acceptor. Experimentally, we could trap the 1:1 CHF3-HCl cyclic complex in an argon matrix, where a blue-shift in the C-H stretching mode of the CHF3 sub-molecule was observed. To understand the nature of the interactions, Atoms in Molecules and Natural Bond Orbital analyses were carried out to unravel the reasons for blue-shifting of the C-H stretching frequency in these complexes. PMID:24979667
Feng, Yong; Liu, Lei; Wang, Jin-Ti; Huang, Hao; Guo, Qing-Xiang
2003-01-01
Composite ab initio CBS-Q and G3 methods were used to calculate the bond dissociation energies (BDEs) of over 200 compounds listed in CRC Handbook of Chemistry and Physics (2002 ed.). It was found that these two methods agree with each other excellently in the calculation of BDEs, and they can predict BDEs within 10 kJ/mol of the experimental values. Using these two methods, it was found that among the examined compounds 161 experimental BDEs are valid because the standard deviation between the experimental and theoretical values for them is only 8.6 kJ/mol. Nevertheless, 40 BDEs listed in the Handbook may be highly inaccurate as the experimental and theoretical values for them differ by over 20 kJ/mol. Furthermore, 11 BDEs listed in the Handbook may be seriously flawed as the experimental and theoretical values for them differ by over 40 kJ/mol. Using the 161 cautiously validated experimental BDEs, we then assessed the performances of the standard density functional (DFT) methods including B3LYP, B3P86, B3PW91, and BH&HLYP in the calculation of BDEs. It was found that the BH&HLYP method performed poorly for the BDE calculations. B3LYP, B3P86, and B3PW91, however, performed reasonably well for the calculation of BDEs with standard deviations of about 12.1-18.0 kJ/mol. Nonetheless, all the DFT methods underestimated the BDEs by 4-17 kJ/mol in average. Sometimes, the underestimation by the DFT methods could be as high as 40-60 kJ/mol. Therefore, the DFT methods were more reliable for relative BDE calculations than for absolute BDE calculations. Finally, it was observed that the basis set effects on the BDEs calculated by the DFT methods were usually small except for the heteroatom-hydrogen BDEs. PMID:14632451
AB INITIO AND CALPHAD THERMODYNAMICS OF MATERIALS
Turchi, P A
2004-04-14
Ab initio electronic structure methods can supplement CALPHAD in two major ways for subsequent applications to stability in complex alloys. The first one is rather immediate and concerns the direct input of ab initio energetics in CALPHAD databases. The other way, more involved, is the assessment of ab initio thermodynamics {acute a} la CALPHAD. It will be shown how these results can be used within CALPHAD to predict the equilibrium properties of multi-component alloys.
NASA Astrophysics Data System (ADS)
Bandekar, Jagdeesh; Klima, Suzanne
1992-10-01
Ab initio quantum chemical computations were carried out on (a) dimethyl ether, (b) N-methyl formamide, (c) dimethyl ether- N-methyl formamide complex, and (d) N-methyl formamide dimer to compute the strengths of hydrogen bonds (H-bonds) between the NH groups and CO and ether COC groups. The basis set used was the 3-21G set of the GAUSSIAN 80 program obtained from QCPE, Bloomington, IN. Variations in the strengths of these two H-bonds with the N . . O distance (where O is either carbonyl or ether group oxygen) were studied and found to be similar in behavior. The strength of the "ether" hydrogen bond is computed to be 10.32 kcal mol -1, which is quite significant compared to the value of 10.11 kcal mol -1 for the more accepted "carbonyl" hydrogen bond. The "ether" hydrogen bond is found to be directional, specific and non-negligible. Work with two more basis sets has indicated that the results so obtained are not dependent on their choice. Possible importance of such a hydrogen bond in polyurethanes, inhalation anesthetics, and depsi-peptides is indicated.
Bankura, Arindam; Chandra, Amalendu
2015-01-28
The dynamics of proton transfer (PT) through hydrogen bonds in a two-dimensional water layer confined between two graphene sheets at room temperature are investigated through ab initio and quantum-classical simulations. The excess proton is found to be mostly solvated as an Eigen cation where the hydronium ion donates three hydrogen bonds to the neighboring water molecules. In the solvation shell of the hydronium ion, the three coordinated water molecules with two donor hydrogen bonds are found to be properly presolvated to accept a proton. Although no hydrogen bond needs to be broken for transfer of a proton to such presolvated water molecules from the hydronium ion, the PT rate is still found to be not as fast as it is for one-dimensional chains. Here, the PT is slowed down as the probability of finding a water with two donor hydrogen bonds in the solvation shell of the hydronium ion is found to be only 25%-30%. The hydroxide ion is found to be solvated mainly as a complex anion where it accepts four H-bonds through its oxygen atom and the hydrogen atom of the hydroxide ion remains free all the time. Here, the presolvation of the hydroxide ion to accept a proton requires that one of its hydrogen bonds is broken and the proton comes from a neighboring water molecule with two acceptor and one donor hydrogen bonds. The coordination number reduction by breaking of a hydrogen bond is a slow process, and also the population of water molecules with two acceptor and one donor hydrogen bonds is only 20%-25% of the total number of water molecules. All these factors together tend to slow down the hydroxide ion migration rate in two-dimensional water layers compared to that in three-dimensional bulk water.
NASA Astrophysics Data System (ADS)
Matar, S. F.; Al Alam, A. F.; Gédéon, D.; Ouaini, N.
2013-11-01
Potential hydrogen storage ternaries Zr3FeH7 and Zr2FeH5, are studied from ab initio with the purpose of identifying changes in electronic structures and bonding properties. Cohesive energy trends: Ecoh. (ZrH2) > Ecoh. (Zr2FeH5) > Ecoh. (Zr3FeH7) > Ecoh. (hypothetic-FeH) indicate a progressive destabilization of the binary hydride ZrH2 through adjoined Fe so that Zr3FeH7 is found less cohesive than Zr2FeH5. From the energy volume equations of states EOS the volume increase upon hydriding the intermetallics leads to higher bulk moduli B0 explained by the Zr/Fe-H bonding. Fe-H bond in Zr2FeH5 leads to annihilate magnetic polarization on Fe whereas Fe magnetic moment develops in Zr3FeH7 identified as ferromagnetic in the ground state. These differences in magnetic behaviors are due to the weakly ferromagnetic Fe largely affected by lattice environment, as opposed to strongly ferromagnetic Co. Hydrogenation of the binary intermetallics weakens the inter-metal bonding and favors the metal-hydrogen bonds leading to more cohesive hydrides as with respect to the pristine binaries. Charge analyses point to covalent like Fe versus ionic Zr and hydrogen charges ranging from covalent H-0.27 to more ionic H-0.5.
NASA Astrophysics Data System (ADS)
Nakamura, Nobuo; Masui, Hirotsugo; Ueda, Takahiro
2000-02-01
Ab initio Hartree-Fock molecular orbital calculations were applied to the crystalline imidazole and its derivatives in order to examine systematically the effect of possible N-H---N type hydrogen bond-ing on the nuclear quadrupole interaction parameters in these materials. The nitrogen quadrupole coupling constant (QCC) and the asymmetry parameter (η) of the electric field gradient (EFG) were found to depend strongly on the size of the molecular clusters, from single molecule, to dimer, trimer and to the infinite molecular chain, i.e., crystalline state, implying that the intermolecular N-H -N hydrogen bond affects significantly the electronic structure of imidazole molecule. A certain correla-tion between the QCC of 14N and the N-H bond distance R was also found and interpreted on the basis of the molecular orbital theory. However, we found that the value of the calculated EFG at the hy-drogen position of the N-H group, or the corresponding QCC value of 2 H, increases drastically as R-3 when R is shorter than about 0.1 nm, due probably to the inapplicability of the Gaussian basis sets to the very short chemical bond as revealed in the actual imidazole derivatives. We suggested that the ob-served N-H distances in imidazole derivatives should be re-examined.
NASA Astrophysics Data System (ADS)
Bock, Charles W.; Trachtman, Mendel; George, Philip
1981-11-01
The CO bond length and the quadratic, cubic and quartic stretching force constants, calculated ab initio using the unscaled 4-31G basis set with full geometry optimization, are reported for three series of monosubstituted carbonyl compounds in which the atom directly bonded to the carbonyl carbon is another carbon, a nitrogen, or an oxygen atom, respectively. The data are analyzed in terms of the In ƒ versus In re relationship, and also the generalized power functions and exponential functions proposed by Herschbach and Laurie. Not only does the atom directly bonded to the carbonyl carbon affect the magnitude of re and the force constants, but the rest of the substituent group is found to be capable of exerting an even greater influence. Within each series of compounds the overall progression from the shortest to the longest CO bonds is tentatively attributed to a diminishing electron density in the bonding region.
Infrared spectroscopy and ab initio study of hydrogen bonded Cl3CD·N(CH3)3 complex in the gas phase
NASA Astrophysics Data System (ADS)
Rutkowski, K. S.; Melikova, S. M.; Linok, O. V.; Czarnik-Matusewicz, B.; Rospenk, M.
2015-02-01
FTIR spectra of the gas phase Cl3CD + TMA mixture have been studied at room temperature in ∼800-4000 cm-1 frequency domain. The formation of the H-bonded Cl3CD…TMA complex has been detected. Spectroscopic parameters of the band ascribed to the complex were evaluated. MP2 frozen core ab initio calculations have been carried out with the Pople-type 6-311++G(d,p) basis set. The equilibrium geometries and harmonic vibrational frequencies of the complex were obtained using CP-corrected gradient techniques. The ''freq = anharm'' option has been tested for Cl3CD monomer and Cl3CD…TMA complex to examine possible anharmonic effects on the vibrations localized on the proton donor. The effects of Darling-Dennison and Fermi resonances on the frequency of the stretching vibration of the CH proton donor were analyzed.
Ab initio calculations on the structure and nature of the hydrogen bonded complex H2S ṡṡṡ HF
NASA Astrophysics Data System (ADS)
Singh, U. Chandra; Kollman, Peter A.
1984-01-01
Ab initio calculations employing an extended basis set and correlation energy estimates at the MP2 (second order Möller-Plesset) level have been used to estimate the binding energy, minimum energy S ...F distance, and minimum energy bisector angle θ (angle between the twofold axis of H2S and the S ...F line) for the hydrogen bonded complex H2S ... HF. The calculated distance and bisector angle θ are in reasonable agreement with experiment; the calculated binding energy can be used to provide a good estimate of the (as yet unknown) experimental value. Morokuma component analyses and simple electrostatic (molecular mechanics) calculations have been carried out on the complex as a function of bisector angle and they demonstrated that, despite suggestions to the contrary, the H2S ... HF bisector angle is predominantly determined by the electrostatic energy.
NASA Astrophysics Data System (ADS)
Yamashita, Takefumi; Takatsuka, Kazuo
2007-02-01
The infrared spectrum of phenol-water cationic cluster, [PhOH•H2O]+, taken by Sawamura et al. [J. Phys. Chem. 100, 8131 (1996)] is puzzling in that the peak due to the stretching mode of the phenolic OH (3657cm-1 for a neutral monomer and 3524cm-1 for PhOH•H2O) seemingly disappears and instead an extremely broad tail extending down to 2900cm-1 is observed. The present authors theoretically ascribe this anomalous spectrum to an inhomogeneous broadening of the OH stretching peak caused by the hydrogen bond, the strength of which has been greatly enhanced by ionization of the phenyl ring. Indeed they estimate that the peak position is at 2300cm-1 and the spectral width can become as wide as 1000cm-1 at the cluster energy of 32kcal/mol. This surprisingly wide broadening can be generic in hydrogen-bond systems, which in turn is useful to study the nature of the hydrogen-bond assisted dynamics in various systems such as those in DNA and proteins. To study the present system quantitatively, the authors have developed an ab initio mixed quantum-classical method, in which the nuclear motions on an adiabatic ab initio potential surface are treated such that only the OH stretching motion is described quantum mechanically, while all the other remaining modes are treated classically with on-the-fly scheme. This method includes the implementation of many numerical methodologies, which enables it to deal with a relatively large molecular system. With this theoretical method, the authors analyze the present anomalous broadening in a great detail. In particular, they suggest that one can extract direct information about the hydrogen-bond dynamics with respect to the clear correlation between the vibrational excitation energy of the OH stretching and intermolecular distance by means of a time-resolved infrared spectroscopy: Reflecting the slow and wide-range variation of the intermolecular distance of the relevant hydrogen bond, the time-resolved spectrum is predicted to vary
Woo, Hin-koon; Wang, Xue B.; Wang, Lai S.; Lau, Kai Chung
2005-12-01
The strength of the low-barrier hydrogen bond in hydrogen maleate in the gas phase was investigated by low-temperature photoelectron spectroscopy and ab initio calculations. Photoelectron spectra of maleic and fumaric acid monoanions (cis-/trans-HO2CCHdCHCO2 -) were obtained at low temperatures and at 193 nm photon energy. Vibrational structure was observed for trans-HO2CCHdCHCO2 - due to the OCO bending modes; however, cis-HO2CCHdCHCO2 - yielded a broad and featureless spectrum. The electron binding energy of cis-HO2CCHdCHCO2 - is about 1 eV blue-shifted relative to trans-HO2CCHdCHCO2 - due to the formation of intramolecular hydrogen bond in the cis-isomer. Theoretical calculations (CCSD(T)/ aug-cc-pVTZ and B3LYP/aug-cc-pVTZ) were carried out to estimate the strength of the intramolecular hydrogen bond in cis-HO2CCHdCHCO2 -. Combining experimental and theoretical calculations yields an estimate of 21.5 ( 2.0 kcal/mol for the intramolecular hydrogen bond strength in hydrogen maleate.
Collective rotation from ab initio theory
NASA Astrophysics Data System (ADS)
Caprio, M. A.; Maris, P.; Vary, J. P.; Smith, R.
2015-08-01
Through ab initio approaches in nuclear theory, we may now seek to quantitatively understand the wealth of nuclear collective phenomena starting from the underlying internucleon interactions. No-core configuration interaction (NCCI) calculations for p-shell nuclei give rise to rotational bands, as evidenced by rotational patterns for excitation energies, electromagnetic moments and electromagnetic transitions. In this review, NCCI calculations of 7-9Be are used to illustrate and explore ab initio rotational structure, and the resulting predictions for rotational band properties are compared with experiment. We highlight the robustness of ab initio rotational predictions across different choices for the internucleon interaction.
Ab initio nuclear structure theory
NASA Astrophysics Data System (ADS)
Negoita, Gianina Alina
Ab initio no core methods have become major tools for understanding the properties of light nuclei based on realistic nucleon-nucleon (NN) and three-nucleon (NNN) interactions. A brief description is provided for the inter-nucleon interactions that fit two-body scattering and bound state data, as well as NNN interactions. Major new progress, including the goal of applying these interactions to solve for properties of nuclei, is limited by convergence issues. That is, with the goal of obtaining high precision solutions of the nuclear many-body Hamiltonian with no core methods (all nucleons treated on the same footing), one needs to proceed to very large basis spaces to achieve a convergence pattern suitable for extrapolation to the exact result. This thesis investigates (1) the similarity renormalization group (SRG) approach to soften the interaction, while preserving its phase shift properties, and (2) adoption of a realistic basis space using Woods-Saxon (WS) single-particle wavefunctions. Both have their advantages and limitations, discussed here. For (1), SRG was demonstrated by applying it to a realistic NN interaction, JISP16, in a harmonic oscillator (HO) representation. The degree of interaction softening achieved through a regulator parameter is examined. For (2), new results are obtained with the realistic JISP16 NN interaction in ab initio calculations of light nuclei 4He, 6He and 12C, using a WS basis optimized to minimize the ground-state energy within the truncated no core shell model. These are numerically-intensive many-body calculations. Finally, to gain insight into the potential for no core investigations of heavier nuclei, an initial investigation was obtained for the odd mass A = 47 - 49 region nuclei straddling 48Ca. The motivation for selecting these nuclei stems from the aim of preparing for nuclear double beta-decay studies of 48Ca. In these heavier systems, phenomenological additions to the realistic NN interaction determined by previous
Piazza, Zachary A.; Li, Wei-Li; Wang, Lai-Sheng E-mail: lai-sheng-wang@brown.edu; Popov, Ivan A.; Boldyrev, Alexander I. E-mail: lai-sheng-wang@brown.edu; Pal, Rhitankar; Cheng Zeng, Xiao
2014-07-21
Photoelectron spectroscopy and ab initio calculations are used to investigate the structures and chemical bonding of the B{sub 25}{sup −} cluster. Global minimum searches reveal a dense potential energy landscape with 13 quasi-planar structures within 10 kcal/mol at the CCSD(T)/6-311+G(d) level of theory. Three quasi-planar isomers (I, II, and III) are lowest in energy and nearly degenerate at the CCSD(T) level of theory, with II and III being 0.8 and 0.9 kcal/mol higher, respectively, whereas at two density functional levels of theory isomer III is the lowest in energy (8.4 kcal/mol more stable than I at PBE0/6-311+G(2df) level). Comparison with experimental photoelectron spectroscopic data shows isomer II to be the major contributor while isomers I and III cannot be ruled out as minor contributors to the observed spectrum. Theoretical analyses reveal similar chemical bonding in I and II, both involving peripheral 2c-2e B−B σ-bonding and delocalized interior σ- and π-bonding. Isomer III has an interesting elongated ribbon-like structure with a π-bonding pattern analogous to those of dibenzopentalene. The high density of low-lying isomers indicates the complexity of the medium-sized boron clusters; the method dependency of predicting relative energies of the low-lying structures for B{sub 25}{sup −} suggests the importance of comparison with experiment in determining the global minima of boron clusters at this size range. The appearance of many low-lying quasi-planar structures containing a hexagonal hole in B{sub 25}{sup −} suggests the importance of this structural feature in maintaining planarity of larger boron clusters.
Ab initio phonon limited transport
NASA Astrophysics Data System (ADS)
Verstraete, Matthieu
We revisit the thermoelectric (TE) transport properties of two champion materials, PbTe and SnSe, using fully first principles methods. In both cases the performance of the material is due to subtle combinations of structural effects, scattering, and phase space reduction. In PbTe anharmonic effects are completely opposite to the predicted quasiharmonic evolution of phonon frequencies and to frequently (and incorrectly) cited extrapolations of experiments. This stabilizes the material at high T, but also tends to enhance its thermal conductivity, in a non linear manner, above 600 Kelvin. This explains why PbTe is in practice limited to room temperature applications. SnSe has recently been shown to be the most efficient TE material in bulk form. This is mainly due to a strongly enhanced carrier concentration and electrical conductivity, after going through a phase transition from 600 to 800 K. We calculate the transport coefficients as well as the defect concentrations ab initio, showing excellent agreement with experiment, and elucidating the origin of the double phase transition as well as the new charge carriers. AH Romero, EKU Gross, MJ Verstraete, and O Hellman PRB 91, 214310 (2015) O. Hellman, IA Abrikosov, and SI Simak, PRB 84 180301 (2011)
Paesani, Francesco; Xantheas, Sotiris S.; Voth, Gregory A.
2009-10-01
A molecular-level description of the unique properties of hydrogen-bond networks is critical for understanding many fundamental physico-chemical processes in aqueous environments. In this article a novel simulation approach, combining an ab-initio based force field for water with a quantum treatment of the nuclear motion, is applied to investigate hydrogen-bond dynamics in liquid water with a specific focus on the relationship of these dynamics to vibrational spectroscopy. Linear and nonlinear infrared (IR) spectra are calculated for liquid water, HOD in D2O and HOD in H2O and discussed in the context of the results obtained using other approaches that have been employed in studies of water dynamics. A comparison between the calculated spectra and the available experimental data yields an overall good agreement, indicating the accuracy of the present simulation approach in describing the properties of liquid water at ambient conditions. Possible improvements on the representation of the underlying water interactions as well as the treatment of the molecular motion at the quantum-mechanical level are also discussed. This research was supported by the Division of Chemical Sciences, Biosciences and Geosciences, US Department of Energy. Battelle operates the Pacific Northwest National Laboratory for the US Department of Energy.
Zabidi, Noriza Ahmad; Kassim, Hasan Abu; Shrivastava, Keshav N.
2008-05-20
Polonium is the only element with a simple cubic (sc) crystal structure. Atoms in solid polonium sit at the corners of a simple cubic unit cell and no where else. Polonium has a valence electron configuration 6s{sup 2}6p{sup 4} (Z = 84). The low temperature {alpha}-phase transforms into the rhombohedral (trigonal) {beta} structure at {approx}348 K. The sc {alpha}-Po unit cell constant is a = 3.345 A. The beta form of polonium ({beta}-Po) has the lattice parameters, a{sub R} = 3.359 A and a rhombohedral angle 98 deg. 13'. We have performed an ab initio electronic structure calculation by using the density functional theory. We have performed the calculation with and without spin-orbit (SO) coupling by using both the LDA and the GGA for the exchange-correlations. The k-points in a simple cubic BZ are determined by R (0.5, 0.5, 0.5), {gamma} (0, 0, 0), X (0.5, 0, 0), M (0.5, 0.5, 0) and {gamma} (0, 0, 0). Other directions of k-points are {gamma} (0, 0, 0), X (0.5, 0, 0), R (0.5, 0.5, 0.5) and {gamma} (0, 0, 0). The SO splittings of p states at the {gamma} point in the GGA+SO scheme for {alpha}-Po are 0.04 eV and 0.02 eV while for the {beta}-Po these are 0.03 eV and 0.97 eV. We have also calculated the vibrational spectra for the unit cells in both the structures. We find that exchanging of a Po atom by Pb atom produces several more bands and destabilizes the {beta} phase.
ERIC Educational Resources Information Center
Halpern, Arthur M.; Glendening, Eric D.
2013-01-01
A project for students in an upper-level course in quantum or computational chemistry is described in which they are introduced to the concepts and applications of a high quality, ab initio treatment of the ground-state potential energy curve (PEC) for H[subscript 2] and D[subscript 2]. Using a commercial computational chemistry application and a…
NASA Astrophysics Data System (ADS)
Klein, Erik; Lukeš, Vladimír
2006-11-01
In this paper, the study of phenol and 37 compounds representing various ortho-, para-, and meta-substituted phenols is presented. Molecules and their radical structures were studied using ab initio methods with inclusion of correlation energy and DFT in order to calculate the O-H bond dissociation enthalpies (BDEs) and vertical ionization potentials (IPs). Calculated BDEs and IPs were compared with available experimental values to ascertain the suitability of used methods, especially for the description of the substituent induced changes in BDE and IP. MP2, MP3, and MP4 methods do not give reliable results, since they significantly underestimate substituent induced changes in BDE and do not reflect distinct effect of substituents related to para and meta position correctly. DFT/B3LYP method reflects the effect of substituents on BDE satisfactorily, though ΔBDEs are in narrower range than experimental values. BDE of phenol was calculated also using CCSD(T) method in various basis sets. Both, DFT and HF methods describe the effect of substituents on IP identically. However, DFT considerably underestimates individual values. HF method gives IPs in very good agreement with experimental data. Obtained results show that dependences of BDEs and IPs on Hammett constants of the substituents are linear. Linearity of DFT BDE vs. IP dependence is even better than the dependences on Hammett constants and obtained equations allow estimating of O-H BDEs of meta- and para-substituted phenols from calculated IPs.
Laitinen, R.S.; Pakkanen, T.A.; Steudel, R.
1987-02-04
Ab initio MO calculations involving the 4-31G* basis set have been used to predict the equilibrium geometries of the hypervalent sulfur hydrides H/sub 2/SS, (HS)/sub 2/SS, H/sub 2/S(SH)/sub 2/, H/sub 2/S(SSH)/sub 2/, and the cyclic H/sub 4/S/sub 4/. The energy changes in their formation from appropriate sulfanes H/sub 2/S/sub n/ (n = 1-4) have been studied with the 6-31G* basis set including the correction for the electron correlation by the second- and third-order Moeller-Plesset perturbation theory. The results are used to discuss the possible pathways in the interconversion reactions between various sulfur compounds containing cumulated SS bonds, for example, the formation of S/sub 7/ from S/sub 8/ for which hypervalent intermediates have been proposed recently. Comparison with experimental evidence is made whenever possible.
NASA Astrophysics Data System (ADS)
Binoy, J.; Prathima, N. B.; Murali Krishna, C.; Santhosh, C.; Hubert Joe, I.; Jayakumar, V. S.
2006-08-01
Acetanilide, a compound of pharmaceutical importance possessing pain-relieving properties due to its blocking the pulse dissipating along the nerve fiber, is subjected to vibrational spectral investigation using NIR FT Raman, FT-IR, and SERS. The geometry, Mulliken charges, and vibrational spectrum of acetanilide have been computed using the Hartree-Fock theory and density functional theory employing the 6-31G (d) basis set. To investigate the influence of intermolecular amide hydrogen bonding, the geometry, charge distribution, and vibrational spectrum of the acetanilide dimer have been computed at the HF/6-31G (d) level. The computed geometries reveal that the acetanilide molecule is planar, while twisting of the secondary amide group with respect to the phenyl ring is found upon hydrogen bonding. The trans isomerism and “amido” form of the secondary amide, hyperconjugation of the C=O group with the adjacent C-C bond, and donor-acceptor interaction have been investigated using computed geometry. The carbonyl stretching band position is found to be influenced by the tendency of the phenyl ring to withdraw nitrogen lone pair, intermolecular hydrogen bonding, conjugation, and hyperconjugation. A decrease in the NH and C=O bond orders and increase in the C-N bond orders due to donor-acceptor interaction can be observed in the vibrational spectra. The SERS spectral analysis reveals that the flat orientation of the molecule on the adsorption plane is preferred.
Phonocatalysis. An ab initio simulation experiment
NASA Astrophysics Data System (ADS)
Kim, Kwangnam; Kaviany, Massoud
2016-06-01
Using simulations, we postulate and show that heterocatalysis on large-bandgap semiconductors can be controlled by substrate phonons, i.e., phonocatalysis. With ab initio calculations, including molecular dynamic simulations, the chemisorbed dissociation of XeF6 on h-BN surface leads to formation of XeF4 and two surface F/h-BN bonds. The reaction pathway and energies are evaluated, and the sorption and reaction emitted/absorbed phonons are identified through spectral analysis of the surface atomic motion. Due to large bandgap, the atomic vibration (phonon) energy transfer channels dominate and among them is the match between the F/h-BN covalent bond stretching and the optical phonons. We show that the chemisorbed dissociation (the pathway activation ascent) requires absorption of large-energy optical phonons. Then using progressively heavier isotopes of B and N atoms, we show that limiting these high-energy optical phonons inhibits the chemisorbed dissociation, i.e., controllable phonocatalysis.
NASA Astrophysics Data System (ADS)
George, Philip; Bock, Charles W.; Trachtman, Mendel
Using the 4-31G basis set augmented with polarization functions on the S-atom all geometrical parameters have been optimized for the planar chain and ring conformers (related by 180° rotation of the OH, SH or NH group) of acrylic acid, the three monothioperformic acids,α-hydroxyacrolein, formyl and thioformylhydroxylamine, glyoxalmonoimine, formimidol and glyoxalmonooxime. The changes in bond lengths and bond angles which accompany the conversion of the chain into the ring conformer are classified under ten headings for a total of eight four-membered, nine five-membered and three six-membered ring systems. Systematic trends are found, indicative of electron transfer, which support the hypothesis that there is a hydrogen-bonding type of interaction in the four-as well as in the five- and six-membered rings. The change in total molecular energy for the conversion reaction is divided up into distortion and bonding energy components, and each partitioned in terms of the expectation energy differencnes ΔEK, ΔVee, ΔVnn and ΔVen. The nature of the distortion and bonding steps, whether attractive-dominant or repulsive-dominant, is correlated with ring size and the presence (or absence) of heavy atoms external to the ring. The internuclear disease OH, O…O and H … O in the OH … O(=C) hydrogen bridge are compared with experimental values for intermolecular hydrogen bonds, and the present findings are discussed in relation to studies in the literature on electron density distributions and chages in Mulliken gross electronic populations.
NASA Astrophysics Data System (ADS)
Bock, Charles W.; Trachtman, Mendel; George, Philip
1982-12-01
The lengths of the terminal double bonds and the central CC single bond, and the cross-coupling constant involving the double bonds, calculated ab initio using the unscaled 4-31G basis set, are reported for 15 planar conformers of conjugated dienes containing the OCC'C″ grouping, 9 planar conformers containing the OCC'N grouping, and 6 planar conformers containing the OCC'O'grouping. The larger values for the CC bond length tend to be associated with the smaller values for the double-bond lengths, and vice versa. A linear relationship holds fairly well between rcc and the sum of the double-bond lengths. The sign of the cross-coupling constant is the same for each type of conjugated diene, negative for foc.c'c″ and foc.c'N' but positive for foc.c'o'. The values of foc.c'c″ and f oc.c'N are larger the smaller rcc' in accord with a linear relationship. For all three types of diene the values for pairs of cis and trans conformers do not differ very much, the ratio fcis/ ftrans being a little greater than 1 in ten cases and a little less than 1 in six cases. This lack of any marked difference between cis and trans conformers suggests that the coupling between the double bonds occurs in the main via the bonded framework of the molecule, and that "through-space" interactions, which would be a special feature of cis conformers, are of relatively little importance.
Den, Takuya S.; Frey, Hans-Martin; Leutwyler, Samuel
2014-11-21
The gas-phase rotational motion of hexafluorobenzene has been measured in real time using femtosecond (fs) time-resolved rotational Raman coherence spectroscopy (RR-RCS) at T = 100 and 295 K. This four-wave mixing method allows to probe the rotation of non-polar gas-phase molecules with fs time resolution over times up to ∼5 ns. The ground state rotational constant of hexafluorobenzene is determined as B{sub 0} = 1029.740(28) MHz (2σ uncertainty) from RR-RCS transients measured in a pulsed seeded supersonic jet, where essentially only the v = 0 state is populated. Using this B{sub 0} value, RR-RCS measurements in a room temperature gas cell give the rotational constants B{sub v} of the five lowest-lying thermally populated vibrationally excited states ν{sub 7/8}, ν{sub 9}, ν{sub 11/12}, ν{sub 13}, and ν{sub 14/15}. Their B{sub v} constants differ from B{sub 0} by between −1.02 MHz and +2.23 MHz. Combining the B{sub 0} with the results of all-electron coupled-cluster CCSD(T) calculations of Demaison et al. [Mol. Phys. 111, 1539 (2013)] and of our own allow to determine the C-C and C-F semi-experimental equilibrium bond lengths r{sub e}(C-C) = 1.3866(3) Å and r{sub e}(C-F) = 1.3244(4) Å. These agree with the CCSD(T)/wCVQZ r{sub e} bond lengths calculated by Demaison et al. within ±0.0005 Å. We also calculate the semi-experimental thermally averaged bond lengths r{sub g}(C-C)=1.3907(3) Å and r{sub g}(C-F)=1.3250(4) Å. These are at least ten times more accurate than two sets of experimental gas-phase electron diffraction r{sub g} bond lengths measured in the 1960s.
NASA Astrophysics Data System (ADS)
Ben Doudou, Bessem; Chen, Jun; Vivet, Alexandre; Poilâne, Christophe
2016-03-01
In this paper, we have investigated the chemical bond interactions between covalently functionalized zigzag (5,0) and (8,0) SWCNT-SWCNT via various covalent linkages. Side-to-side junctions connected via amide, ester and anhydride linkages were particularly studied. The geometries and energy of the forming reaction were investigated using first-principles density functional theory. Furthermore, the band structures and the total density of states (DOS) of the junctions have also been analyzed. Our results show that several promising structures could be obtained by using chemical connection strategy and particularly the junctions formed by coupling amino functionalized SWCNT and carboxylic acid functionalized SWCNT was more favorable.
Li, Junjie; Li, Xiaohu; Iyengar, Srinivasan S
2014-06-10
We discuss a multiconfigurational treatment of the "on-the-fly" electronic structure within the quantum wavepacket ab initio molecular dynamics (QWAIMD) method for coupled treatment of quantum nuclear effects with electronic structural effects. Here, multiple single-particle electronic density matrices are simultaneously propagated with a quantum nuclear wavepacket and other classical nuclear degrees of freedom. The multiple density matrices are coupled through a nonorthogonal configuration interaction (NOCI) procedure to construct the instantaneous potential surface. An adaptive-mesh-guided set of basis functions composed of Gaussian primitives are used to simplify the electronic structure calculations. Specifically, with the replacement of the atom-centered basis functions positioned on the centers of the quantum-mechanically treated nuclei by a mesh-guided band of basis functions, the two-electron integrals used to compute the electronic structure potential surface become independent of the quantum nuclear variable and hence reusable along the entire Cartesian grid representing the quantum nuclear coordinates. This reduces the computational complexity involved in obtaining a potential surface and facilitates the interpretation of the individual density matrices as representative diabatic states. The parametric nuclear position dependence of the diabatic states is evaluated at the initial time-step using a Shannon-entropy-based sampling function that depends on an approximation to the quantum nuclear wavepacket and the potential surface. This development is meant as a precursor to an on-the-fly fully multireference electronic structure procedure embedded, on-the-fly, within a quantum nuclear dynamics formalism. We benchmark the current development by computing structural, dynamic, and spectroscopic features for a series of bihalide hydrogen-bonded systems: FHF(-), ClHCl(-), BrHBr(-), and BrHCl(-). We find that the donor-acceptor structural features are in good
Multiple time step integrators in ab initio molecular dynamics
Luehr, Nathan; Martínez, Todd J.; Markland, Thomas E.
2014-02-28
Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties associated with splitting the ab initio potential into fast and slowly varying components. Here we present two schemes that enable efficient time-scale separation in ab initio calculations: one based on fragment decomposition and the other on range separation of the Coulomb operator in the electronic Hamiltonian. We demonstrate for both water clusters and a solvated hydroxide ion that multiple time-scale molecular dynamics allows for outer time steps of 2.5 fs, which are as large as those obtained when such schemes are applied to empirical potentials, while still allowing for bonds to be broken and reformed throughout the dynamics. This permits computational speedups of up to 4.4x, compared to standard Born-Oppenheimer ab initio molecular dynamics with a 0.5 fs time step, while maintaining the same energy conservation and accuracy.
Ab Initio Studies of Calcium Carbonate Hydration.
Lopez-Berganza, Josue A; Diao, Yijue; Pamidighantam, Sudhakar; Espinosa-Marzal, Rosa M
2015-11-25
Ab initio simulations of large hydrated calcium carbonate clusters are challenging due to the existence of multiple local energy minima. Extensive conformational searches around hydrated calcium carbonate clusters (CaCO3·nH2O for n = 1-18) were performed to find low-energy hydration structures using an efficient combination of Monte Carlo searches, density-functional tight binding (DFTB+) method, and density-functional theory (DFT) at the B3LYP level, or Møller-Plesset perturbation theory at the MP2 level. This multilevel optimization yields several low-energy structures for hydrated calcium carbonate. Structural and energetics analysis of the hydration of these clusters revealed a first hydration shell composed of 12 water molecules. Bond-length and charge densities were also determined for different cluster sizes. The solvation of calcium carbonate in bulk water was investigated by placing the explicitly solvated CaCO3·nH2O clusters in a polarizable continuum model (PCM). The findings of this study provide new insights into the energetics and structure of hydrated calcium carbonate and contribute to the understanding of mechanisms where calcium carbonate formation or dissolution is of relevance. PMID:26505205
Thermochemical data for CVD modeling from ab initio calculations
Ho, P.; Melius, C.F.
1993-12-31
Ab initio electronic-structure calculations are combined with empirical bond-additivity corrections to yield thermochemical properties of gas-phase molecules. A self-consistent set of heats of formation for molecules in the Si-H, Si-H-Cl, Si-H-F, Si-N-H and Si-N-H-F systems is presented, along with preliminary values for some Si-O-C-H species.
NASA Technical Reports Server (NTRS)
Lee, Timothy J.
1989-01-01
HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.
Ab initio infrared and Raman spectra
NASA Astrophysics Data System (ADS)
Fredkin, Donald R.; Komornicki, Andrew; White, Steven R.; Wilson, Kent R.
1983-06-01
We discuss several ways in which molecular absorption and scattering spectra can be computed ab initio, from the fundamental constants of nature. These methods can be divided into two general categories. In the first, or sequential, type of approach, one first solves the electronic part of the Schrödinger equation in the Born-Oppenheimer approximation, mapping out the potential energy, dipole moment vector (for infrared absorption) and polarizability tensor (for Raman scattering) as functions of nuclear coordinates. Having completed the electronic part of the calculation, one then solves the nuclear part of the problem either classically or quantum mechanically. As an example of the sequential ab initio approach, the infrared and Raman rotational and vibrational-rotational spectral band contours for the water molecule are computed in the simplest rigid rotor, normal mode approximation. Quantum techniques are used to calculate the necessary potential energy, dipole moment, and polarizability information at the equilibrium geometry. A new quick, accurate, and easy to program classical technique involving no reference to Euler angles or special functions is developed to compute the infrared and Raman band contours for any rigid rotor, including asymmetric tops. A second, or simultaneous, type of ab initio approach is suggested for large systems, particularly those for which normal mode analysis is inappropriate, such as liquids, clusters, or floppy molecules. Then the curse of dimensionality prevents mapping out in advance the complete potential, dipole moment, and polarizability functions over the whole space of nuclear positions of all atoms, and a solution in which the electronic and nuclear parts of the Born-Oppenheimer approximation are simultaneously solved is needed. A quantum force classical trajectory (QFCT) molecular dynamic method, based on linear response theory, is described, in which the forces, dipole moment, and polarizability are computed quantum
Surface Segregation Energies of BCC Binaries from Ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2003-01-01
We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy method. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameterization. Quantum approximate segregation energies are computed with and without atomistic relaxation. The ab initio calculations are performed without relaxation for the most part, but predicted relaxations from quantum approximate calculations are used in selected cases to compute approximate relaxed ab initio segregation energies. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with other quantum approximate and ab initio theoretical work, and available experimental results.
Collective rotation from ab initio theory
NASA Astrophysics Data System (ADS)
Caprio, Mark A.; Maris, Pieter; Vary, James P.
2015-10-01
The challenge of ab initio nuclear theory is to quantitatively predict the complex and highly-correlated behavior of the nuclear many-body system, starting from the underlying internucleon interactions. We may now seek to understand the wealth of nuclear collective phenomena through ab initio approaches. No-core configuration interaction (NCCI) calculations for p-shell nuclei give rise to rotational bands, as evidenced by rotational patterns for excitation energies, electromagnetic moments, and electromagnetic transitions. In this talk, the intrinsic structure of these bands is discussed, and the predicted rotational bands are compared to experiment. Supported by the US DOE under Award Nos. DE-FG02-95ER-40934, DESC0008485 (SciDAC/NUCLEI), and DE-FG02-87ER40371 and the US NSF under Award No. 0904782. Computational resources provided by NERSC (US DOE Contract No. DE-AC02-05CH11231).
Ab initio infrared and Raman spectra
NASA Technical Reports Server (NTRS)
Fredkin, D. R.; White, S. R.; Wilson, K. R.; Komornicki, A.
1983-01-01
It is pointed out that with increased computer power and improved computational techniques, such as the gradients developed in recent years, it is becoming practical to compute spectra ab initio, from the fundamental constants of nature, for systems of increasing complexity. The present investigation has the objective to explore several possible ab initio approaches to spectra, giving particular attention to infrared and nonresonance Raman. Two approaches are discussed. The sequential approach, in which first the electronic part and then later the nuclear part of the Born-Oppenheimer approximation is solved, is appropriate for small systems. The simultaneous approach, in which the electronic and nuclear parts are solved at the same time, is more appropriate for many-atom systems. A review of the newer quantum gradient techniques is provided, and the infrared and Raman spectral band contours for the water molecule are computed.
NASA Astrophysics Data System (ADS)
Thiessen, P. A.; Treder, H.-J.
Jedes initium wird durch experimenta crucis zum eventus. Jedes theoretisch interpretierbare ex-eventu-Resultat führt auf ein neues Initium. Gerade dies ist die gemeinsame Aussage von Atomistik, Quantenmechanik und Relativitätstheorie.Translated AbstractAb initio vel ex eventu. IIEvery initium becomes an eventus by experimenta crucis. Every theoretically interpretable ex-eventu result leads to a new initium. Right this is the joint assertion of atomism, quantum mechanics, and relativity.
Ab initio Bogoliubov coupled cluster theory
NASA Astrophysics Data System (ADS)
Signoracci, Angelo; Hagen, Gaute; Duguet, Thomas
2014-09-01
Coupled cluster (CC) theory has become a standard method in nuclear theory for realistic ab initio calculations of medium mass nuclei, but remains limited by its requirement of a Slater determinant reference state which reasonably approximates the nuclear system of interest. Extensions of the method, such as equation-of-motion CC, permit the calculation of nuclei with one or two nucleons added or removed from a doubly magic core, yet still only a few dozen nuclei are accessible with modern computational restrictions. In order to extend the applicability of ab initio methods to open-shell systems, the superfluid nature of nuclei must be taken into account. By utilizing Bogoliubov algebra and employing spontaneous symmetry breaking with respect to particle number conservation, superfluid systems can be treated by a single reference state. An ab initio theory to include correlations on top of a Bogoliubov reference state has been developed in the guise of standard CC theory. The formalism and first results of this Bogoliubov coupled cluster theory will be presented to demonstrate the applicability of the method.
M Kondo; T Mates; D Fischer; F Wudl; E Kramer
2011-12-31
Interfaces between phenylacetylene (PA) monolayers and two silicon surfaces, Si(111) and Si(100), are probed by X-ray photoelectron spectroscopy (XPS), ultraviolet photoelectron spectroscopy (UPS), and near-edge X-ray absorption fine structure (NEXAFS) spectroscopy, and the results are analyzed using ab initio molecular orbital calculations. The monolayer systems are prepared via the surface hydrosilylation reaction between PA and hydrogen-terminated silicon surfaces. The following spectral features are obtained for both of the PA-Si(111) and PA-Si(100) systems: a broad {pi}-{pi}* shakeup peak at 292 eV (XPS), a broad first ionization peak at 3.8 eV (UPS), and a low-energy C 1s {yields} {pi}* resonance peak at 284.3 eV (NEXAFS). These findings are ascribed to a styrene-like {pi}-conjugated molecular structure at the PA-Si interface by comparing the experimental data with theoretical analysis results. A conclusion is drawn that the vinyl group can keep its {pi}-conjugation character on the hydrogen-terminated Si(100) [H:Si(100)] surface composed of the dihydride (SiH{sub 2}) groups as well as on hydrogen-terminated Si(111) having the monohydride (SiH) group. The formation mechanism of the PA-Si(100) interface is investigated within cluster ab initio calculations, and the possible structure of the H:Si(100) surface is discussed based on available data.
Heats of Segregation of BCC Binaries from ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2004-01-01
We compare dilute-limit heats of segregation for selected BCC transition metal binaries computed using ab initio and quantum approximate energy methods. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent LMTO-based parameters. Quantum approximate segregation energies are computed with and without atomistic relaxation, while the ab initio calculations are performed without relaxation. Results are discussed within the context of a segregation model driven by strain and bond-breaking effects. We compare our results with full-potential quantum calculations and with available experimental results.
Ab initio no core full configuration approach for light nuclei
NASA Astrophysics Data System (ADS)
Kim, Youngman; Shin, Ik Jae; Maris, Pieter; Vary, James P.; Forssén, Christian; Rotureau, Jimmy
2014-07-01
Comprehensive understanding of the structure and reactions of light nuclei poses theoretical and computational challenges. Still, a number of ab initio approaches have been developed to calculate the properties of atomic nuclei using fundamental interactions among nucleons. Among them, we work with the ab initio no core full configuration (NCFC) method and ab initio no core Gamow Shell Model (GSM). We first review these approaches and present some recent results.
Ab initio no core full configuration approach for light nuclei
NASA Astrophysics Data System (ADS)
Kim, Youngman; Shin, Ik Jae; Maris, Pieter; Vary, James P.; Forssén, Christian; Rotureau, Jimmy
2015-10-01
Comprehensive understanding of the structure and reactions of light nuclei poses theoretical and computational challenges. Still, a number of ab initio approaches have been developed to calculate the properties of atomic nuclei using fundamental interactions among nucleons. Among them, we work with the ab initio no core full configuration (NCFC) method and ab initio no core Gamow Shell Model (GSM). We first review these approaches and present some recent results.
Ab initio melting curve of osmium
NASA Astrophysics Data System (ADS)
Burakovsky, L.; Burakovsky, N.; Preston, D. L.
2015-11-01
The melting curve of osmium up to a pressure P of 500 GPa is obtained from an extensive suite of ab initio quantum molecular dynamics (QMD) simulations using the Z method. The ab initio P =0 melting point of Os is 3370 ±75 K; this range encompasses all of the available data in the literature and corroborates the conclusion of J. W. Arblaster [Platinum Metals Rev. 49, 166 (2005)], 10.1595/147106705X70264 that the melting temperature of pure Os is 3400 ±50 K and that the 3300 K typically quoted in the literature is the melting point of impure Os. The T =0 equation of state (EOS) of Os and the P dependence of the optimized c /a ratio for the hexagonal unit cell, both to pressures ˜900 GPa, are obtained in the ab initio approach as validation of its use. Although excellent agreement with the available experimental data (P ≲80 GPa) is found, it is the third-order Birch-Murnaghan EOS with B0'=5 rather than the more widely accepted B0'=4 that describes the QMD data to higher pressures, in agreement with the more recent experimental EOS by Godwal et al. The theoretical melting curve of Os obtained earlier by Joshi et al. is shown to be inconsistent with our QMD results, and the possible reason for this discrepancy is suggested. Regularities in the melting curves of Os and five other third-row transition metals (Ta, W, Re, Pt, Au) could be used to estimate the currently unknown melting curves of Hf and Ir.
Ab initio quantum chemistry: Methodology and applications
Friesner, Richard A.
2005-01-01
This Perspective provides an overview of state-of-the-art ab initio quantum chemical methodology and applications. The methods that are discussed include coupled cluster theory, localized second-order Moller–Plesset perturbation theory, multireference perturbation approaches, and density functional theory. The accuracy of each approach for key chemical properties is summarized, and the computational performance is analyzed, emphasizing significant advances in algorithms and implementation over the past decade. Incorporation of a condensed-phase environment by means of mixed quantum mechanical/molecular mechanics or self-consistent reaction field techniques, is presented. A wide range of illustrative applications, focusing on materials science and biology, are discussed briefly. PMID:15870212
An improved ab initio structure for fluorine peroxide (FOOF)
NASA Astrophysics Data System (ADS)
Mack, Hans-Georg; Oberhammer, Heinz
1988-03-01
Ab initio calculations with the 6-31G* and Dunning (9s5p/4s2p) basis sets augmented with p and d functions at various levels of theory (RHF, MP2, MP3, and MP4) were carried out on F 2O 2. The best result was obtained at the MP2 level with the Dunning basis plus one set of d functions on fluorine and two sets of d functions on oxygen. These calculations reproduce the experimental bond lengths to within 0.01 Å and the angles to within the experimental uncertainties.
Ab initio non-relativistic spin dynamics
Ding, Feizhi; Goings, Joshua J.; Li, Xiaosong; Frisch, Michael J.
2014-12-07
Many magnetic materials do not conform to the (anti-)ferromagnetic paradigm where all electronic spins are aligned to a global magnetization axis. Unfortunately, most electronic structure methods cannot describe such materials with noncollinear electron spin on account of formally requiring spin alignment. To overcome this limitation, it is necessary to generalize electronic structure methods and allow each electron spin to rotate freely. Here, we report the development of an ab initio time-dependent non-relativistic two-component spinor (TDN2C), which is a generalization of the time-dependent Hartree-Fock equations. Propagating the TDN2C equations in the time domain allows for the first-principles description of spin dynamics. A numerical tool based on the Hirshfeld partitioning scheme is developed to analyze the time-dependent spin magnetization. In this work, we also introduce the coupling between electron spin and a homogenous magnetic field into the TDN2C framework to simulate the response of the electronic spin degrees of freedom to an external magnetic field. This is illustrated for several model systems, including the spin-frustrated Li{sub 3} molecule. Exact agreement is found between numerical and analytic results for Larmor precession of hydrogen and lithium atoms. The TDN2C method paves the way for the ab initio description of molecular spin transport and spintronics in the time domain.
Ab initio non-relativistic spin dynamics
NASA Astrophysics Data System (ADS)
Ding, Feizhi; Goings, Joshua J.; Frisch, Michael J.; Li, Xiaosong
2014-12-01
Many magnetic materials do not conform to the (anti-)ferromagnetic paradigm where all electronic spins are aligned to a global magnetization axis. Unfortunately, most electronic structure methods cannot describe such materials with noncollinear electron spin on account of formally requiring spin alignment. To overcome this limitation, it is necessary to generalize electronic structure methods and allow each electron spin to rotate freely. Here, we report the development of an ab initio time-dependent non-relativistic two-component spinor (TDN2C), which is a generalization of the time-dependent Hartree-Fock equations. Propagating the TDN2C equations in the time domain allows for the first-principles description of spin dynamics. A numerical tool based on the Hirshfeld partitioning scheme is developed to analyze the time-dependent spin magnetization. In this work, we also introduce the coupling between electron spin and a homogenous magnetic field into the TDN2C framework to simulate the response of the electronic spin degrees of freedom to an external magnetic field. This is illustrated for several model systems, including the spin-frustrated Li3 molecule. Exact agreement is found between numerical and analytic results for Larmor precession of hydrogen and lithium atoms. The TDN2C method paves the way for the ab initio description of molecular spin transport and spintronics in the time domain.
Ab Initio Neutron Drops with Chiral Hamiltonians
NASA Astrophysics Data System (ADS)
Potter, Hugh; Maris, Pieter; Vary, James
2015-04-01
Ab initio calculations for neutron drops are of interest for insights into neutron-rich nuclei and neutron star matter, and for examining the neutron-only sector of nucleon-nucleon and 3-nucleon interactions. I present ab initio results calculated using the no-core shell model with 2- and 3-body chiral Hamiltonians for neutron drops up to 20 neutrons confined in a 10 MeV harmonic trap. I discuss ground state energies, internal energies, radii, and evidence for pairing. In addition, excitation energies can be used to investigate the spin-orbit splittings in the p-shell and sd -shell. Prior Green's Function Monte Carlo calculations using the Argonne v8' potential with added 3-nucleon forces serve as a comparison. Supported by DOE Grants DESC0008485 (SciDAC/NUCLEI), DE-FG02-87ER40371, and NSF Grant 0904782; computational resources provided by the Oak Ridge Leadership Computing Facility (DOE Office of Science Contract DE-AC05-00OR22725) under an INCITE award.
Towards Accurate Ab Initio Predictions of the Spectrum of Methane
NASA Technical Reports Server (NTRS)
Schwenke, David W.; Kwak, Dochan (Technical Monitor)
2001-01-01
We have carried out extensive ab initio calculations of the electronic structure of methane, and these results are used to compute vibrational energy levels. We include basis set extrapolations, core-valence correlation, relativistic effects, and Born- Oppenheimer breakdown terms in our calculations. Our ab initio predictions of the lowest lying levels are superb.
Guiding ab initio calculations by alchemical derivatives.
to Baben, M; Achenbach, J O; von Lilienfeld, O A
2016-03-14
We assess the concept of alchemical transformations for predicting how a further and not-tested change in composition would change materials properties. This might help to guide ab initio calculations through multidimensional property-composition spaces. Equilibrium volumes, bulk moduli, and relative lattice stability of fcc and bcc 4d transition metals Zr, Nb, Mo, Tc, Ru, Rh, Pd, and Ag are calculated using density functional theory. Alchemical derivatives predict qualitative trends in lattice stability while equilibrium volumes and bulk moduli are predicted with less than 9% and 28% deviation, respectively. Predicted changes in equilibrium volume and bulk moduli for binary and ternary mixtures of Rh-Pd-Ag are in qualitative agreement even for predicted bulk modulus changes as large as +100% or -50%. Based on these results, it is suggested that alchemical transformations could be meaningful for enhanced sampling in the context of virtual high-throughput materials screening projects. PMID:26979677
Ab initio determination of light hadron masses.
Dürr, S; Fodor, Z; Frison, J; Hoelbling, C; Hoffmann, R; Katz, S D; Krieg, S; Kurth, T; Lellouch, L; Lippert, T; Szabo, K K; Vulvert, G
2008-11-21
More than 99% of the mass of the visible universe is made up of protons and neutrons. Both particles are much heavier than their quark and gluon constituents, and the Standard Model of particle physics should explain this difference. We present a full ab initio calculation of the masses of protons, neutrons, and other light hadrons, using lattice quantum chromodynamics. Pion masses down to 190 mega-electron volts are used to extrapolate to the physical point, with lattice sizes of approximately four times the inverse pion mass. Three lattice spacings are used for a continuum extrapolation. Our results completely agree with experimental observations and represent a quantitative confirmation of this aspect of the Standard Model with fully controlled uncertainties. PMID:19023076
Discovering chemistry with an ab initio nanoreactor
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S.; Martínez, Todd J.
2014-11-02
Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provides detailed physical insight. While theoretical and computational studies have generally focused on specific processes or mechanistic hypotheses, recent methodological and computational advances harken the advent of their principal role in discovery. Here we report the development and application of the ab initio nanoreactor – a highly accelerated, first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor we show new pathways for glycine synthesis frommore » primitive compounds proposed to exist on the early Earth, providing new insight into the classic Urey-Miller experiment. Ultimately, these results highlight the emergence of theoretical and computational chemistry as a tool for discovery in addition to its traditional role of interpreting experimental findings.« less
Ab initio molar volumes and Gaussian radii.
Parsons, Drew F; Ninham, Barry W
2009-02-12
Ab initio molar volumes are calculated and used to derive radii for ions and neutral molecules using a spatially diffuse model of the electron distribution with Gaussian spread. The Gaussian radii obtained can be used for computation of nonelectrostatic ion-ion dispersion forces that underlie Hofmeister specific ion effects. Equivalent hard-sphere radii are also derived, and these are in reasonable agreement with crystalline ionic radii. The Born electrostatic self-energy is derived for a Gaussian model of the electronic charge distribution. It is shown that the ionic volumes used in electrostatic calculations of strongly hydrated cosmotropic ions ought best to include the first hydration shell. Ionic volumes for weakly hydrated chaotropic metal cations should exclude electron overlap (in electrostatic calculations). Spherical radii are calculated as well as nonisotropic ellipsoidal radii for nonspherical ions, via their nonisotropic static polarizability tensors. PMID:19140766
Discovering chemistry with an ab initio nanoreactor
NASA Astrophysics Data System (ADS)
Martinez, Todd
Traditional approaches for modeling chemical reaction networks such as those involved in combustion have focused on identifying individual reactions and using theoretical approaches to explore the underlying mechanisms. Recent advances involving graphical processing units (GPUs), commodity products developed for the videogaming industry, have made it possible to consider a distinct approach wherein one attempts to discover chemical reactions and mechanisms. We provide a brief summary of these developments and then discuss the concept behind the ``ab initio nanoreactor'' which explores the space of possible chemical reactions and molecular species for a given stoichiometry. The nanoreactor concept is exemplified with an example to the Urey-Miller reaction network which has been previously advanced as a potential model for prebiotic chemistry. We briefly discuss some of the future directions envisioned for the development of this nanoreactor concept.
Discovering chemistry with an ab initio nanoreactor
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S.; Martínez, Todd J.
2014-01-01
Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provides detailed physical insight. While theoretical and computational studies have generally focused on specific processes or mechanistic hypotheses, recent methodological and computational advances harken the advent of their principal role in discovery. Here we report the development and application of the ab initio nanoreactor – a highly accelerated, first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, providing new insight into the classic Urey-Miller experiment. These results highlight the emergence of theoretical and computational chemistry as a tool for discovery in addition to its traditional role of interpreting experimental findings. PMID:25411881
Guiding ab initio calculations by alchemical derivatives
NASA Astrophysics Data System (ADS)
to Baben, M.; Achenbach, J. O.; von Lilienfeld, O. A.
2016-03-01
We assess the concept of alchemical transformations for predicting how a further and not-tested change in composition would change materials properties. This might help to guide ab initio calculations through multidimensional property-composition spaces. Equilibrium volumes, bulk moduli, and relative lattice stability of fcc and bcc 4d transition metals Zr, Nb, Mo, Tc, Ru, Rh, Pd, and Ag are calculated using density functional theory. Alchemical derivatives predict qualitative trends in lattice stability while equilibrium volumes and bulk moduli are predicted with less than 9% and 28% deviation, respectively. Predicted changes in equilibrium volume and bulk moduli for binary and ternary mixtures of Rh-Pd-Ag are in qualitative agreement even for predicted bulk modulus changes as large as +100% or -50%. Based on these results, it is suggested that alchemical transformations could be meaningful for enhanced sampling in the context of virtual high-throughput materials screening projects.
Discovering chemistry with an ab initio nanoreactor.
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S; Martínez, Todd J
2014-12-01
Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provide detailed physical insight. Although theoretical and computational studies have generally focused on specific processes or mechanistic hypotheses, recent methodological and computational advances harken the advent of their principal role in discovery. Here we report the development and application of the ab initio nanoreactor--a highly accelerated first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor, we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, which provide new insight into the classic Urey-Miller experiment. These results highlight the emergence of theoretical and computational chemistry as a tool for discovery, in addition to its traditional role of interpreting experimental findings. PMID:25411881
Discovering chemistry with an ab initio nanoreactor
NASA Astrophysics Data System (ADS)
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S.; Martínez, Todd J.
2014-12-01
Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provide detailed physical insight. Although theoretical and computational studies have generally focused on specific processes or mechanistic hypotheses, recent methodological and computational advances harken the advent of their principal role in discovery. Here we report the development and application of the ab initio nanoreactor—a highly accelerated first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor, we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, which provide new insight into the classic Urey-Miller experiment. These results highlight the emergence of theoretical and computational chemistry as a tool for discovery, in addition to its traditional role of interpreting experimental findings.
Ab Initio Modeling of Molecular Radiation
NASA Technical Reports Server (NTRS)
Jaffe, Richard; Schwenke, David
2014-01-01
Radiative emission from excited states of atoms and molecules can comprise a significant fraction of the total heat flux experienced by spacecraft during atmospheric entry at hypersonic speeds. For spacecraft with ablating heat shields, some of this radiative flux can be absorbed by molecular constituents in the boundary layer that are formed by the ablation process. Ab initio quantum mechanical calculations are carried out to predict the strengths of these emission and absorption processes. This talk will describe the methods used in these calculations using, as examples, the 4th positive emission bands of CO and the 1g+ 1u+ absorption in C3. The results of these calculations are being used as input to NASA radiation modeling codes like NeqAir, HARA and HyperRad.
Discovering chemistry with an ab initio nanoreactor
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S.; Martínez, Todd J.
2014-11-02
Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provides detailed physical insight. While theoretical and computational studies have generally focused on specific processes or mechanistic hypotheses, recent methodological and computational advances harken the advent of their principal role in discovery. Here we report the development and application of the ab initio nanoreactor – a highly accelerated, first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, providing new insight into the classic Urey-Miller experiment. Ultimately, these results highlight the emergence of theoretical and computational chemistry as a tool for discovery in addition to its traditional role of interpreting experimental findings.
Ab Initio Calculation of the Hoyle State
Epelbaum, Evgeny; Krebs, Hermann; Lee, Dean; Meissner, Ulf-G.
2011-05-13
The Hoyle state plays a crucial role in the helium burning of stars heavier than our Sun and in the production of carbon and other elements necessary for life. This excited state of the carbon-12 nucleus was postulated by Hoyle as a necessary ingredient for the fusion of three alpha particles to produce carbon at stellar temperatures. Although the Hoyle state was seen experimentally more than a half century ago nuclear theorists have not yet uncovered the nature of this state from first principles. In this Letter we report the first ab initio calculation of the low-lying states of carbon-12 using supercomputer lattice simulations and a theoretical framework known as effective field theory. In addition to the ground state and excited spin-2 state, we find a resonance at -85(3) MeV with all of the properties of the Hoyle state and in agreement with the experimentally observed energy.
Ab initio alpha-alpha scattering
NASA Astrophysics Data System (ADS)
Elhatisari, Serdar; Lee, Dean; Rupak, Gautam; Epelbaum, Evgeny; Krebs, Hermann; Lähde, Timo A.; Luu, Thomas; Meißner, Ulf-G.
2015-12-01
Processes such as the scattering of alpha particles (4He), the triple-alpha reaction, and alpha capture play a major role in stellar nucleosynthesis. In particular, alpha capture on carbon determines the ratio of carbon to oxygen during helium burning, and affects subsequent carbon, neon, oxygen, and silicon burning stages. It also substantially affects models of thermonuclear type Ia supernovae, owing to carbon detonation in accreting carbon-oxygen white-dwarf stars. In these reactions, the accurate calculation of the elastic scattering of alpha particles and alpha-like nuclei—nuclei with even and equal numbers of protons and neutrons—is important for understanding background and resonant scattering contributions. First-principles calculations of processes involving alpha particles and alpha-like nuclei have so far been impractical, owing to the exponential growth of the number of computational operations with the number of particles. Here we describe an ab initio calculation of alpha-alpha scattering that uses lattice Monte Carlo simulations. We use lattice effective field theory to describe the low-energy interactions of protons and neutrons, and apply a technique called the ‘adiabatic projection method’ to reduce the eight-body system to a two-cluster system. We take advantage of the computational efficiency and the more favourable scaling with system size of auxiliary-field Monte Carlo simulations to compute an ab initio effective Hamiltonian for the two clusters. We find promising agreement between lattice results and experimental phase shifts for s-wave and d-wave scattering. The approximately quadratic scaling of computational operations with particle number suggests that it should be possible to compute alpha scattering and capture on carbon and oxygen in the near future. The methods described here can be applied to ultracold atomic few-body systems as well as to hadronic systems using lattice quantum chromodynamics to describe the interactions of
Ab initio alpha-alpha scattering.
Elhatisari, Serdar; Lee, Dean; Rupak, Gautam; Epelbaum, Evgeny; Krebs, Hermann; Lähde, Timo A; Luu, Thomas; Meißner, Ulf-G
2015-12-01
Processes such as the scattering of alpha particles ((4)He), the triple-alpha reaction, and alpha capture play a major role in stellar nucleosynthesis. In particular, alpha capture on carbon determines the ratio of carbon to oxygen during helium burning, and affects subsequent carbon, neon, oxygen, and silicon burning stages. It also substantially affects models of thermonuclear type Ia supernovae, owing to carbon detonation in accreting carbon-oxygen white-dwarf stars. In these reactions, the accurate calculation of the elastic scattering of alpha particles and alpha-like nuclei--nuclei with even and equal numbers of protons and neutrons--is important for understanding background and resonant scattering contributions. First-principles calculations of processes involving alpha particles and alpha-like nuclei have so far been impractical, owing to the exponential growth of the number of computational operations with the number of particles. Here we describe an ab initio calculation of alpha-alpha scattering that uses lattice Monte Carlo simulations. We use lattice effective field theory to describe the low-energy interactions of protons and neutrons, and apply a technique called the 'adiabatic projection method' to reduce the eight-body system to a two-cluster system. We take advantage of the computational efficiency and the more favourable scaling with system size of auxiliary-field Monte Carlo simulations to compute an ab initio effective Hamiltonian for the two clusters. We find promising agreement between lattice results and experimental phase shifts for s-wave and d-wave scattering. The approximately quadratic scaling of computational operations with particle number suggests that it should be possible to compute alpha scattering and capture on carbon and oxygen in the near future. The methods described here can be applied to ultracold atomic few-body systems as well as to hadronic systems using lattice quantum chromodynamics to describe the interactions of
Choudhuri, Jyoti Roy; Chandra, Amalendu
2014-11-21
We have presented a first principles simulation study of the structural and dynamical properties of a liquid-vapor interfacial system of a concentrated (5.3 M) aqueous NaCl solution. We have used ab initio molecular dynamics to examine the structural and dynamical properties of the bulk and interfacial regions. The structural aspects of the system that have been considered here include the inhomogeneous density profiles of ions and water molecules, hydrogen bond distributions, orientational profiles, and also vibrational frequency distributions in the bulk and interfacial regions. It is found that the sodium ions are mostly located in the interior, while the chloride anions occupy a significant portion of the interface of the slab. The water dipoles at the interface prefer to orient parallel to the surface. The dynamical aspects of the interfaces are investigated in terms of diffusion, orientational relaxation, hydrogen bond dynamics, and vibrational spectral diffusion. The results of the interfacial dynamics are compared with those of the corresponding bulk region. It is observed that the interfacial molecules exhibit faster diffusion and orientational relaxation with respect to the bulk. However, the interfacial molecules are found to have longer hydrogen bond lifetimes than those of the bulk. We have also investigated the correlations of hydrogen bond relaxation with the vibrational frequency fluctuations of interfacial water molecules. PMID:25416903
Sakaki, Shigeyoshi; Musashi, Yasuo
1995-03-29
Ab initio MO/MP4, SD-CI, and CCD models were invoked in the consideration of carbon dioxide and ethene insertion into Cu(II)-R bonds. For R = H, CH{sub 3}, and OH, the authors calculated the activation energy and reaction enthalpy for these reactions and presents the related arguments here.
Ab initio derivation of model energy density functionals
NASA Astrophysics Data System (ADS)
Dobaczewski, Jacek
2016-08-01
I propose a simple and manageable method that allows for deriving coupling constants of model energy density functionals (EDFs) directly from ab initio calculations performed for finite fermion systems. A proof-of-principle application allows for linking properties of finite nuclei, determined by using the nuclear nonlocal Gogny functional, to the coupling constants of the quasilocal Skyrme functional. The method does not rely on properties of infinite fermion systems but on the ab initio calculations in finite systems. It also allows for quantifying merits of different model EDFs in describing the ab initio results.
Ab Initio: And a New Era of Airline Pilot Training.
ERIC Educational Resources Information Center
Gesell, Laurence E.
1995-01-01
Expansion of air transportation and decreasing numbers seeking pilot training point to a shortage of qualified pilots. Ab initio training, in which candidates with no flight time are trained to air transport proficiency, could resolve the problem. (SK)
Ab initio two-component Ehrenfest dynamics
NASA Astrophysics Data System (ADS)
Ding, Feizhi; Goings, Joshua J.; Liu, Hongbin; Lingerfelt, David B.; Li, Xiaosong
2015-09-01
We present an ab initio two-component Ehrenfest-based mixed quantum/classical molecular dynamics method to describe the effect of nuclear motion on the electron spin dynamics (and vice versa) in molecular systems. The two-component time-dependent non-collinear density functional theory is used for the propagation of spin-polarized electrons while the nuclei are treated classically. We use a three-time-step algorithm for the numerical integration of the coupled equations of motion, namely, the velocity Verlet for nuclear motion, the nuclear-position-dependent midpoint Fock update, and the modified midpoint and unitary transformation method for electronic propagation. As a test case, the method is applied to the dissociation of H2 and O2. In contrast to conventional Ehrenfest dynamics, this two-component approach provides a first principles description of the dynamics of non-collinear (e.g., spin-frustrated) magnetic materials, as well as the proper description of spin-state crossover, spin-rotation, and spin-flip dynamics by relaxing the constraint on spin configuration. This method also holds potential for applications to spin transport in molecular or even nanoscale magnetic devices.
Ab initio two-component Ehrenfest dynamics
Ding, Feizhi; Goings, Joshua J.; Liu, Hongbin; Lingerfelt, David B.; Li, Xiaosong
2015-09-21
We present an ab initio two-component Ehrenfest-based mixed quantum/classical molecular dynamics method to describe the effect of nuclear motion on the electron spin dynamics (and vice versa) in molecular systems. The two-component time-dependent non-collinear density functional theory is used for the propagation of spin-polarized electrons while the nuclei are treated classically. We use a three-time-step algorithm for the numerical integration of the coupled equations of motion, namely, the velocity Verlet for nuclear motion, the nuclear-position-dependent midpoint Fock update, and the modified midpoint and unitary transformation method for electronic propagation. As a test case, the method is applied to the dissociation of H{sub 2} and O{sub 2}. In contrast to conventional Ehrenfest dynamics, this two-component approach provides a first principles description of the dynamics of non-collinear (e.g., spin-frustrated) magnetic materials, as well as the proper description of spin-state crossover, spin-rotation, and spin-flip dynamics by relaxing the constraint on spin configuration. This method also holds potential for applications to spin transport in molecular or even nanoscale magnetic devices.
Vibrational and ab initio molecular dynamics studies of bradykinin
NASA Astrophysics Data System (ADS)
Święch, Dominika; Kubisiak, Piotr; Andrzejak, Marcin; Borowski, Piotr; Proniewicz, Edyta
2016-07-01
In this study, the comprehensive theoretical and experimental investigations of Raman (RS) and infrared absorption (IR) spectra of bradykinin (BK) are presented. The ab initio Born-Oppenheimer molecular dynamics (BOMD) calculations, in the presence of water molecules that form the first coordination sphere, were used for conformational analysis of the BK structure. Based on the Density Functional Theory (DFT) calculations at the B3LYP/6-311G(d) level the vibrational spectra were interpreted. The calculated frequencies were scaled by means of the effective scaling frequency factor (ESFF) method. The theoretical data, which confirm the compact structure of BK in the presence of the water molecules revealed the remarkable effect of the intermolecular hydrogen bonding on the BK structural properties.
Ab initio H2O in realistic hydrophilic confinement.
Allolio, Christoph; Klameth, Felix; Vogel, Michael; Sebastiani, Daniel
2014-12-15
A protocol for the ab initio construction of a realistic cylindrical pore in amorphous silica, serving as a geometric nanoscale confinement for liquids and solutions, is presented. Upon filling the pore with liquid water at different densities, the structure and dynamics of the liquid inside the confinement can be characterized. At high density, the pore introduces long-range oscillations into the water density profile, which makes the water structure unlike that of the bulk across the entire pore. The tetrahedral structure of water is also affected up to the second solvation shell of the pore wall. Furthermore, the effects of the confinement on hydrogen bonding and diffusion, resulting in a weakening and distortion of the water structure at the pore walls and a slowdown in diffusion, are characterized. PMID:25208765
Ab initio theory of NMR chemical shifts in solids
Louie, S.G. |
1997-12-31
A new formalism for ab initio calculation of the orbital magnetic susceptibility and the NMR chemical shifts in solids and liquids is presented. The approach can be applied to periodic systems such as crystals, surfaces or polymers, and with a supercell technique, to nonperiodic systems such as amorphous materials, liquids, or solids with defects. The formalism is based on the density functional theory in the local density approximation and makes use of a generalized f-sum rule to eliminate the divergent terms that plagued previous theories. Calculations have been successfully carried out for the diamagnetic susceptibility of a number of insulators and for the NMR chemical shifts of a variety of systems including free molecules, ionic crystals, hydrogen-bonded materials and amorphous carbon.
Ab initio potential energy surface and rovibrational states of HBO
NASA Astrophysics Data System (ADS)
Ha, Tae-Kyu; Makarewicz, Jan
1999-01-01
The potential energy surface describing the large-amplitude motion of H around the BO core in the HBO molecule has been determined from ab initio calculations. This surface has been sampled by a set of 170 grid points from a two-dimensional space defined by the stretching and the bending coordinates of the H nucleus. At each grid point, the BO bond length has been optimized using the second-order Møller-Plesset perturbation theory with the basis set aug-cc-pVTZ. The surface has a local minimum for the linear as well as the bent configuration of HBO. A low energy barrier to the linear configuration BOH causes a large-amplitude motion and a strong rovibrational interaction in the molecule. Its rovibrational dynamics is different from the dynamics in bent or quasilinear triatomic molecules.
Heats of Segregation of BCC Binaries from Ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2003-01-01
We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy methods. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameters. Quantum approximate segregation energies are computed with and without atomistic relaxation. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with full-potential quantum calculations and with available experimental results.
NASA Astrophysics Data System (ADS)
George, Philip; Glusker, Jenny P.; Markham, George D.; Trachtman, Mendel; Bock, Charles W.
It is common knowledge that some metal ions prefer to bond to nitrogen atoms, others prefer oxygen and others select sulfur, although the mechanistic origin of this behaviour is not well understood. To provide quantitative data that can illuminate this characteristic difference, we have been performing ab initio molecular orbital calculations on complexes of a wide variety of main group and transition metal ions with simple ligands. In this paper we concentrate on metal ion complexes of NH3 and compare them with the corresponding complexes with H2O (Trachtman et al., 1998, Inorg. Chem., 37, 4221). The results reported here show that for each metal ion the bonding to ammonia in monoammines is intrinsically stronger than that to water in monohydrates, but that the enthalpy of formation of the amine complexes has a greater sensitivity to the nature of the metal ion. For both mono- and divalent transition metal cations the ligand-field stabilization energy (LFSE) in the monoammine complexes is larger than that in the monohydrates, although the larger enthalpies of formation, , for the monoammines are due only in part to the larger LFSE values. The formation of both monoammines and monohydrates is accompanied by a transfer of charge from the ligand hydrogen atoms to the metal ion so that there is relatively little change in the net charge on the oxygen or nitrogen atoms. Hence the increased acidity of the ligand in the metal ion complex is not the result of net electron depletion of the atom that is directly bonded to the metal ion, but rather reflects weakening of the bond of that ligand atom to its proton (O-H or N-H). This characteristic is used by many enzymes to promote catalytic activity.
Ab Initio and Phenomenological Modeling of the Phonon Spectrum of Superhard cp-BC2N
NASA Astrophysics Data System (ADS)
Basalaev, Yu. M.; Kopytov, A. V.; Pavlova, T. Yu.; Poplavnoi, A. S.
2015-11-01
The phonon spectrum of hypothetical superhard cp-BC2N is calculated based on ab initio method of density functional in the center of the Brillouin zone and interpolated over the entire Brillouin zone using the Keating phenomenological model. The interaction parameters are determined by optimization of the IR- and Ramanactive frequencies for a phenomenological model by their comparison with the results of ab initio calculations. Numerical values of short-range interaction constants and charges are in agreement with the characteristics of the chemical bond calculated ab initio. These parameters have transparent physical meaning and chemical nature and can further be used for both qualitative estimations of any physical and physico-chemical quantities and quantitative calculations of the phonon spectra of a number of isostructural compounds. The Keating phenomenological model is used to study the genesis of the phonon spectrum from the spectra of sublattices.
Electron Transport through Polyene Junctions in between Carbon Nanotubes: an Ab Initio Realization
NASA Astrophysics Data System (ADS)
Chen, Yiing-Rei; Chen, Kai-Yu; Dou, Kun-Peng; Tai, Jung-Shen; Lee, Hsin-Han; Kaun, Chao-Cheng
With both ab initio and tight-binding model calculations, we study a system of polyene bridged armchair carbon nanotube electrodes, considering one-polyene and two-polyene cases, to address aspects of quantum transport through junctions with multiple conjugated molecules. The ab initio results of the two-polyene cases not only show the interference effect in transmission, but also the sensitive dependence of such effect on the combination of relative contact sites, which agrees nicely with the tight-binding model. Moreover, we show that the discrepancy mainly brought by ab initio relaxation provides an insight into the influence upon transmission spectra, from the junction's geometry, bonding and effective potential. This work was supported by the Ministry of Science and Technology of the Republic of China under Grant Nos. 99-2112-M-003-012-MY2 and 103-2622-E-002-031, and the National Center for Theoretical Sciences of Taiwan.
Skutterudites under pressure: An ab initio study
Ram, Swetarekha; Kanchana, V.; Valsakumar, M. C.
2014-03-07
Ab initio results on the band structure, density of states, and Fermi surface (FS) properties of LaRu{sub 4}X{sub 12} (X = P, As, Sb) are presented at ambient pressure as well as under compression. The analysis of density of states reveals the major contribution at the Fermi level to be mainly from the Ru-d and X-p states. We have a complicated Fermi surface with both electron and hole characters for all the three compounds which is derived mainly from the Ru-d and X-p states. There is also a simpler FS with hole character derived from the P-p{sub z} orbital for LaRu{sub 4}P{sub 12} and Ru-d{sub z{sup 2}} orbital in the case of As and Sb containing compounds. More interestingly, Fermi surface nesting feature is observed only in the case of the LaRu{sub 4}P{sub 12}. Under compression, we observe the topology of the complicated FS sheet of LaRu{sub 4}As{sub 12} to change around V/V{sub 0} = 0.85, leading to a behaviour similar to that of a multiband superconductor, and in addition, we have two more hole pockets centered around Γ at V/V{sub 0} = 0.8 for the same compound. Apart from this, we find the hole pocket to vanish at V/V{sub 0} = 0.8 in the case of LaRu{sub 4}Sb{sub 12} and the opening of the complicated FS sheet gets reduced. The de Haas van Alphen calculation shows the number of extremal orbits in the complicated sheet to change in As and Sb containing compounds under compression, where we also observe the FS topology to change.
NASA Astrophysics Data System (ADS)
Timko, Jeff; Kuyucak, Serdar
2012-11-01
Polarization is an important component of molecular interactions and is expected to play a particularly significant role in inhomogeneous environments such as pores and interfaces. Here we investigate the effects of polarization in the gramicidin A ion channel by performing quantum mechanics/molecular mechanics molecular dynamics (MD) simulations and comparing the results with those obtained from classical MD simulations with non-polarizable force fields. We consider the dipole moments of backbone carbonyl groups and channel water molecules as well as a number of structural quantities of interest. The ab initio results show that the dipole moments of the carbonyl groups and water molecules are highly sensitive to the hydrogen bonds (H-bonds) they participate in. In the absence of a K+ ion, water molecules in the channel are quite mobile, making the H-bond network highly dynamic. A central K+ ion acts as an anchor for the channel waters, stabilizing the H-bond network and thereby increasing their average dipole moments. In contrast, the K+ ion has little effect on the dipole moments of the neighboring carbonyl groups. The weakness of the ion-peptide interactions helps to explain the near diffusion-rate conductance of K+ ions through the channel. We also address the sampling issue in relatively short ab initio MD simulations. Results obtained from a continuous 20 ps ab initio MD simulation are compared with those generated by sampling ten windows from a much longer classical MD simulation and running each window for 2 ps with ab initio MD. Both methods yield similar results for a number of quantities of interest, indicating that fluctuations are fast enough to justify the short ab initio MD simulations.
Ab Initio Studies of Stratospheric Ozone Depletion Chemistry
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Head-Gordon, Martin; Langhoff, Stephen R. (Technical Monitor)
1995-01-01
An overview of the current understanding of ozone depletion chemistry, particularly with regards the formation of the so-called Antarctic ozone hole, will be presented together with an outline as to how ab initio quantum chemistry can be used to further our understanding of stratospheric chemistry. The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of fluorine, chlorine, bromine and nitrogen oxide species will be demonstrated by presentation of some example studies. The ab initio results will be shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the theoretical results are shown to fill in the gaps and to resolve experimental controversies. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of halogen oxide species will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of experimental studies.
THERMODYNAMICS OF MATERIALS: FROM AB INITIO TO PHENOMENOLOGY
Turchi, P A
2004-09-24
Quantum mechanical-based (or ab initio) methods are used to predict the stability properties of materials although their application is limited to relatively simple systems in terms of structures and number of alloy components. However thermodynamics of complex multi-component alloys requires a more versatile approach afforded within the CALPHAD formalism. Despite its success, the lack of experimental data very often prevents the design of robust thermodynamic databases. After a brief survey of ab initio methodologies and CALPHAD, it will be shown how ab initio electronic structure methods can supplement in two ways CALPHAD for subsequent applications. The first one is rather immediate and concerns the direct input of ab initio energetics in CALPHAD databases. The other way, more involved, is the assessment of ab initio thermodynamics '{acute a} la CALPHAD'. It will be shown how these results can be used within CALPHAD to predict the equilibrium properties of multi-component alloys. Finally, comments will be made on challenges and future prospects.
Barrett, B R; Navratil, P; Vary, J P
2011-04-11
A long-standing goal of nuclear theory is to determine the properties of atomic nuclei based on the fundamental interactions among the protons and neutrons (i.e., nucleons). By adopting nucleon-nucleon (NN), three-nucleon (NNN) and higher-nucleon interactions determined from either meson-exchange theory or QCD, with couplings fixed by few-body systems, we preserve the predictive power of nuclear theory. This foundation enables tests of nature's fundamental symmetries and offers new vistas for the full range of complex nuclear phenomena. Basic questions that drive our quest for a microscopic predictive theory of nuclear phenomena include: (1) What controls nuclear saturation; (2) How the nuclear shell model emerges from the underlying theory; (3) What are the properties of nuclei with extreme neutron/proton ratios; (4) Can we predict useful cross sections that cannot be measured; (5) Can nuclei provide precision tests of the fundamental laws of nature; and (6) Under what conditions do we need QCD to describe nuclear structure, among others. Along with other ab initio nuclear theory groups, we have pursued these questions with meson-theoretical NN interactions, such as CD-Bonn and Argonne V18, that were tuned to provide high-quality descriptions of the NN scattering phase shifts and deuteron properties. We then add meson-theoretic NNN interactions such as the Tucson-Melbourne or Urbana IX interactions. More recently, we have adopted realistic NN and NNN interactions with ties to QCD. Chiral perturbation theory within effective field theory ({chi}EFT) provides us with a promising bridge between QCD and hadronic systems. In this approach one works consistently with systems of increasing nucleon number and makes use of the explicit and spontaneous breaking of chiral symmetry to expand the strong interaction in terms of a dimensionless constant, the ratio of a generic small momentum divided by the chiral symmetry breaking scale taken to be about 1 GeV/c. The resulting NN
Cirtog, M; Asselin, P; Soulard, P; Tremblay, B; Madebène, B; Alikhani, M E; Georges, R; Moudens, A; Goubet, M; Huet, T R; Pirali, O; Roy, P
2011-03-31
A series of hydrogen bonded complexes involving oxirane and water molecules have been studied. In this paper we report on the vibrational study of the oxirane-water complex (CH(2))(2)O-H(2)O. Neon matrix experiments and ab initio anharmonic vibrational calculations have been performed, providing a consistent set of vibrational frequencies and anharmonic coupling constants. The implementation of a new large flow supersonic jet coupled to the Bruker IFS 125 HR spectrometer at the infrared AILES beamline of the French synchrotron SOLEIL (Jet-AILES) enabled us to record first jet-cooled Fourier transform infrared spectra of oxirane-water complexes at different resolutions down to 0.2 cm(-1). Rovibrational parameters and a lower bound of the predissociation lifetime of 25 ps for the v(OH)(b) = 1 state have been derived from the rovibrational analysis of the ν(OH)(b) band contour recorded at respective rotational temperatures of 12 K (Jet-AILES) and 35 K (LADIR jet). PMID:21381647
Ab initio simulation of transport phenomena in rarefied gases.
Sharipov, Felix; Strapasson, José L
2012-09-01
Ab initio potentials are implemented into the direct simulation Monte Carlo (DSMC) method. Such an implementation allows us to model transport phenomena in rarefied gases without any fitting parameter of intermolecular collisions usually extracted from experimental data. Applying the method proposed by Sharipov and Strapasson [Phys. Fluids 24, 011703 (2012)], the use of ab initio potentials in the DSMC requires the same computational efforts as the widely used potentials such as hard spheres, variable hard sphere, variable soft spheres, etc. At the same time, the ab initio potentials provide more reliable results than any other one. As an example, the transport coefficients of a binary mixture He-Ar, viz., viscosity, thermal conductivity, and thermal diffusion factor, have been calculated for several values of the mole fraction. PMID:23030889
Microsolvation of methyl hydrogen peroxide: Ab initio quantum chemical approach
NASA Astrophysics Data System (ADS)
Kulkarni, Anant D.; Rai, Dhurba; Bartolotti, Libero J.; Pathak, Rajeev K.
2009-08-01
Methyl hydrogen peroxide (MHP), one of the simplest organic hydroperoxides, is a strong oxidant, with enhanced activity in aqueous ambience. The present study investigates, at the molecular level, the role of hydrogen bonding that is conducive to cluster formation of MHP with water molecules from its peroxide end, with the methyl group remaining hydrophobic for up to five water molecules. Ab initio quantum chemical computations on MHP⋯(H2O)n, [n =1-5] are performed at second order Møller-Plesset (MP2) perturbation theory employing the basis sets 6-31G(d,p) and 6-311++G(2d,2p) to study the cluster formation of MHP with water molecules from its peroxide end and hydrophobic hydration due to the methyl group. Successive addition of water molecules alters the hydrogen bonding pattern, which leads to changes in overall cluster geometry and in turn to IR vibrational frequency shifts. Molecular co-operativity in these clusters is gauged directly through a detailed many-body interaction energy analysis. Molecular electrostatic potential maps are shown to have a bearing on predicting further growth of these clusters, which is duly corroborated through sample calculations for MHP⋯(H2O)8. Further, a continuum solvation model calculation for energetically stable clusters suggests that this study should serve as a precursor for pathways to aqueous solvation of MHP.
Ab initio study of the phenol-water cation radical
NASA Astrophysics Data System (ADS)
Hobza, Pavel; Burcl, Rudolf; Špirko, Vladimír; Dopfer, Otto; Müller-Dethlefs, Klaus; Schlag, Edward W.
1994-07-01
The phenol-water cation radical has been investigated by ab initio theory using the spin-restricted open-shell Hartree-Fock and spin-restricted open-shell second-order Møller-Plesset theories with 3-21G*(O) and 6-31G* basis sets. The full geometrical optimization was performed for several hydrogen-bonded structures and one hemibonded structure. Clearly, the most stable structure has been found for Cs symmetry with the linear hydrogen bond between the proton of the OH group of the phenol cation radical and the oxygen of the water, and the water hydrogens pointing away from the phenyl ring. For this structure harmonic (and for some intermolecular modes anharmonic) vibrational frequencies have been computed for various isotopic complexes. The computed shifts of phenol-localized intramolecular modes on complexation and on deuteration as well as the calculated intermolecular frequencies of the different isotopic complexes allow for an assignment of vibrational frequencies observed in the experimental zero-kinetic-energy (ZEKE) photoelectron spectra. Five out of a possible six intermolecular vibrations and several intramolecular modes have been assigned, including the 18b vibration which shows a strong blue shift in frequency upon complexation. Structure and properties of the phenol-water cation radical are compared with those of the corresponding neutral complex.
Ab initio calculations and analysis of chemical bonding in SrTiO3 and SrZrO3 cubic crystals
NASA Astrophysics Data System (ADS)
Evarestov, R. A.; Tupitsyn, I. I.; Bandura, A. V.; Alexandrov, V. E.
The possibility of the different first-principles methods to describe the chemical bonding in SrTiO3 and SrZrO3 cubic crystals is investigated. The local properties of the electronic structure (atomic charges, bond orders, atomic delocalization indexes, and polarization fractions) were calculated with different methods: traditional Mulliken population analysis in LCAO calculations, two projection techniques in plane-wave (PW) calculations, population analysis based on Wannier-type atomic orbitals, and chemical bonding analysis based on the localized Wannier functions for occupied (valence band) LCAO states. All the techniques considered except the traditional Mulliken analysis demonstrate that the ionicity of chemical bonding in SrZrO3 is larger than in SrTiO3, in agreement with the Zr and Ti electronegativities relation and the relative bandgaps observed.
NASA Astrophysics Data System (ADS)
Ojha, Animesh K.; Srivastava, Sunil K.; Koster, J.; Shukla, M. K.; Leszczynski, J.; Asthana, B. P.; Kiefer, W.
2004-02-01
Raman spectra of neat formamide (HCONH 2) and its binary mixture (HCONH 2+CH 3OH) with hydrogen donor solvent, methanol (CH 3OH) were investigated using a sensitive scanning multichannel detection scheme, which is simultaneously more precise also, especially when the observed Raman line profile has multiple component bands. The spectra in the two regions, namely 1200-1500 and 1500-1800 cm -1 were recorded with varying mole fractions of the reference molecule, HCONH 2, from 0.1 to 0.9. The spectra in the region 1200-1500 cm -1 show a broad band at ˜1312 cm -1, which shows a peculiar concentration dependence, and a relatively sharp peak at ˜1392 cm -1, whose peak position is not influenced by concentration. The spectra in the region 1500-1800 cm -1 also show two peaks, one at ˜1593 cm -1 and the other one at ˜1668 cm -1 which are assigned to NH 2 bending and ν(CO) stretching vibrations, respectively. Both these Raman bands show an appreciable upshift of ˜15-20 cm -1 and the one at ˜1668 cm -1 has also a distinct asymmetry towards higher wavenumber. The optimized geometries and vibrational wavenumbers of various normal modes for neat formamide as well as its hydrogen-bonded complexes were also calculated using ab initio theory at the MP2 level. The results have been used to understand and explain the concentration dependent changes in the spectral features in terms of hydrogen bonding and self-association.
Ab initio computations of photodissociation products of CFC alternatives
Tai, S.; Illinger, K.H.; Kenny, J.E.
1995-12-31
Ab initio computations, have already been used to examine the energetics of the photodissociation of stratospheric chlorofluorocarbons. Our awn research has investigated the ab initio computation of vibrational frequencies and infrared intensities of CF{sub 3}CH{sub 2}F, CF{sub 3}CF{sub 2}H, and CF{sub 3}CH{sub 3}; continuing research will attempt to expand these computations to the energetics of the photodissociation of these molecules, since sane of the most common types of chlorofluorocarbon substitutes are hydrofluoroethanes.
Ab Initio Structure Analysis Using Laboratory Powder Diffraction Data
NASA Astrophysics Data System (ADS)
Sasaki, Akito
Today, laboratory X-ray diffractometers are seeing increasingly wide use in the ab initio crystal structure analysis of organic powder samples. This is because optics and optical devices have been improved, making it possible to obtain precise integrated intensities of reflections in high 2-theta ranges. Another reason is that one can use direct-space methods, which do not require “high-resolution diffraction data”, much more easily than before. Described here are some key points to remember when performig ab initio crystal structure analysis using powder diffraction data from organic compounds.
Esrafili, Mehdi D
2012-05-01
We report geometries, stabilization energies, symmetry adapted perturbation theory (SAPT) and quantum theory of atoms in molecules (QTAIM) analyses of a series of carbene-BX(3) complexes, where X = H, OH, NH(2), CH(3), CN, NC, F, Cl, and Br. The stabilization energies were calculated at HF, B3LYP, MP2, MP4 and CCSD(T)/aug-cc-pVDZ levels of theory using optimized geometries of all the complexes obtained from B3LYP/aug-cc-pVTZ. Quantitatively, all the complexes indicate the presence of B-C(carbene) interaction due to the short B-C(carbene) distances. Inspection of stabilization energies reveals that the interaction energies increase in the order NH(2) > OH > CH(3) > F > H > Cl > Br > NC > CN, which is the opposite trend shown in the binding distances. Considering the SAPT results, it is found that electrostatic effects account for about 50% of the overall attraction of the studied complexes. By comparison, the induction components of these interactions represent about 40% of the total attractive forces. Despite falling in a region of charge depletion with nabla(2)ρ(BCP) >0, the B-C(carbene) bond critical points (BCPs) are characterized by a reasonably large value of the electron density (ρ(BCP)) and H(BCP) <0, indicating that the potential energy overcomes the kinetic energy density at BCP and the B-C(carbene) bond is a polar covalent bond. PMID:21877151
Zaleski, Daniel P; Mullaney, John C; Bittner, Dror M; Tew, David P; Walker, Nicholas R; Legon, Anthony C
2015-10-28
Strongly bound complexes (CH2)3⋯MCl (M = Cu or Ag), formed by non-covalent interaction of cyclopropane and either cuprous chloride or argentous chloride, have been generated in the gas phase by means of the laser ablation of either copper or silver metal in the presence of supersonically expanded pulses of a gas mixture containing small amounts of cyclopropane and carbon tetrachloride in a large excess of argon. The rotational spectra of the complexes so formed were detected with a chirped-pulse, Fourier transform microwave spectrometer and analysed to give rotational constants and Cu and Cl nuclear quadrupole coupling constants for eight isotopologues of each of (CH2)3⋯CuCl and (CH2)3⋯AgCl. The geometry of each of these complexes was established unambiguously to have C(2v) symmetry, with the three C atoms coplanar, and with the MCl molecule lying along a median of the cyclopropane C3 triangle. This median coincides with the principal inertia axis a in each of the two complexes (CH2)3⋯MCl. The M atom interacts with the pseudo-π bond linking the pair of equivalent carbon atoms (F)C (F = front) nearest to it, so that M forms a non-covalent bond to one C-C edge of the cyclopropane molecule. The (CH2)3⋯MCl complexes have similar angular geometries to those of the hydrogen- and halogen-bonded analogues (CH2)3⋯HCl and (CH2)3⋯ClF, respectively. Quantitative details of the geometries were determined by interpretation of the observed rotational constants and gave results in good agreement with those from ab initio calculations carried out at the CCSD(T)(F12*)/aug-cc-pVTZ-F12 level of theory. Interesting geometrical features are the lengthening of the (F)C-(F)C bond and the shrinkage of the two equivalent (B)C-(F)C (B = back) bonds relative to the C-C bond in cyclopropane itself. The expansions of the (F)C-(F)C bond are 0.1024(9) Å and 0.0727(17) Å in (CH2)3⋯CuCl and (CH2)3⋯AgCl, respectively, according to the determined r0 geometries. The C-C bond
Zaleski, Daniel P.; Mullaney, John C.; Bittner, Dror M.; Walker, Nicholas R. E-mail: nick.walker@newcastle.ac.uk; Tew, David P.; Legon, Anthony C. E-mail: nick.walker@newcastle.ac.uk
2015-10-28
Strongly bound complexes (CH{sub 2}){sub 3}⋯MCl (M = Cu or Ag), formed by non-covalent interaction of cyclopropane and either cuprous chloride or argentous chloride, have been generated in the gas phase by means of the laser ablation of either copper or silver metal in the presence of supersonically expanded pulses of a gas mixture containing small amounts of cyclopropane and carbon tetrachloride in a large excess of argon. The rotational spectra of the complexes so formed were detected with a chirped-pulse, Fourier transform microwave spectrometer and analysed to give rotational constants and Cu and Cl nuclear quadrupole coupling constants for eight isotopologues of each of (CH{sub 2}){sub 3}⋯CuCl and (CH{sub 2}){sub 3}⋯AgCl. The geometry of each of these complexes was established unambiguously to have C{sub 2v} symmetry, with the three C atoms coplanar, and with the MCl molecule lying along a median of the cyclopropane C{sub 3} triangle. This median coincides with the principal inertia axis a in each of the two complexes (CH{sub 2}){sub 3}⋯MCl. The M atom interacts with the pseudo-π bond linking the pair of equivalent carbon atoms {sup F}C (F = front) nearest to it, so that M forms a non-covalent bond to one C—C edge of the cyclopropane molecule. The (CH{sub 2}){sub 3}⋯MCl complexes have similar angular geometries to those of the hydrogen- and halogen-bonded analogues (CH{sub 2}){sub 3}⋯HCl and (CH{sub 2}){sub 3}⋯ClF, respectively. Quantitative details of the geometries were determined by interpretation of the observed rotational constants and gave results in good agreement with those from ab initio calculations carried out at the CCSD(T)(F12*)/aug-cc-pVTZ-F12 level of theory. Interesting geometrical features are the lengthening of the {sup F}C—{sup F}C bond and the shrinkage of the two equivalent {sup B}C—{sup F}C (B = back) bonds relative to the C—C bond in cyclopropane itself. The expansions of the {sup F}C—{sup F}C bond are 0
NASA Astrophysics Data System (ADS)
Zaleski, Daniel P.; Mullaney, John C.; Bittner, Dror M.; Tew, David P.; Walker, Nicholas R.; Legon, Anthony C.
2015-10-01
Strongly bound complexes (CH2)3⋯MCl (M = Cu or Ag), formed by non-covalent interaction of cyclopropane and either cuprous chloride or argentous chloride, have been generated in the gas phase by means of the laser ablation of either copper or silver metal in the presence of supersonically expanded pulses of a gas mixture containing small amounts of cyclopropane and carbon tetrachloride in a large excess of argon. The rotational spectra of the complexes so formed were detected with a chirped-pulse, Fourier transform microwave spectrometer and analysed to give rotational constants and Cu and Cl nuclear quadrupole coupling constants for eight isotopologues of each of (CH2)3⋯CuCl and (CH2)3⋯AgCl. The geometry of each of these complexes was established unambiguously to have C2v symmetry, with the three C atoms coplanar, and with the MCl molecule lying along a median of the cyclopropane C3 triangle. This median coincides with the principal inertia axis a in each of the two complexes (CH2)3⋯MCl. The M atom interacts with the pseudo-π bond linking the pair of equivalent carbon atoms FC (F = front) nearest to it, so that M forms a non-covalent bond to one C—C edge of the cyclopropane molecule. The (CH2)3⋯MCl complexes have similar angular geometries to those of the hydrogen- and halogen-bonded analogues (CH2)3⋯HCl and (CH2)3⋯ClF, respectively. Quantitative details of the geometries were determined by interpretation of the observed rotational constants and gave results in good agreement with those from ab initio calculations carried out at the CCSD(T)(F12*)/aug-cc-pVTZ-F12 level of theory. Interesting geometrical features are the lengthening of the FC—FC bond and the shrinkage of the two equivalent BC—FC (B = back) bonds relative to the C—C bond in cyclopropane itself. The expansions of the FC—FC bond are 0.1024(9) Å and 0.0727(17) Å in (CH2)3⋯CuCl and (CH2)3⋯AgCl, respectively, according to the determined r0
Motif based Hessian matrixfor ab initio geometry optimization ofnanostructures
Zhao, Zhengji; Wang, Lin-Wang; Meza, Juan
2006-04-05
A simple method to estimate the atomic degree Hessian matrixof a nanosystem is presented. The estimated Hessian matrix, based on themotif decomposition of the nanosystem, can be used to accelerate abinitio atomic relaxations with speedups of 2 to 4 depending on the sizeof the system. In addition, the programing implementation for using thismethod in a standard ab initio package is trivial.
Towards SiC Surface Functionalization: An Ab Initio Study
Cicero, G; Catellani, A
2005-01-28
We present a microscopic model of the interaction and adsorption mechanism of simple organic molecules on SiC surfaces as obtained from ab initio molecular dynamics simulations. Our results open the way to functionalization of silicon carbide, a leading candidate material for bio-compatible devices.
Ab initio pseudopotential band calculation of organic conductors
Ishibashi, Shoji; Kohyama, Masanori
1999-12-01
The authors have calculated the band structures of organic conductors TTF-TCNQ and {beta}-(BEDT-TTF){sub 2}I{sub 3} using the ab initio plane-wave pseudopotential method within the local-density approximation (LDA). The Fermi-surface shape and the origin of bands near the Fermi level are investigated for each compound.
Takahasi, S. ); Curtiss, L.A.; Gosztola, D.; Koura, N. ); Loong, C.K.; Saboungi, M.L. . Materials Science Div.)
1993-04-01
The Raman and neutron scattering spectra of 46 mol% AlCl[sub 3] -54 mol% 1-ethyl-3-methyl imidazolium chloride (EMIC) and 67 mol% AlCl[sub 3] - 33 mol% EMIC melts are presented. Ab initio molecular orbital calculations have been carried out on structures of chloroaluminate anion and EMI cation and the interaction between anion and cation.
Haldar, Susanta; Gnanasekaran, Ramachandran; Hobza, Pavel
2015-10-28
Dissociation energies (D0) of 11 H-bonded and 11 dispersion-bound complexes were calculated as the sum of interaction energies and the change of zero-point vibrational energies (ΔZPVE). The structures of H-bonded complexes were optimized at the RI-MP2/cc-pVTZ level, at which deformation and harmonic ΔZPVE energies were also calculated. The structures of dispersion-bound complexes were optimized at the DFT-D3 level, and harmonic ΔZPVE energies were determined at the same level as well. For comparison, CCSD(T)/CBS D0 energies were also evaluated for both types of complexes. The CCSD(T)/CBS interaction energy was constructed as the sum of MP2/CBS interaction energy, extrapolated from aug-cc-pVTZ and aug-cc-pVQZ basis sets, and ΔCCSD(T) correction, determined with the aug-cc-pVDZ basis set. The ΔZPVE energies were determined for all complexes at the harmonic level and for selected complexes, these energies were also calculated using second-order vibration perturbation (VPT2) theory. For H-bonded complexes, the harmonic CCSD(T)/CBS D0 energies were in better agreement with the experimental values (with a mean relative error (MRE) of 6.2%) than the RI-MP2/cc-pVTZ D0 (a MRE of 12.3%). The same trend was found for dispersion-bound complexes (6.2% (MRE) at CCSD(T)/CBS and 7.7% (MRE) at the DFT-D3 level). When the anharmonic ΔZPVE term was included instead of harmonic one, the agreement between theoretical and experimental D0 deteriorated for H-bonded as well as dispersion-bound complexes. Finally, the applicability of "diagonal approximation" for determining the anharmonic ΔZPVE was shown. For the phenolH2O complex, the ΔZPVE energy calculated at the VPT2 level and on the basis of "diagonal approximation" differed by less than 0.1 kcal mol(-1). PMID:26392236
Ab initio molecular dynamics of solvation effects on reactivity at electrified interfaces.
Herron, Jeffrey A; Morikawa, Yoshitada; Mavrikakis, Manos
2016-08-23
Using ab initio molecular dynamics as implemented in periodic, self-consistent (generalized gradient approximation Perdew-Burke-Ernzerhof) density functional theory, we investigated the mechanism of methanol electrooxidation on Pt(111). We investigated the role of water solvation and electrode potential on the energetics of the first proton transfer step, methanol electrooxidation to methoxy (CH3O) or hydroxymethyl (CH2OH). The results show that solvation weakens the adsorption of methoxy to uncharged Pt(111), whereas the binding energies of methanol and hydroxymethyl are not significantly affected. The free energies of activation for breaking the C-H and O-H bonds in methanol were calculated through a Blue Moon Ensemble using constrained ab initio molecular dynamics. Calculated barriers for these elementary steps on unsolvated, uncharged Pt(111) are similar to results for climbing-image nudged elastic band calculations from the literature. Water solvation reduces the barriers for both C-H and O-H bond activation steps with respect to their vapor-phase values, although the effect is more pronounced for C-H bond activation, due to less disruption of the hydrogen bond network. The calculated activation energy barriers show that breaking the C-H bond of methanol is more facile than the O-H bond on solvated negatively biased or uncharged Pt(111). However, with positive bias, O-H bond activation is enhanced, becoming slightly more facile than C-H bond activation. PMID:27503889
Ab initio investigation of grain boundary cohesion in Al alloys
NASA Astrophysics Data System (ADS)
Zhang, Shengjun; Kontsevoi, Oleg Y.; Freeman, A. J.; Olson, G. B.
2010-03-01
Strength and hardness of aluminum alloys can be substantially increased by alloying with Mg, Zn, Cu, Si, and other elements. The main drawback of Al alloys is their susceptibility to stress corrosion cracking, which is caused by alloying impurities segregated at grain boundaries. We investigated the embrittling and cohesion-enhancing effects of impurities on a σ5(012)[100] grain boundary in Al by means of the full-potential linearized augmented plane-wave (FLAPW) method within the framework of the Rice-Wang thermodynamic model and within the ab initio tensile test approach. We calculated segregation energies, analyzed local atomic configurations, electronic structures and spatial charge density distributions around segregated impurities, and identified the roles of atomic size and the bonding behavior of the impurity with the surrounding Al atoms. The results show that He, H and Na are strong embrittlers, Zn is a weak embrittler, while Sc, B, Cu and Mg are cohesion enhancers. We further evaluated the effect of co-alloying with two or more elements on grain boundary strength. This work provides a fundamental basis for the design of high strength Al alloys.
FTIR, Raman spectra and ab initio calculations of 2-mercaptobenzothiazole.
Rai, Amareshwar K; Singh, Rachana; Singh, K N; Singh, V B
2006-02-01
FTIR and Raman spectra of a rubber vulcanization accelerator, 2-mercaptobenzothiazole (MBT), were recorded in the solid phase. The harmonic vibrational wavenumbers, for both the toutomeric forms of MBT, as well as for its dimeric complex, have been calculated, using ab initio RHF and density functional B3LYP methods invoking different basis sets upto RHF/6-31G** and B3LYP/6-31G** and the results were compared with the experimental values. Conformational studies have been also carried out regarding its toutomeric monomer forms and its dimer form. With all the basis sets the thione form of MBT (II) is predicted to be more stable than thiol form (I) and dimeric conformation (III) is predicted to be more stable with monomeric conformations (I) and (II). Vibrational assignments have been made, and it has been found that the calculated normal mode frequencies of dimeric conformation (III) are required for the analysis of IR and Raman bands of the MBT. The predicted shift in NH- stretching vibration towards the lower wave number side with the B3LYP/6-31G** calculations for the most stable dimer form (III), is in better agreement with experimental results. The intermolecular sulfur-nitrogen distance in N-H...S hydrogen bond was found to be 3.35 angstroms from these calculations, is also in agreement to the experimental value. PMID:16098806
Ab initio Raman spectroscopy of water under extreme conditions
NASA Astrophysics Data System (ADS)
Rozsa, Viktor; Pan, Ding; Wan, Quan; Galli, Giulia
Water exhibits one of the most complex phase diagrams of any binary compound. Despite extensive studies, the melting lines of high-pressure ice phases remain very controversial, with reports differing by hundreds of Kelvin. The boundary between ice VII and liquid phase is particularly disputed, with recent work exploring plasticity and amorphization mediating the transition. Raman measurements are often used to fingerprint melting, yet their interpretation is difficult without atomistic modeling. Here, we report a study of high P/T water where we computed Raman spectra using a method combining ab initio molecular dynamics and density functional perturbation theory, as implemented in the Qbox code. Spectra were computed for the liquid at 10 and 20 GPa, both at 1000 K, and for solid ice VII (20 GPa, 500 K). Decomposing the spectra into inter and intra molecular contributions provided insight into the dynamics of the hydrogen-bonded network at extreme conditions. The relevance of our simulation results for models of water in Earth, Uranus, and Neptune will be discussed, and an interpretation of existing experiments at high pressure will be presented.
Miller, A.E.
1993-01-01
This thesis is composed of two studies in catalysis. The first is an exploration, using computational techniques, of the mechanism for the first carbon-carbon bond formation in the methanol to gasoline (MTG) reaction. The second is a study of the factors important to the understanding of ligand effects on metal-metal bonds, and in particular, to metal clusters. Three possibilities were considered as ways to activate a carbon in the MTG process prior to formation of C2 or higher hydrocarbons. These were a free radical mechanism, a surface ylide mechanism, and a possible defect site which might lead to steric crowding of CH[sub 2] groups. Although the free radical mechanism was found to be thermodynamically within the parameters of the MTG process, it contained a high transition state. Consideration of the molecules available prior to hydrocarbon build-up and their specific electronic structure, led to the view that the available carbon atoms (methanol, dimethyl ether, etc.) were unlikely to be activated by a free radical intermediate. The surface-stabilized ylide which has been proposed as an intermediate by many was studied to determine if in fact the ylide was stabilized. The total energy of the ylide was compared to that of the naked site on the zeolite and free methylene. Free methylene ranged, depending on the geometry of the ylide, between 50 and 80 kcal more stable. These numbers are qualitatively correct, but more electron correlation would have to be incorporated in the calculation to get an accurate value for the destabilization. Starting from a defect site, two CH[sub 2] groups were each attached to two oxygen atoms. It was thought that two CH[sub 2] groups would take up considerably more space than either than original Al atom or the four hydrogens. Molecular Mechanics calculations showed the zeolite to be sufficiently flexible to prevent crowding of the CH[sub 2]'s.
Ab initio molecular crystal structures, spectra, and phase diagrams.
Hirata, So; Gilliard, Kandis; He, Xiao; Li, Jinjin; Sode, Olaseni
2014-09-16
Conspectus Molecular crystals are chemists' solids in the sense that their structures and properties can be understood in terms of those of the constituent molecules merely perturbed by a crystalline environment. They form a large and important class of solids including ices of atmospheric species, drugs, explosives, and even some organic optoelectronic materials and supramolecular assemblies. Recently, surprisingly simple yet extremely efficient, versatile, easily implemented, and systematically accurate electronic structure methods for molecular crystals have been developed. The methods, collectively referred to as the embedded-fragment scheme, divide a crystal into monomers and overlapping dimers and apply modern molecular electronic structure methods and software to these fragments of the crystal that are embedded in a self-consistently determined crystalline electrostatic field. They enable facile applications of accurate but otherwise prohibitively expensive ab initio molecular orbital theories such as Møller-Plesset perturbation and coupled-cluster theories to a broad range of properties of solids such as internal energies, enthalpies, structures, equation of state, phonon dispersion curves and density of states, infrared and Raman spectra (including band intensities and sometimes anharmonic effects), inelastic neutron scattering spectra, heat capacities, Gibbs energies, and phase diagrams, while accounting for many-body electrostatic (namely, induction or polarization) effects as well as two-body exchange and dispersion interactions from first principles. They can fundamentally alter the role of computing in the studies of molecular crystals in the same way ab initio molecular orbital theories have transformed research practices in gas-phase physical chemistry and synthetic chemistry in the last half century. In this Account, after a brief summary of formalisms and algorithms, we discuss applications of these methods performed in our group as compelling
NASA Astrophysics Data System (ADS)
Saini, Jyoti; Viswanathan, K. S.
2016-08-01
When hydrogen bonded complexes are formed with more than one contact, the question arises if these multiple contacts operate synergistically. Propargyl alcohol-H2O complex presents a good case study to address this question, which is discussed in this work. Complexes of propargyl alcohol (PA) and H2O were studied experimentally using matrix isolation infrared spectroscopy, which was supported by quantum chemical computations performed at the M06-2X and MP2 level of theories, using 6-311++G (d,p) and aug-cc-pVDZ basis sets. A 1:1 PA-H2O complex was identified in the experiments and corroborated by our computations, where the PA was in the gauche conformation. This complex, which was a global minimum, showed dual interactions, one of which was an n-σ interaction between the O-H group of PA and the O of H2O, while the second was a H···π contact between the O-H group of H2O and the π system of PA. We explored if the two interactions in the 1:1 complex exhibited synergism. We finally argue that the two interactions showed antagonism rather than synergism. Our computations indicated three other local minima for the 1:1 complexes; though these local minima were not identified in our experiments. Atoms-in-molecules and energy decomposition analysis executed through LMO-EDA were also performed to understand the nature of intermolecular interactions in the PA-H2O complexes. We have also revisited the problem of conformations of PA, with a view to understanding the reasons for gauche conformational preferences in PA.
Ab-initio calculations on melting of thorium
NASA Astrophysics Data System (ADS)
Mukherjee, D.; Sahoo, B. D.; Joshi, K. D.; Kaushik, T. C.; Gupta, Satish C.
2016-05-01
Ab-initio molecular dynamics study has been performed on face centered cubic structured thorium to determine its melting temperature at room pressure. The ion-electron interaction potential energy calculated as a function of temperature for three volumes (a0)3 and (1.02a0)3 and (1.04a0)3 increases gradually with temperature and undergoes a sharp jump at ~2200 K, ~2100 K and ~1800 K, respectively. Here, a0 = 5.043 Å is the equilibrium lattice parameter at 0 K obtained from ab-initio calculations. These jumps in interaction energy are treated as due to the onset of melting and corresponding temperatures as melting point. The melting point of 2100 K is close to the experimental value of 2023K. Further, the same has been verified by plotting the atomic arrangement evolved at various temperatures and corresponding pair correlation functions.
Ab initio molecular dynamics: concepts, recent developments, and future trends.
Iftimie, Radu; Minary, Peter; Tuckerman, Mark E
2005-05-10
The methodology of ab initio molecular dynamics, wherein finite-temperature dynamical trajectories are generated by using forces computed "on the fly" from electronic structure calculations, has had a profound influence in modern theoretical research. Ab initio molecular dynamics allows chemical processes in condensed phases to be studied in an accurate and unbiased manner, leading to new paradigms in the elucidation of microscopic mechanisms, rationalization of experimental data, and testable predictions of new phenomena. The purpose of this work is to give a brief introduction to the technique and to review several important recent developments in the field. Several illustrative examples showing the power of the technique have been chosen. Perspectives on future directions in the field also will be given. PMID:15870204
Recent progress in ab initio density matrix renormalization group methodology
NASA Astrophysics Data System (ADS)
Hachmann, Johannes; Dorando, Jonathan J.; Kin-Lic Chan, Garnet
2008-03-01
We present some recent developments in the ab initio density matrix renormalization group (DMRG) method for quantum chemical problems, in particular our local, quadratic scaling algorithm [1] for low dimensional systems. This method is particularly suited for the description of strong nondynamic correlation, and allows us to compute numerically exact (FCI) correlated energies for large active spaces, up to one order of magnitude larger then can be done by conventional CASCI techniques. Other features of this method are its inherent multireference nature, compactness, variational results, size-consistency and size-extensivity. In addition we will review the problems (predominantly organic electronic materials) on which we applied the ab initio DMRG: 1) metal-insulator transition in hydrogen chains [1] 2) all-trans polyacetylene [1] 3) acenes [2] 4) polydiacetylenes [3]. References [1] Hachmann, Cardoen, Chan, JCP 125 (2006), 144101. [2] Hachmann, Dorando, Avil'es, Chan, JCP 127 (2007), 134309. [3] unpublished.
Understanding phonon transport in thermoelectric materials using ab initio approaches
NASA Astrophysics Data System (ADS)
Broido, David
Good thermoelectric materials have low phonon thermal conductivity, kph. Accurate theories to describe kph are important components in developing predictive models of thermoelectric efficiency that can help guide synthesis and measurement efforts. We have developed ab initio approaches to calculate kph, in which phonon modes and phonon scattering rates are computed using interatomic force constants determined from density functional theory, and a full solution of the Boltzmann transport equation for phonons is implemented. A recent approach to calculate interatomic force constants using ab initio molecular dynamics has yielded a good description of the thermal properties of Bi2Te3. But, the complexity of new promising candidate thermoelectric materials introduces computational challenges in assessing their thermal properties. An example is germanane, a germanium based hydrogen-terminated layered semiconductor, which we will discuss in this talk.
Ab Initio Calculations Of Light-Ion Reactions
Navratil, P; Quaglioni, S; Roth, R; Horiuchi, W
2012-03-12
The exact treatment of nuclei starting from the constituent nucleons and the fundamental interactions among them has been a long-standing goal in nuclear physics. In addition to the complex nature of nuclear forces, one faces the quantum-mechanical many-nucleon problem governed by an interplay between bound and continuum states. In recent years, significant progress has been made in ab initio nuclear structure and reaction calculations based on input from QCD employing Hamiltonians constructed within chiral effective field theory. In this contribution, we present one of such promising techniques capable of describing simultaneously both bound and scattering states in light nuclei. By combining the resonating-group method (RGM) with the ab initio no-core shell model (NCSM), we complement a microscopic cluster approach with the use of realistic interactions and a microscopic and consistent description of the clusters. We discuss applications to light nuclei scattering, radiative capture and fusion reactions.
Spin-orbit decomposition of ab initio nuclear wave functions
NASA Astrophysics Data System (ADS)
Johnson, Calvin W.
2015-03-01
Although the modern shell-model picture of atomic nuclei is built from single-particle orbits with good total angular momentum j , leading to j -j coupling, decades ago phenomenological models suggested that a simpler picture for 0 p -shell nuclides can be realized via coupling of the total spin S and total orbital angular momentum L . I revisit this idea with large-basis, no-core shell-model calculations using modern ab initio two-body interactions and dissect the resulting wave functions into their component L - and S -components. Remarkably, there is broad agreement with calculations using the phenomenological Cohen-Kurath forces, despite a gap of nearly 50 years and six orders of magnitude in basis dimensions. I suggest that L -S decomposition may be a useful tool for analyzing ab initio wave functions of light nuclei, for example, in the case of rotational bands.
Exploring complex chemical reactions by ab-initio simulation
NASA Astrophysics Data System (ADS)
Parrinello, Michele
1998-03-01
Recent progress in the ab-initio molecular dynamics method and the power of parallel computing, allow the detailed study of complex chemical reaction of great industrial relevance. We illustrate this unprecedented capability by investigating the second generation Ziegler-Natta catalytic process. In this inhomogeneous catalyst, a polymerization reaction is induced by TiCl4 molecules deposited on an MgCl2 solid support. A density functional based ab-initio molecular dynamics calculation conducted with a minimum of initial assumption allows to understand the nature of the catalytic center and to determine the reaction path with the associated free energy barrier. Furthermore our calculation can explain in a nontrivial way the stereo-selectivity of the process.
Ab Initio Electronic Relaxation Times and Transport in Noble Metals
NASA Astrophysics Data System (ADS)
Mustafa, Jamal I.; Bernardi, Marco; Neaton, Jeffrey B.; Louie, Steven G.
Relaxation times employed to study electron transport in metals are typically assumed to be constants and obtained empirically using the Drude model. Here, we employ ab initio calculations to compute the electron-phonon relaxation times of Cu, Ag, and Au, and find that they vary significantly on the Fermi surface, spanning ~15 -45 fs. We compute room temperature resistivities in excellent agreement with experiment by combining GW bandstructures, Wannier-interpolated band velocities, and ab initio relaxation times. Our calculations are compared to other approximations used for the relaxation times. Additionally, an importance sampling scheme is introduced to speed up the convergence of resistivity and transport calculations by sampling directly points on the Fermi surface. This work was supported by NSF Grant No. DMR15-1508412 and U.S. DOE under Contract No. DE-AC02-05CH11231. Computational resources have been provided by DOE at LBNL's NERSC facility.
Ab initio Monte Carlo investigation of small lithium clusters.
Srinivas, S.
1999-06-16
Structural and thermal properties of small lithium clusters are studied using ab initio-based Monte Carlo simulations. The ab initio scheme uses a Hartree-Fock/density functional treatment of the electronic structure combined with a jump-walking Monte Carlo sampling of nuclear configurations. Structural forms of Li{sub 8} and Li{sub 9}{sup +} clusters are obtained and their thermal properties analyzed in terms of probability distributions of the cluster potential energy, average potential energy and configurational heat capacity all considered as a function of the cluster temperature. Details of the gradual evolution with temperature of the structural forms sampled are examined. Temperatures characterizing the onset of structural changes and isomer coexistence are identified for both clusters.
Towards AB Initio Calculation of the Circular Dichroism of Peptides
NASA Astrophysics Data System (ADS)
Molteni, E.; Onida, G.; Tiana, G.
2012-08-01
In this work we plan to use ab initio spectroscopy calculations to compute circular dichroism (CD) spectra of peptides. CD provides information on protein secondary structure content; peptides, instead, remain difficult to address, due to their tendency to adopt multiple conformations in equilibrium. Therefore peptides are an interesting test-case for ab initio calculation of CD spectra. As a first application, we focus on the (83-92) fragment of HIV-1 protease, which is known to be involved in the folding and dimerization of this protein. As a preliminary step, we performed classical molecular dynamics (MD) simulations, in order to obtain a set of representative conformers of the peptide. Then, on some of the obtained conformations, we calculated absorption spectra at the independent particle, RPA and TDLDA levels, showing the presence of charge transfer excitations, and their influence on spectral features.
Song, Lingchun; Han, Jaebeom; Lin, Yen-lin; Xie, Wangshen; Gao, Jiali
2010-01-01
The explicit polarization (X-Pol) method has been examined using ab initio molecular orbital theory and density functional theory. The X-Pol potential was designed to provide a novel theoretical framework for developing next-generation force fields for biomolecular simulations. Importantly, the X-Pol potential is a general method, which can be employed with any level of electronic structure theory. The present study illustrates the implementation of the X-Pol method using ab initio Hartree—Fock theory and hybrid density functional theory. The computational results are illustrated by considering a set of bimolecular complexes of small organic molecules and ions with water. The computed interaction energies and hydrogen bond geometries are in good accord with CCSD(T) calculations and B3LYP/aug-cc-pVDZ optimizations. PMID:19618944
GAUSSIAN 76: An ab initio Molecular Orbital Program
DOE R&D Accomplishments Database
Binkley, J. S.; Whiteside, R.; Hariharan, P. C.; Seeger, R.; Hehre, W. J.; Lathan, W. A.; Newton, M. D.; Ditchfield, R.; Pople, J. A.
1978-01-01
Gaussian 76 is a general-purpose computer program for ab initio Hartree-Fock molecular orbital calculations. It can handle basis sets involving s, p and d-type Gaussian functions. Certain standard sets (STO-3G, 4-31G, 6-31G*, etc.) are stored internally for easy use. Closed shell (RHF) or unrestricted open shell (UHF) wave functions can be obtained. Facilities are provided for geometry optimization to potential minima and for limited potential surface scans.
The study of molecular spectroscopy by ab initio methods
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.
1991-01-01
This review illustrates the potential of theory for solving spectroscopic problems. The accuracy of approximate techniques for including electron correlation have been calibrated by comparison with full configuration-interaction calculations. Examples of the application of ab initio calculations to vibrational, rotational, and electronic spectroscopy are given. It is shown that the state-averaged, complete active space self-consistent field, multireference configuration-interaction procedure provides a good approach for treating several electronic states accurately in a common molecular orbital basis.
Ab Initio Nuclear Structure and Reaction Calculations for Rare Isotopes
Draayer, Jerry P.
2014-09-28
We have developed a novel ab initio symmetry-adapted no-core shell model (SA-NCSM), which has opened the intermediate-mass region for ab initio investigations, thereby providing an opportunity for first-principle symmetry-guided applications to nuclear structure and reactions for nuclear isotopes from the lightest p-shell systems to intermediate-mass nuclei. This includes short-lived proton-rich nuclei on the path of X-ray burst nucleosynthesis and rare neutron-rich isotopes to be produced by the Facility for Rare Isotope Beams (FRIB). We have provided ab initio descriptions of high accuracy for low-lying (including collectivity-driven) states of isotopes of Li, He, Be, C, O, Ne, Mg, Al, and Si, and studied related strong- and weak-interaction driven reactions that are important, in astrophysics, for further understanding stellar evolution, X-ray bursts and triggering of s, p, and rp processes, and in applied physics, for electron and neutrino-nucleus scattering experiments as well as for fusion ignition at the National Ignition Facility (NIF).
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics
Makhov, Dmitry V.; Shalashilin, Dmitrii V.; Glover, William J.; Martinez, Todd J.
2014-08-07
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as “cloning,” in analogy to the “spawning” procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, “trains,” as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions.
Ab Initio Study of Defect Properties in YPO4
Gao, Fei; Xiao, Haiyan Y.; Zhou, Yungang; Devanathan, Ramaswami; Hu, Shenyang Y.; Li, Yulan; Sun, Xin; Khaleel, Mohammad A.
2012-03-01
Ab initio methods based on density functional theory have been used to calculate the formation energies of intrinsic defects, including vacancies, interstitials, antisites and Frenkel pairs in YPO4 under the O-rich and Y2O3-rich, and the O-rich and Y-rich conditions. The larger size of the yttrium atom may give rise to higher formation energy of the phosphorus antisite defect. In general, the formation energies of anion interstitials are much smaller than those of cation interstitials for both conditions considered. It is of greatly interest to find that the relative stabilities among the same types of interstitials are independent of the reference states. The most stable configuration for oxygen interstitials is an O-O split interstitial near the Ta site, while the most stable configuration for cation interstitials is a tetrahedral interstitial near the Ta site. The cation split interstitials are unfavorable in YPO4, with much higher formation energies. Furthermore, the properties of Frenkel pairs are compared with those calculated using empirical potentials. The results reveal that both ab initio and empirical potential calculations show a similar trend in the formation energies of Frenkel pairs, but the formation energies obtained by empirical potentials are much larger than those calculated by ab initio method.
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics
NASA Astrophysics Data System (ADS)
Makhov, Dmitry V.; Glover, William J.; Martinez, Todd J.; Shalashilin, Dmitrii V.
2014-08-01
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as "cloning," in analogy to the "spawning" procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions.
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics.
Makhov, Dmitry V; Glover, William J; Martinez, Todd J; Shalashilin, Dmitrii V
2014-08-01
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as "cloning," in analogy to the "spawning" procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions. PMID:25106573
A Complete and Accurate Ab Initio Repeat Finding Algorithm.
Lian, Shuaibin; Chen, Xinwu; Wang, Peng; Zhang, Xiaoli; Dai, Xianhua
2016-03-01
It has become clear that repetitive sequences have played multiple roles in eukaryotic genome evolution including increasing genetic diversity through mutation, changes in gene expression and facilitating generation of novel genes. However, identification of repetitive elements can be difficult in the ab initio manner. Currently, some classical ab initio tools of finding repeats have already presented and compared. The completeness and accuracy of detecting repeats of them are little pool. To this end, we proposed a new ab initio repeat finding tool, named HashRepeatFinder, which is based on hash index and word counting. Furthermore, we assessed the performances of HashRepeatFinder with other two famous tools, such as RepeatScout and Repeatfinder, in human genome data hg19. The results indicated the following three conclusions: (1) The completeness of HashRepeatFinder is the best one among these three compared tools in almost all chromosomes, especially in chr9 (8 times of RepeatScout, 10 times of Repeatfinder); (2) in terms of detecting large repeats, HashRepeatFinder also performed best in all chromosomes, especially in chr3 (24 times of RepeatScout and 250 times of Repeatfinder) and chr19 (12 times of RepeatScout and 60 times of Repeatfinder); (3) in terms of accuracy, HashRepeatFinder can merge the abundant repeats with high accuracy. PMID:26272474
Ab initio calculations of reactions with light nuclei
NASA Astrophysics Data System (ADS)
Quaglioni, Sofia; Hupin, Guillaume; Calci, Angelo; Navrátil, Petr; Roth, Robert
2016-03-01
An ab initio (i.e., from first principles) theoretical framework capable of providing a unified description of the structure and low-energy reaction properties of light nuclei is desirable to further our understanding of the fundamental interactions among nucleons, and provide accurate predictions of crucial reaction rates for nuclear astrophysics, fusion-energy research, and other applications. In this contribution we review ab initio calculations for nucleon and deuterium scattering on light nuclei starting from chiral two- and three-body Hamiltonians, obtained within the framework of the ab initio no-core shell model with continuum. This is a unified approach to nuclear bound and scattering states, in which square-integrable energy eigenstates of the A-nucleon system are coupled to (A-a)+a target-plus-projectile wave functions in the spirit of the resonating group method to obtain an efficient description of the many-body nuclear dynamics both at short and medium distances and at long ranges.
Chiroptical properties of unsubstituted carbohydrates: Ab initio and semiempirical studies
NASA Astrophysics Data System (ADS)
Parra C., Alejandro
Ab initio calculations support assignment of the vacuum ultraviolet circular dichroism (CD) of simple saccharides to 11A 1 --> 21B1 and 11A 1 --> 11A2 transitions centered on the oxygen atoms of the acetal group treated as two weakly coupled ether chromophores. The calculations are consistent with assignments previously made on the basis of a deconvolution of CD spectra. Estimates of the oxygen centered contributions to magnetic transition dipole moments were made. Semiempirical calculations were performed to model the NaD molar optical rotation of 1,6- and 3,6- anhydrosugars. For 1,6-anhydrosugars, current parameters produce reasonable agreement with experimental values. For 3,6-anhydrosugars, modifications to the ether parameters had to be introduced. The most relevant included a reorientation of the bond-centered s-->s* transition dipole charges in the ether chromophore to a C2v orientation, and a shift from prolate polarizability ellipsoids to general ellipsoids. These changes result in good agreement with experimental Na D molar rotations for 3,6-anhydrosugars. A low energy CD band arises in 3,6- and 1,6-anhydrosugars when agreement with the experimental NaD molar rotations is achieved. It is proposed that this band is a real feature in the spectrum. The origin of the band is primarily the interaction between b1 symmetry- oriented transition dipoles in the COC groups with other transition dipoles in the molecule. Comparison with experimental spectra leads to an assignment of this band to 11A1 --> 21B1 transitions centered on the COC groups.
Ab initio yield curve dynamics [rapid communication
NASA Astrophysics Data System (ADS)
Hawkins, Raymond J.; Roy Frieden, B.; D'Anna, Joseph L.
2005-09-01
We derive an equation of motion for interest-rate yield curves by applying a minimum Fisher information variational approach to the implied probability density. By construction, solutions to the equation of motion recover observed bond prices. More significantly, the form of the resulting equation explains the success of the Nelson Siegel approach to fitting static yield curves and the empirically observed modal structure of yield curves. A practical numerical implementation of this equation of motion is found by using the Karhunen Lòeve expansion and Galerkin's method to formulate a reduced-order model of yield curve dynamics.
Stress reduction of Cu-doped diamond-like carbon films from ab initio calculations
Li, Xiaowei; Ke, Peiling; Wang, Aiying
2015-01-15
Structure and properties of Cu-doped diamond-like carbon films (DLC) were investigated using ab initio calculations. The effect of Cu concentrations (1.56∼7.81 at.%) on atomic bond structure was mainly analyzed to clarify the residual stress reduction mechanism. Results showed that with introducing Cu into DLC films, the residual compressive stress decreased firstly and then increased for each case with the obvious deterioration of mechanical properties, which was in agreement with the experimental results. Structural analysis revealed that the weak Cu-C bond and the relaxation of both the distorted bond angles and bond lengths accounted for the significant reduction of residual compressive stress, while at the higher Cu concentration the increase of residual stress attributed to the existence of distorted Cu-C structures and the increased fraction of distorted C-C bond lengths.
NASA Astrophysics Data System (ADS)
Benassi, R.; Bertarini, C.; Hilfert, L.; Kempter, G.; Kleinpeter, E.; Spindler, J.; Taddei, F.; Thomas, S.
2000-03-01
The structure of a number of 2- exo-methylene substituted quinazolines and benzodiazepines, respectively, 1, 3a, b, 4( X=-CN, -COOEt ) and their 2-cyanoimino substituted analogues 2, 3c, d( X=-CN, -SO 2C 6H 4-Me (p) was completely assigned by the whole arsenal of 1D and 2D NMR spectroscopic methods. The E/ Z isomerism at the exo-cyclic double bond was determined by both NMR spectroscopy and confirmed by ab initio quantum chemical calculations; the Z isomer is the preferred one, its amount proved dependent on steric hindrance. Due to the push-pull effect in this part of the molecules the restricted rotation about the partial C 2,C 11 and C 2,N 11 double bonds, could also be studied and the barrier to rotation measured by dynamic NMR spectroscopy. The free energies of activation of this dynamic process proved very similar along the compounds studied but being dependent on the polarity of the solvent. Quantum chemical calculations at the ab initio level were employed to prove the stereochemistry at the exo-cyclic partial double bonds of 1- 4, to calculate the barriers to rotation but also to discuss in detail both the ground and the transition state of the latter dynamic process in order to better understand electronic, inter- and intramolecular effects on the barrier to rotation which could be determined experimentally. In the cyanoimino substituted compounds 2, 3c, d, the MO ab initio calculations evidence the isomer interconversion to be better described by the internal rotation process than by the lateral shift mechanism.
The Crystal Structure of Impurity Centers Tm^{2+} and Eu^{2+} in SrCl2 : Ab Initio Calculations
NASA Astrophysics Data System (ADS)
Chernyshev, V. A.; Serdcev, A. V.; Petrov, V. P.; Nikiforov, A. E.
2016-01-01
Ab initio calculations of the impurity centers Tm^{2+} thulium and europium Eu^{2+} in SrCl2 and MeF2 (Me = Ca, Sr, Ba) were carried out at low (zero) temperature. The crystal structure of impurity centers was investigated. Charge density maps show that the bonds formed by the rare-earth ions have an ionic character. The crystal structures, lattice dynamics, and band structures of MeF2 and SrCl2 were calculated at low temperature. Ab initio calculations were performed in periodic CRYSTAL code within the framework of the MO LCAO approach by using hybrid DFT functionals.
NASA Astrophysics Data System (ADS)
Resat, Marianne Sowa; Smolanoff, Jason N.; Goldman, Ilyse B.; Anderson, Scott L.
1994-06-01
We report a combined experimental and theoretical study of the reaction of small carbon cluster cations with N2O aimed at understanding the reaction mechanism and how it is affected by the electronic and geometric structure of the C+n reactants. Cross sections for reaction of C+n (n=3-12) with N2O were measured over a collision energy range from 0.1-10 eV, using a guided ion beam tandem mass spectrometer. Ab initio calculations were used to examine the structure and energetics of reactant and product species. Small clusters, which are linear, react with no activation barrier, resulting in either oxide or nitride formation. The branching between oxide and nitride channels shows a strong even-odd alternation, with even clusters preferentially forming nitrides. This appears to be correlated with an even/odd alternation in the ionization potential of the CnN. The larger, monocyclic C+n have activation barriers for reaction, and a completely different product distribution. Secondary reactions of the primary oxide and nitride products were studied at high N2O pressures. Products containing two O or two N atoms are not observed, but it is possible to add one of each. Possible reaction mechanisms are discussed and supported by thermochemistry derived from spin restricted ab initio calculations.
Rana, Malay Kumar; Chandra, Amalendu
2013-05-28
The behavior of water near a graphene sheet is investigated by means of ab initio and classical molecular dynamics simulations. The wetting of the graphene sheet by ab initio water and the relation of such behavior to the strength of classical dispersion interaction between surface atoms and water are explored. The first principles simulations reveal a layered solvation structure around the graphene sheet with a significant water density in the interfacial region implying no drying or cavitation effect. It is found that the ab initio results of water density at interfaces can be reproduced reasonably well by classical simulations with a tuned dispersion potential between the surface and water molecules. Calculations of vibrational power spectrum from ab initio simulations reveal a shift of the intramolecular stretch modes to higher frequencies for interfacial water molecules when compared with those of the second solvation later or bulk-like water due to the presence of free OH modes near the graphene sheet. Also, a weakening of the water-water hydrogen bonds in the vicinity of the graphene surface is found in our ab initio simulations as reflected in the shift of intermolecular vibrational modes to lower frequencies for interfacial water molecules. The first principles calculations also reveal that the residence and orientational dynamics of interfacial water are somewhat slower than those of the second layer or bulk-like molecules. However, the lateral diffusion and hydrogen bond relaxation of interfacial water molecules are found to occur at a somewhat faster rate than that of the bulk-like water molecules. The classical molecular dynamics simulations with tuned Lennard-Jones surface-water interaction are found to produce dynamical results that are qualitatively similar to those of ab initio molecular dynamics simulations. PMID:23742495
Borisov, Yurii A.; Garrett, Bruce C.; Kobanovskii, Y. A.; Bilera, I. V.; Buravtsev, N. N.
2003-08-07
In this study, we established a correlation between the enthalpy of cleavage of the C=C bond in fluorine-substituted olefins giving rise to two carbenes in the electronic ground state and the distribution of the electron density on this bond.
Macromolecular ab initio phasing enforcing secondary and tertiary structure
Millán, Claudia; Sammito, Massimo; Usón, Isabel
2015-01-01
Ab initio phasing of macromolecular structures, from the native intensities alone with no experimental phase information or previous particular structural knowledge, has been the object of a long quest, limited by two main barriers: structure size and resolution of the data. Current approaches to extend the scope of ab initio phasing include use of the Patterson function, density modification and data extrapolation. The authors’ approach relies on the combination of locating model fragments such as polyalanine α-helices with the program PHASER and density modification with the program SHELXE. Given the difficulties in discriminating correct small substructures, many putative groups of fragments have to be tested in parallel; thus calculations are performed in a grid or supercomputer. The method has been named after the Italian painter Arcimboldo, who used to compose portraits out of fruit and vegetables. With ARCIMBOLDO, most collections of fragments remain a ‘still-life’, but some are correct enough for density modification and main-chain tracing to reveal the protein’s true portrait. Beyond α-helices, other fragments can be exploited in an analogous way: libraries of helices with modelled side chains, β-strands, predictable fragments such as DNA-binding folds or fragments selected from distant homologues up to libraries of small local folds that are used to enforce nonspecific tertiary structure; thus restoring the ab initio nature of the method. Using these methods, a number of unknown macromolecules with a few thousand atoms and resolutions around 2 Å have been solved. In the 2014 release, use of the program has been simplified. The software mediates the use of massive computing to automate the grid access required in difficult cases but may also run on a single multicore workstation (http://chango.ibmb.csic.es/ARCIMBOLDO_LITE) to solve straightforward cases. PMID:25610631
Ab Initio Computation of the Energies of Circular Quantum Dots
Lohne, M. Pedersen; Hagen, Gaute; Hjorth-Jensen, M.; Kvaal, S.; Pederiva, F.
2011-01-01
We perform coupled-cluster and diffusion Monte Carlo calculations of the energies of circular quantum dots up to 20 electrons. The coupled-cluster calculations include triples corrections and a renormalized Coulomb interaction defined for a given number of low-lying oscillator shells. Using such a renormalized Coulomb interaction brings the coupled-cluster calculations with triples correlations in excellent agreement with the diffusion Monte Carlo calculations. This opens up perspectives for doing ab initio calculations for much larger systems of electrons.
The implementation of ab initio quantum chemistry calculations on transporters.
Cooper, M D; Hillier, I H
1991-06-01
The RHF and geometry optimization sections of the ab initio quantum chemistry code, GAMESS, have been optimized for a network of parallel microprocessors, Inmos T800-20 transputers, using both indirect and direct SCF techniques. The results indicate great scope for implementation of such codes on small parallel computer systems, very high efficiencies having been achieved, particularly in the cases of direct SCF and geometry optimization with large basis sets. The work, although performed upon one particular parallel system, the Meiko Computing Surface, is applicable to a wide range of parallel systems with both shared and distributed memory. PMID:1919615
Ab initio electronic properties of dual phosphorus monolayers in silicon
2014-01-01
In the midst of the epitaxial circuitry revolution in silicon technology, we look ahead to the next paradigm shift: effective use of the third dimension - in particular, its combination with epitaxial technology. We perform ab initio calculations of atomically thin epitaxial bilayers in silicon, investigating the fundamental electronic properties of monolayer pairs. Quantitative band splittings and the electronic density are presented, along with effects of the layers’ relative alignment and comments on disordered systems, and for the first time, the effective electronic widths of such device components are calculated. PMID:25246862
Ab initio vibrational and dielectric properties of Y V O
NASA Astrophysics Data System (ADS)
Vali, R.
2009-10-01
For the yttrium orthovanadate Y V O with a tetragonal zircon-type structure, the first complete set of Raman-active and IR-active phonon modes has been calculated using ab initio density functional perturbation theory. The calculated IR reflectivity spectra are in good agreement with available experimental data. We report the calculated frequencies of three Raman-active modes that could not be detected experimentally and a new assignment of the experimental Raman data. The contributions of each IR-active phonon modes to static dielectric tensor have been determined.
Potential Dependence of Electrochemical Barriers from ab Initio Calculations.
Chan, Karen; Nørskov, Jens K
2016-05-01
We present a simple and computationally efficient method to determine the potential dependence of the activation energies for proton-electron transfer from a single ab initio barrier calculation. We show that the potential dependence of the activation energy is given by the partial charge transferred at the transition state. The method is evaluated against the potential dependence determined explicitly through multiple calculations at varying potential. We show that the transfer coefficient is given by the charge transferred from the initial to transition state, which has significant implications for electrochemical kinetics. PMID:27088442
Ab-Initio Shell Model with a Core
Lisetskiy, A F; Barrett, B R; Kruse, M; Navratil, P; Stetcu, I; Vary, J P
2008-06-04
We construct effective 2- and 3-body Hamiltonians for the p-shell by performing 12{h_bar}{Omega} ab initio no-core shell model (NCSM) calculations for A=6 and 7 nuclei and explicitly projecting the many-body Hamiltonians onto the 0{h_bar}{Omega} space. We then separate these effective Hamiltonians into 0-, 1- and 2-body contributions (also 3-body for A=7) and analyze the systematic behavior of these different parts as a function of the mass number A and size of the NCSM basis space. The role of effective 3- and higher-body interactions for A > 6 is investigated and discussed.
Ab initio Study of He Stability in hcp-Ti
Dai, Yunya; Yang, Li; Peng, SM; Long, XG; Gao, Fei; Zu, Xiaotao T.
2010-12-20
The stability of He in hcp-Ti was studied using ab initio method based on density functional theory. The results indicate that a single He atom prefers to occupy the tetrahedral site rather than the octahedral site. The interaction of He defects with Ti atoms has been used to explain the relative stabilities of He point defects in hcp-Ti. The relative stability of He defects in hcp-Ti is useful for He clustering and bubble nucleation in metal tritides, which provides the basis for development of improved atomistic models.
Ab Initio Calculations Applied to Problems in Metal Ion Chemistry
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.; Partridge, Harry; Arnold, James O. (Technical Monitor)
1994-01-01
Electronic structure calculations can provide accurate spectroscopic data (such as molecular structures) vibrational frequencies, binding energies, etc.) that have been very useful in explaining trends in experimental data and in identifying incorrect experimental measurements. In addition, ab initio calculations. have given considerable insight into the many interactions that make the chemistry of transition metal systems so diverse. In this review we focus on cases where calculations and experiment have been used to solve interesting chemical problems involving metal ions. The examples include cases where theory was used to differentiate between disparate experimental values and cases where theory was used to explain unexpected experimental results.
Ab initio study of hydrogen on beryllium surfaces
NASA Astrophysics Data System (ADS)
Bachurin, D. V.; Vladimirov, P. V.
2015-11-01
Static ab initio calculations were performed for five principal hexagonal close-packed beryllium surfaces: basal, prismatic (type I and II) and pyramidal (type I and II). The basal plane was found to be the most energetically favorable, while the energies of the prismatic (type I) and pyramidal (type I) planes were slightly higher followed by the type II planes. Beryllium is known to show extreme interlayer distance relaxation near the surface. Up to five outermost atomic layers were involved in surface relaxation. The presence of hydrogen on the beryllium surfaces led to a noticeable reduction of the surface energy.
Ab initio study of neutron drops with chiral Hamiltonians
NASA Astrophysics Data System (ADS)
Potter, H. D.; Fischer, S.; Maris, P.; Vary, J. P.; Binder, S.; Calci, A.; Langhammer, J.; Roth, R.
2014-12-01
We report ab initio calculations for neutron drops in a 10 MeV external harmonic-oscillator trap using chiral nucleon-nucleon plus three-nucleon interactions. We present total binding energies, internal energies, radii and odd-even energy differences for neutron numbers N = 2- 18 using the no-core shell model with and without importance truncation. Furthermore, we present total binding energies for N = 8 , 16 , 20 , 28 , 40 , 50 obtained in a coupled-cluster approach. Comparisons with quantum Monte Carlo results, where available, using Argonne v8‧ with three-nucleon interactions reveal important dependences on the chosen Hamiltonian.
Ab-initio study of transition metal hydrides
Sharma, Ramesh; Shukla, Seema Dwivedi, Shalini Sharma, Yamini
2014-04-24
We have performed ab initio self consistent calculations based on Full potential linearized augmented plane wave (FP-LAPW) method to investigate the optical and thermal properties of yttrium hydrides. From the band structure and density of states, the optical absorption spectra and specific heats have been calculated. The band structure of Yttrium metal changes dramatically due to hybridization of Y sp orbitals with H s orbitals and there is a net charge transfer from metal to hydrogen site. The electrical resistivity and specific heats of yttrium hydrides are lowered but the thermal conductivity is slightly enhanced due to increase in scattering from hydrogen sites.
Debela, T. T.; Wang, X. D.; Cao, Q. P.; Zhang, D. X.; Wang, S. Y.; Wang, Cai-Zhuang; Jiang, J. Z.
2013-12-12
Atomic structure transitions of liquid niobium during solidification, at different temperatures from 3200 to 1500 K, were studied by using ab initio molecular dynamics simulations. The local atomic structure variations with temperature are investigated by using the pair-correlation function, the structure factor, the bond-angle distribution function, the Honeycutt–Anderson index, Voronoi tessellation and the cluster alignment methods. Our results clearly show that, upon quenching, the icosahedral short-range order dominates in the stable liquid and supercooled liquid states before the system transforms to crystalline body-center cubic phase at a temperature of about 1830 K.
Ab initio study on electronically excited states of lithium isocyanide, LiNC
NASA Astrophysics Data System (ADS)
Yasumatsu, Hisato; Jeung, Gwang-Hi
2014-01-01
The electronically excited states of the lithium isocyanide molecule, LiNC, were studied by means of ab initio calculations. The bonding nature of LiNC up to ˜10 eV is discussed on the basis of the potential energy surfaces according to the interaction between the ion-pair and covalent states. The ion-pair states are described by Coulomb attractive interaction in the long distance range, while the covalent ones are almost repulsive or bound with a very shallow potential dent. These two states interact each other to form adiabatic potential energy surfaces with non-monotonic change in the potential energy with the internuclear distance.
NASA Astrophysics Data System (ADS)
Pietrucci, Fabio; Andreoni, Wanda
2011-08-01
Social permutation invariant coordinates are introduced describing the bond network around a given atom. They originate from the largest eigenvalue and the corresponding eigenvector of the contact matrix, are invariant under permutation of identical atoms, and bear a clear signature of an order-disorder transition. Once combined with ab initio metadynamics, these coordinates are shown to be a powerful tool for the discovery of low-energy isomers of molecules and nanoclusters as well as for a blind exploration of isomerization, association, and dissociation reactions.
The role of Metals in Amyloid Aggregation: A Test Case for ab initio Simulations
Minicozzi, V.; Rossi, G. C.; Stellato, F.; Morante, S.
2007-12-26
First principle ab initio molecular dynamics simulations of the Car-Parrinello type have proved to be of invaluable help in understanding the microscopic mechanisms of chemical bonding both in solid state physics and in structural biophysics. In this work we present as test cases the study of the Cu coordination mode in two especially important examples: Prion protein and {beta}-amyloids. Using medium size PC-clusters as well as larger parallel platforms, we are able to deal with systems comprising 300 to 500 atoms and 1000 to 1500 electrons for as long as 2-3 ps. We present structural results which confirm indications coming from NMR and XAS data.
Enhancing mechanical toughness of aluminum surfaces by nano-boron implantation: An ab initio study
NASA Astrophysics Data System (ADS)
Zhu, Zhen; Kwon, Dae-Gyeon; Kwon, Young-Kyun; Tománek, David
2015-01-01
Searching for ways to enhance surface hardness of aluminum, we study the equilibrium structure, stability, elastic properties and formation dynamics of a boron-enriched surface using ab initio density functional calculations. We used molecular dynamics simulations to model the implantation of energetic boron nanoparticles in Al and identify structural arrangements that optimize the formation of strong covalent Bsbnd Al bonds. Nano-indentation simulations based on constrained optimization suggest that presence of boron nanostructures in the subsurface region enhances significantly the mechanical hardness of aluminum surfaces.
NASA Astrophysics Data System (ADS)
Garrone, E.; Ugliengo, P.; Ghiotti, G.; Borello, E.; Saunders, V. R.
1993-08-01
High-quality ab initio computational results for the interaction of N 2O with the silanol molecule H 3SiOH, taken as representative of the isolated hydroxyl at the silica surface, are compared with IR data concerning the real system measured at 195 K. Besides a non-specific interaction, features typical of hydrogen bonding are observed: two types of complexes are formed simultaneously, differing in their 2 v1 mode, probably one nitrogen-bonded, the other oxygen-bonded. Ab initio results are in excellent agreement as far as binding energies are concerned: the two complexes show comparable stabilities, once electron correlation has been taken into account; however, SCF-calculated IR features do not allow a definite vibrational assignment of the two complexes.
Ab initio study of hot electrons in GaAs.
Bernardi, Marco; Vigil-Fowler, Derek; Ong, Chin Shen; Neaton, Jeffrey B; Louie, Steven G
2015-04-28
Hot carrier dynamics critically impacts the performance of electronic, optoelectronic, photovoltaic, and plasmonic devices. Hot carriers lose energy over nanometer lengths and picosecond timescales and thus are challenging to study experimentally, whereas calculations of hot carrier dynamics are cumbersome and dominated by empirical approaches. In this work, we present ab initio calculations of hot electrons in gallium arsenide (GaAs) using density functional theory and many-body perturbation theory. Our computed electron-phonon relaxation times at the onset of the Γ, L, and X valleys are in excellent agreement with ultrafast optical experiments and show that the ultrafast (tens of femtoseconds) hot electron decay times observed experimentally arise from electron-phonon scattering. This result is an important advance to resolve a controversy on hot electron cooling in GaAs. We further find that, contrary to common notions, all optical and acoustic modes contribute substantially to electron-phonon scattering, with a dominant contribution from transverse acoustic modes. This work provides definitive microscopic insight into hot electrons in GaAs and enables accurate ab initio computation of hot carriers in advanced materials. PMID:25870287
Ab initio dynamics of the cytochrome P450 hydroxylation reaction
Elenewski, Justin E.; Hackett, John C
2015-01-01
The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis. PMID:25681906
Three-cluster dynamics within an ab initio framework
Quaglioni, Sofia; Romero-Redondo, Carolina; Navratil, Petr
2013-09-26
In this study, we introduce a fully antisymmetrized treatment of three-cluster dynamics within the ab initio framework of the no-core shell model/resonating-group method. Energy-independent nonlocal interactions among the three nuclear fragments are obtained from realistic nucleon-nucleon interactions and consistent ab initio many-body wave functions of the clusters. The three-cluster Schrödinger equation is solved with bound-state boundary conditions by means of the hyperspherical-harmonic method on a Lagrange mesh. We discuss the formalism in detail and give algebraic expressions for systems of two single nucleons plus a nucleus. Using a soft similarity-renormalization-group evolved chiral nucleon-nucleon potential, we apply the method to amore » 4He+n+n description of 6He and compare the results to experiment and to a six-body diagonalization of the Hamiltonian performed within the harmonic-oscillator expansions of the no-core shell model. Differences between the two calculations provide a measure of core (4He) polarization effects.« less
Oxidation of GaN: An ab initio thermodynamic approach
NASA Astrophysics Data System (ADS)
Jackson, Adam J.; Walsh, Aron
2013-10-01
GaN is a wide-band-gap semiconductor used in high-efficiency light-emitting diodes and solar cells. The solid is produced industrially at high chemical purities by deposition from a vapor phase, and oxygen may be included at this stage. Oxidation represents a potential path for tuning its properties without introducing more exotic elements or extreme processing conditions. In this work, ab initio computational methods are used to examine the energy potentials and electronic properties of different extents of oxidation in GaN. Solid-state vibrational properties of Ga, GaN, Ga2O3, and a single substitutional oxygen defect have been studied using the harmonic approximation with supercells. A thermodynamic model is outlined which combines the results of ab initio calculations with data from experimental literature. This model allows free energies to be predicted for arbitrary reaction conditions within a wide process envelope. It is shown that complete oxidation is favorable for all industrially relevant conditions, while the formation of defects can be opposed by the use of high temperatures and a high N2:O2 ratio.
Ab initio thermodynamic model for magnesium carbonates and hydrates.
Chaka, Anne M; Felmy, Andrew R
2014-09-01
An ab initio thermodynamic framework for predicting properties of hydrated magnesium carbonate minerals has been developed using density-functional theory linked to macroscopic thermodynamics through the experimental chemical potentials for MgO, water, and CO2. Including semiempirical dispersion via the Grimme method and small corrections to the generalized gradient approximation of Perdew, Burke, and Ernzerhof for the heat of formation yields a model with quantitative agreement for the benchmark minerals brucite, magnesite, nesquehonite, and hydromagnesite. The model shows how small differences in experimental conditions determine whether nesquehonite, hydromagnesite, or magnesite is the result of laboratory synthesis from carbonation of brucite, and what transformations are expected to occur on geological time scales. Because of the reliance on parameter-free first-principles methods, the model is reliably extensible to experimental conditions not readily accessible to experiment and to any mineral composition for which the structure is known or can be hypothesized, including structures containing defects, substitutions, or transitional structures during solid state transformations induced by temperature changes or processes such as water, CO2, or O2 diffusion. Demonstrated applications of the ab initio thermodynamic framework include an independent means to evaluate differences in thermodynamic data for lansfordite, predicting the properties of Mg analogues of Ca-based hydrated carbonates monohydrocalcite and ikaite, which have not been observed in nature, and an estimation of the thermodynamics of barringtonite from the stoichiometry and a single experimental observation. PMID:24679248
Ab Initio Thermodynamic Model for Magnesium Carbonates and Hydrates
Chaka, Anne M.; Felmy, Andrew R.
2014-03-28
An ab initio thermodynamic framework for predicting properties of hydrated magnesium carbonate minerals has been developed using density-functional theory linked to macroscopic thermodynamics through the experimental chemical potentials for MgO, water, and CO2. Including semiempirical dispersion via the Grimme method and small corrections to the generalized gradient approximation of Perdew, Burke, and Ernzerhof for the heat of formation yields a model with quantitative agreement for the benchmark minerals brucite, magnesite, nesquehonite, and hydromagnesite. The model shows how small differences in experimental conditions determine whether nesquehonite, hydromagnesite, or magnesite is the result of laboratory synthesis from carbonation of brucite, and what transformations are expected to occur on geological time scales. Because of the reliance on parameter-free first principles methods, the model is reliably extensible to experimental conditions not readily accessible to experiment and to any mineral composition for which the structure is known or can be hypothesized, including structures containing defects, substitutions, or transitional structures during solid state transformations induced by temperature changes or processes such as water, CO2, or O2 diffusion. Demonstrated applications of the ab initio thermodynamic framework include an independent means to evaluate differences in thermodynamic data for lansfordite, predicting the properties of Mg analogs of Ca-based hydrated carbonates monohydrocalcite and ikaite which have not been observed in nature, and an estimation of the thermodynamics of barringtonite from the stoichiometry and a single experimental observation.
Ab initio prediction of the critical thickness of a precipitate
NASA Astrophysics Data System (ADS)
Sampath, S.; Janisch, R.
2013-09-01
Segregation and precipitation of second phases in metals and metallic alloys is an important phenomenon that has a strong influence on the mechanical properties of the material. Models exist that describe the growth of coherent, semi-coherent and incoherent precipitates. One important parameter of these models is the energy of the interface between matrix and precipitate. In this work we apply ab initio density functional theory calculations to obtain this parameter and to understand how it depends on chemical composition and mechanical strain at the interface. Our example is a metastable Mo-C phase, the body-centred tetragonal structure, which exists as a semi-coherent precipitate in body-centred cubic molybdenum. The interface of this precipitate is supposed to change from coherent to semi-coherent during the growth of the precipitate. We predict the critical thickness of the precipitate by calculating the different contributions to a semi-coherent interface energy by means of ab initio density functional theory calculations. The parameters in our model include the elastic strain energy stored in the precipitate, as well as a misfit dislocation energy that depends on the dislocation core width and the dislocation spacing. Our predicted critical thickness agrees well with experimental observations.
Ab initio dynamics of the cytochrome P450 hydroxylation reaction
NASA Astrophysics Data System (ADS)
Elenewski, Justin E.; Hackett, John C.
2015-02-01
The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis.
Ab initio dynamics of the cytochrome P450 hydroxylation reaction
Elenewski, Justin E.; Hackett, John C
2015-02-14
The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis.
AB INITIO SIMULATIONS FOR MATERIAL PROPERTIES ALONG THE JUPITER ADIABAT
French, Martin; Becker, Andreas; Lorenzen, Winfried; Nettelmann, Nadine; Bethkenhagen, Mandy; Redmer, Ronald; Wicht, Johannes
2012-09-15
We determine basic thermodynamic and transport properties of hydrogen-helium-water mixtures for the extreme conditions along Jupiter's adiabat via ab initio simulations, which are compiled in an accurate and consistent data set. In particular, we calculate the electrical and thermal conductivity, the shear and longitudinal viscosity, and diffusion coefficients of the nuclei. We present results for associated quantities like the magnetic and thermal diffusivity and the kinematic shear viscosity along an adiabat that is taken from a state-of-the-art interior structure model. Furthermore, the heat capacities, the thermal expansion coefficient, the isothermal compressibility, the Grueneisen parameter, and the speed of sound are calculated. We find that the onset of dissociation and ionization of hydrogen at about 0.9 Jupiter radii marks a region where the material properties change drastically. In the deep interior, where the electrons are degenerate, many of the material properties remain relatively constant. Our ab initio data will serve as a robust foundation for applications that require accurate knowledge of the material properties in Jupiter's interior, e.g., models for the dynamo generation.
Entropy of Liquid Water from Ab Initio Molecular Dynamics
NASA Astrophysics Data System (ADS)
Spanu, Leonardo; Zhang, Cui; Galli, Giulia
2012-02-01
The debate on the structural properties of water has been mostly based on the calculation of pair correlation functions. However, the simulation of thermodynamic and spectroscopic quantities may be of great relevance for the characterization of liquid water properties. We have computed the entropy of liquid water using a two-phase thermodynamic model and trajectories generated by ab initio molecular dynamics simulations [1]. In an attempt to better understand the performance of several density functionals in simulating liquid water, we have performed ab initio molecular dynamics using semilocal, hybrid [2] and van der Waals density functionals [3]. We show that in all cases, at the experimental equilibrium density and at temperatures in the vicinity of 300 K, the computed entropies are underestimated, with respect to experiment, and the liquid exhibits a degree of tetrahedral order higher than in experiments. We also discuss computational strategies to simulate spectroscopic properties of water, including infrared and Raman spectra.[4pt] [1] C.Zhang, L.Spanu and G.Galli, J.Phys.Chem. B 2011 (in press)[0pt] [2] C.Zhang, D.Donadio, F.Gygi and G.Galli, J. Chem. Theory Comput. 7, 1443 (2011)[0pt] [3] C.Zhang, J.Wu, G.Galli and F.Gygi, J. Chem. Theory Comput. 7, 3061 (2011)
Unified ab initio approaches to nuclear structure and reactions
NASA Astrophysics Data System (ADS)
Navrátil, Petr; Quaglioni, Sofia; Hupin, Guillaume; Romero-Redondo, Carolina; Calci, Angelo
2016-05-01
The description of nuclei starting from the constituent nucleons and the realistic interactions among them has been a long-standing goal in nuclear physics. In addition to the complex nature of the nuclear forces, with two-, three- and possibly higher many-nucleon components, one faces the quantum-mechanical many-nucleon problem governed by an interplay between bound and continuum states. In recent years, significant progress has been made in ab initio nuclear structure and reaction calculations based on input from QCD-employing Hamiltonians constructed within chiral effective field theory. After a brief overview of the field, we focus on ab initio many-body approaches—built upon the no-core shell model—that are capable of simultaneously describing both bound and scattering nuclear states, and present results for resonances in light nuclei, reactions important for astrophysics and fusion research. In particular, we review recent calculations of resonances in the 6He halo nucleus, of five- and six-nucleon scattering, and an investigation of the role of chiral three-nucleon interactions in the structure of 9Be. Further, we discuss applications to the 7Be {({{p}},γ )}8{{B}} radiative capture. Finally, we highlight our efforts to describe transfer reactions including the 3H{({{d}},{{n}})}4He fusion.
Ab initio tight-binding Hamiltonian for transition metal dichalcogenides
NASA Astrophysics Data System (ADS)
Fang, Shiang; Kuate Defo, Rodrick; Shirodkar, Sharmila N.; Lieu, Simon; Tritsaris, Georgios A.; Kaxiras, Efthimios
2015-11-01
We present an accurate ab initio tight-binding Hamiltonian for the transition metal dichalcogenides, MoS2, MoSe2, WS2, WSe2, with a minimal basis (the d orbitals for the metal atoms and p orbitals for the chalcogen atoms) based on a transformation of the Kohn-Sham density functional theory Hamiltonian to a basis of maximally localized Wannier functions. The truncated tight-binding Hamiltonian, with only on-site, first, and partial second neighbor interactions, including spin-orbit coupling, provides a simple physical picture and the symmetry of the main band-structure features. Interlayer interactions between adjacent layers are modeled by transferable hopping terms between the chalcogen p orbitals. The full-range tight-binding Hamiltonian can be reduced to hybrid-orbital k .p effective Hamiltonians near the band extrema that capture important low-energy excitations. These ab initio Hamiltonians can serve as the starting point for applications to interacting many-body physics including optical transitions and Berry curvature of bands, of which we give some examples.
Decarboxylation of furfural on Pd(111): Ab initio molecular dynamics simulations
NASA Astrophysics Data System (ADS)
Xue, Wenhua; Dang, Hongli; Shields, Darwin; Liu, Yingdi; Jentoft, Friederike; Resasco, Daniel; Wang, Sanwu
2013-03-01
Furfural conversion over metal catalysts plays an important role in the studies of biomass-derived feedstocks. We report ab initio molecular dynamics simulations for the decarboxylation process of furfural on the palladium surface at finite temperatures. We observed and analyzed the atomic-scale dynamics of furfural on the Pd(111) surface and the fluctuations of the bondlengths between the atoms in furfural. We found that the dominant bonding structure is the parallel structure in which the furfural plane, while slightly distorted, is parallel to the Pd surface. Analysis of the bondlength fluctuations indicates that the C-H bond is the aldehyde group of a furfural molecule is likely to be broken first, while the C =O bond has a tendency to be isolated as CO. Our results show that the reaction of decarbonylation dominates, consistent with the experimental measurements. Supported by DOE (DE-SC0004600). Simulations and calculations were performed on XSEDE's and NERSC's supercomputers.
{ital Ab initio} study of silicon in the R8 phase
Pfrommer, B.G.; Cote, M.; Louie, S.G.; Cohen, M.L.
1997-09-01
We present a detailed {ital ab initio} study of the electronic and structural properties of the recently discovered R8 phase of silicon. Within the framework of density-functional theory in the local-density approximation and using pseudopotentials with a plane-wave basis, we study the energetics of the R8 phase compared to the other tetrahedrally bonded diamond, {beta}-Sn, and BC8 phases. The bonding properties and the pressure dependence of the bond lengths of the BC8 and R8 phases are investigated. An analysis of the band structure reveals that R8 silicon could be a semimetal or semiconductor with a small, indirect band gap. Our computed density of states of R8 silicon shows a sharpening of the valence-band edge similar to the one observed for amorphous silicon. {copyright} {ital 1997} {ital The American Physical Society}
Development of Novel Analytical Method for Ab Initio Powder Structural Analysis
NASA Astrophysics Data System (ADS)
Sakata, Makoto; Nishibori, Eiji; Sawa, Hiroshi
Genetic Algorithm (GA) applied to ab initio structure determination from synchrotron powder diffraction is described. It seems to have an advantage over other real space methods for ab initio structure determination because of the existence of schema theorem. As an example, the case of Prednisolone Succinate is shown in some detail. Future development of GA in crystallography is briefly described.
Ab Initio Potential Energy Surfaces and the Calculation of Accurate Vibrational Frequencies
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Dateo, Christopher E.; Martin, Jan M. L.; Taylor, Peter R.; Langhoff, Stephen R. (Technical Monitor)
1995-01-01
Due to advances in quantum mechanical methods over the last few years, it is now possible to determine ab initio potential energy surfaces in which fundamental vibrational frequencies are accurate to within plus or minus 8 cm(exp -1) on average, and molecular bond distances are accurate to within plus or minus 0.001-0.003 Angstroms, depending on the nature of the bond. That is, the potential energy surfaces have not been scaled or empirically adjusted in any way, showing that theoretical methods have progressed to the point of being useful in analyzing spectra that are not from a tightly controlled laboratory environment, such as vibrational spectra from the interstellar medium. Some recent examples demonstrating this accuracy will be presented and discussed. These include the HNO, CH4, C2H4, and ClCN molecules. The HNO molecule is interesting due to the very large H-N anharmonicity, while ClCN has a very large Fermi resonance. The ab initio studies for the CH4 and C2H4 molecules present the first accurate full quartic force fields of any kind (i.e., whether theoretical or empirical) for a five-atom and six-atom system, respectively.
An Ab Initio Approach Towards Engineering Fischer-Tropsch Surface Chemistry
Matthew Neurock
2005-06-13
As petroleum prices continue to rise and the United States seeks to reduce its dependency on foreign oil, there is a renewed interest in the research and development of more efficient and alternative energy sources, such as fuel cells. One approach is to utilize processes that can produce long-chain hydrocarbons from other sources. One such reaction is Fischer-Tropsch synthesis. Fischer-Tropsch synthesis is a process by which syngas (CO and H{sub 2}) is converted to higher molecular weight hydrocarbons. The reaction involves a complex set of bond-breaking and bond-making reactions, such as CO and H{sub 2} activation, hydrocarbon hydrogenation reactions, and hydrocarbon coupling reactions. This report details our initial construction of an ab initio based kinetic Monte Carlo code that can be used to begin to simulate Fischer-Tropsch synthesis over model Co(0001) surfaces. The code is based on a stochastic kinetic formalism that allows us to explicitly track the transformation of all reactants, intermediates and products. The intrinsic kinetics for the simulations were derived from the ab initio results that we reported in previous year summaries.
NASA Astrophysics Data System (ADS)
Jensen, Per; Buenker, Robert J.; Hirsch, Gerhard; Rai, Sachchida N.
We have calculated ab initio the three-dimensional potential-energy surface of the NH2 molecule at 145 nuclear geometries spanning energy ranges of about 18 000 cm-1 for the NH stretch and 12 000 cm-1 for the bend. The ab initio configuration-interaction calculations were done using the multireference MRD-CI method. The calculated equilibrium configuration has NH bond length re = 1·0207 Å and bond angle α = 103·1°. The rotational-vibrational energies for 14NH2, 14NHD and 14ND2 were calculated variationally using the Morse-oscillator rigid-bender internal-dynamics Hamiltonian. For 14NH2 we calculate that υ1 = 3267 (3219) cm-1, υ2 = 1462 (1497) cm-1 and υ3 = 3283 (3301) cm-1, where experimental values are given in parentheses.
NASA Astrophysics Data System (ADS)
Philipp, Dean Michael
Methodology is discussed for mixed ab initio quantum mechanics/molecular mechanics modeling of systems where the quantum mechanics (QM) and molecular mechanics (MM) regions are within the same molecule. The ab initio QM calculations are at the restricted Hartree-Fock level using the pseudospectral method of the Jaguar program while the MM part is treated with the OPLS force fields implemented in the IMPACT program. The interface between the QM and MM regions, in particular, is elaborated upon, as it is dealt with by ``breaking'' bonds at the boundaries and using Boys-localized orbitals found from model molecules in place of the bonds. These orbitals are kept frozen during QM calculations. The mixed modeling presented here can be used for single point energy calculations and geometry optimizations. Results from tests of the method to find relative conformational energies and geometries of alanine tetrapeptides are presented along with comparisons to pure QM and pure MM calculations.
The ab-initio density matrix renormalization group in practice
Olivares-Amaya, Roberto; Hu, Weifeng; Sharma, Sandeep; Yang, Jun; Chan, Garnet Kin-Lic; Nakatani, Naoki
2015-01-21
The ab-initio density matrix renormalization group (DMRG) is a tool that can be applied to a wide variety of interesting problems in quantum chemistry. Here, we examine the density matrix renormalization group from the vantage point of the quantum chemistry user. What kinds of problems is the DMRG well-suited to? What are the largest systems that can be treated at practical cost? What sort of accuracies can be obtained, and how do we reason about the computational difficulty in different molecules? By examining a diverse benchmark set of molecules: π-electron systems, benchmark main-group and transition metal dimers, and the Mn-oxo-salen and Fe-porphine organometallic compounds, we provide some answers to these questions, and show how the density matrix renormalization group is used in practice.
Ab initio engineering of materials with stacked hexagonal tin frameworks.
Shao, Junping; Beaufils, Clément; Kolmogorov, Aleksey N
2016-01-01
The group-IV tin has been hypothesized to possess intriguing electronic properties in an atom-thick hexagonal form. An attractive pathway of producing sizable 2D crystallites of tin is based on deintercalation of bulk compounds with suitable tin frameworks. Here, we have identified a new synthesizable metal distannide, NaSn2, with a 3D stacking of flat hexagonal layers and examined a known compound, BaSn2, with buckled hexagonal layers. Our ab initio results illustrate that despite being an exception to the 8-electron rule, NaSn2 should form under pressures easily achievable in multi-anvil cells and remain (meta)stable under ambient conditions. Based on calculated Z2 invariants, the predicted NaSn2 may display topologically non-trivial behavior and the known BaSn2 could be a strong topological insulator. PMID:27387140
Ab Initio Calculations of Excited Carrier Dynamics in Gallium Nitride
NASA Astrophysics Data System (ADS)
Jhalani, Vatsal; Bernardi, Marco
Bulk wurtzite GaN is the primary material for blue light-emission technology. The radiative processes in GaN are regulated by the dynamics of excited (or so-called ``hot'') carriers, through microscopic processes not yet completely understood. We present ab initio calculations of electron-phonon (e-ph) scattering rates for hot carriers in GaN. Our work combines density functional theory to compute the electronic states, and density functional perturbation theory to obtain the phonon dispersions and e-ph coupling matrix elements. These quantities are interpolated on fine Brillouin zone grids with maximally localized Wannier functions, to converge the e-ph scattering rates within 5 eV of the band edges. We resolve the contribution of the different phonon modes to the total scattering rate, and study the impact on the relaxation times of the long-range Fröhlich interaction due to the longitudinal-optical phonon modes.
Highly anisotropic thermal conductivity of arsenene: An ab initio study
NASA Astrophysics Data System (ADS)
Zeraati, Majid; Vaez Allaei, S. Mehdi; Abdolhosseini Sarsari, I.; Pourfath, Mahdi; Donadio, Davide
2016-02-01
Elemental two-dimensional (2D) materials exhibit intriguing heat transport and phononic properties. Here we have investigated the lattice thermal conductivity of newly proposed arsenene, the 2D honeycomb structure of arsenic, using ab initio calculations. Solving the Boltzmann transport equation for phonons, we predict a highly anisotropic thermal conductivity of 30.4 and 7.8 W/mK along the zigzag and armchair directions, respectively, at room temperature. Our calculations reveal that phonons with mean free paths between 20 nm and 1 μ m provide the main contribution to the large thermal conductivity in the zigzag direction; mean free paths of phonons contributing to heat transport in the armchair directions range between 20 and 100 nm. The obtained anisotropic thermal conductivity and feasibility of synthesis, in addition to high electron mobility reported elsewhere, make arsenene a promising material for nanoelectronic applications and thermal management.
Ab initio study of II-(VI)2 dichalcogenides.
Olsson, P; Vidal, J; Lincot, D
2011-10-12
The structural stabilities of the (Zn,Cd)(S,Se,Te)(2) dichalcogenides have been determined ab initio. These compounds are shown to be stable in the pyrite phase, in agreement with available experiments. Structural parameters for the ZnTe(2) pyrite semiconductor compound proposed here are presented. The opto-electronic properties of these dichalcogenide compounds have been calculated using quasiparticle GW theory. Bandgaps, band structures and effective masses are proposed as well as absorption coefficients and refraction indices. The compounds are all indirect semiconductors with very flat conduction band dispersion and high absorption coefficients. The work functions and surface properties are predicted. The Te and Se based compounds could be of interest as absorber materials in photovoltaic applications. PMID:21937783
Ab initio study of guanine damage by hydroxyl radical.
Chaban, Galina M; Wang, Dunyou; Huo, Winifred M
2015-01-15
Multiconfigurational ab initio methods are used in this study to examine two initial reactions that take place during the OH radical attack of the DNA base guanine: a ring opening reaction and a hydrogen transfer reaction. The same reactions are also studied in the presence of a single water molecule. The ring opening reaction has a moderate barrier height of ∼20-25 kcal/mol that is relatively insensitive to the presence of water. The barrier of the H-transfer reaction, on the other hand, is lowered from ∼50 to ∼22 kcal/mol when one water molecule is added, thus becoming comparable to the barrier height of the ring opening reaction. PMID:25517252
Ab initio water pair potential with flexible monomers.
Jankowski, Piotr; Murdachaew, Garold; Bukowski, Robert; Akin-Ojo, Omololu; Leforestier, Claude; Szalewicz, Krzysztof
2015-03-26
A potential energy surface for the water dimer with explicit dependence on monomer coordinates is presented. The surface was fitted to a set of previously published interaction energies computed on a grid of over a quarter million points in the 12-dimensional configurational space using symmetry-adapted perturbation theory and coupled-cluster methods. The present fit removes small errors in published fits, and its accuracy is critically evaluated. The minimum and saddle-point structures of the potential surface were found to be very close to predictions from direct ab initio optimizations. The computed second virial coefficients agreed well with experimental values. At low temperatures, the effects of monomer flexibility in the virial coefficients were found to be much smaller than the quantum effects. PMID:25687650
Reactive Monte Carlo sampling with an ab initio potential
NASA Astrophysics Data System (ADS)
Leiding, Jeff; Coe, Joshua D.
2016-05-01
We present the first application of reactive Monte Carlo in a first-principles context. The algorithm samples in a modified NVT ensemble in which the volume, temperature, and total number of atoms of a given type are held fixed, but molecular composition is allowed to evolve through stochastic variation of chemical connectivity. We discuss general features of the method, as well as techniques needed to enhance the efficiency of Boltzmann sampling. Finally, we compare the results of simulation of NH3 to those of ab initio molecular dynamics (AIMD). We find that there are regions of state space for which RxMC sampling is much more efficient than AIMD due to the "rare-event" character of chemical reactions.
Efficient Ab initio Modeling of Random Multicomponent Alloys
NASA Astrophysics Data System (ADS)
Jiang, Chao; Uberuaga, Blas P.
2016-03-01
We present in this Letter a novel small set of ordered structures (SSOS) method that allows extremely efficient ab initio modeling of random multicomponent alloys. Using inverse II-III spinel oxides and equiatomic quinary bcc (so-called high entropy) alloys as examples, we demonstrate that a SSOS can achieve the same accuracy as a large supercell or a well-converged cluster expansion, but with significantly reduced computational cost. In particular, because of this efficiency, a large number of quinary alloy compositions can be quickly screened, leading to the identification of several new possible high-entropy alloy chemistries. The SSOS method developed here can be broadly useful for the rapid computational design of multicomponent materials, especially those with a large number of alloying elements, a challenging problem for other approaches.
High-throughput ab-initio dilute solute diffusion database
Wu, Henry; Mayeshiba, Tam; Morgan, Dane
2016-01-01
We demonstrate automated generation of diffusion databases from high-throughput density functional theory (DFT) calculations. A total of more than 230 dilute solute diffusion systems in Mg, Al, Cu, Ni, Pd, and Pt host lattices have been determined using multi-frequency diffusion models. We apply a correction method for solute diffusion in alloys using experimental and simulated values of host self-diffusivity. We find good agreement with experimental solute diffusion data, obtaining a weighted activation barrier RMS error of 0.176 eV when excluding magnetic solutes in non-magnetic alloys. The compiled database is the largest collection of consistently calculated ab-initio solute diffusion data in the world. PMID:27434308
Ab initio engineering of materials with stacked hexagonal tin frameworks
Shao, Junping; Beaufils, Clément; Kolmogorov, Aleksey N.
2016-01-01
The group-IV tin has been hypothesized to possess intriguing electronic properties in an atom-thick hexagonal form. An attractive pathway of producing sizable 2D crystallites of tin is based on deintercalation of bulk compounds with suitable tin frameworks. Here, we have identified a new synthesizable metal distannide, NaSn2, with a 3D stacking of flat hexagonal layers and examined a known compound, BaSn2, with buckled hexagonal layers. Our ab initio results illustrate that despite being an exception to the 8-electron rule, NaSn2 should form under pressures easily achievable in multi-anvil cells and remain (meta)stable under ambient conditions. Based on calculated Z2 invariants, the predicted NaSn2 may display topologically non-trivial behavior and the known BaSn2 could be a strong topological insulator. PMID:27387140
Ab initio calculation of the shock Hugoniot of bulk silicon
NASA Astrophysics Data System (ADS)
Strickson, Oliver; Artacho, Emilio
2016-03-01
We describe how ab initio molecular dynamics can be used to determine the Hugoniot locus (states accessible by a shock wave) for materials with a number of stable phases, and with an approximate treatment of plasticity and yield, without having to simulate these phenomena directly. We consider the case of bulk silicon, with forces from density-functional theory, up to 70 GPa. The fact that shock waves can split into multiple waves due to phase transitions or yielding is taken into account here by specifying the strength of any preceding waves explicitly based on their yield strain. Points corresponding to uniaxial elastic compression along three crystal axes and a number of postshock phases are given, including a plastically yielded state, approximated by an isotropic stress configuration following an elastic wave of predetermined strength. The results compare well to existing experimental data for shocked silicon.
Interatomic Coulombic decay widths of helium trimer: Ab initio calculations
Kolorenč, Přemysl; Sisourat, Nicolas
2015-12-14
We report on an extensive study of interatomic Coulombic decay (ICD) widths in helium trimer computed using a fully ab initio method based on the Fano theory of resonances. Algebraic diagrammatic construction for one-particle Green’s function is utilized for the solution of the many-electron problem. An advanced and universal approach to partitioning of the configuration space into discrete states and continuum subspaces is described and employed. Total decay widths are presented for all ICD-active states of the trimer characterized by one-site ionization and additional excitation of an electron into the second shell. Selected partial decay widths are analyzed in detail, showing how three-body effects can qualitatively change the character of certain relaxation transitions. Previously unreported type of three-electron decay processes is identified in one class of the metastable states.
Interatomic Coulombic decay widths of helium trimer: Ab initio calculations.
Kolorenč, Přemysl; Sisourat, Nicolas
2015-12-14
We report on an extensive study of interatomic Coulombic decay (ICD) widths in helium trimer computed using a fully ab initio method based on the Fano theory of resonances. Algebraic diagrammatic construction for one-particle Green's function is utilized for the solution of the many-electron problem. An advanced and universal approach to partitioning of the configuration space into discrete states and continuum subspaces is described and employed. Total decay widths are presented for all ICD-active states of the trimer characterized by one-site ionization and additional excitation of an electron into the second shell. Selected partial decay widths are analyzed in detail, showing how three-body effects can qualitatively change the character of certain relaxation transitions. Previously unreported type of three-electron decay processes is identified in one class of the metastable states. PMID:26671378
Ab initio calculation of thermodynamic properties of silicon
NASA Astrophysics Data System (ADS)
Wei, Siqing; Li, Changlin; Chou, M. Y.
1994-11-01
We present a fully ab initio calculation of the thermodynamic properties for silicon within the quasiharmonic approximation, making use of volume-dependent phonon frequencies obtained from pseudopotential local-density calculations. The temperature dependence of the thermal-expansion coefficient, specific heat (at constant volume), and other related quantities are studied. We confirm that the thermal-expansion coefficient behaves differently in three temperature regions: positive for temperature below 15 K, negative between 15 and 125 K, and positive again above 125 K. This finding agrees with experiment. The abnormal (negative) thermal-expansion coefficient at low temperatures is explained through a detailed study of mode Grüneisen parameters. Both specific-heat and thermal-expansion-coefficient values calculated are in excellent agreement with experiment up to a few hundred kelvin.
Ab initio Potential Energy Surface for H-H2
NASA Technical Reports Server (NTRS)
Partridge, Harry; Bauschlicher, Charles W., Jr.; Stallcop, James R.; Levin, Eugene
1993-01-01
Ab initio calculations employing large basis sets are performed to determine an accurate potential energy surface for H-H2 interactions for a broad range of separation distances. At large distances, the spherically averaged potential determined from the calculated energies agrees well with the corresponding results determined from dispersion coefficients; the van der Waals well depth is predicted to be 75 +/- (mu)E(sub h). Large basis sets have also been applied to reexamine the accuracy of theoretical repulsive potential energy surfaces. Multipolar expansions of the computed H-H2 potential energy surface are reported for four internuclear separation distances (1.2, 1.401, 1.449, and 1.7a(sub 0) of the hydrogen molecule. The differential elastic scattering cross section calculated from the present results is compared with the measurements from a crossed beam experiment.
Ab initio calculations of grain boundaries in bcc metals
NASA Astrophysics Data System (ADS)
Scheiber, Daniel; Pippan, Reinhard; Puschnig, Peter; Romaner, Lorenz
2016-03-01
In this study, we compute grain boundary (GB) properties for a large set of GBs in bcc transition metals with a special focus on W, Mo and Fe using ab initio density functional theory (DFT) and semi-empirical second nearest neighbour modified embedded atom method (2NN-MEAM) potentials. The GB properties include GB energies, surface energies, GB excess volume and work of separation, which we analyse and then compare to experimental data. We find that the used 2NN-MEAM potentials can predict general trends of GB properties, but do not always reproduce the GB ground state structure and energy found with DFT. In particular, our results explain the experimental finding that W and Mo prefer intergranular fracture, while other bcc metals prefer transgranular cleavage.
Transport coefficients in diamond from ab-initio calculations
NASA Astrophysics Data System (ADS)
Löfâs, Henrik; Grigoriev, Anton; Isberg, Jan; Ahuja, Rajeev
2013-03-01
By combining the Boltzmann transport equation with ab-initio electronic structure calculations, we obtain transport coefficients for boron-doped diamond. We find the temperature dependence of the resistivity and the hall coefficients in good agreement with experimental measurements. Doping in the samples is treated via the rigid band approximation and scattering is treated in the relaxation time approximation. In contrast to previous results, the acoustic phonon scattering is the dominating scattering mechanism for the considered doping range. At room temperature, we find the thermopower, S, in the range 1-1.6 mV/K and the power factor, S2σ, in the range 0.004-0.16 μW /cm K2.
Ab initio engineering of materials with stacked hexagonal tin frameworks
NASA Astrophysics Data System (ADS)
Shao, Junping; Beaufils, Clément; Kolmogorov, Aleksey N.
2016-07-01
The group-IV tin has been hypothesized to possess intriguing electronic properties in an atom-thick hexagonal form. An attractive pathway of producing sizable 2D crystallites of tin is based on deintercalation of bulk compounds with suitable tin frameworks. Here, we have identified a new synthesizable metal distannide, NaSn2, with a 3D stacking of flat hexagonal layers and examined a known compound, BaSn2, with buckled hexagonal layers. Our ab initio results illustrate that despite being an exception to the 8-electron rule, NaSn2 should form under pressures easily achievable in multi-anvil cells and remain (meta)stable under ambient conditions. Based on calculated Z2 invariants, the predicted NaSn2 may display topologically non-trivial behavior and the known BaSn2 could be a strong topological insulator.
Ab initio electronic stopping power of protons in bulk materials
NASA Astrophysics Data System (ADS)
Shukri, Abdullah Atef; Bruneval, Fabien; Reining, Lucia
2016-01-01
The electronic stopping power is a crucial quantity for ion irradiation: it governs the deposited heat, the damage profile, and the implantation depth. Whereas experimental data are readily available for elemental solids, the data are much more scarce for compounds. Here we develop a fully ab initio computational scheme based on linear response time-dependent density-functional theory to predict the random electronic stopping power (RESP) of materials without any empirical fitting. We show that the calculated RESP compares well with experimental data, when at full convergence, with the inclusion of the core states and of the exchange correlation. We evaluate the unexpectedly limited magnitude of the nonlinear terms in the RESP by comparing with other approaches based on the time propagation of time-dependent density-functional theory. Finally, we check the validity of a few empirical rules of thumbs that are commonly used to estimate the electronic stopping power.
Isofulminic acid, HONC: Ab initio theory and microwave spectroscopy.
Mladenović, Mirjana; Lewerenz, Marius; McCarthy, Michael C; Thaddeus, Patrick
2009-11-01
Isofulminic acid, HONC, the most energetic stable isomer of isocyanic acid HNCO, higher in energy by 84 kcal/mol, has been detected spectroscopically by rotational spectroscopy supported by coupled cluster electronic structure calculations. The fundamental rotational transitions of the normal, carbon-13, oxygen-18, and deuterium isotopic species have been detected in the centimeter band in a molecular beam by Fourier transform microwave spectroscopy, and rotational constants and nitrogen and deuterium quadrupole coupling constants have been derived. The measured constants agree well with those predicted by ab initio calculations. A number of other electronic and spectroscopic parameters of isofulminic acid, including the dipole moment, vibrational frequencies, infrared intensities, and centrifugal distortion constants have been calculated at a high level of theory. Isofulminic acid is a good candidate for astronomical detection with radio telescopes because it is highly polar and its more stable isomers (HNCO, HOCN, and HCNO) have all been identified in space. PMID:19895013
Ab initio study of helium behavior in titanium tritides
Liang, J. H.; Dai, Yunya; Yang, Li; Peng, SM; Fan, K. M.; Long, XG; Zhou, X. S.; Zu, Xiaotao; Gao, Fei
2013-03-01
Ab initio calculations based on density functional theory have been performed to investigate the relative stability of titanium tritides and the helium behavior in stable titanium tritides. The results show that the β-phase TiT1.5 without two tritium along the [100] direction (TiT1.5[100]) is more stable than other possible structures. The stability of titanium tritides decrease with the increased generation of helium in TiT1.5[100]. In addition, helium generated by tritium decay prefers locating at a tetrahedral site, and favorably migrates between two neighbor vacant tetrahedral sites through an intermediate octahedral site in titanium tritides, with a migration energy of 0.23 eV. Furthermore, helium is easily accumulated on a (100) plane in β-phase TiT1.5[100].
High-throughput ab-initio dilute solute diffusion database.
Wu, Henry; Mayeshiba, Tam; Morgan, Dane
2016-01-01
We demonstrate automated generation of diffusion databases from high-throughput density functional theory (DFT) calculations. A total of more than 230 dilute solute diffusion systems in Mg, Al, Cu, Ni, Pd, and Pt host lattices have been determined using multi-frequency diffusion models. We apply a correction method for solute diffusion in alloys using experimental and simulated values of host self-diffusivity. We find good agreement with experimental solute diffusion data, obtaining a weighted activation barrier RMS error of 0.176 eV when excluding magnetic solutes in non-magnetic alloys. The compiled database is the largest collection of consistently calculated ab-initio solute diffusion data in the world. PMID:27434308
Ab initio quantum dynamics using coupled-cluster.
Kvaal, Simen
2012-05-21
The curse of dimensionality (COD) limits the current state-of-the-art ab initio propagation methods for non-relativistic quantum mechanics to relatively few particles. For stationary structure calculations, the coupled-cluster (CC) method overcomes the COD in the sense that the method scales polynomially with the number of particles while still being size-consistent and extensive. We generalize the CC method to the time domain while allowing the single-particle functions to vary in an adaptive fashion as well, thereby creating a highly flexible, polynomially scaling approximation to the time-dependent Schrödinger equation. The method inherits size-consistency and extensivity from the CC method. The method is dubbed orbital-adaptive time-dependent coupled-cluster, and is a hierarchy of approximations to the now standard multi-configurational time-dependent Hartree method for fermions. A numerical experiment is also given. PMID:22612082
Efficient Ab initio Modeling of Random Multicomponent Alloys.
Jiang, Chao; Uberuaga, Blas P
2016-03-11
We present in this Letter a novel small set of ordered structures (SSOS) method that allows extremely efficient ab initio modeling of random multicomponent alloys. Using inverse II-III spinel oxides and equiatomic quinary bcc (so-called high entropy) alloys as examples, we demonstrate that a SSOS can achieve the same accuracy as a large supercell or a well-converged cluster expansion, but with significantly reduced computational cost. In particular, because of this efficiency, a large number of quinary alloy compositions can be quickly screened, leading to the identification of several new possible high-entropy alloy chemistries. The SSOS method developed here can be broadly useful for the rapid computational design of multicomponent materials, especially those with a large number of alloying elements, a challenging problem for other approaches. PMID:27015491
Ab initio correlated calculations of rare-gas dimer quadrupoles
NASA Astrophysics Data System (ADS)
Donchev, Alexander G.
2007-10-01
This paper reports ab initio calculations of rare gas ( RG=Kr , Ar, Ne, and He) dimer quadrupoles at the second order of Møller-Plesset perturbation theory (MP2). The study reveals the crucial role of the dispersion contribution to the RG2 quadrupole in the neighborhood of the equilibrium dimer separation. The magnitude of the dispersion quadrupole is found to be much larger than that predicted by the approximate model of Hunt. As a result, the total MP2 quadrupole moment is significantly smaller than was assumed in virtually all previous related studies. An analytical model for the distance dependence of the RG2 quadrupole is proposed. The model is based on the effective-electron approach of Jansen, but replaces the original Gaussian approximation to the electron density in an RG atom by an exponential one. The role of the nonadditive contribution in RG3 quadrupoles is discussed.
Ab Initio Screening of CO2-philic Groups.
Tian, Ziqi; Saito, Tomonori; Jiang, De-En
2015-04-23
Ab initio calculations were used to identify CO2-philic groups. Over 55 neutral molecules were screened for CO2 affinity via binding energetics. It is found that poly(ethylene oxide)s (PEO) oligomers with more than three repeating units are good CO2-binding groups, consistent with the high-performance of PEO-based materials for CO2/N2 separation. More interestingly, two triazole groups linked with a methylene chain are also excellent for CO2 binding with a favorable interaction of more than 28 kJ/mol, indicating that polymers or covalent-organic frameworks (COFs) with triazoles may be utilized for CO2 capture. This work provides a useful guide to introduce promising organic groups into polymeric membranes and COFs for CO2/N2 separation media. PMID:25825811
Ab initio methods for nuclear properties - a computational physics approach
NASA Astrophysics Data System (ADS)
Maris, Pieter
2011-04-01
A microscopic theory for the structure and reactions of light nuclei poses formidable challenges for high-performance computing. Several ab-initio methods have now emerged that provide nearly exact solutions for some nuclear properties. The ab-initio no-core full configuration (NCFC) approach is based on basis space expansion methods and uses Slater determinants of single-nucleon basis functions to express the nuclear wave function. In this approach, the quantum many-particle problem becomes a large sparse matrix eigenvalue problem. The eigenvalues of this matrix give us the binding energies, and the corresponding eigenvectors the nuclear wave functions. These wave functions can be employed to evaluate experimental quantities. In order to reach numerical convergence for fundamental problems of interest, the matrix dimension often exceeds 1 billion, and the number of nonzero matrix elements may saturate available storage on present-day leadership class facilities. I discuss different strategies for distributing and solving this large sparse matrix on current multicore computer architectures, including methods to deal with with memory bottleneck. Several of these strategies have been implemented in the code MFDn, which is a parallel fortran code for nuclear structure calculations. I will show scaling behavior and compare the performance of the pure MPI version with the hybrid MPI/OpenMP code on Cray XT4 and XT5 platforms. For large core counts (typically 5,000 and above), the hybrid version is more efficient than pure MPI. With this code, we have been able to predict properties of the unstable nucleus 14F, which have since been confirmed by experiments. I will also give an overview of other recent results for nuclei in the A = 6 to 16 range with 2- and 3-body interactions. Supported in part by US DOE Grant DE-FC02-09ER41582.
n-type conductivity in Si-doped amorphous AlN: an ab initio investigation
NASA Astrophysics Data System (ADS)
Durandurdu, Murat
2016-04-01
We report the electronic structure and topology of a heavily Si-doped amorphous aluminium nitride (Al37.5Si12.5N50) using ab initio simulations. The amorphous Al37.5Si12.5N50 system is found to be structurally similar to pure amorphous aluminium nitride. It has an average coordination number of about 3.9 and exhibits a small amount of Si-Si homopolar bonds. The formation of Si-Al bonds is not very favourable. Electronic structure calculations reveal that the Si doping has a negligible effect on the band gap width but causes delocalization of the valence band tail states and a shift of the Fermi level towards the conduction band. Thus, amorphous Al37.5Si12.5N50 alloys show n-type conductivity.
Ab initio Study on Ionization Energies of 3-Amino-1-propanol
NASA Astrophysics Data System (ADS)
Wang, Ke-dong; Jia, Ying-bin; Lai, Zhen-jiang; Liu, Yu-fang
2011-06-01
Fourteen conformers of 3-amino-1-propanol as the minima on the potential energy surface are examined at the MP2/6-311++G** level. Their relative energies calculated at B3LYP, MP3 and MP4 levels of theory indicated that two most stable conformers display the intramolecular OH···N hydrogen bonds. The vertical ionization energies of these conformers calculated with ab initio electron propagator theory in the P3/aug-cc-pVTZ approximation are in agreement with experimental data from photoelectron spectroscopy. Natural bond orbital analyses were used to explain the differences of IEs of the highest occupied molecular ortibal of conformers. Combined with statistical mechanics principles, conformational distributions at various temperatures are obtained and the temperature dependence of photoelectron spectra is interpreted.
Quantum fluctuations and isotope effects in ab initio descriptions of water
Wang, Lu; Markland, Thomas E.; Ceriotti, Michele
2014-09-14
Isotope substitution is extensively used to investigate the microscopic behavior of hydrogen bonded systems such as liquid water. The changes in structure and stability of these systems upon isotope substitution arise entirely from the quantum mechanical nature of the nuclei. Here, we provide a fully ab initio determination of the isotope exchange free energy and fractionation ratio of hydrogen and deuterium in water treating exactly nuclear quantum effects and explicitly modeling the quantum nature of the electrons. This allows us to assess how quantum effects in water manifest as isotope effects, and unravel how the interplay between electronic exchange and correlation and nuclear quantum fluctuations determine the structure of the hydrogen bond in water.
Ab initio atomic recombination reaction energetics on model heat shield surfaces
NASA Technical Reports Server (NTRS)
Senese, Fredrick; Ake, Robert
1992-01-01
Ab initio quantum mechanical calculations on small hydration complexes involving the nitrate anion are reported. The self-consistent field method with accurate basis sets has been applied to compute completely optimized equilibrium geometries, vibrational frequencies, thermochemical parameters, and stable site labilities of complexes involving 1, 2, and 3 waters. The most stable geometries in the first hydration shell involve in-plane waters bridging pairs of nitrate oxygens with two equal and bent hydrogen bonds. A second extremely labile local minimum involves out-of-plane waters with a single hydrogen bond and lies about 2 kcal/mol higher. The potential in the region of the second minimum is extremely flat and qualitatively sensitive to changes in the basis set; it does not correspond to a true equilibrium structure.
3,5-Difluorobenzonitrile: ab initio calculations, FTIR and Raman spectra.
Rastogi, V K; Alcolea Palafox, M; Tanwar, R P; Mittal, Lalit
2002-07-01
Geometry, vibrational wavenumbers and several thermodynamic parameters were calculated using ab initio quantum chemical methods for the 3,5-difluorobenzonitrile molecule. The results were compared with the experimental values. With the help of three specific scaling procedures, the observed vibrational wavenumbers in FTIR and Raman spectra were analysed and assigned to different normal modes of the molecule. Most of the modes have wavenumbers in the expected range and the error obtained was in general very low. Using PEDs the contributions were determined for the different modes to each wavenumber. From the PED, it is apparent that the frequency corresponding to C[triple bond]N stretching contains 87% contribution from the C[triple bond]N stretching force constant and it mixes with the C-CN stretching mode 13 to the extent of 12%. Other general conclusions were also deduced. PMID:12164497
reaxFF Reactive Force Field for Disulfide Mechanochemistry, Fitted to Multireference ab Initio Data.
Müller, Julian; Hartke, Bernd
2016-08-01
Mechanochemistry, in particular in the form of single-molecule atomic force microscopy experiments, is difficult to model theoretically, for two reasons: Covalent bond breaking is not captured accurately by single-determinant, single-reference quantum chemistry methods, and experimental times of milliseconds or longer are hard to simulate with any approach. Reactive force fields have the potential to alleviate both problems, as demonstrated in this work: Using nondeterministic global parameter optimization by evolutionary algorithms, we have fitted a reaxFF force field to high-level multireference ab initio data for disulfides. The resulting force field can be used to reliably model large, multifunctional mechanochemistry units with disulfide bonds as designed breaking points. Explorative calculations show that a significant part of the time scale gap between AFM experiments and dynamical simulations can be bridged with this approach. PMID:27415976
An Ab Initio Approach Towards Engineering Fischer-Tropsch Surface Chemistry
Matthew Neurock; David A. Walthall
2006-05-07
One of the greatest societal challenges over the next decade is the production of cheap, renewable energy for the 10 billion people that inhabit the earth. This will require the development of various different energy sources potentially including fuels derived from methane, coal, and biomass and alternatives sources such as solar, wind and nuclear energy. One approach will be to synthesize gasoline and other fuels from simpler hydrocarbons such as CO derived from methane or other U.S. based sources such as coal. Syngas (CO and H{sub 2}) can be readily converted into higher molecular weight hydrocarbons through Fischer-Tropsch synthesis. Fischer-Tropsch synthesis involves the initiation or activation of CO and H{sub 2} bonds, the subsequent propagation steps including hydrogenation and carbon-carbon coupling, followed by chain termination reactions. Commercially viable catalysts include supported Co and Co-alloys. Over the first two years of this project we have used ab initio methods to determine the adsorption energies for all reactants, intermediates, and products along with the overall reaction energies and their corresponding activation barriers over the Co(0001) surface. Over the third year of the project we developed and advanced an ab initio-based kinetic Monte Carlo simulation code to simulate Fischer Tropsch synthesis. This report details our work over the last year which has focused on the derivation of kinetic parameters for the elementary steps involved in FT synthesis from ab initio density functional theoretical calculations and the application of the kinetic Monte Carlo algorithm to simulate the initial rates of reaction for FT over the ideal Co(0001) surface. The results from our simulations over Co(0001) indicate the importance of stepped surfaces for the activation of adsorbed CO. In addition, they demonstrate that the dominant CH{sub x}* surface intermediate under steady state conditions is CH*. This strongly suggests that hydrocarbon coupling
4He Thermophysical Properties: New Ab Initio Calculations
Hurly, John J.; Mehl, James B.
2007-01-01
Since 2000, atomic physicists have reduced the uncertainty of the helium-helium “ab initio” potential; for example, from approximately 0.6 % to 0.1 % at 4 bohr, and from 0.8 % to 0.1 % at 5.6 bohr. These results led us to: (1) construct a new inter-atomic potential ϕ07, (2) recalculate values of the second virial coefficient, the viscosity, and the thermal conductivity of 4He from 1 K to 10,000 K, and (3), analyze the uncertainties of the thermophysical properties that propagate from the uncertainty of ϕ07 and from the Born-Oppenheimer approximation of the electron-nucleon quantum mechanical system. We correct minor errors in a previous publication [J. J. Hurly and M. R. Moldover, J. Res. Nat. Inst. Standards Technol. 105, 667 (2000)] and compare our results with selected data published after 2000. The ab initio results tabulated here can serve as standards for the measurement of thermophysical properties. PMID:27110456
Infrared signature of micro-hydration in the organophosphate sarin: An ab initio study
Alam, Todd M.; Pearce, Charles Joseph
2015-06-28
The infrared (IR) spectra of micro-hydrated Sarin•(H2O)n clusters containing between one and four explicit waters have been studied using ab initio density functional theory (DFT) methods. The phosphate group P=O bond vibration region (~1270 to 1290 cm–1) revealed the largest frequency variation with hydration, with a frequency red shift reflecting the direct hydrogen bond formation between the P=O of Sarin and water. Small variations to the P-F stretch (~810 to 815 cm–1) and the C-O-P vibrational modes (~995 to 1004 cm–1) showed that the water interactions with these functional groups were minor, and that the structures of Sarin were notmore » extensively perturbed in the hydrated complexes. Increasing the number of explicit hydration waters produced only small vibrational changes in the lowest free energy complexes. These minor changes were consistent with a single water-phosphate hydrogen bond being the dominant structure, though a second water-phosphate hydrogen bond was observed in some complexes and was identified by an additional red shift of the P=O bond vibration. As a result, the H2O•H2O vibrational modes (~3450 to 3660 cm–1) increased in complexity with higher hydration levels and reflect the extended hydrogen bonding networks formed between the explicit waters in the hydrated Sarin clusters.« less
NASA Technical Reports Server (NTRS)
Jackels, C. F.
1985-01-01
Ab initio quantum chemical techniques are used to investigate covalently-bonded and hydrogen-bonded species that may be important intermediates in the reaction of hydroxyl and hydroperoxyl radicals. Stable structures of both types are identified. Basis sets of polarized double zeta quality and large scale configuration interaction wave functions are utilized. Based on electronic energies, the covalently bonded HOOOH species is 26.4 kcal/mol more stable than the OH and HO2 radicals. Similarly, the hydrogen bonded HO---HO2 species has an electronic energy 4.7 kcal/mol below that of the component radicals, after correction is made for the basis set superposition error. The hydrogen bonded form is planar, possesses one relatively normal hydrogen bond, and has the lowest energy 3A' and 1A' states that are essentially degenerate. The 1A" and 3A" excited states produced by rotation of the unpaired OH electron into the molecular plane are very slightly bound.
A combined photoelectron spectroscopy and ab initio study of the quasi-planar B24- cluster
NASA Astrophysics Data System (ADS)
Popov, Ivan A.; Piazza, Zachary A.; Li, Wei-Li; Wang, Lai-Sheng; Boldyrev, Alexander I.
2013-10-01
The structure and chemical bonding of the 24-atom boron cluster are investigated using photoelectron spectroscopy and ab initio calculations. The joint experimental and theoretical investigation shows that B24- possesses a quasi-planar structure containing fifteen outer and nine inner atoms with six of the inner atoms forming a filled pentagonal moiety. The central atom of the pentagonal moiety is puckered out of plane by 0.9 Å, reminiscent of the six-atom pentagonal caps of the well-known B12 icosahedral unit. The next closest isomer at the ROCCSD(T) level of theory has a tubular double-ring structure. Comparison of the simulated spectra with the experimental data shows that the global minimum quasi-planar B24- isomer is the major contributor to the observed photoelectron spectrum, while the tubular isomer has no contribution to the experiment. Chemical bonding analyses reveal that the periphery of the quasi-planar B24 constitutes 15 classical 2c-2e B-B σ-bonds, whereas delocalized σ- and π-bonds are found in the interior of the cluster with one unique 6c-2e π-bond responsible for bonding in the B-centered pentagon. The current work suggests that the 24-atom boron cluster continues to be quasi-2D, albeit the tendency to form filled pentagonal units, characteristic of 3D cage-like structures of bulk boron, is observed.
Sum, A.K.; Sandler, S.I.
2000-02-17
The results of ab initio calculations for cyclic clusters of methanol, ethanol, 1-propanol, and methanethiol are presented. Dimer, trimer, and tetramer clusters of all four compounds are studied, as are pentamer and hexamer clusters of methanol. From optimized clusters at HG/6--31G**, total energies and binding energies were calculated with both the HF and MP2 theories using the aug-cc-pVDZ basis set. Accurate binding energies were also calculated for the dimer and trimer of methanol using symmetry-adapted perturbation theory with the same basis set. Intermolecular and intramolecular distances, charge distribution of binding sites, binding energies, and equilibrium constants were computed to determine the hydrogen bond cooperativity effect for each species. The cooperativity effect, exclusive to hydrogen bonding systems, results form specific forces among the molecules, in particular charge-transfer processes and the greater importance of interactions between molecules not directly hydrogen bonded because of the longer range of the interactions. The ratios of equilibrium constants for forming multimer hydrogen bonds to that for dimer hydrogen bond formation increase rapidly with the cluster size, in contrast to the constant value commonly used in thermodynamic models for hydrogen bonding liquids.
Operator evolution for ab initio theory of light nuclei
NASA Astrophysics Data System (ADS)
Schuster, Micah; Quaglioni, Sofia; Johnson, Calvin; Jurgenson, Eric; Navrátil, Petr
2014-09-01
The past two decades have seen a revolution in ab initio calculations of nuclear properties. One key element has been the development of a rigorous effective interaction theory, applying unitary transformations to soften the nuclear Hamiltonian and hence accelerate the convergence as a function of the model space size. For consistency, however, one ought to apply the same transformation to other operators when calculating transitions and mean values from the eigenstates of the renormalized Hamiltonian. Working in a translationally invariant harmonic oscillator basis for the two- and three-nucleon systems, we evolve the Hamiltonian, square radius, and total dipole strength operators by the similarity renormalization group (SRG). The inclusion of up to three-body matrix elements in the 4He nucleus all but completely restores the invariance of the expectation values under the transformation. We also consider a Gaussian operator with adjustable range; short ranges have the largest absolute renormalization when including two- and three-body induced terms, while at long ranges the induced three-body contribution takes on increased relative importance. The past two decades have seen a revolution in ab initio calculations of nuclear properties. One key element has been the development of a rigorous effective interaction theory, applying unitary transformations to soften the nuclear Hamiltonian and hence accelerate the convergence as a function of the model space size. For consistency, however, one ought to apply the same transformation to other operators when calculating transitions and mean values from the eigenstates of the renormalized Hamiltonian. Working in a translationally invariant harmonic oscillator basis for the two- and three-nucleon systems, we evolve the Hamiltonian, square radius, and total dipole strength operators by the similarity renormalization group (SRG). The inclusion of up to three-body matrix elements in the 4He nucleus all but completely restores
Ab initio calculation and anharmonic force field of hypochlorous acid, HOCl
NASA Astrophysics Data System (ADS)
Halonen, L.; Ha, T.-K.
1988-03-01
Ab initio calculations on HOCl have been performed at the third-order Møller-Plesset perturbation theory level to determine the equilibrium structure and the anharmonic force field. An empirical anharmonic force field based on the ab initio results is obtained using available experimental vibration-rotation data. Four of the six harmonic and six of the ten cubic force constants have been determined experimentally, the remaining values being fixed at the ab initio values. A good fit to the experimental vibration-rotation data of four isotopic species is obtained.
Ab initio study of weakly bound halogen complexes: RX⋯PH3.
Georg, Herbert C; Fileti, Eudes E; Malaspina, Thaciana
2013-01-01
Ab initio calculations were employed to study the role of ipso carbon hybridization in halogenated compounds RX (R=methyl, phenyl, acetyl, H and X=F, Cl, Br and I) and its interaction with a phosphorus atom, as occurs in the halogen bonded complex type RX⋯PH3. The analysis was performed using ab initio MP2, MP4 and CCSD(T) methods. Systematic energy analysis found that the interaction energies are in the range -4.14 to -11.92 kJ mol(-1) (at MP2 level without ZPE correction). Effects of electronic correlation levels were evaluated at MP4 and CCSD(T) levels and a reduction of up to 27% in interaction energy obtained in MP2 was observed. Analysis of the electrostatic maps confirms that the PhCl⋯PH3 and all MeX⋯PH3 complexes are unstable. NBO analysis suggested that the charge transfer between the moieties is bigger when using iodine than bromine and chlorine. The electrical properties of these complexes (dipole and polarizability) were determined and the most important observed aspect was the systematic increase at the dipole polarizability, given by the interaction polarizability. This increase is in the range of 0.7-6.7 u.a. (about 3-7%). PMID:22895850