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.
NASA Astrophysics Data System (ADS)
Buck, Henk
In an effort to overcome a significant difference between high-level ab initio calculations and X-ray data of DNA duplexes, Fonseca Guerra et al. 9-11 studied a number of model systems for A-T and G-C basepairs at various levels of nonlocal Density Functional Theory. There was an excellent agreement with the gas-phase experimental bond enthalpies for the A-T and G-C basepairs. On the other hand the hydrogen bond lengths between the bases differ from the X-ray results. After introduction of a molecular environment as local water and Na+ ions, the agreement between theory and experiment was excellent. However, careful analysis shows that this picture is far from correct. In fact, the model was constructed as a backbone-modified DNA duplex in which the nonbonding oxygens of the phosphate linkages are completely shielded by proton addition. Experimental results with respect to backbone-modified DNAs clearly show that changes in the backbone focused on phosphate shielding result in DNA duplexes with a variety in conformational behavior. In addition to an analysis of the aforementioned contradiction, we also give molecular mechanics calculations which show that the A-T and G-C bond enthalpies are of the same order as the corresponding results of Fonseca Guerra et al. under the condition of complete anionic shielding of the nonbonded oxygens in the phosphate linkages.
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.
NASA Astrophysics Data System (ADS)
Liu, Haiyan; Lu, Zhenyu; Cisneros, G. Andrés; Yang, Weitao
2004-07-01
The determination of reaction paths for enzyme systems remains a great challenge for current computational methods. In this paper we present an efficient method for the determination of minimum energy reaction paths with the ab initio quantum mechanical/molecular mechanical approach. Our method is based on an adaptation of the path optimization procedure by Ayala and Schlegel for small molecules in gas phase, the iterative quantum mechanical/molecular mechanical (QM/MM) optimization method developed earlier in our laboratory and the introduction of a new metric defining the distance between different structures in the configuration space. In this method we represent the reaction path by a discrete set of structures. For each structure we partition the atoms into a core set that usually includes the QM subsystem and an environment set that usually includes the MM subsystem. These two sets are optimized iteratively: the core set is optimized to approximate the reaction path while the environment set is optimized to the corresponding energy minimum. In the optimization of the core set of atoms for the reaction path, we introduce a new metric to define the distances between the points on the reaction path, which excludes the soft degrees of freedom from the environment set and includes extra weights on coordinates describing chemical changes. Because the reaction path is represented by discrete structures and the optimization for each can be performed individually with very limited coupling, our method can be executed in a natural and efficient parallelization, with each processor handling one of the structures. We demonstrate the applicability and efficiency of our method by testing it on two systems previously studied by our group, triosephosphate isomerase and 4-oxalocrotonate tautomerase. In both cases the minimum energy paths for both enzymes agree with the previously reported paths.
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
Comparison of DFT and ab initio QM/MM methods for modelling reaction in chorismate synthase
NASA Astrophysics Data System (ADS)
Lawan, Narin; Ranaghan, Kara E.; Manby, Frederick R.; Mulholland, Adrian J.
2014-07-01
Quantum mechanics/molecular mechanics (QM/MM) methods are a popular tool in the investigation of enzyme reactions. Here, we compare B3LYP density functional theory (DFT) and ab initio QM/MM methods for modelling the conversion of 5-enolpyruvylshikimate-3-phosphate to chorismate in chorismate synthase. Good agreement with experimental data is only obtained at the SCS-MP2/CHARMM27 level for a reaction mechanism in which phosphate elimination precedes proton transfer. B3LYP predicts reaction energetics that are qualitatively wrong, stressing the need for ab initio QM/MM methods, and caution in interpretation of DFT results for this enzyme.
Ab initio modelling: Genesis of crystal structures
NASA Astrophysics Data System (ADS)
van de Walle, Axel
2005-05-01
Genetic algorithms prove useful to distil a complex quantum mechanical calculation of interatomic interactions down to its simplest mathematical expression. This makes it possible to predict the structure of new compounds from first principles.
Ab Initio Modelling of Steady Rotating Stars
NASA Astrophysics Data System (ADS)
Rieutord, Michel; Espinosa Lara, Francisco
Modelling isolated rotating stars at any rotation rate is a challenge for the next generation of stellar models. These models will couple dynamical aspects of rotating stars, like angular momentum and chemicals transport, with classical chemical evolution, gravitational contraction or mass-loss. Such modelling needs to be achieved in two dimensions, combining the calculation of the structure of the star, its mean flows and the time-evolution of the whole. We present here a first step in this challenging programme. It leads to the first self-consistent two-dimensional models of rotating stars in a steady state generated by the ESTER code. In these models the structure (pressure, density and temperature) and the flow fields are computed in a self-consistent way allowing the prediction of the differential rotation and the associated meridian circulation of the stars. After a presentation of the physical properties of such models and the numerical methods at work, we give the first grid of such models describing massive and intermediate-mass stars for a selection of rotation rates up to 90 % of the breakup angular velocity.
Student Modeling and Ab Initio Language Learning.
ERIC Educational Resources Information Center
Heift, Trude; Schulze, Mathias
2003-01-01
Provides examples of student modeling techniques that have been employed in computer-assisted language learning over the past decade. Describes two systems for learning German: "German Tutor" and "Geroline." Shows how a student model can support computerized adaptive language testing for diagnostic purposes in a Web-based language learning…
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 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.
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.
NASA Astrophysics Data System (ADS)
Pandey, Prasenjit; Chakraborty, Tanmoy; Mukherjee, Asok K.
2013-10-01
Ab initio theory at the HF/6-311G(d,p) level has been used to compute the hydrogen bonding thermodynamics in bulk liquid ethanol. Inter-cluster hydrogen bonding is assumed to mimic the H-bonding in bulk ethanol. Rotation of the clusters has been neglected, but translational and vibrational motions are taken into account for calculating bulk thermodynamic parameters. Results are well in agreement with an earlier report [J. Chem. Phys. 116, 4212 (2002)]. For a more accurate dipole moment of monomer, MP2/6-311++G(d,p) calculation was done. Use of the computed thermodynamic data in a statistical model yields the Kirkwood-Frohlich correlation factor and the dielectric constant of ethanol (21.0) close to the experimental value, 24.3 at 298 K.
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)
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.
TOPICAL REVIEW: Ab initio symplectic no-core shell model
NASA Astrophysics Data System (ADS)
Dytrych, T.; Sviratcheva, K. D.; Draayer, J. P.; Bahri, C.; Vary, J. P.
2008-12-01
The no-core shell model (NCSM) is a prominent ab initio method that yields a good description of the low-lying states in few-nucleon systems as well as in more complex p-shell nuclei. Nevertheless, its applicability is limited by the rapid growth of the many-body basis with larger model spaces and increasing number of nucleons. The symplectic no-core shell model (Sp-NCSM) aspires to extend the scope of the NCSM beyond the p-shell region by augmenting the conventional spherical harmonic oscillator basis with the physically relevant symplectic \\SpR{3} symmetry-adapted configurations of the symplectic shell model that describe naturally the monopole-quadrupole vibrational and rotational modes, and also partially incorporate α-cluster correlations. In this review, the models underpinning the Sp-NCSM approach, namely, the NCSM, the Elliott SU(3) model and the symplectic shell model, are discussed. Following this, a prescription for constructing translationally invariant symplectic configurations in the spherical harmonic oscillator basis is given. This prescription is utilized to unveil the extent to which symplectic configurations enter into low-lying states in 12C and 16O nuclei calculated within the framework of the NCSM with the JISP16 realistic nucleon-nucleon interaction. The outcomes of this proof-of-principle study are presented in detail.
Ab initio modelling of methane hydrate thermophysical properties.
Jendi, Z M; Servio, P; Rey, A D
2016-04-21
The key thermophysical properties of methane hydrate were determined using ab initio modelling. Using density functional theory, the second-order elastic constants, heat capacity, compressibility, and thermal expansion coefficient were calculated. A wide and relevant range of pressure-temperature conditions were considered, and the structures were assessed for stability using the mean square displacement and radial distribution functions. Methane hydrate was found to be elastically isotropic with a linear dependence of the bulk modulus on pressure. Equally significant, multi-body interactions were found to be important in hydrates, and water-water interactions appear to strongly influence compressibility like in ice Ih. While the heat capacity of hydrate was found to be higher than that of ice, the thermal expansion coefficient was significantly lower, most likely due to the lower rigidity of hydrates. The mean square displacement gave important insight into stability, heat capacity, and elastic moduli, and the radial distribution functions further confirmed stability. The presented results provide a much needed atomistic thermoelastic characterization of methane hydrates and are essential input for the large-scale applications of hydrate detection and production. PMID:27019976
AN AB INITIO MODEL FOR COSMIC-RAY MODULATION
Engelbrecht, N. E.; Burger, R. A.
2013-07-20
A proper understanding of the effects of turbulence on the diffusion and drift of cosmic rays (CRs) is of vital importance for a better understanding of CR modulation in the heliosphere. This study presents an ab initio model for CR modulation, incorporating for the first time the results yielded by a two-component turbulence transport model. This model is solved for solar minimum heliospheric conditions, utilizing boundary values chosen so that model results are in reasonable agreement with spacecraft observations of turbulence quantities in the solar ecliptic plane and along the out-of-ecliptic trajectory of the Ulysses spacecraft. These results are employed as inputs for modeled slab and two-dimensional (2D) turbulence energy spectra. The modeled 2D spectrum is chosen based on physical considerations, with a drop-off at the very lowest wavenumbers. There currently exist no models or observations for the wavenumber where this drop-off occurs, and it is considered to be the only free parameter in this study. The modeled spectra are used as inputs for parallel mean free path expressions based on those derived from quasi-linear theory and perpendicular mean free paths from extended nonlinear guiding center theory. Furthermore, the effects of turbulence on CR drifts are modeled in a self-consistent way, also employing a recently developed model for wavy current sheet drift. The resulting diffusion and drift coefficients are applied to the study of galactic CR protons and antiprotons using a 3D, steady-state CR modulation code, and sample solutions in fair to good agreement with multiple spacecraft observations are presented.
An investigation of ab initio shell-model interactions derived by no-core shell model
NASA Astrophysics Data System (ADS)
Wang, XiaoBao; Dong, GuoXiang; Li, QingFeng; Shen, CaiWan; Yu, ShaoYing
2016-09-01
The microscopic shell-model effective interactions are mainly based on the many-body perturbation theory (MBPT), the first work of which can be traced to Brown and Kuo's first attempt in 1966, derived from the Hamada-Johnston nucleon-nucleon potential. However, the convergence of the MBPT is still unclear. On the other hand, ab initio theories, such as Green's function Monte Carlo (GFMC), no-core shell model (NCSM), and coupled-cluster theory with single and double excitations (CCSD), have made many progress in recent years. However, due to the increasing demanding of computing resources, these ab initio applications are usually limited to nuclei with mass up to A = 16. Recently, people have realized the ab initio construction of valence-space effective interactions, which is obtained through a second-time renormalization, or to be more exactly, projecting the full-manybody Hamiltonian into core, one-body, and two-body cluster parts. In this paper, we present the investigation of such ab initio shell-model interactions, by the recent derived sd-shell effective interactions based on effective J-matrix Inverse Scattering Potential (JISP) and chiral effective-field theory (EFT) through NCSM. In this work, we have seen the similarity between the ab initio shellmodel interactions and the interactions obtained by MBPT or by empirical fitting. Without the inclusion of three-body (3-bd) force, the ab initio shell-model interactions still share similar defects with the microscopic interactions by MBPT, i.e., T = 1 channel is more attractive while T = 0 channel is more repulsive than empirical interactions. The progress to include more many-body correlations and 3-bd force is still badly needed, to see whether such efforts of ab initio shell-model interactions can reach similar precision as the interactions fitted to experimental data.
NASA Astrophysics Data System (ADS)
Faghaninia, Alireza; Ager, Joel W.; Lo, Cynthia S.
2015-06-01
Accurate models of carrier transport are essential for describing the electronic properties of semiconductor materials. To the best of our knowledge, the current models following the framework of the Boltzmann transport equation (BTE) either rely heavily on experimental data (i.e., semiempirical), or utilize simplifying assumptions, such as the constant relaxation time approximation (BTE-cRTA). While these models offer valuable physical insights and accurate calculations of transport properties in some cases, they often lack sufficient accuracy—particularly in capturing the correct trends with temperature and carrier concentration. We present here a transport model for calculating low-field electrical drift mobility and Seebeck coefficient of n -type semiconductors, by explicitly considering relevant physical phenomena (i.e., elastic and inelastic scattering mechanisms). We first rewrite expressions for the rates of elastic scattering mechanisms, in terms of ab initio properties, such as the band structure, density of states, and polar optical phonon frequency. We then solve the linear BTE to obtain the perturbation to the electron distribution—resulting from the dominant scattering mechanisms—and use this to calculate the overall mobility and Seebeck coefficient. Therefore, we have developed an ab initio model for calculating mobility and Seebeck coefficient using the Boltzmann transport (aMoBT) equation. Using aMoBT, we accurately calculate electrical transport properties of the compound n -type semiconductors, GaAs and InN, over various ranges of temperature and carrier concentration. aMoBT is fully predictive and provides high accuracy when compared to experimental measurements on both GaAs and InN, and vastly outperforms both semiempirical models and the BTE-cRTA. Therefore, we assert that this approach represents a first step towards a fully ab initio carrier transport model that is valid in all compound semiconductors.
Ab initio no-core shell model with continuum
NASA Astrophysics Data System (ADS)
Navratil, Petr
2008-04-01
The ab initio no-core shell model (NCSM) is a many-body approach to nuclear structure of light nuclei. The NCSM adopts an effective interaction theory to transform fundamental inter-nucleon interactions into effective interactions for a specified nucleus in a selected harmonic oscillator basis space [1]. The method is capable of predicting nuclear structure from inter-nucleon forces derived from quantum chromodynamics by means of chiral effective field theory [2]. NCSM extensions to the microscopic description of nuclear reactions are now under development. In my talk, I will first discuss our recent calculations of the ^4He total photo-absorption cross section using two- and three-nucleon interactions from chiral effective field theory [3]. I will then outline our effort to augment the NCSM by the resonating group method (RGM) technique to develop a new method capable of describing simultaneously both bound states and nuclear reactions on light nuclei [4]. This approach, which preserves translational symmetry and the Pauli principle, will allow us to calculate cross sections of reactions important for astrophysics and describe weakly-bound systems from first principles. I will present our first phase shift results for neutron scattering off ^3H, ^4He and ^7Li and proton scattering off ^3He, ^4He and ^7Be using realistic nucleon-nucleon potentials. 3mm [1] P. Navr'atil, J. P. Vary and B. R. Barrett, Phys. Rev. C 62, 054311 (2000). [2] P. Navr'atil and V. G. Gueorguiev and J. P. Vary, W. E. Ormand and A. Nogga, Phys. Rev. Lett. 99, 042501 (2007). [3] S. Quaglioni and P. Navr'atil, Phys. Lett. B 652, 370 (2007). [4] S. Quaglioni and P. Navr'atil, arXiv:0712.0855.
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.
Ab Initio Quantum Mechanical/Molecular Mechanical Studies of Histone Modifying Enzymes
NASA Astrophysics Data System (ADS)
Zhang, Yingkai
Histone proteins that form the nucleosome core are subject to a variety of post-translational transformations. These histone modifications make up the histone code which extends the information in the genetic code and is emerging as an essential mechanism to regulate gene expression. In spite of a current flurry of significant advances in experimental studies, there has been little theoretical understanding regarding how enzymes generate or remove these modifications. Very recently, we have made excellent progresses in investigating two such important histone-modifying enzyme families: zinc-dependent histone deacetylases (HDACs) and histone lysine methyltransferases (HKMTs). Our studies on a histonedeacetylase- like protein HDLP suggested a novel catalytic mechanism. The simulations on HKMT SET7/9 have characterized the histone lysine methylation reaction and elucidated the origin of enzyme catalysis. Our computational approaches centered on the pseudobond ab initio quantum mechanical/molecular mechanical (QM/MM) method, which allows for accurate modeling of the chemistry at the reaction active site while properly including the effects of the protein environment
Keegan, Ronan M.; Bibby, Jaclyn; Thomas, Jens; Xu, Dong; Zhang, Yang; Mayans, Olga; Winn, Martyn D.; Rigden, Daniel J.
2015-02-01
Two ab initio modelling programs solve complementary sets of targets, enhancing the success of AMPLE with small proteins. AMPLE clusters and truncates ab initio protein structure predictions, producing search models for molecular replacement. Here, an interesting degree of complementarity is shown between targets solved using the different ab initio modelling programs QUARK and ROSETTA. Search models derived from either program collectively solve almost all of the all-helical targets in the test set. Initial solutions produced by Phaser after only 5 min perform surprisingly well, improving the prospects for in situ structure solution by AMPLE during synchrotron visits. Taken together, the results show the potential for AMPLE to run more quickly and successfully solve more targets than previously suspected.
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.
Ab initio charge-carrier mobility model for amorphous molecular semiconductors
NASA Astrophysics Data System (ADS)
Massé, Andrea; Friederich, Pascal; Symalla, Franz; Liu, Feilong; Nitsche, Robert; Coehoorn, Reinder; Wenzel, Wolfgang; Bobbert, Peter A.
2016-05-01
Accurate charge-carrier mobility models of amorphous organic molecular semiconductors are essential to describe the electrical properties of devices based on these materials. The disordered nature of these semiconductors leads to percolative charge transport with a large characteristic length scale, posing a challenge to the development of such models from ab initio simulations. Here, we develop an ab initio mobility model using a four-step procedure. First, the amorphous morphology together with its energy disorder and intermolecular charge-transfer integrals are obtained from ab initio simulations in a small box. Next, the ab initio information is used to set up a stochastic model for the morphology and transfer integrals. This stochastic model is then employed to generate a large simulation box with modeled morphology and transfer integrals, which can fully capture the percolative charge transport. Finally, the charge-carrier mobility in this simulation box is calculated by solving a master equation, yielding a mobility function depending on temperature, carrier concentration, and electric field. We demonstrate the procedure for hole transport in two important molecular semiconductors, α -NPD and TCTA. In contrast to a previous study, we conclude that spatial correlations in the energy disorder are unimportant for α -NPD. We apply our mobility model to two types of hole-only α -NPD devices and find that the experimental temperature-dependent current density-voltage characteristics of all devices can be well described by only slightly decreasing the simulated energy disorder strength.
NASA Astrophysics Data System (ADS)
He, Yuping
2015-03-01
We present calculations of the thermal transport coefficients of Si-based clathrates and solar perovskites, as obtained from ab initio calculations and models, where all input parameters derived from first principles. We elucidated the physical mechanisms responsible for the measured low thermal conductivity in Si-based clatherates and predicted their electronic properties and mobilities, which were later confirmed experimentally. We also predicted that by appropriately tuning the carrier concentration, the thermoelectric figure of merit of Sn and Pb based perovskites may reach values ranging between 1 and 2, which could possibly be further increased by optimizing the lattice thermal conductivity through engineering perovskite superlattices. Work done in collaboration with Prof. G. Galli, and supported by DOE/BES Grant No. DE-FG0206ER46262.
Ab initio modeling of the motional Stark effect on MAST
De Bock, M. F. M.; Conway, N. J.; Walsh, M. J.; Carolan, P. G.; Hawkes, N. C.
2008-10-15
A multichord motional Stark effect (MSE) system has recently been built on the MAST tokamak. In MAST the {pi} and {sigma} lines of the MSE spectrum overlap due to the low magnetic field typical for present day spherical tokamaks. Also, the field curvature results in a large change in the pitch angle over the observation volume. The measured polarization angle does not relate to one local pitch angle but to an integration over all pitch angles in the observation volume. The velocity distribution of the neutral beam further complicates the measurement. To take into account volume effects and velocity distribution, an ab initio code was written that simulates the MSE spectrum on MAST. The code is modular and can easily be adjusted for other tokamaks. The code returns the intensity, polarized fraction, and polarization angle as a function of wavelength. Results of the code are presented, showing the effect on depolarization and wavelength dependence of the polarization angle. The code is used to optimize the design and calibration of the MSE diagnostic.
Dane Morgan
2010-06-10
The project began March 13, 2006, allocated for three years, and received a one year extension from March 13, 2009 to March 12, 2010. It has now completed 48 of 48 total months. The project was focused on using ab initio methods to gain insights into radiation induced segregation (RIS) in Ni-Fe-Cr alloys. The project had the following key accomplishments • Development of a large database of ab initio energetics that can be used by many researchers in the future for increased understanding of this system. For example, we have the first calculations showing a dramatic stabilization effect of Cr-Cr interstitial dumbbells in Ni. • Prediction of both vacancy and interstitial diffusion constants for Ni-Cr and Ni-Fe for dilute Cr and Fe. This work included generalization of widely used multifrequency models to make use of ab initio derived energetics and thermodynamics. • Prediction of qualitative trends of RIS from vacancy and interstitial mechanisms, suggesting the two types of defect fluxes drive Cr RIS in opposite directions. • Detailed kinetic Monte Carlo modeling of diffusion by vacancy mechanism in Ni-Cr as a function of Cr concentration. The results demonstrate that Cr content can have a significant effect on RIS. • Development of a quantitative RIS transport model, including models for thermodynamic factors and boundary conditions.
Giese, Timothy J; York, Darrin M
2016-06-14
A new approach for performing Particle Mesh Ewald in ab initio quantum mechanical/molecular mechanical (QM/MM) simulations with extended atomic orbital basis sets is presented. The new approach, the Ambient-Potential Composite Ewald (CEw) method, does not perform the QM/MM interaction with Mulliken charges nor electrostatically fit charges. Instead the nuclei and electron density interact directly with the MM environment, but in a manner that avoids the use of dense Fourier transform grids. By performing the electrostatics with the underlying QM density, the CEw method avoids self-consistent field instabilities that have been encountered with simple charge mapping procedures. Potential of mean force (PMF) profiles of the p-nitrophenyl phosphate dissociation reaction in explicit solvent are computed from PBE0/6-31G* QM/MM molecular dynamics simulations with various electrostatic protocols. The CEw profiles are shown to be stable with respect to real-space Ewald cutoff, whereas the PMFs computed from truncated and switched electrostatics produce artifacts. PBE0/6-311G**, AM1/d-PhoT, and DFTB2 QM/MM simulations are performed to generate two-dimensional PMF profiles of the phosphoryl transesterification reactions with ethoxide and phenoxide leaving groups. The semiempirical models incorrectly produce a concerted ethoxide mechanism, whereas PBE0 correctly produces a stepwise mechanism. The ab initio reaction barriers agree more closely to experiment than the semiempirical models. The failure of Mulliken-charge QM/MM-Ewald is analyzed. PMID:27171914
An ab initio molecular dynamics study of the roaming mechanism of the H2+HOC+ reaction
NASA Astrophysics Data System (ADS)
Yu, Hua-Gen
2011-08-01
We report here a direct ab initio molecular dynamics study of the p-/o-H2+HOC+ reaction on the basis of the accurate SAC-MP2 potential energy surface. The quasi-classical trajectory method was employed. This work largely focuses on the study of reaction mechanisms. A roaming mechanism was identified for this molecular ion-molecule reaction. The driving forces behind the roaming mechanism were thoroughly investigated by using a trajectory dynamics approach. In addition, the thermal rate coefficients of the H2+HOC+ reaction were calculated in the temperature range [25, 300] K and are in good agreement with experiments.
Relating Ab Initio Mechanical Behavior of Intergranular Glassy Films in Γ-Si3N4 to Continuum Scales
NASA Astrophysics Data System (ADS)
Ouyang, L.; Chen, J.; Ching, W.; Misra, A.
2006-05-01
Nanometer thin intergranular glassy films (IGFs) form in polycrystalline ceramics during sintering at high temperatures. The structure and properties of these IGFs are significantly changed by doping with rare earth elements. We have performed highly accurate large-scale ab initio calculations of the mechanical properties of both undoped and Yittria doped (Y-IGF) model by theoretical uniaxial tensile experiments. Uniaxial strain was applied by incrementally stretching the super cell in one direction, while the other two dimensions were kept constant. At each strain, all atoms in the model were fully relaxed using Vienna Ab initio Simulation Package VASP. The relaxed model at a given strain serves as the starting position for the next increment of strain. This process is carried on until the total energy (TE) and stress data show that the "sample" is fully fractured. Interesting differences are seen between the stress-strain response of undoped and Y-doped models. For the undoped model, the stress-strain behavior indicates that the initial atomic structure of the IGF is such that there is negligible coupling between the x- and the y-z directions. However, once the behavior becomes non- linear the lateral stresses increase, indicating that the atomic structure evolves with loading [1]. To relate the ab initio calculations to the continuum scales we analyze the atomic-scale deformation field under this uniaxial loading [1]. The applied strain in the x-direction is mostly accommodated by the IGF part of the model and the crystalline part experiences almost negligible strain. As the overall strain on the sample is incrementally increased, the local strain field evolves such that locations proximal to the softer spots attract higher strains. As the load progresses, the strain concentration spots coalesce and eventually form persistent strain localization zone across the IGF. The deformation pattern obtained through ab initio calculations indicates that it is possible to
Temperature dependent mechanical properties of Mo-Si-B compounds via ab initio molecular dynamics
NASA Astrophysics Data System (ADS)
Dharmawardhana, C. C.; Sakidja, R.; Aryal, S.; Ching, W. Y.
2013-07-01
A new method was proposed to obtain high temperature mechanical properties with a combination of ab initio molecular dynamics and stress-strain analyses. It was applied to compounds in the Mo-Si-B ternary system, namely, T1 (Mo5Si3) and T2 (Mo5SiB2) phases. The calculated coefficient of thermal expansion, thermal expansion anisotropy, and elastic constants agree well with those from the available experiments. The method enables us to theoretically access these properties up to 2000 K.
Stoddard, Nathan; Pichler, Peter; Duscher, Gerd J M; Windl, Wolfgang
2005-01-01
In this Letter, we present ab initio results identifying a new diffusion path for the nitrogen pair complex in silicon, resulting in an effective diffusivity of 67exp(-2.38 eV/kT) cm{sup 2}/s. This nudged elastic band result is compared with other nitrogen diffusion paths and mechanisms, and is determined to have unmatched agreement with experimental results. It is also shown that careful consideration of total energy corrections and use of a fully temperature-dependent diffusion prefactor have modest but important effects on the calculation of diffusivity for paired and for interstitial nitrogen.
Knockout reactions from p-shell nuclei : tests of ab initio structure models.
Grinyer, G. F.; Bazin, D.; Gade, A.; Tostevin, J. A.; Adrich, P.; Bowen, M. D.; Brown, B. A.; Campbell, C. M.; Cook, J. M.; Glasmacher, T.; McDaniel, S.; Navratil, P.; Obertelli, A.; Quaglioni, S.; Siwek, K.; Terry, J. R.; Weisshaar, D.; Wiringa, R. B.
2011-04-22
Absolute cross sections have been determined following single neutron knockout reactions from {sup 10}Be and {sup 10}C at intermediate energy. Nucleon density distributions and bound-state wave function overlaps obtained from both variational Monte Carlo (VMC) and no core shell model (NCSM) ab initio calculations have been incorporated into the theoretical description of knockout reactions. Comparison to experimental cross sections demonstrates that the VMC approach, with the inclusion of 3-body forces, provides the best overall agreement while the NCSM and conventional shell-model calculations both overpredict the cross sections by 20% to 30% for {sup 10}Be and by 40% to 50% for {sup 10}C, respectively. This study gains new insight into the importance of 3-body forces and continuum effects in light nuclei and provides a sensitive technique to assess the accuracy of ab initio calculations for describing these effects.
Knockout Reactions from p-Shell Nuclei: Tests of Ab Initio Structure Models
Grinyer, G. F.; Bazin, D.; Adrich, P.; Obertelli, A.; Weisshaar, D.; Gade, A.; Bowen, M. D.; Brown, B. A.; Campbell, C. M.; Cook, J. M.; Glasmacher, T.; McDaniel, S.; Siwek, K.; Terry, J. R.; Tostevin, J. A.; Navratil, P.; Quaglioni, S.; Wiringa, R. B.
2011-04-22
Absolute cross sections have been determined following single neutron knockout reactions from {sup 10}Be and {sup 10}C at intermediate energy. Nucleon density distributions and bound-state wave function overlaps obtained from both variational Monte Carlo (VMC) and no core shell model (NCSM) ab initio calculations have been incorporated into the theoretical description of knockout reactions. Comparison to experimental cross sections demonstrates that the VMC approach, with the inclusion of 3-body forces, provides the best overall agreement while the NCSM and conventional shell-model calculations both overpredict the cross sections by 20% to 30% for {sup 10}Be and by 40% to 50% for {sup 10}C, respectively. This study gains new insight into the importance of 3-body forces and continuum effects in light nuclei and provides a sensitive technique to assess the accuracy of ab initio calculations for describing these effects.
Point defect modeling in materials: Coupling ab initio and elasticity approaches
NASA Astrophysics Data System (ADS)
Varvenne, Céline; Bruneval, Fabien; Marinica, Mihai-Cosmin; Clouet, Emmanuel
2013-10-01
Modeling point defects at an atomic scale requires careful treatment of the long-range atomic relaxations. This elastic field can strongly affect point defect properties calculated in atomistic simulations because of the finite size of the system under study. This is an important restriction for ab initio methods which are limited to a few hundred atoms. We propose an original approach coupling ab initio calculations and linear elasticity theory to obtain the properties of an isolated point defect for reduced supercell sizes. The reliability and benefit of our approach are demonstrated for three problematic cases: the self-interstitial in zirconium, clusters of self-interstitials in iron, and the neutral vacancy in silicon.
Ab Initio No-Core Shell Model Calculations Using Realistic Two- and Three-Body Interactions
Navratil, P; Ormand, W E; Forssen, C; Caurier, E
2004-11-30
There has been significant progress in the ab initio approaches to the structure of light nuclei. One such method is the ab initio no-core shell model (NCSM). Starting from realistic two- and three-nucleon interactions this method can predict low-lying levels in p-shell nuclei. In this contribution, we present a brief overview of the NCSM with examples of recent applications. We highlight our study of the parity inversion in {sup 11}Be, for which calculations were performed in basis spaces up to 9{Dirac_h}{Omega} (dimensions reaching 7 x 10{sup 8}). We also present our latest results for the p-shell nuclei using the Tucson-Melbourne TM three-nucleon interaction with several proposed parameter sets.
Density-matrix based determination of low-energy model Hamiltonians from ab initio wavefunctions
Changlani, Hitesh J.; Zheng, Huihuo; Wagner, Lucas K.
2015-09-14
We propose a way of obtaining effective low energy Hubbard-like model Hamiltonians from ab initio quantum Monte Carlo calculations for molecular and extended systems. The Hamiltonian parameters are fit to best match the ab initio two-body density matrices and energies of the ground and excited states, and thus we refer to the method as ab initio density matrix based downfolding. For benzene (a finite system), we find good agreement with experimentally available energy gaps without using any experimental inputs. For graphene, a two dimensional solid (extended system) with periodic boundary conditions, we find the effective on-site Hubbard U{sup ∗}/t to be 1.3 ± 0.2, comparable to a recent estimate based on the constrained random phase approximation. For molecules, such parameterizations enable calculation of excited states that are usually not accessible within ground state approaches. For solids, the effective Hamiltonian enables large-scale calculations using techniques designed for lattice models.
NASA Astrophysics Data System (ADS)
Tapia, O.; Andres, J.; Aullo, J. M.; Bränden, C.-I.
1985-11-01
The electronic mechanisms of a model hydride transfer reaction are theoretically studied with ab inito RHF and UHF SCF MO procedures at the 4-31G basis set level and analytical gradient methods. The model system describes the reduction of cyclopropenyl cation to cyclopropene by the oxidation of lithium hydride to lithium cation. The molecular fragments corresponding to the asymptotic reactive channels characterizing the stepwise mechanisms currently discussed in the literature have been characterized. The binding energy between the fragments is estimated within a simple electrostatic approximate scheme. The results show that a hydride-ion mechanism is a likely pathway for this particular system. The system is thereafter thoroughly studied from the supermolecule approach. Reaction paths for the ground and first triplet electronic states have been calculated. The hypersurface is explored from a geometrical disposition of the reactants that mimics the one found in several dehydrogenases (perpendicular configuration). A hydride ion is found to be the particle transferred on the unconstrained as well as the constrained reaction pathways in the ground electronic state. In the triplet state (perpendicular configuration) the mechanism is stepwise: electron transfer followed by a hydrogen atom transfer. It has been noticed that the perpendicular geometrical disposition of the reactants plays an important role by polarizing the susceptible cyclopropene C-H bond in the sense of increasing the electronic density at the hydrogen nucleus. This provides a clue to rationalize several dehydrogenase's active site structure and mechanism. The reactant molecular complex found in the inverted potential energy curves, namely the LiH---Cp+ association has an electronic distribution which can be described as a hydride ion cementing two electron deficient centers corresponding to the cyclopropenyl and the lithium cations. Direct CI calculations confirm the overall picture obtained above.
Monte Carlo Shell Model for ab initio nuclear structure
NASA Astrophysics Data System (ADS)
Abe, T.; Maris, P.; Otsuka, T.; Shimizu, N.; Utsuno, Y.; Vary, J. P.
2014-03-01
We report on our recent application of the Monte Carlo Shell Model to no-core calculations. At the initial stage of the application, we have performed benchmark calculations in the p-shell region. Results are compared with those in the Full Configuration Interaction and No-Core Full Configuration methods. These are found to be consistent with each other within quoted uncertainties when they could be quantified. The preliminary results in Nshell = 5 reveal the onset of systematic convergence pattern.
Ab Initio QM/MM Study of the Ester-hydrolysis Reaction Mechanism in Haloalkane Dehalogenase
NASA Astrophysics Data System (ADS)
Zhang, Yiming; Zhou, Yu; Nayak, Saroj; Garcia, Angel
2006-03-01
Ab Initio QM/MM calculations are used to investigate the ester-hydrolysis step of dichloroethane hydrolysis catalyzed by haloalkane dehalogenase. Amino acids around the active site (which includes ASP124, HIS289, ASP260, TRP125, TRP175), dichoroethane and water are treated by QM at a level of HF/6-31G(d,p). The remainder of the protein and solvent are treated classically. Two scenarios of hydrolysis mechanism for the alkyl-enzyme intermediate have been considered. In one, the HIS289-catalyzed water oxygen could be incorporated in the carboxylate group of ASP124, leading the cleavage of one of the original carbonyl bonds on ASP124. In the other, the ASP124 and HIS289 as general base, activate water as the nucleophilic agent, which attacks the alkyl carbon in substrate. The reaction paths and potential energy profiles are compared for both mechanisms.
NASA Astrophysics Data System (ADS)
Zeng, Xiancheng; Hu, Hao; Hu, Xiangqian; Cohen, Aron J.; Yang, Weitao
2008-03-01
Electron transfer (ET) reactions are one of the most important processes in chemistry and biology. Because of the quantum nature of the processes and the complicated roles of the solvent, theoretical study of ET processes is challenging. To simulate ET processes at the electronic level, we have developed an efficient density functional theory (DFT) quantum mechanical (QM)/molecular mechanical (MM) approach that uses the fractional number of electrons as the order parameter to calculate the redox free energy of ET reactions in solution. We applied this method to study the ET reactions of the aqueous metal complexes Fe(H2O)62+/3+ and Ru(H2O)62+/3+. The calculated oxidation potentials, 5.82 eV for Fe(II/III) and 5.14 eV for Ru(II/III), agree well with the experimental data, 5.50 and 4.96 eV, for iron and ruthenium, respectively. Furthermore, we have constructed the diabatic free energy surfaces from histogram analysis based on the molecular dynamics trajectories. The resulting reorganization energy and the diabatic activation energy also show good agreement with experimental data. Our calculations show that using the fractional number of electrons (FNE) as the order parameter in the thermodynamic integration process leads to efficient sampling and validate the ab initio QM/MM approach in the calculation of redox free energies.
Ab initio modeling of 2D layered organohalide lead perovskites.
Fraccarollo, Alberto; Cantatore, Valentina; Boschetto, Gabriele; Marchese, Leonardo; Cossi, Maurizio
2016-04-28
A number of 2D layered perovskites A2PbI4 and BPbI4, with A and B mono- and divalent ammonium and imidazolium cations, have been modeled with different theoretical methods. The periodic structures have been optimized (both in monoclinic and in triclinic systems, corresponding to eclipsed and staggered arrangements of the inorganic layers) at the DFT level, with hybrid functionals, Gaussian-type orbitals and dispersion energy corrections. With the same methods, the various contributions to the solid stabilization energy have been discussed, separating electrostatic and dispersion energies, organic-organic intralayer interactions and H-bonding effects, when applicable. Then the electronic band gaps have been computed with plane waves, at the DFT level with scalar and full relativistic potentials, and including the correlation energy through the GW approximation. Spin orbit coupling and GW effects have been combined in an additive scheme, validated by comparing the computed gap with well known experimental and theoretical results for a model system. Finally, various contributions to the computed band gaps have been discussed on some of the studied systems, by varying some geometrical parameters and by substituting one cation in another's place. PMID:27131557
Ab initio modeling of 2D layered organohalide lead perovskites
NASA Astrophysics Data System (ADS)
Fraccarollo, Alberto; Cantatore, Valentina; Boschetto, Gabriele; Marchese, Leonardo; Cossi, Maurizio
2016-04-01
A number of 2D layered perovskites A2PbI4 and BPbI4, with A and B mono- and divalent ammonium and imidazolium cations, have been modeled with different theoretical methods. The periodic structures have been optimized (both in monoclinic and in triclinic systems, corresponding to eclipsed and staggered arrangements of the inorganic layers) at the DFT level, with hybrid functionals, Gaussian-type orbitals and dispersion energy corrections. With the same methods, the various contributions to the solid stabilization energy have been discussed, separating electrostatic and dispersion energies, organic-organic intralayer interactions and H-bonding effects, when applicable. Then the electronic band gaps have been computed with plane waves, at the DFT level with scalar and full relativistic potentials, and including the correlation energy through the GW approximation. Spin orbit coupling and GW effects have been combined in an additive scheme, validated by comparing the computed gap with well known experimental and theoretical results for a model system. Finally, various contributions to the computed band gaps have been discussed on some of the studied systems, by varying some geometrical parameters and by substituting one cation in another's place.
Structure and mechanical properties of cement and intermetallic compounds via ab-initio simulations
NASA Astrophysics Data System (ADS)
Dharmawardhana, Chamila Chathuranga
Calcium silicate hydrates comprise a class of minerals formed synthetically during Portland cement hydration or naturally through various geological processes. The importance of these minerals is immense since they are the primary binding phases for Portland cement derived construction materials. Efforts spanning centuries have been devoted to understand the structural aspects of cohesion in these minerals. In recent years, the focus has progressively turned to atomic level comprehension. Structurally these minerals can range from crystalline to highly disordered amorphous phases. This thesis focuses upon unraveling the nature of chemical bonding in a large subset of calcium silicate hydrate (CSH) crystals. Thus their electronic structure was calculated and bonding mechanisms were investigated quantitatively. Results highlight a wide range of contributions from each type of bonding (Si-O, Ca-O, O-H and hydrogen bond) with respect to silicate polymerization, crystal symmetry, water and OH content. Consequently, total bond order density (TBOD) was designated as the overall single criterion for characterizing crystal cohesion. The TBOD categorization indicates that a rarely known orthorhombic phase Suolunite is closest to the ideal composition and structure of cement. Present work finds the relationship of partial bond order density (PBOD) of each bond species, especially HBs to the mechanical properties of CSH crystals. This can be used as a basis to validate existing C-S-H models and to build improved ones. This work goes further and validates the recently proposed models (2014) for C-S-H (I) phase on the same basis of proposed electronic structure parameters. Then the respective Calcium aluminosilicate hydrates C-A-S-H (I) phase models are proposed. Finally, these results lead to improved interpretations and construction of realistic atomistic models of cement hydrates. Ab initio molecular dynamics (AIMD) could be vital to solve critical problems in complex
Converging sequences in the ab initio no-core shell model
Forssen, C.; Vary, J. P.; Caurier, E.; Navratil, P.
2008-02-15
We demonstrate the existence of multiple converging sequences in the ab initio no-core shell model. By examining the underlying theory of effective operators, we expose the physical foundations for the alternative pathways to convergence. This leads us to propose a revised strategy for evaluating effective interactions for A-body calculations in restricted model spaces. We suggest that this strategy is particularly useful for applications to nuclear processes in which states of both parities are used simultaneously, such as for transition rates. We demonstrate the utility of our strategy with large-scale calculations in light nuclei.
NASA Astrophysics Data System (ADS)
Meisel, David D.; Szasz, Csilla; Kero, Johan
2008-06-01
The Arecibo UHF radar is able to detect the head-echos of micron-sized meteoroids up to velocities of 75 km/s over a height range of 80 140 km. Because of their small size there are many uncertainties involved in calculating their above atmosphere properties as needed for orbit determination. An ab initio model of meteor ablation has been devised that should work over the mass range 10-16 kg to 10-7 kg, but the faint end of this range cannot be observed by any other method and so direct verification is not possible. On the other hand, the EISCAT UHF radar system detects micrometeors in the high mass part of this range and its observations can be fit to a “standard” ablation model and calibrated to optical observations (Szasz et al. 2007). In this paper, we present a preliminary comparison of the two models, one observationally confirmable. Among the features of the ab initio model that are different from the “standard” model are: (1) uses the experimentally based low pressure vaporization theory of O’Hanlon (A users’s guide to vacuum technology, 2003) for ablation, (2) uses velocity dependent functions fit from experimental data on heat transfer, luminosity and ionization efficiencies measured by Friichtenicht and Becker (NASA Special Publication 319: 53, 1973) for micron sized particles, (3) assumes a density and temperature dependence of the micrometeoroids and ablation product specific heats, (4) assumes a density and size dependent value for the thermal emissivity and (5) uses a unified synthesis of experimental data for the most important meteoroid elements and their oxides through least square fits (as functions of temperature, density, and/or melting point) of the tables of thermodynamic parameters given in Weast (CRC Handbook of Physics and Chemistry, 1984), Gray (American Institute of Physics Handbook, 1972), and Cox (Allen’s Astrophysical Quantities 2000). This utilization of mostly experimentally determined data is the main reason for
NASA Astrophysics Data System (ADS)
Dieterich, Johannes M.; Werner, Hans-Joachim; Mata, Ricardo A.; Metz, Sebastian; Thiel, Walter
2010-01-01
Energy and free energy barriers for acetaldehyde conversion in aldehyde oxidoreductase are determined for three reaction pathways using quantum mechanical/molecular mechanical (QM/MM) calculations on the solvated enzyme. Ab initio single-point QM/MM energies are obtained at the stationary points optimized at the DFT(B3LYP)/MM level. These ab initio calculations employ local correlation treatments [LMP2 and LCCSD(T0)] in combination with augmented triple- and quadruple-zeta basis sets, and the final coupled cluster results include MP2-based corrections for basis set incompleteness and for the domain approximation. Free energy perturbation (FEP) theory is used to generate free energy profiles at the DFT(B3LYP)/MM level for the most important reaction steps by sampling along the corresponding reaction paths using molecular dynamics. The ab initio and FEP QM/MM results are combined to derive improved estimates of the free energy barriers, which differ from the corresponding DFT(B3LYP)/MM energy barriers by about 3 kcal mol-1. The present results confirm the qualitative mechanistic conclusions from a previous DFT(B3LYP)/MM study. Most favorable is a three-step Lewis base catalyzed mechanism with an initial proton transfer from the cofactor to the Glu869 residue, a subsequent nucleophilic attack that yields a tetrahedral intermediate (IM2), and a final rate-limiting hydride transfer. The competing metal center activated pathway has the same final step but needs to overcome a higher barrier in the initial step on the route to IM2. The concerted mechanism has the highest free energy barrier and can be ruled out. While confirming the qualitative mechanistic scenario proposed previously on the basis of DFT(B3LYP)/MM energy profiles, the present ab initio and FEP QM/MM calculations provide corrections to the barriers that are important when aiming at high accuracy.
Symmetry-Adapted Ab Initio Shell Model for Nuclear Structure Calculations
NASA Astrophysics Data System (ADS)
Draayer, J. P.; Dytrych, T.; Launey, K. D.; Langr, D.
2012-05-01
An innovative concept, the symmetry-adapted ab initio shell model, that capitalizes on partial as well as exact symmetries that underpin the structure of nuclei, is discussed. This framework is expected to inform the leading features of nuclear structure and reaction data for light and medium mass nuclei, which are currently inaccessible by theory and experiment and for which predictions of modern phenomenological models often diverge. We use powerful computational and group-theoretical algorithms to perform ab initio CI (configuration-interaction) calculations in a model space spanned by SU(3) symmetry-adapted many-body configurations with the JISP16 nucleon-nucleon interaction. We demonstrate that the results for the ground states of light nuclei up through A = 16 exhibit a strong dominance of low-spin and high-deformation configurations together with an evident symplectic structure. This, in turn, points to the importance of using a symmetry-adapted framework, one based on an LS coupling scheme with the associated spatial configurations organized according to deformation.
NASA Technical Reports Server (NTRS)
Lawson, John W.; Bauschlicher, Charles W.; Daw, Murray
2011-01-01
Refractory materials such as metallic borides, often considered as ultra high temperature ceramics (UHTC), are characterized by high melting point, high hardness, and good chemical inertness. These materials have many applications which require high temperature materials that can operate with no or limited oxidation. Ab initio, first principles methods are the most accurate modeling approaches available and represent a parameter free description of the material based on the quantum mechanical equations. Using these methods, many of the intrinsic properties of these material can be obtained. We performed ab initio calculations based on density functional theory for the UHTC materials ZrB2 and HfB2. Computational results are presented for structural information (lattice constants, bond lengths, etc), electronic structure (bonding motifs, densities of states, band structure, etc), thermal quantities (phonon spectra, phonon densities of states, specific heat), as well as information about point defects such as vacancy and antisite formation energies.
Nonadiabatic ab initio dynamics of a model protonated Schiff base of 9-cis retinal.
Chung, Wilfredo Credo; Nanbu, Shinkoh; Ishida, Toshimasa
2010-08-19
The dynamics of the photoisomerization of a model protonated Schiff base of 9-cis retinal in isorhodopsin is investigated using nonadiabatic molecular dynamics simulation combined with ab initio quantum chemical calculations on-the-fly. The quantum chemical part is treated at the complete-active space self-consistent field level for six electrons in six active pi orbitals with the 6-31G basis set (CASSCF(6,6)/6-31G). The probabilities of nonadiabatic transitions between the S(1) ((1)pipi*) and S(0) states are estimated in light of the Zhu-Nakamura theory. The photoinduced cis-trans isomerization of 9-cis retinal proceeds slower than that of its 11-cis analogue and at a lower quantum yield, confirming experimental observations. An energetic barrier in the excited state impedes the elongation and twist of the C(9)=C(10) stretch and torsion coordinates, respectively, resulting in the trapping of trajectories before transition. Consequently, the isomerization takes longer time and the transition more often occurs at smaller twist angle of =C(8)-C(9)=C(10)-C(11)=, which leads to regeneration of the 9-cis reactant. Thus, neither the smaller twist observed in the X-ray crystal nor the slower movement of nuclei in the transition region would be the main reason for the longer reaction time and lower yield. A well-known space-saving asynchronous bicycle pedal or crankshaft photoisomerization mechanism is found to be operational in 9-cis retinal. The simulation in vacuo suggests that the excited-state barrier and the photoisomerization itself are intrinsic properties of the visual chromophore and not triggered mainly by the protein environment that surrounds the chromophore. PMID:20666503
AB INITIO Modeling of Thermomechanical Properties of Mo-Based Alloys for Fossil Energy Conversion
Ching, Wai-Yim
2013-12-31
In this final scientific/technical report covering the period of 3.5 years started on July 1, 2011, we report the accomplishments on the study of thermo-mechanical properties of Mo-based intermetallic compounds under NETL support. These include computational method development, physical properties investigation of Mo-based compounds and alloys. The main focus is on the mechanical and thermo mechanical properties at high temperature since these are the most crucial properties for their potential applications. In particular, recent development of applying ab initio molecular dynamic (AIMD) simulations to the T1 (Mo{sub 5}Si{sub 3}) and T2 (Mo{sub 5}SiB{sub 2}) phases are highlighted for alloy design in further improving their properties.
NASA Astrophysics Data System (ADS)
Liu, Shi-Yu; Liu, Shiyang; Li, De-Jun; Wang, Sanwu; Guo, Jing; Shen, Yaogen
2015-02-01
Utilizing a combination of ab initio density-functional theory and thermodynamics formalism, we have established the microscopic mechanisms for oxidation of the binary and ternary alloy surfaces and provided a clear explanation for the experimental results of the oxidation. We construct three-dimensional surface phase diagrams (SPDs) for oxygen adsorption on three different Nb-X(110) (X = Ti, Al or Si) binary alloy surfaces. On the basis of the obtained SPDs, we conclude a general microscopic mechanism for the thermodynamic oxidation, that is, under O-rich conditions, a uniform single-phase SPD (type I) and a nonuniform double-phase SPD (type II) correspond to the sustained complete selective oxidation and the non-sustained partial selective oxidation by adding the X element, respectively. Furthermore, by revealing the framework of thermodynamics for the oxidation mechanism of ternary alloys through the comparison of the surface energies of two separated binary alloys, we provide an understanding for the selective oxidation behavior of the Nb ternary alloy surfaces. Using these general microscopic mechanisms, one could predict the oxidation behavior of any binary and multi-component alloy surfaces based on thermodynamics considerations.
Liu, Shi-Yu; Liu, Shiyang; Li, De-Jun; Wang, Sanwu; Guo, Jing; Shen, Yaogen
2015-02-14
Utilizing a combination of ab initio density-functional theory and thermodynamics formalism, we have established the microscopic mechanisms for oxidation of the binary and ternary alloy surfaces and provided a clear explanation for the experimental results of the oxidation. We construct three-dimensional surface phase diagrams (SPDs) for oxygen adsorption on three different Nb-X(110) (X = Ti, Al or Si) binary alloy surfaces. On the basis of the obtained SPDs, we conclude a general microscopic mechanism for the thermodynamic oxidation, that is, under O-rich conditions, a uniform single-phase SPD (type I) and a nonuniform double-phase SPD (type II) correspond to the sustained complete selective oxidation and the non-sustained partial selective oxidation by adding the X element, respectively. Furthermore, by revealing the framework of thermodynamics for the oxidation mechanism of ternary alloys through the comparison of the surface energies of two separated binary alloys, we provide an understanding for the selective oxidation behavior of the Nb ternary alloy surfaces. Using these general microscopic mechanisms, one could predict the oxidation behavior of any binary and multi-component alloy surfaces based on thermodynamics considerations.
Mechanism of alkane dehydrogenation catalyzed by acidic zeolites: Ab initio transition path sampling
NASA Astrophysics Data System (ADS)
Bučko, Tomáš; Benco, Lubomir; Dubay, Orest; Dellago, Christoph; Hafner, Jürgen
2009-12-01
The dehydrogenation of propane over acidic chabazite has been studied using ab initio density-functional simulations in combination with static transition-state searches and dynamic transition path sampling (TPS) methods at elevated temperatures. The acidic zeolite has been modeled both using a small cluster and a large periodic model consisting of two unit cells, the TPS simulations allow to account for the effect of temperature and entropy. In agreement with experimental observations we find propene as the dominant reaction product and that the barrier for the dehydrogenation of a methyl group is higher than that for a methylene group. However, whereas all studies based on small cluster models (including the present one) conclude that the reaction proceeds via the formation of an alkoxy intermediate, our TPS studies based on a large periodic model lead to the conclusion that propene formation occurs via the formation of various forms of propyl cations stabilized by entropy, while the formation of an alkoxy species is a relatively rare event. It was observed only in 15% of the reactive trajectories for methyl dehydrogenation and even in only 8% of the methylene dehydrogenation reactions. Our studies demonstrate the importance of entropic effects and the need to account for the structure and flexibility of the zeolitic framework by using large periodic models.
NASA Astrophysics Data System (ADS)
Milowska, Karolina; Birowska, Magdalena; Majewski, Jacek A.
2012-02-01
We present results of extensive theoretical studies of mechanical, electric, and magnetic properties of functionalized carbon nanotubes (CNTs). Our studies are based on the ab initio calculations in the framework of the density functional theory. We have performed calculations for various metallic and semiconductor single wall CNTs, functionalized with simple organic molecules such as OH, COOH, NHn, CHn and metals, Al, Fe, Ni, Cu, Zn, and Pd. We have determined the stability of the functionalized CNTs, their elastic moduli, conductance, and magnetic moments (in the case of CNTs decorated with magnetic ions). These studies shed light on physical mechanisms governing the binding of the adsorbed molecules and also provide valuable quantitative predictions that are of importance for design of novel composite materials and functional devices. In particular, we find out that the Young's modulus of functionalized CNTs is smaller than in the case of bare CNTs, however it is large enough to provide a strong enforcement of composites. The functionalization with molecules leads also to the metallization of semiconducting CNTs, being relevant in the context of CNT interconnects, whereas the functionalization with metals might be used to cut CNTs into ribbons.
An ab initio model for the modulation of galactic cosmic-ray electrons
Engelbrecht, N. E.; Burger, R. A.
2013-12-20
The modulation of galactic cosmic-ray electrons is studied using an ab initio three-dimensional steady state cosmic-ray modulation code in which the effects of turbulence on both the diffusion and drift of these cosmic-rays are treated as self-consistently as possible. A significant refinement is that a recent two-component turbulence transport model is used. This model yields results in reasonable agreement with observations of turbulence quantities throughout the heliosphere. The sensitivity of computed galactic electron intensities to choices of various turbulence parameters pertaining to the dissipation range of the slab turbulence spectrum, and to the choice of model of dynamical turbulence, is demonstrated using diffusion coefficients derived from the quasi-linear and extended nonlinear guiding center theories. Computed electron intensities and latitude gradients are also compared with spacecraft observations.
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.
NASA Astrophysics Data System (ADS)
Olsson, P. A. T.; Kese, K.; Kroon, M.; Alvarez Holston, A.-M.
2015-06-01
In this work we report the results of an ab initio study of the transgranular fracture toughness and cleavage of brittle zirconium hydrides. We use the Griffith-Irwin relation to assess the fracture toughness using calculated surface energy and estimated isotropic Voigt-Reuss-Hill averages of the elastic constants. The calculated fracture toughness values are found to concur well with experimental data, which implies that fracture is dominated by cleavage failure. To investigate the cleavage energetics, we model the decohesion process. To describe the interplanar interaction we adopt Rose’s universal binding energy relation, which is found to reproduce the behaviour accurately. The modelling shows that the work of fracture and ductility decreases with increasing hydrogen content.
Large basis ab initio shell model investigation of {sup 9}Be and {sup 11}Be
Forssen, C.; Navratil, P.; Ormand, W.E.; Caurier, E.
2005-04-01
We present the first ab initio structure investigation of the loosely bound {sup 11}Be nucleus, together with a study of the lighter isotope {sup 9}Be. The nuclear structure of these isotopes is particularly interesting because of the appearance of a parity-inverted ground state in {sup 11}Be. Our study is performed in the framework of the ab initio no-core shell model. Results obtained using four different, high-precision two-nucleon interactions, in model spaces up to 9({Dirac_h}/2{pi}){omega}, are shown. For both nuclei, and all potentials, we reach convergence in the level ordering of positive- and negative-parity spectra separately. Concerning their relative position, the positive-parity states are always too high in excitation energy, but a fast drop with respect to the negative-parity spectrum is observed when the model space is increased. This behavior is most dramatic for {sup 11}Be. In the largest model space we were able to reach, the 1/2{sup +} level has dropped down to become either the first or the second excited state, depending on which interaction we use. We also observe a contrasting behavior in the convergence patterns for different two-nucleon potentials and argue that a three-nucleon interaction is needed to explain the parity inversion. Furthermore, large-basis calculations of {sup 13}C and {sup 11}B are performed. This allows us to study the systematics of the position of the first unnatural-parity state in the N=7 isotone and the A=11 isobar. The {sup 11}B run in the 9({Dirac_h}/2{pi}){omega} model space involves a matrix with dimension exceeding 1.1x10{sup 9}, and is our largest calculation so far. We present results on binding energies, excitation spectra, level configurations, radii, electromagnetic observables, and {sup 10}Be+n overlap functions.
Baima, Jacopo; Zelferino, Alessandro; Olivero, Paolo; Erba, Alessandro; Dovesi, Roberto
2016-01-21
Quantum-mechanical ab initio calculations are performed to elucidate the vibrational spectroscopic features of a common irradiation-induced defect in diamond, i.e. the neutral vacancy. Raman spectra are computed analytically through a Coupled-Perturbed-Hartree-Fock/Kohn-Sham approach as a function of both different defect spin states and defect concentration. The experimental Raman features of defective diamond located in the 400-1300 cm(-1) spectral range, i.e. below the first-order line of pristine diamond at 1332 cm(-1), are well reproduced, thus corroborating the picture according to which, at low damage densities, this spectral region is mostly affected by non-graphitic sp(3) defects. No peaks above 1332 cm(-1) are found, thus ruling out previous tentative assignments of different spectral features (at 1450 and 1490 cm(-1)) to the neutral vacancy. The perturbation introduced by the vacancy to the thermal nuclear motion of carbon atoms in the defective lattice is discussed in terms of atomic anisotropic displacement parameters (ADPs), computed from converged lattice dynamics calculations. PMID:26686374
Ab initio calculations of the mechanical and electronic properties of strained Si nanowires
NASA Astrophysics Data System (ADS)
Leu, Paul W.; Svizhenko, Alexei; Cho, Kyeongjae
2008-06-01
This paper reports a systematic study of the mechanical and electronic properties of strained small diameter (0.7-2.6 nm) silicon nanowires (Si NWs) using ab initio density functional theory calculations. The values of Young’s modulus, Poisson ratio, band gap, effective mass, work function, and deformation potentials are calculated for ⟨110⟩ and ⟨111⟩ Si NWs. We find that quantum confinement in ⟨110⟩ Si NWs splits conduction band valleys and decreases transport effective mass compared to the bulk case. Consequently, additional tensile strain should not lead to further significant electron mobility improvement. An interesting finding we report in this paper is that under compressive strain, there is a dramatic decrease in deformation potentials of ⟨110⟩ Si NWs, which may result in a strong increase in electron mobilities, despite a concurrent increase in effective mass. We also observe a similar strain-induced counterplay of hole deformation potentials and effective masses for both ⟨110⟩ and ⟨111⟩ Si NWs. Finally, we do not see any significant effect of tensile or compressive strain on electron effective masses and deformation potentials in ⟨111⟩ Si NWs. The sudden changes in effective mass and deformation potentials are concurrent with a change in the conduction and valence band edge states. In ⟨110⟩ NWs, this change corresponds to a transition from direct-to-indirect band gap under strain.
NASA Astrophysics Data System (ADS)
Matsuda, Yuki
This dissertation focuses on ab-initio quantum mechanical calculations of nanoelectronics in three research topics: contact resistance properties of carbon nanotubes and graphenes (Chapters 1 through 3), electrical properties of carbon nanotubes (Chapter 4) and silicon nanowires (Chapter 5). Through all the chapters, the aim of the research is to provide useful guidelines for experimentalists. Chapter 1 presents the contact resistance of metal electrode-carbon nanotube and metal electrode-graphene interfaces for various deposited metals, based on first-principles quantum mechanical density functional and matrix Green's function methods. Chapters 2 and 3 describe inventive ways to enhance contact resistance properties as well as mechanical stabilities using "molecular anchors" (Chapter 2) or using "end-contacted" (or end-on) electrodes (Chapter 3). Chapters 1 through 3 also provide useful guidelines for nanotube assembly process which is one of the main obstacles in nanoelectronics. Chapter 4 shows accurate and detailed band structure properties of single-walled carbon nanotubes using B3LYP hybrid functional, which are critical parameters in determining the electronic properties such as small band gaps (˜0.1 eV) and effective masses. Chapter 5 details both structural and electronic properties of silicon nanowires. These results lead to the findings controlling the diameter and surface coverage by adsorbates (e.g., hydrogen) of silicon nanowires can be effectively used to optimize their properties for various applications. All the theoretical results are compared with other theoretical studies and experimental data. Notably, electronic studies using B3LYP show excellent agreement with experimental studies quantitatively, which previous quantum mechanical calculations had failed. These studies show how quantum mechanical predictions of complex phenomena can be effectively investigated computationally in nanomaterials and nanodevices. Given the difficulty, expense
NASA Astrophysics Data System (ADS)
Bacca, Sonia
2016-04-01
A brief review of models to describe nuclear structure and reactions properties is presented, starting from the historical shell model picture and encompassing modern ab initio approaches. A selection of recent theoretical results on observables for exotic light and medium-mass nuclei is shown. Emphasis is given to the comparison with experiment and to what can be learned about three-body forces and continuum properties.
NASA Astrophysics Data System (ADS)
Hagiwara, Yohsuke; Tateno, Masaru
2010-10-01
We review the recent research on the functional mechanisms of biological macromolecules using theoretical methodologies coupled to ab initio quantum mechanical (QM) treatments of reaction centers in proteins and nucleic acids. Since in most cases such biological molecules are large, the computational costs of performing ab initio calculations for the entire structures are prohibitive. Instead, simulations that are jointed with molecular mechanics (MM) calculations are crucial to evaluate the long-range electrostatic interactions, which significantly affect the electronic structures of biological macromolecules. Thus, we focus our attention on the methodologies/schemes and applications of jointed QM/MM calculations, and discuss the critical issues to be elucidated in biological macromolecular systems.
Řezáč, Jan; Hobza, Pavel
2014-08-12
Hydrogen fluoride dimer is a perfect model system for studying hydrogen bonding. Its size makes it possible to apply the most advanced theoretical methods available, yet it is a full-featured complex of molecules with nontrivial electronic structure and dynamic properties. Moreover, the dissociation energy of the HF dimer has been measured experimentally with an unparalleled accuracy of ±1 cm(-1)(Bohac et al. J. Chem. Phys. 1992, 9, 6681). In this work, we attempt to reproduce it by purely ab initio means, using advanced quantum-mechanical computational methods free of any empiricism. The purpose of this study is to demonstrate the capabilities of today's computational chemistry and to point out its limitations by identifying the contributions that introduce the largest uncertainty into the result. The dissociation energy is calculated using a composite scheme including large basis set CCSD(T) calculations, contributions of higher excitations up to CCSDTQ, relativistic and diagonal Born-Oppenheimer corrections and anharmonic vibrational calculations. The error of the calculated dissociation energy is 0.07 kcal/mol (25 cm(-1), 2.5%) when compared to the experiment. The major part of this error can be attributed to the inaccuracy of the calculations of the zero-point vibrational energy. PMID:26588277
Liu, Hanchao; Wang, Yimin; Bowman, Joel M.
2015-05-21
The calculation and characterization of the IR spectrum of liquid water have remained a challenge for theory. In this paper, we address this challenge using a combination of ab initio approaches, namely, a quantum treatment of IR spectrum using the ab initio WHBB water potential energy surface and a refined ab initio dipole moment surface. The quantum treatment is based on the embedded local monomer method, in which the three intramolecular modes of each embedded H{sub 2}O monomer are fully coupled and also coupled singly to each of six intermolecular modes. The new dipole moment surface consists of a previous spectroscopically accurate 1-body dipole moment surface and a newly fitted ab initio intrinsic 2-body dipole moment. A detailed analysis of the new dipole moment surface in terms of the coordinate dependence of the effective atomic charges is done along with tests of it for the water dimer and prism hexamer double-harmonic spectra against direct ab initio calculations. The liquid configurations are taken from previous molecular dynamics calculations of Skinner and co-workers, using the TIP4P plus E3B rigid monomer water potential. The IR spectrum of water at 300 K in the range of 0–4000 cm{sup −1} is calculated and compared with experiment, using the ab initio WHBB potential and new ab initio dipole moment, the q-TIP4P/F potential, which has a fixed-charged description of the dipole moment, and the TTM3-F potential and dipole moment surfaces. The newly calculated ab initio spectrum is in very good agreement with experiment throughout the above spectral range, both in band positions and intensities. This contrasts to results with the other potentials and dipole moments, especially the fixed-charge q-TIP4P/F model, which gives unrealistic intensities. The calculated ab initio spectrum is analyzed by examining the contribution of various transitions to each band.
NASA Astrophysics Data System (ADS)
Liu, Hanchao; Wang, Yimin; Bowman, Joel M.
2015-05-01
The calculation and characterization of the IR spectrum of liquid water have remained a challenge for theory. In this paper, we address this challenge using a combination of ab initio approaches, namely, a quantum treatment of IR spectrum using the ab initio WHBB water potential energy surface and a refined ab initio dipole moment surface. The quantum treatment is based on the embedded local monomer method, in which the three intramolecular modes of each embedded H2O monomer are fully coupled and also coupled singly to each of six intermolecular modes. The new dipole moment surface consists of a previous spectroscopically accurate 1-body dipole moment surface and a newly fitted ab initio intrinsic 2-body dipole moment. A detailed analysis of the new dipole moment surface in terms of the coordinate dependence of the effective atomic charges is done along with tests of it for the water dimer and prism hexamer double-harmonic spectra against direct ab initio calculations. The liquid configurations are taken from previous molecular dynamics calculations of Skinner and co-workers, using the TIP4P plus E3B rigid monomer water potential. The IR spectrum of water at 300 K in the range of 0-4000 cm-1 is calculated and compared with experiment, using the ab initio WHBB potential and new ab initio dipole moment, the q-TIP4P/F potential, which has a fixed-charged description of the dipole moment, and the TTM3-F potential and dipole moment surfaces. The newly calculated ab initio spectrum is in very good agreement with experiment throughout the above spectral range, both in band positions and intensities. This contrasts to results with the other potentials and dipole moments, especially the fixed-charge q-TIP4P/F model, which gives unrealistic intensities. The calculated ab initio spectrum is analyzed by examining the contribution of various transitions to each band.
Ab initio diffuse-interface model for lithiated electrode interface evolution
NASA Astrophysics Data System (ADS)
Stournara, Maria E.; Kumar, Ravi; Qi, Yue; Sheldon, Brian W.
2016-07-01
The study of chemical segregation at interfaces, and in particular the ability to predict the thickness of segregated layers via analytical expressions or computational modeling, is a fundamentally challenging topic in the design of novel heterostructured materials. This issue is particularly relevant for the phase-field (PF) methodology, which has become a prominent tool for describing phase transitions. These models rely on phenomenological parameters that pertain to the interfacial energy and thickness, quantities that cannot be experimentally measured. Instead of back-calculating these parameters from experimental data, here we combine a set of analytical expressions based on the Cahn-Hilliard approach with ab initio calculations to compute the gradient energy parameter κ and the thickness λ of the segregated Li layer at the LixSi-Cu interface. With this bottom-up approach we calculate the thickness λ of the Li diffuse interface to be on the order of a few nm, in agreement with prior experimental secondary ion mass spectrometry observations. Our analysis indicates that Li segregation is primarily driven by solution thermodynamics, while the strain contribution in this system is relatively small. This combined scheme provides an essential first step in the systematic evaluation of the thermodynamic parameters of the PF methodology, and we believe that it can serve as a framework for the development of quantitative interface models in the field of Li-ion batteries.
Ab initio diffuse-interface model for lithiated electrode interface evolution.
Stournara, Maria E; Kumar, Ravi; Qi, Yue; Sheldon, Brian W
2016-07-01
The study of chemical segregation at interfaces, and in particular the ability to predict the thickness of segregated layers via analytical expressions or computational modeling, is a fundamentally challenging topic in the design of novel heterostructured materials. This issue is particularly relevant for the phase-field (PF) methodology, which has become a prominent tool for describing phase transitions. These models rely on phenomenological parameters that pertain to the interfacial energy and thickness, quantities that cannot be experimentally measured. Instead of back-calculating these parameters from experimental data, here we combine a set of analytical expressions based on the Cahn-Hilliard approach with ab initio calculations to compute the gradient energy parameter κ and the thickness λ of the segregated Li layer at the Li_{x}Si-Cu interface. With this bottom-up approach we calculate the thickness λ of the Li diffuse interface to be on the order of a few nm, in agreement with prior experimental secondary ion mass spectrometry observations. Our analysis indicates that Li segregation is primarily driven by solution thermodynamics, while the strain contribution in this system is relatively small. This combined scheme provides an essential first step in the systematic evaluation of the thermodynamic parameters of the PF methodology, and we believe that it can serve as a framework for the development of quantitative interface models in the field of Li-ion batteries. PMID:27575197
Ab initio derivation of multi-orbital extended Hubbard model for molecular crystals
NASA Astrophysics Data System (ADS)
Tsuchiizu, Masahisa; Omori, Yukiko; Suzumura, Yoshikazu; Bonnet, Marie-Laure; Robert, Vincent
2012-01-01
From configuration interaction (CI) ab initio calculations, we derive an effective two-orbital extended Hubbard model based on the gerade (g) and ungerade (u) molecular orbitals (MOs) of the charge-transfer molecular conductor (TTM-TTP)I3 and the single-component molecular conductor [Au(tmdt)2]. First, by focusing on the isolated molecule, we determine the parameters for the model Hamiltonian so as to reproduce the CI Hamiltonian matrix. Next, we extend the analysis to two neighboring molecule pairs in the crystal and we perform similar calculations to evaluate the inter-molecular interactions. From the resulting tight-binding parameters, we analyze the band structure to confirm that two bands overlap and mix in together, supporting the multi-band feature. Furthermore, using a fragment decomposition, we derive the effective model based on the fragment MOs and show that the staking TTM-TTP molecules can be described by the zig-zag two-leg ladder with the inter-molecular transfer integral being larger than the intra-fragment transfer integral within the molecule. The inter-site interactions between the fragments follow a Coulomb law, supporting the fragment decomposition strategy.
Ab Initio Enhanced calphad Modeling of Actinide-Rich Nuclear Fuels
Morgan, Dane; Yang, Yong Austin
2013-10-28
The process of fuel recycling is central to the Advanced Fuel Cycle Initiative (AFCI), where plutonium and the minor actinides (MA) Am, Np, and Cm are extracted from spent fuel and fabricated into new fuel for a fast reactor. Metallic alloys of U-Pu-Zr-MA are leading candidates for fast reactor fuels and are the current basis for fast spectrum metal fuels in a fully recycled closed fuel cycle. Safe and optimal use of these fuels will require knowledge of their multicomponent phase stability and thermodynamics (Gibbs free energies). In additional to their use as nuclear fuels, U-Pu-Zr-MA contain elements and alloy phases that pose fundamental questions about electronic structure and energetics at the forefront of modern many-body electron theory. This project will validate state-of-the-art electronic structure approaches for these alloys and use the resulting energetics to model U-Pu-Zr-MA phase stability. In order to keep the work scope practical, researchers will focus on only U-Pu-Zr-{Np,Am}, leaving Cm for later study. The overall objectives of this project are to: Provide a thermodynamic model for U-Pu-Zr-MA for improving and controlling reactor fuels; and, Develop and validate an ab initio approach for predicting actinide alloy energetics for thermodynamic modeling.
Ab initio kinetics and thermal decomposition mechanism of mononitrobiuret and 1,5-dinitrobiuret
NASA Astrophysics Data System (ADS)
Sun, Hongyan; Vaghjiani, Ghanshyam L.
2015-05-01
Mononitrobiuret (MNB) and 1,5-dinitrobiuret (DNB) are tetrazole-free, nitrogen-rich, energetic compounds. For the first time, a comprehensive ab initio kinetics study on the thermal decomposition mechanisms of MNB and DNB is reported here. In particular, the intramolecular interactions of amine H-atom with electronegative nitro O-atom and carbonyl O-atom have been analyzed for biuret, MNB, and DNB at the M06-2X/aug-cc-pVTZ level of theory. The results show that the MNB and DNB molecules are stabilized through six-member-ring moieties via intramolecular H-bonding with interatomic distances between 1.8 and 2.0 Å, due to electrostatic as well as polarization and dispersion interactions. Furthermore, it was found that the stable molecules in the solid state have the smallest dipole moment amongst all the conformers in the nitrobiuret series of compounds, thus revealing a simple way for evaluating reactivity of fuel conformers. The potential energy surface for thermal decomposition of MNB was characterized by spin restricted coupled cluster theory at the RCCSD(T)/cc-pV∞ Z//M06-2X/aug-cc-pVTZ level. It was found that the thermal decomposition of MNB is initiated by the elimination of HNCO and HNN(O)OH intermediates. Intramolecular transfer of a H-atom, respectively, from the terminal NH2 group to the adjacent carbonyl O-atom via a six-member-ring transition state eliminates HNCO with an energy barrier of 35 kcal/mol and from the central NH group to the adjacent nitro O-atom eliminates HNN(O)OH with an energy barrier of 34 kcal/mol. Elimination of HNN(O)OH is also the primary process involved in the thermal decomposition of DNB, which processes C2v symmetry. The rate coefficients for the primary decomposition channels for MNB and DNB were quantified as functions of temperature and pressure. In addition, the thermal decomposition of HNN(O)OH was analyzed via Rice-Ramsperger-Kassel-Marcus/multi-well master equation simulations, the results of which reveal the formation
Ab initio kinetics and thermal decomposition mechanism of mononitrobiuret and 1,5-dinitrobiuret.
Sun, Hongyan; Vaghjiani, Ghanshyam L
2015-05-28
Mononitrobiuret (MNB) and 1,5-dinitrobiuret (DNB) are tetrazole-free, nitrogen-rich, energetic compounds. For the first time, a comprehensive ab initio kinetics study on the thermal decomposition mechanisms of MNB and DNB is reported here. In particular, the intramolecular interactions of amine H-atom with electronegative nitro O-atom and carbonyl O-atom have been analyzed for biuret, MNB, and DNB at the M06-2X/aug-cc-pVTZ level of theory. The results show that the MNB and DNB molecules are stabilized through six-member-ring moieties via intramolecular H-bonding with interatomic distances between 1.8 and 2.0 Å, due to electrostatic as well as polarization and dispersion interactions. Furthermore, it was found that the stable molecules in the solid state have the smallest dipole moment amongst all the conformers in the nitrobiuret series of compounds, thus revealing a simple way for evaluating reactivity of fuel conformers. The potential energy surface for thermal decomposition of MNB was characterized by spin restricted coupled cluster theory at the RCCSD(T)/cc-pV∞ Z//M06-2X/aug-cc-pVTZ level. It was found that the thermal decomposition of MNB is initiated by the elimination of HNCO and HNN(O)OH intermediates. Intramolecular transfer of a H-atom, respectively, from the terminal NH2 group to the adjacent carbonyl O-atom via a six-member-ring transition state eliminates HNCO with an energy barrier of 35 kcal/mol and from the central NH group to the adjacent nitro O-atom eliminates HNN(O)OH with an energy barrier of 34 kcal/mol. Elimination of HNN(O)OH is also the primary process involved in the thermal decomposition of DNB, which processes C2v symmetry. The rate coefficients for the primary decomposition channels for MNB and DNB were quantified as functions of temperature and pressure. In addition, the thermal decomposition of HNN(O)OH was analyzed via Rice-Ramsperger-Kassel-Marcus/multi-well master equation simulations, the results of which reveal the formation
Ab Initio Kinetics and Thermal Decomposition Mechanism of Mononitrobiuret and 1,5- Dinitrobiuret
Sun, Hongyan; Vaghjiani, Ghanshyam G.
2015-05-26
Mononitrobiuret (MNB) and 1,5-dinitrobiuret (DNB) are tetrazole-free, nitrogen-rich, energetic compounds. For the first time, a comprehensive ab initio kinetics study on the thermal decomposition mechanisms of MNB and DNB is reported here. In particular, the intramolecular interactions of amine H-atom with electronegative nitro O-atom and carbonyl O-atom have been analyzed for biuret, MNB, and DNB at the M06-2X/aug-cc-pVTZ level of theory. The results show that the MNB and DNB molecules are stabilized through six-member-ring moieties via intramolecular H-bonding with interatomic distances between 1.8 and 2.0 Å, due to electrostatic as well as polarization and dispersion interactions. Furthermore, it was found that the stable molecules in the solid state have the smallest dipole moment amongst all the conformers in the nitrobiuret series of compounds, thus revealing a simple way for evaluating reactivity of fuel conformers. The potential energy surface for thermal decomposition of MNB was characterized by spin restricted coupled cluster theory at the RCCSD(T)/cc-pV∞ Z//M06-2X/aug-cc-pVTZ level. It was found that the thermal decomposition of MNB is initiated by the elimination of HNCO and HNN(O)OH intermediates. Intramolecular transfer of a H-atom, respectively, from the terminal NH2 group to the adjacent carbonyl O-atom via a six-member-ring transition state eliminates HNCO with an energy barrier of 35 kcal/mol and from the central NH group to the adjacent nitro O-atom eliminates HNN(O)OH with an energy barrier of 34 kcal/mol. Elimination of HNN(O)OH is also the primary process involved in the thermal decomposition of DNB, which processes C2v symmetry. The rate coefficients for the primary decomposition channels for MNB and DNB were quantified as functions of temperature and pressure. In addition, the thermal decomposition of HNN(O)OH was analyzed via Rice–Ramsperger–Kassel–Marcus/multi-well master equation simulations, the results of which reveal the
Köster, Andreas; Spura, Thomas; Rutkai, Gábor; Kessler, Jan; Wiebeler, Hendrik; Vrabec, Jadran; Kühne, Thomas D
2016-07-15
The accuracy of water models derived from ab initio molecular dynamics simulations by means on an improved force-matching scheme is assessed for various thermodynamic, transport, and structural properties. It is found that although the resulting force-matched water models are typically less accurate than fully empirical force fields in predicting thermodynamic properties, they are nevertheless much more accurate than generally appreciated in reproducing the structure of liquid water and in fact superseding most of the commonly used empirical water models. This development demonstrates the feasibility to routinely parametrize computationally efficient yet predictive potential energy functions based on accurate ab initio molecular dynamics simulations for a large variety of different systems. © 2016 Wiley Periodicals, Inc. PMID:27232117
NASA Astrophysics Data System (ADS)
Neukirch, Amanda; Nei, Wanyi; Pedesseau, Laurent; Even, Jacky; Katan, Claudine; Mohite, Aditya; Tretiak, Segrei
2015-03-01
The need for an inexpensive, clean, and plentiful source of energy has generated large amounts of research in an assortment of solution processed organic and hybrid organic-inorganic solar cells. A relative newcomer to the field of solution processed photovoltaics is the lead halide perovskite solar cell. In the past 5 years, the efficiencies of devices made from this material have increased from 3.5% to nearly 20%. Despite the rapid development of organic-inorganic perovskite solar cells, a thorough understanding of the fundamental photophysical processes driving the high performance of these devices is not well understood. I am using state-of-the-art ab initio computational techniques in order to characterize the properties at the interface of perovskite devices in order to aide in materials design and device engineering. I will present an in-depth analysis of the electronic and optical properties of bulk and surface states of pure and mixed halide systems. The high-level static quantum mechanical calculations, including spin-orbit-coupling and the many body GW approach, identify the key electronic states involved in photoinduced dynamics. This knowledge provides important information on how the optical properties change with variations to the system. Supported by the DOE, the LANL LDRD program XW11, and CNLS.
A Simple ab Initio Model for the Hydrated Electron That Matches Experiment.
Kumar, Anil; Walker, Jonathan A; Bartels, David M; Sevilla, Michael D
2015-08-27
Since its discovery over 50 years ago, the "structure" and properties of the hydrated electron have been a subject for wonderment and also fierce debate. In the present work we seriously explore a minimal model for the aqueous electron, consisting of a small water anion cluster embedded in a polarized continuum, using several levels of ab initio calculation and basis set. The minimum energy "zero Kelvin" structure found for any 4-water (or larger) anion cluster, at any post-Hartree–Fock theory level, is very similar to a recently reported embedded-DFT-in-classical-water-MD simulation (Uhlig, Marsalek, and Jungwirth, J. Phys. Chem. Lett. 2012, 3, 3071−3075), with four OH bonds oriented toward the maximum charge density in a small central "void". The minimum calculation with just four water molecules does a remarkably good job of reproducing the resonance Raman properties, the radius of gyration derived from the optical spectrum, the vertical detachment energy, and the hydration free energy. For the first time we also successfully calculate the EPR g-factor and (low temperature ice) hyperfine couplings. The simple tetrahedral anion cluster model conforms very well to experiment, suggesting it does in fact represent the dominant structural motif of the hydrated electron. PMID:26275103
Tunneling of electrons via rotor-stator molecular interfaces: Combined ab initio and model study
NASA Astrophysics Data System (ADS)
Petreska, Irina; Ohanesjan, Vladimir; Pejov, Ljupčo; Kocarev, Ljupčo
2016-07-01
Tunneling of electrons through rotor-stator anthracene aldehyde molecular interfaces is studied with a combined ab initio and model approach. Molecular electronic structure calculated from first principles is utilized to model different shapes of tunneling barriers. Together with a rectangular barrier, we also consider a sinusoidal shape that captures the effects of the molecular internal structure more realistically. Quasiclassical approach with the Simmons' formula for current density is implemented. Special attention is paid on conformational dependence of the tunneling current. Our results confirm that the presence of the side aldehyde group enhances the interesting electronic properties of the pure anthracene molecule, making it a bistable system with geometry dependent transport properties. We also investigate the transition voltage and we show that conformation-dependent field emission could be observed in these molecular interfaces at realistically low voltages. The present study accompanies our previous work where we investigated the coherent transport via strongly coupled delocalized orbital by application of Non-equilibrium Green's Function Formalism.
Determination of a silane intermolecular force field potential model from an ab initio calculation
Li, Arvin Huang-Te; Chao, Sheng D.; Chang, Chien-Cheng
2010-12-15
Intermolecular interaction potentials of the silane dimer in 12 orientations have been calculated by using the Hartree-Fock (HF) self-consistent theory and the second-order Moeller-Plesset (MP2) perturbation theory. We employed basis sets from Pople's medium-size basis sets [up to 6-311++G(3df, 3pd)] and Dunning's correlation consistent basis sets (up to the triply augmented correlation-consistent polarized valence quadruple-zeta basis set). We found that the minimum energy orientations were the G and H conformers. We have suggested that the Si-H attractions, the central silicon atom size, and electronegativity play essential roles in weakly binding of a silane dimer. The calculated MP2 potential data were employed to parametrize a five-site force field for molecular simulations. The Si-Si, Si-H, and H-H interaction parameters in a pairwise-additive, site-site potential model for silane molecules were regressed from the ab initio energies.
Ab-initio molecular modeling of interfaces in tantalum-carbon system
Balani, Kantesh; Mungole, Tarang; Bakshi, Srinivasa Rao; Agarwal, Arvind
2012-03-15
Processing of ultrahigh temperature TaC ceramic material with sintering additives of B{sub 4}C and reinforcement of carbon nanotubes (CNTs) gives rise to possible formation of several interfaces (Ta{sub 2}C-TaC, TaC-CNT, Ta{sub 2}C-CNT, TaB{sub 2}-TaC, and TaB{sub 2}-CNT) that could influence the resultant properties. Current work focuses on interfaces developed during spark plasma sintering of TaC-system and performing ab initio molecular modeling of the interfaces generated during processing of TaC-B{sub 4}C and TaC-CNT composites. The energy of the various interfaces has been evaluated and compared with TaC-Ta{sub 2}C interface. The iso-surface electronic contours are extracted from the calculations eliciting the enhanced stability of TaC-CNT interface by 72.2%. CNTs form stable interfaces with Ta{sub 2}C and TaB{sub 2} phases with a reduction in the energy by 35.8% and 40.4%, respectively. The computed Ta-C-B interfaces are also compared with experimentally observed interfaces in high resolution TEM images.
Beyond Born-Mayer: Improved Models for Short-Range Repulsion in ab Initio Force Fields.
Van Vleet, Mary J; Misquitta, Alston J; Stone, Anthony J; Schmidt, J R
2016-08-01
Short-range repulsion within intermolecular force fields is conventionally described by either Lennard-Jones (A/r(12)) or Born-Mayer (A exp(-Br)) forms. Despite their widespread use, these simple functional forms are often unable to describe the interaction energy accurately over a broad range of intermolecular distances, thus creating challenges in the development of ab initio force fields and potentially leading to decreased accuracy and transferability. Herein, we derive a novel short-range functional form based on a simple Slater-like model of overlapping atomic densities and an iterated stockholder atom (ISA) partitioning of the molecular electron density. We demonstrate that this Slater-ISA methodology yields a more accurate, transferable, and robust description of the short-range interactions at minimal additional computational cost compared to standard Lennard-Jones or Born-Mayer approaches. Finally, we show how this methodology can be adapted to yield the standard Born-Mayer functional form while still retaining many of the advantages of the Slater-ISA approach. PMID:27337546
NASA Astrophysics Data System (ADS)
Vitali, Ettore; Shi, Hao; Qin, Mingpu; Zhang, Shiwei
The possibility of calculating dynamical correlation functions from first principles provides a unique opportunity to explore the manifold of the excited states of a quantum many-body system. Such calculations allow us to predict interesting physical properties like spectral functions, excitation spectra and charge and spin gaps, which are more difficult to access from usual equilibrium calculations. We address the ab-initio calculation of dynamical Green functions and two-body correlation functions in the Auxiliary-field Quantum Monte Carlo method, using the two-dimensional Hubbard model as an example. When the sign problem is not present, an unbiased estimation of imaginary time correlation functions is obtained. We discuss in detail the complexity and the stability of the calculations. Moreover, we propose a new approach which is expected to be very useful when dealing with dilute systems, e.g. for cold gases, allowing calculations with a very favorable complexity in the system size. Supported by NSF, DOE SciDAC, and Simons Foundation.
Mei, J.; Cooper, B.R.; Hao, Y.G.; Scoy, F.L. Van
1994-12-31
Molecular dynamics simulations have been performed to study thermal expansions of Ni-rich (fcc structure) Ni/Cr alloys (which serve as the basis for practical superalloy systems). This has been done using ab initio interatomic potentials with no experimental input. The coefficient of thermal expansion (CTE) as a function of temperature has been calculated. By admixing Re and Me atoms into fee Ni and the fee alloy system Ni/Cr, additive effects on the thermal expansion have been predicted. While addition of Cr lowers the CTE of Ni, and moderate addition of Mo lowers the CTE of Ni over a wide temperature range, moderate addition of Re raises the CTE of both Ni and Ni/Cr alloys over a significant temperature range. An explanation for the contrasting effect of additive Re on the CTE, based on a one-dimensional atomic chain model, is that the trade-off, between atomic volume effects increasing the CTE over that of pure Ni and pair-potential effects (exemplified by the Grueneisen parameter) decreasing the CTE from that of pure nickel, changes for Re compared to Cr and Mo.
Ab initio molecular dynamics simulation of proton hopping in a model polymer membrane.
Devanathan, Ram; Idupulapati, Nagesh; Baer, Marcel D; Mundy, Christopher J; Dupuis, Michel
2013-12-27
We report the results of ab initio molecular dynamics simulations of a model Nafion polymer membrane initially equilibrated using classical molecular dynamics simulations. We studied three hydration levels (λ) of 3, 9, and 15 H2O/SO3(-) corresponding to dry, hydrated, and saturated fuel cell membrane, respectively. The barrier for proton transfer from the SO3(-)-H3O(+) contact ion pair to a solvent-separated ion pair decreased from 2.3 kcal/mol for λ = 3 to 0.8 kcal/mol for λ = 15. The barrier for proton transfer between two water molecules was in the range from 0.7 to 0.8 kcal/mol for the λ values studied. The number of proton shuttling events between a pair of water molecules is an order of magnitude more than the number of proton hops across three distinct water molecules. The proton diffusion coefficient at λ = 15 is about 0.9 × 10(-5) cm(2)/s, which is in good agreement with experiment and our previous quantum hopping molecular dynamics simulations. PMID:24320080
Beyond Born-Mayer: Improved models for short-range repulsion in ab initio force fields
Van Vleet, Mary J.; Misquitta, Alston J.; Stone, Anthony J.; Schmidt, Jordan R.
2016-06-23
Short-range repulsion within inter-molecular force fields is conventionally described by either Lennard-Jones or Born-Mayer forms. Despite their widespread use, these simple functional forms are often unable to describe the interaction energy accurately over a broad range of inter-molecular distances, thus creating challenges in the development of ab initio force fields and potentially leading to decreased accuracy and transferability. Herein, we derive a novel short-range functional form based on a simple Slater-like model of overlapping atomic densities and an iterated stockholder atom (ISA) partitioning of the molecular electron density. We demonstrate that this Slater-ISA methodology yields a more accurate, transferable, andmore » robust description of the short-range interactions at minimal additional computational cost compared to standard Lennard-Jones or Born-Mayer approaches. Lastly, we show how this methodology can be adapted to yield the standard Born-Mayer functional form while still retaining many of the advantages of the Slater-ISA approach.« less
A nonlocal, ab initio model of dissociative electron attachment and vibrational excitation of NO
Trevisan, Cynthia S.; Houfek, Karel; Zhang, Zhiyong; Orel, Ann E.; McCurdy, C. William; Rescigno, Thomas N.
2005-02-01
We present the results of an ab initio study of elastic scattering and vibrational excitation of NO by electron impact in the low-energy (0-2 eV) region where the cross sections are dominated by resonance contributions. The 3Sigma-, 1Delta and 1Sigma+ NO- resonance lifetimes are taken from our earlier study [Phys. Rev. A 69, 062711 (2004)], but the resonance energies used here are obtained from new configuration-interaction studies. Here we employ a more elaborate nonlocal treatment of the nuclear dynamics, which is found to remedy the principal deficiencies of the local complex potential model we employed in our earlier study, and gives cross sections in better agreement with the most recent experiments. We also present cross sections for dissociative electron attachment to NO leading to groundstate products. The calculations show that, while the peak cross sections starting from NO in its ground vibrational state are very small, the cross sections are extremely sensitive to vibrational excitation of the target and should be readily observable for target NO molecules excited to v = 10 and above.
Cluster form factor calculation in the ab initio no-core shell model
Navratil, Petr
2004-11-01
We derive expressions for cluster overlap integrals or channel cluster form factors for ab initio no-core shell model (NCSM) wave functions. These are used to obtain the spectroscopic factors and can serve as a starting point for the description of low-energy nuclear reactions. We consider the composite system and the target nucleus to be described in the Slater determinant (SD) harmonic oscillator (HO) basis while the projectile eigenstate to be expanded in the Jacobi coordinate HO basis. This is the most practical case. The spurious center of mass components present in the SD bases are removed exactly. The calculated cluster overlap integrals are translationally invariant. As an illustration, we present results of cluster form factor calculations for <{sup 5}He vertical bar{sup 4}He+n>, <{sup 5}He vertical bar{sup 3}H+d>, <{sup 6}Li vertical bar{sup 4}He+d>, <{sup 6}Be vertical bar{sup 3}He+{sup 3}He>, <{sup 7}Li vertical bar{sup 4}He+{sup 3}H>, <{sup 7}Li vertical bar{sup 6}Li+n>, <{sup 8}Be vertical bar{sup 6}Li+d>, <{sup 8}Be vertical bar{sup 7}Li+p>, <{sup 9}Li vertical bar{sup 8}Li+n>, and <{sup 13}C vertical bar{sup 12}C+n>, with all the nuclei described by multi-({Dirac_h}/2{pi}){omega} NCSM wave functions.
JUPITER MODELS WITH IMPROVED AB INITIO HYDROGEN EQUATION OF STATE (H-REOS.2)
Nettelmann, N.; Becker, A.; Redmer, R.; Holst, B.
2012-05-01
The amount and distribution of heavy elements in Jupiter gives indications on the process of its formation and evolution. Core mass and metallicity predictions, however, depend on the equations of state (EOSs) used and on model assumptions. We present an improved ab initio hydrogen EOS, H-REOS.2, and compute the internal structure and thermal evolution of Jupiter within the standard three-layer approach. The advance over our previous Jupiter models with H-REOS.1 by Nettelmann et al. is that the new models are also consistent with the observed {approx}> 2 times solar heavy element abundances in Jupiter's atmosphere. Such models have a rock core mass M{sub c} = 0-8 M{sub Circled-Plus }, total mass of heavy elements M{sub Z} = 28-32 M{sub Circled-Plus }, a deep internal layer boundary at {>=}4 Mbar, and a cooling time of 4.4-5.0 Gyr when assuming homogeneous evolution. We also calculate two-layer models in the manner of Militzer et al. and find a comparable large core of 16-21 M{sub Circled-Plus }, out of which {approx}11 M{sub Circled-Plus} is helium, but a significantly higher envelope metallicity of 4.5 times solar. According to our preferred three-layer models, neither the characteristic frequency ({nu}{sub 0} {approx} 156 {mu}Hz) nor the normalized moment of inertia ({lambda} {approx}0.276) is sensitive to the core mass but accurate measurements could well help to rule out some classes of models.
NASA Astrophysics Data System (ADS)
Aryal, Sita Ram
The alumino-silicate solid solution series (Al 4+2xSi2-2 xO10-x) is an important class of ceramics. Except for the end member (x=0), Al2 SiO5 the crystal structures of the other phases, called mullite, have partially occupied sites. Stoichiometric supercell models for the four mullite phases 3Al2O 3 · 2SiO2 · 2Al 2O3 · SiO2, 4 Al2O3· SiO 2, 9Al2O3 · SiO2, and iota-Al2 O3 (iota-alumina) are constructed starting from experimentally reported crystal structures. A large number of models were built for each phase and relaxed using the Vienna ab initio simulation package (VASP) program. The model with the lowest total energy for a given x was chosen as the representative structure for that phase. Electronic structure and mechanical properties of mullite phases were studied via first-principles calculations. Of the various phases of transition alumina, iota-Al 2O3 is the least well known. In addition structural details have not, until now, been available. It is the end member of the aluminosilicate solid solution series with x=1. Based on a high alumina content mullite phase, a structural model for iota- Al2O3 is constructed. The simulated x-ray diffraction (XRD) pattern of this model agrees well with a measured XRD pattern. The iota-Al2 O3 is a highly disordered ultra-low-density phase of alumina with a theoretical density of 2854kg/m3. Using this theoretically constructed model, elastic, thermodynamic, electronic, and spectroscopic properties of iota-Al2 O3 have been calculated and compared it with those of alpha- Al2O3 and gamma- Al2O3. Boron carbide (B4C) undergoes an amorphization under high velocity impacts. The mechanism of amorphization is not clear. Ab initio methods are used to carry out large-scale uniaxial compression simulations on two polytypes of stoichiometric boron carbide (B4C), B 11C-CBC, and B12- CCC where B11C or B12 is the 12-atom icosahedron and CBC or CCC is the three-atom chain. The simulations were performed on large supercells of 180 atoms
The ab initio model potential method. Second series transition metal elements
NASA Astrophysics Data System (ADS)
Barandiarán, Zoila; Seijo, Luis; Huzinaga, Sigeru
1990-10-01
The ab initio core method potential model (AIMP) has already been presented in its nonrelativistic version and applied to the main group and first series transition metal elements [J. Chem. Phys. 86, 2132 (1987); 91, 7011 (1989)]. In this paper we extend the AIMP method to include relativistic effects within the Cowan-Griffin approximation and we present relativistic Zn-like core model potentials and valence basis sets, as well as their nonrelativistic Zn-like core and Kr-like core counterparts. The pilot molecular calculations on YO, TcO, AgO, and AgH reveal that the 4p orbital is indeed a core orbital only at the end part of the series, whereas the 4s orbital can be safely frozen from Y to Cd. The all-electron and model potential results agree in 0.01-0.02 Å in Re and 25-50 cm-1 in ν¯e if the same type of valence part of the basis set is used. The comparison of the relativistic results on AgH with those of the all-electron Dirac-Fock calculations by Lee and McLean is satisfactory: the absolute value of Re is reproduced within the 0.01 Å margin and the relativistic contraction of 0.077 Å is also very well reproduced (0.075 Å). Finally, the relative magnitude of the effects of the core orbital change, mass-velocity potential, and Darwin potential on the net relativistic effects are analyzed in the four molecules studied.
Ab initio kinetics and thermal decomposition mechanism of mononitrobiuret and 1,5-dinitrobiuret
Sun, Hongyan E-mail: ghanshyam.vaghjiani@us.af.mil; Vaghjiani, Ghanshyam L. E-mail: ghanshyam.vaghjiani@us.af.mil
2015-05-28
Mononitrobiuret (MNB) and 1,5-dinitrobiuret (DNB) are tetrazole-free, nitrogen-rich, energetic compounds. For the first time, a comprehensive ab initio kinetics study on the thermal decomposition mechanisms of MNB and DNB is reported here. In particular, the intramolecular interactions of amine H-atom with electronegative nitro O-atom and carbonyl O-atom have been analyzed for biuret, MNB, and DNB at the M06-2X/aug-cc-pVTZ level of theory. The results show that the MNB and DNB molecules are stabilized through six-member-ring moieties via intramolecular H-bonding with interatomic distances between 1.8 and 2.0 Å, due to electrostatic as well as polarization and dispersion interactions. Furthermore, it was found that the stable molecules in the solid state have the smallest dipole moment amongst all the conformers in the nitrobiuret series of compounds, thus revealing a simple way for evaluating reactivity of fuel conformers. The potential energy surface for thermal decomposition of MNB was characterized by spin restricted coupled cluster theory at the RCCSD(T)/cc-pV∞ Z//M06-2X/aug-cc-pVTZ level. It was found that the thermal decomposition of MNB is initiated by the elimination of HNCO and HNN(O)OH intermediates. Intramolecular transfer of a H-atom, respectively, from the terminal NH{sub 2} group to the adjacent carbonyl O-atom via a six-member-ring transition state eliminates HNCO with an energy barrier of 35 kcal/mol and from the central NH group to the adjacent nitro O-atom eliminates HNN(O)OH with an energy barrier of 34 kcal/mol. Elimination of HNN(O)OH is also the primary process involved in the thermal decomposition of DNB, which processes C{sub 2v} symmetry. The rate coefficients for the primary decomposition channels for MNB and DNB were quantified as functions of temperature and pressure. In addition, the thermal decomposition of HNN(O)OH was analyzed via Rice–Ramsperger–Kassel–Marcus/multi-well master equation simulations, the results of which
NASA Astrophysics Data System (ADS)
Akin-Ojo, Omololu; Song, Yang; Wang, Feng
2008-08-01
A new method called adaptive force matching (AFM) has been developed that is capable of producing high quality force fields for condensed phase simulations. This procedure involves the parametrization of force fields to reproduce ab initio forces obtained from condensed phase quantum-mechanics/molecular-mechanics (QM/MM) calculations. During the procedure, the MM part of the QM/MM is iteratively improved so as to approach ab initio quality. In this work, the AFM method has been tested to parametrize force fields for liquid water so that the resulting force fields reproduce forces calculated using the ab initio MP2 and the Kohn-Sham density functional theory with the Becke-Lee-Yang-Parr (BLYP) and Becke three-parameter LYP (B3LYP) exchange correlation functionals. The AFM force fields generated in this work are very simple to evaluate and are supported by most molecular dynamics (MD) codes. At the same time, the quality of the forces predicted by the AFM force fields rivals that of very expensive ab initio calculations and are found to successfully reproduce many experimental properties. The site-site radial distribution functions (RDFs) obtained from MD simulations using the force field generated from the BLYP functional through AFM compare favorably with the previously published RDFs from Car-Parrinello MD simulations with the same functional. Technical aspects of AFM such as the optimal QM cluster size, optimal basis set, and optimal QM method to be used with the AFM procedure are discussed in this paper.
Ab initio model of salicylate adsorbed onto Al{sub 2}O{sub 3} and illite clay
Kubicki, J.D.; Itoh, M.J.; Apitz, S.E.
1996-10-01
Organic-mineral surface chemistry plays a significant role in numerous geochemical processes such as global carbon cycling, weathering, and contaminant fate and transport. Knowledge of bonding mechanisms between naturally-occurring organic matter (NOM) and minerals is necessary in environmental science. This research examines surface complexation of salicylic acid (which is often used as an analog for NOM) adsorbed onto Al{sub 2}O{sub 3} and illite. ATR-FTIR spectra of the adsorbed complexes were measured and compared to theoretical vibrational spectra of possible surface configurations derived form molecular orbital (MO) calculations. A variety of Al- and Si-salicylate complexes were modeled with ab initio MO calculations. The theoretical vibrational spectrum that best fits the observed spectra corresponds to a salicylate anion bonded to an octahedral Al{sup 3+} ion via a C-O-Al ester-type linkage. These results support the configuration proposed in Biber and Stumm for salicylate adsorbed onto Al{sub 2}O{sub 3}.
Ab initio joint density-functional theory of solvated electrodes, with model and explicit solvation
NASA Astrophysics Data System (ADS)
Arias, Tomas
2015-03-01
First-principles guided design of improved electrochemical systems has the potential for great societal impact by making non-fossil-fuel systems economically viable. Potential applications include improvements in fuel-cells, solar-fuel systems (``artificial photosynthesis''), supercapacitors and batteries. Economical fuel-cell systems would enable zero-carbon footprint transportation, solar-fuel systems would directly convert sunlight and water into hydrogen fuel for such fuel-cell vehicles, supercapacitors would enable nearly full recovery of energy lost during vehicle braking thus extending electric vehicle range and acceptance, and economical high-capacity batteries would be central to mitigating the indeterminacy of renewable resources such as wind and solar. Central to the operation of all of the above electrochemical systems is the electrode-electrolyte interface, whose underlying physics is quite rich, yet remains remarkably poorly understood. The essential underlying technical challenge to the first principles studies which could explore this physics is the need to properly represent simultaneously both the interaction between electron-transfer events at the electrode, which demand a quantum mechanical description, and multiscale phenomena in the liquid environment such as the electrochemical double layer (ECDL) and its associated shielding, which demand a statistical description. A direct ab initio approach to this challenge would, in principle, require statistical sampling and thousands of repetitions of already computationally demanding quantum mechanical calculations. This talk will begin with a brief review of a recent advance, joint density-functional theory (JDFT), which allows for a fully rigorous and, in principle, exact representation of the thermodynamic equilibrium between a system described at the quantum-mechanical level and a liquid environment, but without the need for costly sampling. We then shall demonstrate how this approach applies in
Timoshenko, J.; Shivhare, A.; Scott, R. W.; Lu, D.; Frenkel, A. I.
2016-06-30
We adopted ab-initio X-ray Absorption Near Edge Structure (XANES) modelling for structural refinement of local environments around metal impurities in a large variety of materials. Our method enables both direct modelling, where the candidate structures are known, and the inverse modelling, where the unknown structural motifs are deciphered from the experimental spectra. We present also estimates of systematic errors, and their influence on the stability and accuracy of the obtained results. We illustrate our approach by following the evolution of local environment of palladium atoms in palladium-doped gold thiolate clusters upon chemical and thermal treatments.
NASA Astrophysics Data System (ADS)
Liang, Wenkel
This dissertation consists of two general parts: (I) developments of optimization algorithms (both nuclear and electronic degrees of freedom) for time-independent molecules and (II) novel methods, first-principle theories and applications in time dependent molecular structure modeling. In the first part, we discuss in specific two new algorithms for static geometry optimization, the eigenspace update (ESU) method in nonredundant internal coordinate that exhibits an enhanced performace with up to a factor of 3 savings in computational cost for large-sized molecular systems; the Car-Parrinello density matrix search (CP-DMS) method that enables direct minimization of the SCF energy as an effective alternative to conventional diagonalization approach. For the second part, we consider the time dependence and first presents two nonadiabatic dynamic studies that model laser controlled molecular photo-dissociation for qualitative understandings of intense laser-molecule interaction, using ab initio direct Ehrenfest dynamics scheme implemented with real-time time-dependent density functional theory (RT-TDDFT) approach developed in our group. Furthermore, we place our special interest on the nonadiabatic electronic dynamics in the ultrafast time scale, and presents (1) a novel technique that can not only obtain energies but also the electron densities of doubly excited states within a single determinant framework, by combining methods of CP-DMS with RT-TDDFT; (2) a solvated first-principles electronic dynamics method by incorporating the polarizable continuum solvation model (PCM) to RT-TDDFT, which is found to be very effective in describing the dynamical solvation effect in the charge transfer process and yields a consistent absorption spectrum in comparison to the conventional linear response results in solution. (3) applications of the PCM-RT-TDDFT method to study the intramolecular charge-transfer (CT) dynamics in a C60 derivative. Such work provides insights into the
NASA Astrophysics Data System (ADS)
Ran, Shi-Ju
2016-05-01
In this work, a simple and fundamental numeric scheme dubbed as ab initio optimization principle (AOP) is proposed for the ground states of translational invariant strongly correlated quantum lattice models. The idea is to transform a nondeterministic-polynomial-hard ground-state simulation with infinite degrees of freedom into a single optimization problem of a local function with finite number of physical and ancillary degrees of freedom. This work contributes mainly in the following aspects: (1) AOP provides a simple and efficient scheme to simulate the ground state by solving a local optimization problem. Its solution contains two kinds of boundary states, one of which play the role of the entanglement bath that mimics the interactions between a supercell and the infinite environment, and the other gives the ground state in a tensor network (TN) form. (2) In the sense of TN, a novel decomposition named as tensor ring decomposition (TRD) is proposed to implement AOP. Instead of following the contraction-truncation scheme used by many existing TN-based algorithms, TRD solves the contraction of a uniform TN in an opposite way by encoding the contraction in a set of self-consistent equations that automatically reconstruct the whole TN, making the simulation simple and unified; (3) AOP inherits and develops the ideas of different well-established methods, including the density matrix renormalization group (DMRG), infinite time-evolving block decimation (iTEBD), network contractor dynamics, density matrix embedding theory, etc., providing a unified perspective that is previously missing in this fields. (4) AOP as well as TRD give novel implications to existing TN-based algorithms: A modified iTEBD is suggested and the two-dimensional (2D) AOP is argued to be an intrinsic 2D extension of DMRG that is based on infinite projected entangled pair state. This paper is focused on one-dimensional quantum models to present AOP. The benchmark is given on a transverse Ising
Combined electron beam imaging and ab initio modeling of T{sub 1} precipitates in Al-Li-Cu alloys
Dwyer, C.; Weyland, M.; Chang, L. Y.; Muddle, B. C.
2011-05-16
Among the many considerable challenges faced in developing a rational basis for advanced alloy design, establishing accurate atomistic models is one of the most fundamental. Here we demonstrate how advanced imaging techniques in a double-aberration-corrected transmission electron microscope, combined with ab initio modeling, have been used to determine the atomic structure of embedded 1 nm thick T{sub 1} precipitates in precipitation-hardened Al-Li-Cu aerospace alloys. The results provide an accurate determination of the controversial T{sub 1} structure, and demonstrate how next-generation techniques permit the characterization of embedded nanostructures in alloys and other nanostructured materials.
NASA Astrophysics Data System (ADS)
Ansari, R.; Faghihnasiri, M.; Malakpour, S.; Sahmani, S.
2015-07-01
In the current investigation, ab initio calculations are performed to explore the influence of electric field on the mechanical properties of bilayer boron nitride with AB stacking order (AB-2LBN). To accomplish this, density functional theory (DFT) within the framework of generalized gradient approximation (GGA) is implemented. It is demonstrated that the electric field has significant effects on Young's modulus and Poisson's ratio of AB-2LBN when its magnitude is small. With increasing the magnitude of electric field, these effects diminish so that the mechanical properties with and without considering the electric field become approximately identical. Also, it is shown that the equilibrium strain energy decreases linearly by increasing the magnitude of applied electric field.
Plotnikov, Nikolay; Kamerlin, Shina Caroline Lynn; Warshel, Arieh
2011-01-01
Recent years have seen tremendous effort in the development of approaches with which to obtain quantum mechanics/molecular mechanics (QM/MM) free energies for reactions in the condensed phase. Nevertheless, there remain significant challenges to address, particularly the high computational cost involved in performing proper configurational sampling and in particular in obtaining ab initio QM/MM (QM(ai)/MM) free energy surfaces. One increasingly popular approach that seems to offer an ideal way to progress in this direction is the elegant metadynamics (MTD) approach. However, in the current work we point out the subtle efficiency problems associated with this approach, and illustrate that we have at hand what is arguably a more powerful approach. More specifically, we demonstrate the effectiveness of an updated version of our original idea of using a classical reference potential for QM(ai)/MM calculations [J. Phys. Chem. B. 102 (1998), 2293)], which we refer to as “paradynamics” (PD). This approach is based on the use of an empirical valence bond (EVB) reference potential, which is already similar to the real ab initio potential. The reference potential is fitted to the ab initio potential by an iterative and, to a great degree, automated refinement procedure. The corresponding free energy profile is then constructed using the refined EVB potential, and the linear response approximation (LRA) is used to evaluate the QM(ai)/MM activation free energy barrier. The automated refinement of the EVB surface (and thus the reduction of the difference between the reference and ab initio potentials) is a key factor in accelerating the convergence of the LRA approach. We apply our PD approach to a test reaction, namely the SN2 reaction between chloride ion and methyl chloride, and demonstrate that, at present, this approach is far more powerful and cost effective than the metadynamics approach (at least in its current implementation). We also discuss the general features
Electronic and mechanical properties of ZnX (X = S, Se and Te)--An ab initio study
Verma, Ajay Singh; Sharma, Sheetal; Jindal, Vijay Kumar; Sarkar, Bimal Kumar
2011-12-12
Zinc chalcogenides (ZnX, X = S, Se and Te) have been increasing attention as wide and direct band gap semiconductor for blue and ultraviolet optical devices. This paper analyzes electronic and mechanical properties of these materials by ab initio pseudo-potential method that uses non conserving pseudopotentials in fully nonlocal form, as implemented in SIESTA code. In this approach the local density approximation (LDA) is used for the exchange-correlation (XC) potential. The calculations are given for band gap, elastic constants (C{sub 11}, C{sub 12} and C{sub 44}), shear modulus, and Young's modulus. The results are in very good agreement with previous theoretical calculations and available experimental data.
An ab initio study of the size-dependent mechanical behavior of single-walled AlN nanotubes
NASA Astrophysics Data System (ADS)
Hao, Jun-Hua; Wang, Yu-Fang; Yin, Yu-Hua; Jiang, Run; Wang, Yun-Feng; Jin, Qing-Hua
2015-07-01
Employing ab initio electronic structure calculations combined with the linear combination of atomic orbitals (LCAO) we have investigated a size dependence of mechanical behavior in single-walled AlN nanotubes with armchair and zigzag forms. A simple procedure of nanotubes construction based on the wurtzite (0 0 1) slab with monolayer rolling and subsequent cylindrical coordinate system introduction is suggested. The present calculations indicate that the Young's modulus and electronic band gap of these tubes are increased monotonically as the radius increases, but decreases with the Al-N bond length. In addition, the amount of charge transfer calculated by the Mulliken's population analysis is introduced to explain clearly the strength of bonding between Al and N atoms in single-walled AlN nanotubes.
Ab-initio calculation study on the formation mechanism of boron-oxygen complexes in c-Si
Yu, Xuegong; Chen, Peng; Chen, Xianzi; Liu, Yong; Yang, Deren
2015-07-15
Boron-oxygen (B-O) complex in crystalline silicon (c-Si) solar cells is responsible for the light-induced efficiency degradation of solar cell. However, the formation mechanism of B-O complex is not clear yet. By Ab-initio calculation, it is found that the stagger-type oxygen dimer (O{sub 2i}{sup st}) should be the component of B-O complex, whose movement occurs through its structure reconfiguration at low temperature, instead of its long-distance diffusion. The O{sub 2i}{sup st} can form two stable “latent centers” with the B{sub s}, which are recombination-inactive. The latent centers can be evolved into the metastable recombination centers via their structure transformation in the presence of excess carriers. These results can well explain the formation behaviors of B-O complexes in c-Si.
Ab initio-predicted micro-mechanical performance of refractory high-entropy alloys
Li, Xiaoqing; Tian, Fuyang; Schönecker, Stephan; Zhao, Jijun; Vitos, Levente
2015-01-01
Recently developed high-entropy alloys (HEAs) consisting of multiple principal elements represent a new field of metallurgy and have demonstrated appealing properties for a wide range of applications. Using ab initio alloy theory, we reveal the alloying effect on the elastic properties and the ideal tensile strength (ITS) in the [001] direction of four body-centered cubic (bcc) refractory HEAs based on Zr, V, Ti, Nb, and Hf. We find that these HEAs show high elastic anisotropy and large positive Cauchy pressure, suggesting good extrinsic ductility. Starting from ZrNbHf, it is found that the ITS decreases with equimolar Ti addition. On the other hand, if both Ti and V are added to ZrNbHf, the ITS is enhanced by about 42%. An even more captivating effect is the ITS increase by about 170%, if Ti and V are substituted for Hf. The alloying effect on the ITS is explained by the d-band filling. An intrinsic brittle-to-ductile transition is found in terms of the failure mode under uniaxial tension. These investigations suggest that intrinsically ductile HEAs with high ideal strength can be achieved by controlling the proportion of group four elements to group five elements. PMID:26199145
Ab initio-predicted micro-mechanical performance of refractory high-entropy alloys.
Li, Xiaoqing; Tian, Fuyang; Schönecker, Stephan; Zhao, Jijun; Vitos, Levente
2015-01-01
Recently developed high-entropy alloys (HEAs) consisting of multiple principal elements represent a new field of metallurgy and have demonstrated appealing properties for a wide range of applications. Using ab initio alloy theory, we reveal the alloying effect on the elastic properties and the ideal tensile strength (ITS) in the [001] direction of four body-centered cubic (bcc) refractory HEAs based on Zr, V, Ti, Nb, and Hf. We find that these HEAs show high elastic anisotropy and large positive Cauchy pressure, suggesting good extrinsic ductility. Starting from ZrNbHf, it is found that the ITS decreases with equimolar Ti addition. On the other hand, if both Ti and V are added to ZrNbHf, the ITS is enhanced by about 42%. An even more captivating effect is the ITS increase by about 170%, if Ti and V are substituted for Hf. The alloying effect on the ITS is explained by the d-band filling. An intrinsic brittle-to-ductile transition is found in terms of the failure mode under uniaxial tension. These investigations suggest that intrinsically ductile HEAs with high ideal strength can be achieved by controlling the proportion of group four elements to group five elements. PMID:26199145
He, Rongxing; Li, Lei; Zhong, Jie; Zhu, Chongqin; Francisco, Joseph S; Zeng, Xiao Cheng
2016-04-26
Solar emission produces copious nitrosonium ions (NO(+)) in the D layer of the ionosphere, 60 to 90 km above the Earth's surface. NO(+) is believed to transfer its charge to water clusters in that region, leading to the formation of gaseous nitrous acid (HONO) and protonated water cluster. The dynamics of this reaction at the ionospheric temperature (200-220 K) and the associated mechanistic details are largely unknown. Using ab initio molecular dynamics (AIMD) simulations and transition-state search, key structures of the water hydrates-tetrahydrate NO(+)(H2O)4 and pentahydrate NO(+)(H2O)5-are identified and shown to be responsible for HONO formation in the ionosphere. The critical tetrahydrate NO(+)(H2O)4 exhibits a chain-like structure through which all of the lowest-energy isomers must go. However, most lowest-energy isomers of pentahydrate NO(+)(H2O)5 can be converted to the HONO-containing product, encountering very low barriers, via a chain-like or a three-armed, star-like structure. Although these structures are not the global minima, at 220 K, most lowest-energy NO(+)(H2O)4 and NO(+)(H2O)5 isomers tend to channel through these highly populated isomers toward HONO formation. PMID:27071120
NASA Astrophysics Data System (ADS)
Ribeiro, Mauro C. C.; Almeida, Luiz C. J.
2000-09-01
Ab initio calculations have been performed in order to investigate a recently proposed polarizable model [M. C. C. Ribeiro, Phys. Rev. B 61, 3297 (2000)] for molecular dynamics (MD) simulation of the molten salt Ca0.4K0.6(NO3)1.4. On the basis of the electronegativity equalization method, polarization effects in the MD simulations have been introduced by a fluctuating charge (FC) model for the nitrate ion. Partial charges in the nitrate ion are obtained by ab initio calculations at several levels of theory, and compared with previously proposed models for MD simulations of nitrate melts. Charge fluctuation is achieved in the ab initio calculations by using positive probe charges placed around a nitrate ion. The parameters of the FC model are corroborated by comparison of the ab initio partial charges with the ones obtained directly by the electronegativity equalization method. Simulated annealing of a cluster including two double-charged cations and two nitrate ions shows that very different structures are obtained depending on whether the FC model or its nonpolarizable counterpart is considered. Ab initio calculations show that the structure of this cluster is strongly dependent on polarization effects in the nitrate ions.
Wong, Kin-Yiu; Yuqing, Xu; York, Darrin M.
2014-01-01
Detailed understandings of the reaction mechanisms of RNA catalysis in various environments can have profound importance for many applications, ranging from the design of new biotechnologies to the unraveling of the evolutionary origin of life. An integral step in the nucleolytic RNA catalysis is self-cleavage of RNA strands by 2′-O-transphosphorylation. Key to elucidating a reaction mechanism is determining the molecular structure and bonding characteristics of transition state. A direct and powerful probe of transition state is measuring isotope effects on biochemical reactions, particularly if we can reproduce isotope effect values from quantum calculations. This paper significantly extends the scope of our previous joint experimental and theoretical work in examining isotope effects on enzymatic and non-enzymatic 2′-O-transphosphorylation reaction models that mimic reactions catalyzed by RNA enzymes (ribozymes), and protein enzymes such as ribonuclease A (RNase A). Native reactions are studied, as well as reactions with thio substitutions representing chemical modifications often used in experiments to probe mechanism. Here, we report and compare results from eight levels of electronic-structure calculations for constructing the potential energy surfaces in kinetic and equilibrium isotope effects (KIE and EIE) computations, including a “gold-standard” coupled-cluster level of theory [CCSD(T)]. In addition to the widely-used Bigeleisen equation for estimating KIE and EIE values, internuclear anharmonicity and quantum tunneling effects were also computed using our recently-developed ab initio path-integral method, i.e., automated integration-free path-integral (AIF-PI) method. The results of this work establish an important set of benchmarks that serve to guide calculations of KIE and EIE for RNA catalysis. PMID:24841935
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.
Time-domain ab initio modeling of photoinduced dynamics at nanoscale interfaces.
Wang, Linjun; Long, Run; Prezhdo, Oleg V
2015-04-01
Nonequilibrium processes involving electronic and vibrational degrees of freedom in nanoscale materials are under active experimental investigation. Corresponding theoretical studies are much scarcer. The review starts with the basics of time-dependent density functional theory, recent developments in nonadiabatic molecular dynamics, and the fusion of the two techniques. Ab initio simulations of this kind allow us to directly mimic a great variety of time-resolved experiments performed with pump-probe laser spectroscopies. The focus is on the ultrafast photoinduced charge and exciton dynamics at interfaces formed by two complementary materials. We consider purely inorganic materials, inorganic-organic hybrids, and all organic interfaces, involving bulk semiconductors, metallic and semiconducting nanoclusters, graphene, carbon nanotubes, fullerenes, polymers, molecular crystals, molecules, and solvent. The detailed atomistic insights available from time-domain ab initio studies provide a unique description and a comprehensive understanding of the competition between electron transfer, thermal relaxation, energy transfer, and charge recombination processes. These advances now make it possible to directly guide the development of organic and hybrid solar cells, as well as photocatalytic, electronic, spintronic, and other devices relying on complex interfacial dynamics. PMID:25622188
NASA Technical Reports Server (NTRS)
Smith, Grant D.; Jaffe, R. L.; Yoon, D. Y.; Arnold, James O. (Technical Monitor)
1994-01-01
Conformational energy contours of perfluoroalkanes, determined from ab initio calculations, confirm the well-known spitting of trans states into two minima at plus or minus 17 degrees but also show that the gauche states split as well, with minima at plus or minus 124 degrees and plus or minus 84 in order to relieve steric crowding. The directions of such split distortions from the perfectly staggered states are strongly coupled for adjacent pairs of bonds in a manner identical to the intradyad pair for poly (isobutylene) chains. These conformational characteristics are fully represented by a six-state rotational isomeric state (RIS) model for PTFE comprised of t(+), t(-), g(sup +)+, g(sup +)-, g(sup -) + and g(sup -)-states, located at the split energy minima. The resultant 6 x 6 statistical weight matrix is described by first-order interaction parameters for the g+(+) (ca. 0.6 kcal/mol) and g+- (ca. 2.0 kcal/mol) states, and second order parameters for the g(sup +)+g(sup +)+ (ca 0.6 kcal/mol) and g(sup +)+g(sup -)+ (ca. 1.0 kcal/mol) states. This six-state RIS model, without adjustment of the geometric or energy parameters as determined from the ab initio calculations, predicts the unperturbed chain dimensions and the fraction of gauche bonds as a function of temperature for PTFE in good agreement with available experimental values.
NASA Astrophysics Data System (ADS)
Xia, W. S.; Lin, M. C.
2001-03-01
The mechanism for the OH+HNO3 reaction has been studied by ab initio molecular orbital calculations at the G2M(cc3) level of theory. Four complexes and four transition states have been found and confirmed by intrinsic reaction coordinate analyses. The commonly assumed six-membered ring complex formed by hydrogen bonding of the OH radical with HNO3, -ON(O)OH…OH-, was found to be stable by 8.1 kcal/mol; its decomposition producing NO3+H2O was predicted to have a barrier of 11.6 kcal/mol. A five-membered ring complex, -ON(O)OH…O(H)-, with the H atom of the OH radical placed out of the ring plane, was found to have a stability of 5.3 kcal/mol; it fragments to form NO3+H2O with a barrier of 6.6 kcal/mol. Two additional complexes, which are the mirror image of each other with a 7.4 kcal/mol binding energy, were found to be related to the OH exchange reaction with a 13.3 kcal/mol barrier above the complexes. The direct abstraction process producing H2O2 and NO2 was predicted to have a large barrier of 24.4 kcal/mol, insignificant to atmospheric chemistry. The rate constant has been calculated at 200-1500 K and 0-760 Torr. The results show that the reaction has strong pressure and tunneling effects below room temperature. In addition, the rate constants for the decay of OH and OD (in OD+DNO3) evaluated by kinetic modeling compare reasonably well with experimental data below room temperature. The unusually pronounced kinetic isotope effect observed experimentally, kH/kD⩾10, could be accounted for by the combination of the greater tunneling rate in the H system and the large redissociation rate of stabilized complexes in the D system. The rate constant predicted for the production of H2O and NO3 in the temperature range 750-1500 K can be effectively represented by the expression k=1.45×10-23 T3.5 exp(+839/T) cm3/s.
NASA Astrophysics Data System (ADS)
Kim, Heung-Sik; Kin-Ho Lee, Eric; Kim, Yong Baek
2015-12-01
The recently discovered three-dimensional hyperhoneycomb iridate, β-Li2IrO3, has raised hopes for the realization of the dominant Kitaev interaction between spin-orbit entangled local moments due to its near-ideal lattice structure. If true, this material may lie close to the sought-after quantum spin-liquid phase in three dimensions. Utilizing ab initio electronic structure calculations, we first show that the spin-orbit entangled basis, j\\text{eff} = 1/2 , correctly captures the low-energy electronic structure. The effective spin model derived in the strong-coupling limit supplemented by the ab initio results is shown to be dominated by the Kitaev interaction. We demonstrated that the possible range of parameters is consistent with a non-coplanar spiral magnetic order found in a recent experiment. All of these analyses suggest that β-Li2IrO3 may be the closest among known materials to the Kitaev spin-liquid regime.
Ab initio approach to model x-ray diffraction in warm dense matter.
Vorberger, J; Gericke, D O
2015-03-01
It is demonstrated how the static electron-electron structure factor in warm dense matter can be obtained from density functional theory in combination with quantum Monte Carlo data. In contrast to theories assuming well-separated bound and free states, this ab initio approach yields also valid results for systems close to the Mott transition (pressure ionization), where bound states are strongly modified and merge with the continuum. The approach is applied to x-ray Thomson scattering and compared to predictions of the Chihara formula whereby we use the ion-ion and electron-ion structure from the same simulations. The results show significant deviations of the screening cloud from the often applied Debye-like form. PMID:25871229
Lister, C.J.; McCutchan, E.A.
2014-06-15
A new generation of ab-initio calculations, based on realistic two- and three-body forces, is having a profound impact on our view of how nuclei work. To improve the numerical methods, and the parameterization of 3-body forces, new precise data are needed. Electromagnetic transitions are very sensitive to the dynamics which drive mixing between configurations. We have made a series of precise (< 3%) measurements of electromagnetic transitions in the A=10 nuclei {sup 10}C and {sup 10}Be by using the Doppler Shift Attenuation method carefully. Many interesting features can be reproduced including the strong α clustering. New measurements on {sup 8}Be and {sup 12}Be highlight the interplay between the alpha clusters and their valence neutrons.
Liu, Lihong; Liu, Jian; Martinez, Todd J.
2015-12-17
Here, we investigate the photoisomerization of a model retinal protonated Schiff base (trans-PSB3) using ab initio multiple spawning (AIMS) based on multi-state second order perturbation theory (MSPT2). Discrepancies between the photodynamical mechanism computed with three-root state-averaged complete active space self-consistent field (SA-3-CASSCF, which does not include dynamic electron correlation effects) and MSPT2 show that dynamic correlation is critical in this photoisomerization reaction. Furthermore, we show that the photodynamics of trans-PSB3 is not well described by predictions based on minimum energy conical intersections (MECIs) or minimum energy conical intersection (CI) seam paths. Instead, most of the CIs involved in the photoisomerizationmore » are far from MECIs and minimum energy CI seam paths. Thus, both dynamical nuclear effects and dynamic electron correlation are critical to understanding the photochemical mechanism.« less
Liu, Lihong; Liu, Jian; Martinez, Todd J.
2015-12-17
Here, we investigate the photoisomerization of a model retinal protonated Schiff base (trans-PSB3) using ab initio multiple spawning (AIMS) based on multi-state second order perturbation theory (MSPT2). Discrepancies between the photodynamical mechanism computed with three-root state-averaged complete active space self-consistent field (SA-3-CASSCF, which does not include dynamic electron correlation effects) and MSPT2 show that dynamic correlation is critical in this photoisomerization reaction. Furthermore, we show that the photodynamics of trans-PSB3 is not well described by predictions based on minimum energy conical intersections (MECIs) or minimum energy conical intersection (CI) seam paths. Instead, most of the CIs involved in the photoisomerization are far from MECIs and minimum energy CI seam paths. Thus, both dynamical nuclear effects and dynamic electron correlation are critical to understanding the photochemical mechanism.
Moin, Syed Tarique; Hofer, Thomas S.; Weiss, Alexander K. H.; Rode, Bernd M.
2013-07-07
Ab initio quantum mechanical charge field molecular dynamics (QMCF-MD) were successfully applied to Cu(II) embedded in water to elucidate structure and to understand dynamics of ligand exchange mechanism. From the simulation studies, it was found that using an extended large quantum mechanical region including two shells of hydration is required for a better description of the dynamics of exchanging water molecules. The structural features characterized by radial distribution function, angular distribution function and other analytical parameters were consistent with experimental data. The major outcome of this study was the dynamics of exchange mechanism and reactions in the first hydration shell that could not be studied so far. The dynamical data such as mean residence time of the first shell water molecules and other relevant data from the simulations are close to the results determined experimentally. Another major characteristic of hydrated Cu(II) is the Jahn-Teller distortion which was also successfully reproduced, leading to the final conclusion that the dominating aqua complex is a 6-coordinated species. The ab initio QMCF-MD formalism proved again its capabilities of unraveling even ambiguous properties of hydrated species that are far difficult to explore by any conventional quantum mechanics/molecular mechanics (QM/MM) approach or experiment.
NASA Astrophysics Data System (ADS)
Vincent, E.; Becquart, C. S.; Domain, C.
2007-02-01
The embrittlement of pressure vessel steels under radiation has been long ago correlated with the presence of Cu solutes. Other solutes such as Ni, Mn and Si are now suspected to contribute also to the embrittlement. The interactions of these solutes with radiation induced point defects thus need to be characterized properly in order to understand the elementary mechanisms behind the formation of the clusters formed upon radiation. Ab initio calculations based on the density functional theory have been performed to determine the interactions of point defects with solute atoms in dilute FeX alloys (X = Cu, Mn, Ni or Si) in order to build a database used to parameterise an atomic kinetic Monte Carlo model. Some results of irradiation damage in dilute Fe-CuNiMnSi alloys obtained with this model are presented.
Wu, Linmin; Zhang, Jing
2015-12-14
The mechanical properties of Li{sub x}CoO{sub 2} under various Li concentrations and associated anisotropy have been systematically studied using the first principles method. During the lithium intercalation process, the Young's modulus, bulk modulus, shear modulus, and ultimate strength increase with increasing lithium concentration. Strong anisotropy of mechanical properties between a-axis and c-axis in Li{sub x}CoO{sub 2} is identified at low lithium concentrations, and the anisotropy decreases with increasing lithium concentration. The observed lithium concentration dependence and anisotropy are explained by analyzing the charge transfer using Bader charge analysis, bond order analysis, and bond strength by investigating partial density of states and charge density difference. With the decrease of Li concentration, the charge depletion in the bonding regions increases, indicating a weaker Co-O bond strength. Additionally, the Young's modulus, bulk modulus, shear modulus, and toughness are obtained by simulating ab initio tensile tests. From the simulated stress-strain curves, Li{sub x}CoO{sub 2} shows the highest toughness, which is in contraction with Pugh criterion prediction based on elastic properties only.
NASA Astrophysics Data System (ADS)
Wu, Linmin; Zhang, Jing
2015-12-01
The mechanical properties of LixCoO2 under various Li concentrations and associated anisotropy have been systematically studied using the first principles method. During the lithium intercalation process, the Young's modulus, bulk modulus, shear modulus, and ultimate strength increase with increasing lithium concentration. Strong anisotropy of mechanical properties between a-axis and c-axis in LixCoO2 is identified at low lithium concentrations, and the anisotropy decreases with increasing lithium concentration. The observed lithium concentration dependence and anisotropy are explained by analyzing the charge transfer using Bader charge analysis, bond order analysis, and bond strength by investigating partial density of states and charge density difference. With the decrease of Li concentration, the charge depletion in the bonding regions increases, indicating a weaker Co-O bond strength. Additionally, the Young's modulus, bulk modulus, shear modulus, and toughness are obtained by simulating ab initio tensile tests. From the simulated stress-strain curves, LixCoO2 shows the highest toughness, which is in contraction with Pugh criterion prediction based on elastic properties only.
Forssen, C; Caurier, E; Navratil, P
2008-12-23
Recently, charge radii and ground-state electromagnetic moments of Li and Be isotopes were measured precisely. We have performed large-scale ab initio no-core shell model calculations for these isotopes using high-precision nucleon-nucleon potentials. Our computed charge radii, quadrupole and magnetic-dipole moments are in a good agreement with the measurements with the exception of the {sup 11}Li charge radius. The overall trends of all observables are well reproduced. The magnetic moments are in particular well described. Also, we are able to reproduce the small magnitude of the {sup 6}Li quadrupole moment and with the CD-Bonn NN potential also its correct sign.
Forssen, C.; Caurier, E.; Navratil, P.
2009-02-15
Recently, charge radii and ground-state electromagnetic moments of Li and Be isotopes were measured precisely. We have performed large-scale ab initio no-core shell model calculations for these isotopes using high-precision nucleon-nucleon potentials. The isotopic trends of our computed charge radii and quadrupole and magnetic-dipole moments are in good agreement with experimental results with the exception of the {sup 11}Li charge radius. The magnetic moments are in particular well described, whereas the absolute magnitudes of the quadrupole moments are about 10% too small. The small magnitude of the {sup 6}Li quadrupole moment is reproduced, and with the CD-Bonn NN potential, also its correct sign.
Pigozzi, Giancarlo; Janczak-Rusch, Jolanta; Passerone, Daniele; Antonio Pignedoli, Carlo; Patscheider, Joerg; Jeurgens, Lars P. H.; Antusek, Andrej; Parlinska-Wojtan, Magdalena; Bissig, Vinzenz
2012-10-29
Nano-sized Ag-Cu{sub 8nm}/AlN{sub 10nm} multilayers were deposited by reactive DC sputtering on {alpha}-Al{sub 2}O{sub 3}(0001) substrates. Investigation of the phase constitution and interface structure of the multilayers evidences a phase separation of the alloy sublayers into nanosized grains of Ag and Cu. The interfaces between the Ag grains and the quasi-single-crystalline AlN sublayers are semi-coherent, whereas the corresponding Cu/AlN interfaces are incoherent. The orientation relationship between Ag and AlN is constant throughout the entire multilayer stack. These observations are consistent with atomistic models of the interfaces as obtained by ab initio calculations.
Surface electron density models for accurate ab initio molecular dynamics with electronic friction
NASA Astrophysics Data System (ADS)
Novko, D.; Blanco-Rey, M.; Alducin, M.; Juaristi, J. I.
2016-06-01
Ab initio molecular dynamics with electronic friction (AIMDEF) is a valuable methodology to study the interaction of atomic particles with metal surfaces. This method, in which the effect of low-energy electron-hole (e-h) pair excitations is treated within the local density friction approximation (LDFA) [Juaristi et al., Phys. Rev. Lett. 100, 116102 (2008), 10.1103/PhysRevLett.100.116102], can provide an accurate description of both e-h pair and phonon excitations. In practice, its applicability becomes a complicated task in those situations of substantial surface atoms displacements because the LDFA requires the knowledge at each integration step of the bare surface electron density. In this work, we propose three different methods of calculating on-the-fly the electron density of the distorted surface and we discuss their suitability under typical surface distortions. The investigated methods are used in AIMDEF simulations for three illustrative adsorption cases, namely, dissociated H2 on Pd(100), N on Ag(111), and N2 on Fe(110). Our AIMDEF calculations performed with the three approaches highlight the importance of going beyond the frozen surface density to accurately describe the energy released into e-h pair excitations in case of large surface atom displacements.
NASA Astrophysics Data System (ADS)
Chu, Xingli; Zhang, Yanxing; Li, Shasha; Yang, Zongxian
2014-04-01
The first-principles method based on density functional theory (DFT) is used to investigate the reaction mechanism for the adsorption of H2S on the oxygen-enriched yttria-stabilized zirconia (YSZ + O) (111) surface. It is found that the H2S dissociation processes have low energy barriers (< 0.5 eV) and high exothermicities (2.5 eV), and the dissociative S atoms may result in the poisoning of the YSZ + O surface by forming the SO and the hyposulfite (SO22 -) species with very strong bonds to the surface. In addition, using the ab initio atomistic thermodynamics method, the surface regeneration or de-sulfurization process of a sulfur-poisoned (i.e. sulfur-covered) YSZ + O(111) surface is studied. According to the phase diagram, the adsorbed atomic sulfur can be oxidized to SO2 and removed from the YSZ + O surface by introducing oxidizing reagents, e.g. O2 and H2O.
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
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.
Mechanical and dynamical stability of TiAsTe compound from ab initio calculations
NASA Astrophysics Data System (ADS)
Deligoz, E.; Ozisik, H.
2015-07-01
The first-principles calculations are employed to provide a fundamental understanding of the structural features and relative thermodynamical, mechanical and phonon stability of TiAsTe compound. The calculated lattice parameters are in good agreement with available experimental results. We have computed elastic constants, its derived moduli and ratios that characterize mechanical properties for the first time. The calculated elastic constants indicate that these materials are mechanically stable at ambient condition. The minimum thermal conductivities of TiAsTe are calculated using both Clarke's model and Cahill's model. Furthermore, the elastic anisotropy has been visualized in detail by plotting the directional dependence of compressibility, Young's modulus and shear modulus. Our results suggest strong elastic anisotropy for this compound. Additionally, the phonon spectra and phonon density of states are also obtained and discussed. The full phonon dispersion calculations confirm the dynamic stability of TiAsTe.
Diffusion Mechanism of Hydrogen in Amorphous Silicon:
Su, Y.-S.; Pantelides, Sokrates T
2002-01-01
The mechanism of H migration in amorphous Si has remained an unresolved problem. The main issue is the small activation energy (1.5 eV) relative to the known strength of Si-H bonds (2--3.5 eV). We report first-principles finite-temperature simulations which demonstrate vividly that H is not released spontaneously, as proposed by most models, but awaits the arrival of a floating bond (FB). The 'migrating species' is an FB-H complex, with H jumping from Si to Si and the FB literally floating around it. Migration stops when the FB veers away.
NASA Astrophysics Data System (ADS)
Churakov, S. V.; Parrinello, M.
2003-04-01
The carbonation of brucite (Mg(OH)2) has been considered as a potential technology for cleaning industrial carbon dioxide waste. The kinetics of the reaction Mg(OH)2 + CO2 -> MgCO3 + H2O have been studied experimentally at 573°C by Bearat at al. [1]. Their experiments suggest that the carbonation of magnesium hydroxide proceeds by the reaction Mg(OH)2 -> MgO + H2O followed by the adsorption of CO2 molecules on the dehydrated brucite surface. Due to the large difference in volumes between Mg(OH)2 and MgO, dehydration causes the formation of dislocations and cracks, allowing water molecules to leave the brucite surface and facilitating the advance of the carbonation front in the bulk solid. The detailed mechanism of this process is however unknown. We used the Car-Parrinello ab initio molecular dynamics method to study the structure and dynamics of the (0001), (1-100) and (11-20) surfaces of brucit and calculated the enthalpy and activation barrier of H2O nucleation and dehydration on different surfaces. The results obtained are in agreement with previous studies of brucite dehydration by Masini and Bernasconi [2]. The reactive Car-Parrinello molecular dynamics method [3] has been applied to investigate the detailed mechanism of the dehydration-carbonation reaction at the (1-100) interface of brucite with the gas phase. Based on the results of our MD simulations and the calculated enthalpy of CO2 adsorption on the dehydrated brucite surfaces we propose a mechanism for the dehydration/carbonation reaction. [1] Bearat H, McKelvy MJ, Chizmeshya AVG, Sharma R, Carpenter RW (2002) J. Amer. Ceram. Soc. 85(4):742 [2] Masini P and Bernasconi M (2001) J. Phys. Cond. Mat. 13: 1-12 [3] Iannuzzi M, Laio A and Parrinello M (2003) Phys. Rev. Lett. (submitted)
Zhou, Yanzi; Wang, Shenglong; Zhang, Yingkai
2010-01-01
Acetylcholinesterase (AChE) is a remarkably efficient serine hydrolase responsible for the termination of impulse signaling at cholinergic synapses. By employing Born-Oppenheimer molecular dynamics simulations with B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized its complete catalytic reaction mechanism for hydrolyzing neurotransmitter acetylcholine (ACh) and determined its multi-step free energy reaction profiles for the first time. In both acylation and deacylation reaction stages, the first step involves the nucleophilic attack to the carbonyl carbon with the triad His447 serving as the general base, and leads to a tetrahedral covalent intermediate stabilized by the oxyanion hole. From the intermediate to the product, the orientation of His447 ring needs to be adjusted very slightly, and then the proton transfers from His447 to the product and the break of the scissile bond happen spontaneously. For the three-pronged oxyanion hole, it only makes two hydrogen bonds with the carbonyl oxygen at either the initial reactant or the final product state, but the third hydrogen bond is formed and stable at all transition and intermediate states during the catalytic process. At the intermediate state of the acylation reaction, a short and low-barrier hydrogen bond (LBHB) is found to be formed between two catalytic triad residues His447 and Glu334, and the spontaneous proton transfer between two residues has been observed. However, it is only about 1 ~ 2 kcal/mol stronger than the normal hydrogen bond. In comparison with previous theoretical investigations of the AChE catalytic mechanism, our current study clearly demonstrates the power and advantages of employing Born-Oppenheimer ab initio QM/MM MD simulations in characterizing enzyme reaction mechanisms. PMID:20550161
Ab initio bubble-driven denaturation of double-stranded DNA: Self-mechanical theory.
Kuetche, Victor K
2016-07-21
Among the different theoretical models of the open-site-driven DNA-denaturation found in the literature, very few interests are actually paid to the fundamental unzipping process of the double-stranded DNA within the vicinity of its ground state condensate. In this paper, we address an alternative to better understand the process of denaturation of such a macromolecule by investigating the onset of its dynamics around its equilibrium state. We show that from the initiation of the transcription bubble by the promoter to the termination state, the open-states of the strands evolve dynamically while generating some localized waveguide channels with elastic scattering properties. We properly discuss the nonlinear dynamics of these structures within the viewpoint of the self-mechanical theory while inferring to the physical structure of the findings and their potential issues. PMID:27113786
Thermodynamic and mechanical properties of TiC from ab initio calculation
Dang, D. Y.; Fan, J. L.; Gong, H. R.
2014-07-21
The temperature-dependent thermodynamic and mechanical properties of TiC are systematically investigated by means of a combination of density-functional theory, quasi-harmonic approximation, and thermal electronic excitation. It is found that the quasi-harmonic Debye model should be pertinent to reflect thermodynamic properties of TiC, and the elastic properties of TiC decease almost linearly with the increase of temperature. Calculations also reveal that TiC possesses a pronounced directional pseudogap across the Fermi level, mainly due to the strong hybridization of Ti 3d and C 2p states. Moreover, the strong covalent bonding of TiC would be enhanced (reduced) with the decrease (increase) of temperature, while the change of volume (temperature) should have negligible effect on density of states at the Fermi level. The calculated results agree well with experimental observations in the literature.
NASA Astrophysics Data System (ADS)
Xu, Li-Hong; Lees, Ronald M.; Hougen, Jon T.
2013-06-01
The ν_{2}, ν_{3} and ν_{9} CH stretching modes of methanol in the 3μm region exhibit a significant amount of torsion-vibration interaction, as illustrated for ν_{9} by the facts that: (i) the three hydrogen atoms each pass through a plane of symmetry of the molecule twice during the course of one full internal rotation motion, once at a minimum and once at a maximum in the three-fold potential energy curve, (ii) the H atom in the plane of symmetry is nearly motionless for the ν_{9} mode, and therefore (iii) the property of remaining motionless must be transferred from one H to another six times during one full internal rotation motion. In this talk we examine quantitatively the general phenomenon of torsion-vibration interaction in the methyl top stretching modes in two ways. First, we present plots of normal modes produced in Gaussian projected frequency calculations that are expressed either in terms of several sets of internal coordinates, or in terms of Cartesian displacement vectors for the methyl hydrogen atoms. Some of these plots display a nearly three-fold sine or cosine behavior, where the sine or cosine behavior is dictated by group-theoretical symmetry arguments. Other plots display stunning features ranging from loss of simple three-fold oscillatory pattern to cusp-like peaks or dips. Somewhat surprisingly, none of our ab initio plots for methanol exhibit a sign change after a 2π internal rotation of the methyl top. Second, we present a relatively simple model for the three CH stretching motions, characterized by three parameters associated with: (i) a vibrational A/E energy difference, (ii) a Jahn-Teller-like torsion-vibration interaction term within the vibrational E state, and (iii) a Renner-Teller-like torsion-vibration interaction term within the E state. This model gives nearly quantitative agreement with both the regular and irregular features of the ab initio plots. The good agreement suggests that various aspects of the physics of the
NASA Astrophysics Data System (ADS)
Jackson, Bret; Nattino, Francesco; Kroes, Geert-Jan
2014-08-01
The dissociative chemisorption of methane on metal surfaces is of great practical and fundamental importance. Not only is it the rate-limiting step in the steam reforming of natural gas, the reaction exhibits interesting mode-selective behavior and a strong dependence on the temperature of the metal. We present a quantum model for this reaction on Ni(100) and Ni(111) surfaces based on the reaction path Hamiltonian. The dissociative sticking probabilities computed using this model agree well with available experimental data with regard to variation with incident energy, substrate temperature, and the vibrational state of the incident molecule. We significantly expand the vibrational basis set relative to earlier studies, which allows reaction probabilities to be calculated for doubly excited initial vibrational states, though it does not lead to appreciable changes in the reaction probabilities for singly excited initial states. Sudden models used to treat the center of mass motion parallel to the surface are compared with results from ab initio molecular dynamics and found to be reasonable. Similar comparisons for molecular rotation suggest that our rotationally adiabatic model is incorrect, and that sudden behavior is closer to reality. Such a model is proposed and tested. A model for predicting mode-selective behavior is tested, with mixed results, though we find it is consistent with experimental studies of normal vs. total (kinetic) energy scaling. Models for energy transfer into lattice vibrations are also examined.
Jackson, Bret; Nattino, Francesco; Kroes, Geert-Jan
2014-08-07
The dissociative chemisorption of methane on metal surfaces is of great practical and fundamental importance. Not only is it the rate-limiting step in the steam reforming of natural gas, the reaction exhibits interesting mode-selective behavior and a strong dependence on the temperature of the metal. We present a quantum model for this reaction on Ni(100) and Ni(111) surfaces based on the reaction path Hamiltonian. The dissociative sticking probabilities computed using this model agree well with available experimental data with regard to variation with incident energy, substrate temperature, and the vibrational state of the incident molecule. We significantly expand the vibrational basis set relative to earlier studies, which allows reaction probabilities to be calculated for doubly excited initial vibrational states, though it does not lead to appreciable changes in the reaction probabilities for singly excited initial states. Sudden models used to treat the center of mass motion parallel to the surface are compared with results from ab initio molecular dynamics and found to be reasonable. Similar comparisons for molecular rotation suggest that our rotationally adiabatic model is incorrect, and that sudden behavior is closer to reality. Such a model is proposed and tested. A model for predicting mode-selective behavior is tested, with mixed results, though we find it is consistent with experimental studies of normal vs. total (kinetic) energy scaling. Models for energy transfer into lattice vibrations are also examined.
Jackson, Bret; Nattino, Francesco; Kroes, Geert-Jan
2014-08-01
The dissociative chemisorption of methane on metal surfaces is of great practical and fundamental importance. Not only is it the rate-limiting step in the steam reforming of natural gas, the reaction exhibits interesting mode-selective behavior and a strong dependence on the temperature of the metal. We present a quantum model for this reaction on Ni(100) and Ni(111) surfaces based on the reaction path Hamiltonian. The dissociative sticking probabilities computed using this model agree well with available experimental data with regard to variation with incident energy, substrate temperature, and the vibrational state of the incident molecule. We significantly expand the vibrational basis set relative to earlier studies, which allows reaction probabilities to be calculated for doubly excited initial vibrational states, though it does not lead to appreciable changes in the reaction probabilities for singly excited initial states. Sudden models used to treat the center of mass motion parallel to the surface are compared with results from ab initio molecular dynamics and found to be reasonable. Similar comparisons for molecular rotation suggest that our rotationally adiabatic model is incorrect, and that sudden behavior is closer to reality. Such a model is proposed and tested. A model for predicting mode-selective behavior is tested, with mixed results, though we find it is consistent with experimental studies of normal vs. total (kinetic) energy scaling. Models for energy transfer into lattice vibrations are also examined. PMID:25106565
Molt, Robert W; Bartlett, Rodney J; Watson, Thomas; Bazanté, Alexandre P
2012-12-13
We have identified the major conformers of CL-20 explosive, otherwise known as 2,4,6,8,10,12-hexanitrohexaazaisowurtzitane, more formally known as 2,4,6,8,10,12-hexanitrohexaazatetracyclo[5.5.0.0]-dodecane, via Monte Carlo search in conformational space through molecular mechanics and subsequent quantum mechanical refinement using perturbation theory. Our search produced enough conformers to account for all of the various forms of CL-20 found in crystals. This suggests that our methodology will be useful in studying the conformational landscape of other nitramines. The energy levels of the conformers found are all within 0.25 eV of one another based on MBPT(2)/6-311G(d,p); consequently, without further refinement from a method such as coupled cluster theory, all conformers may reasonably be populated at STP in the gas phase. We also report the harmonic vibrational frequencies of conformers, including the implications on the mechanism of detonation. In particular, we establish that the weakest N-N nitramine of CL-20 is the cyclohexane equatorial nitramine. This preliminary mapping of the conformers of CL-20 makes it possible to study the mechanism of detonation of this explosive rigorously in future work. PMID:23136867
Ab initio calculations on the compensation mechanisms in InP
NASA Astrophysics Data System (ADS)
Schmidt, T. M.; Miwa, R. H.; Fazzio, A.; Mota, R.
2001-01-01
Compensation mechanisms which become semi-insulating InP from intrinsic n-type materials are investigated by first principles total energy calculations. Associated to the possible mechanisms two intrinsic defects are shown to play an important role in the transition, induced by annealing, from the n-type to the semi-insulating character. The single defect V In3- results from the In-vacancy-H-complex, which, before the annealing, was identified by infra-red spectra as responsible for the n-type conductivity, and the P In+V P complex defect comes from the diffusion of the neutral In vacancy. Both defects are shown to present localized levels inside the band gap which compensate for the free electron.
Enhancing mechanical properties of calcite by Mg substitutions: An ab initio study
NASA Astrophysics Data System (ADS)
Elstnerova, Pavlina; Friak, Martin; Hickel, Tilmann; Fabritius, Helge Otto; Lymperakis, Liverios; Petrov, Michal; Raabe, Dierk; Neugebauer, Joerg; Nikolov, Svetoslav; Zigler, Andreas; Hild, Sabine
2011-03-01
Arthropoda representing a majority of all known animal species are protected by an exoskeleton formed by their cuticle. The cuticle represents a hierarchically structured multifunctional bio-composite based on chitin and proteins. Some groups like Crustacea reinforce the load-bearing parts of their cuticle with calcite. As the calcite sometimes contains Mg it was speculated that Mg may have a stiffening impact on the mechanical properties of the cuticle. We present a theoretical parameter-free quantum-mechanical study of thermodynamic, structural and elastic properties of Mg-substituted calcite. Our results show that substituting Ca by Mg causes an almost linear decrease in the crystal volume with Mg concentration and of substituted crystals. As a consequence the calcite crystals become stiffer giving rise e.g. to substantially increased bulk moduli.
NASA Astrophysics Data System (ADS)
Wu, Hai-Ying; Chen, Ya-Hong; Zhou, Ping; Han, Xiang-Yu; Liu, Zi-Jiang
2014-09-01
The structural, electronic, and mechanical stability properties of magnesium sulfide in different phases are presented using the plane wave pseudopotential method within the generalized gradient approximation. Eight different phases such as rocksalt (B1), zincblende (B3), wurtzite (B4), nickel arsenide (B8), cesium chloride (B2), PH4I-type (B11), FeSi-type (B28), and MnP-type (B31) are considered in great detail. The calculated ground-state properties of these phases are consistent with available experimental and theoretical data. It is found that MgS in the B1 and B8 phases are indirect band gap materials, the B3, B4, B11, B28, and B31 phases are all direct gap materials, while the B2 phase displays the metallic character. The B1, B3, B4, B8, B28, and B31 phases are mechanically stable at ambient conditions, but the B2 and B11 phases are mechanically unstable under zero pressure and zero temperature
NASA Astrophysics Data System (ADS)
Yakunin, Konstantin N.; Mezzacappa, Anthony; Marronetti, Pedro; Yoshida, Shin'ichirou; Bruenn, Stephen W.; Hix, W. Raphael; Lentz, Eric J.; Bronson Messer, O. E.; Harris, J. Austin; Endeve, Eirik; Blondin, John M.; Lingerfelt, Eric J.
2015-10-01
We present the gravitational waveforms computed in ab initio two-dimensional core collapse supernova models evolved with the chimera code for progenitor masses between 12 and 25 M⊙ . All models employ multifrequency neutrino transport in the ray-by-ray approximation, state-of-the-art weak interaction physics, relativistic transport corrections such as the gravitational redshift of neutrinos, two-dimensional hydrodynamics with the commensurate relativistic corrections, Newtonian self-gravity with a general-relativistic monopole correction, and the Lattimer-Swesty equation of state with 220 MeV compressibility, and begin with the most recent Woosley-Heger nonrotating progenitors in this mass range. All of our models exhibit robust explosions. Therefore, our waveforms capture all stages of supernova development: 1) a relatively short and weak prompt signal, 2) a quiescent stage, 3) a strong signal due to convection and standing accretion shock instability activity, 4) termination of active accretion onto the proto-neutron star, and 5) a slowly increasing tail that reaches a saturation value. Fourier decomposition shows that the gravitational wave signals we predict should be observable by AdvLIGO for Galactic events across the range of progenitors considered here. The fundamental limitation of these models is in their imposition of axisymmetry. Further progress will require counterpart three-dimensional models, which are underway.
Milowska, Karolina Z.; Birowska, Magdalena; Majewski, Jacek A.
2013-12-04
We present exemplary results of extensive studies of structural, mechanical and electronic properties of covalent functionalization of carbon nanotubes (CNTs). We report new results for metallic (9,0), and semiconducting (10,0) single-wall carbon nanotubes (CNT) functionalized with -COOH, -OH, and both groups with concentration up to 12.5%. Our studies are performed in the framework of the density functional theory (DFT). We discuss here the stability, local and global changes in structure, elastic moduli (Young's, Shear, and Bulk), electronic structure and resulting band gaps, as a function of the density of the adsorbed molecules.
Mechanisms of branching reactions in melanin formation - Ab initio quantum engineering approach -
NASA Astrophysics Data System (ADS)
Kishida, Ryo; Menez Aspera, Susan; Kasai, Hideaki
Melanin, a pigment found in animals, consists of two types of oligomeric unit: eumelanin and pheomelanin. The color of the skin, the hair, and the eyes is controlled by the ratio of eumelanin/pheomelanin production. Especially, dopachrome and dopaquinone are the precursor molecules of melanin which directly affect the composition of melanin through their branching reactions. Dopachrome is converted into two possible monomers of eumelanin. Dopaquinone can undergo both eumelanin and pheomelanin synthesis. To understand the mechanisms and controlling factors that govern the conversions, reactions of the two molecules are investigated using density functional theory-based first-principles calculations. Our results deepen mechanistic understanding of the reactions and open possibilities to design properties and functions of melanin. In this talk, we will discuss about the competitions of the branching reactions.
Ab initio study of mechanism of forming a Spiro-Si-heterocyclic ring compound
NASA Astrophysics Data System (ADS)
Lu, Xiuhui; Ji, Hua
2013-04-01
The mechanism of the cycloaddition reaction of forming a spiro-Si-heterocyclic ring compound between singlet dichloroalkylidenesilylene (Cl2C=Si :) and ethene has been investigated with CCSD(T)//MP2/6-31G* method. From the potential energy profile, it can be predicted that the reaction has one dominant reaction pathway. The presented rule of this reaction is that the 3 p unoccupied orbital of Si in dichloroalkylidene and the π orbital of ethene forming the π → p donor-acceptor bond, resulting in the formation of a three-membered ring intermediate. Ring-enlargement effect make the three-membered ring intermediate isomerizes to a four-membered ring silylene. Due to sp 3 hybridization of Si atom in the four-membered ring silylene, the four-membered ring silylene further combines with ethene to form a spiro-Si-heterocyclic ring compound.
Guo, Jing-hua; Zhang, Hong; Tang, Yongjian; Cheng, Xinlu
2013-02-28
The hydrogen spillover mechanism, including the H chemisorption, diffusion, and H(2) associative desorption on the surface of COFs and H atoms migration from metal catalyst to COFs, have been studied via density functional theory (DFT) calculation. The results described herein show that each sp(2) C atom on COFs' surface can adsorb one H atom with the bond length d(C-H) between 1.11 and 1.14 Å, and the up-down arrangement of the adsorbed H atoms is the most stable configuration. By counting the chemisorption binding sites for these COFs, we can predict the saturation storage densities. High hydrogen storage densities show that the gravimetric uptakes of COFs are in the range of 5.13-6.06 wt%. The CI-NEB calculations reveal that one H atom diffusing along the C-C path on HHTP surface should overcome the 1.41-2.16 eV energy barrier. We chose tetrahedral Pt(4) cluster and HHTP as the representative catalyst and substrate, respectively, to study the H migration from metal cluster to COFs. At most, two H atoms can migrate from Pt(4) cluster to HHTP substrate. The migration reaction is an endothermic process, undergoing an activation barrier of 1.87 eV and 0.57 eV for the first and second H migration process, respectively. Three types of H(2) associative desorption from hydrogenated COFs were studied: (I) the two H adatoms recombining to one H(2) molecule with a recombination barrier of 4.28 eV, (II) the abstraction of adsorbed H atoms by gas-phase hydrogen atoms through ER type recombination reactions with a recombination barrier of 1.05 eV, (III) the H(2) desorption through the reverse spillover mechanism with an energy barrier of 2.90 eV. PMID:23338125
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.
Ab initio calculation of mechanical and thermal properties of U2Mo intermetallic
NASA Astrophysics Data System (ADS)
Jaroszewicz, S.; Losada, E. L.; Garcés, J. E.; Mosca, H. O.
2013-10-01
We present a study of structural, elastic and thermodynamic properties of tetragonal (C11b) U2Mo by means of density-functional theory based calculations using full-potential linearized augmented plane wave method. In this approach the generalized gradient approximation were used for the exchange-correlation potential calculation. The optimized lattice parameters are in excellent agreement with the experimental data. Through the Debye-Grüneisen model the temperature and pressure dependence of equation of state, bulk modulus, thermal expansion and specific heat have been obtained and discussed in the range of pressure 0-20 GPa and the temperature 0-800 K.
Uranus and Neptune structure models with ab initio EOS data for CH4, NH3, and H2O
NASA Astrophysics Data System (ADS)
Nettelmann, Nadine; Fortney, Jonathan; Hamel, Sebastien; Bethkenhagen, Mandy; Redmer, Ronald
2014-05-01
Uranus and Neptune are supposed to be rich in ices in their deep interiors as their mean density closely resembles that of liquid water. Moreover, highly super-solar abundances of CH4 and CO, indicative of internal water, have been observed in their atmospheres. We here compare ab initio equations of state for CH4, NH3, and H2O and apply them to compute ice-rich, adiabatic internal structure models of Uranus and Neptune. The explicit consideration of the light ices CH4 and NH3 allows us to put tighter constraints on the minimum H/He abundance in their deep interior, which was found to be non-zero in all previous Uranus and in most of the Neptune models that were based on water as a proxy for ices. In particular, we investigate if hydrogen in the deep interior can solely be a result of assumed Carbon sedimentation (diamond rain), as an alternative scenario to the early accretion of H/He containing material during the formation of the planets. We conclude by discussing the deep internal H/He abundance in light of rock-rich and warmer-than-adiabatic interiors, which has been suggested to explain Uranus' low intrinsic luminosity. Our models serve to better understand the formation and bulk composition of Neptune-sized planets.
Miller, J.; Miaskiewicz, K.; Osman, R.
1993-12-01
Studies of ring-saturated pyrimidine base lesions are used to illustrate an integrated modeling approach that combines quantum-chemical calculations with molecular dynamics simulation. Electronic-structure calculations on the lesions in Isolation reveal strong conformational preferences due to interactions between equatorial substituents to the pyrimidine ring. Large distortions of DNA should result when these interactions force the methyl group of thymine to assume an axial orientation, as is the case for thymine glycol but not for dihydrothymine. Molecular dynamics simulations of the dodecamer d(CGCGAATTCGCG){sub 2} with and without a ring-saturated thymine lesion at position T7 support this conclusion. Implications of these studies for recognition of thymine lesions by endonuclease III are also discussed.
7Be(p,(gamma))8B S-factor from Ab Initio No-Core Shell Model Wave Functions
Navratil, P; Bertulani, C A; Caurier, E
2005-12-02
Nuclear structure of {sup 7}Be, {sup 8}B and {sup 7,8}Li is studied within the ab initio no-core shell model (NCSM). Starting from high-precision nucleon-nucleon (NN) interactions, wave functions of {sup 7}Be and {sup 8}B bound states are obtained in basis spaces up to 10 h bar{Omega} and used to calculate channel cluster form factors (overlap integrals) of the {sup 8}B ground state with {sup 7}Be+p. Due to the use of the harmonic oscillator (HO) basis, the overlap integrals have incorrect asymptotic properties. We fix this problem in two alternative ways. First, by a Woods-Saxon (WS) potential solution fit to the interior of the NCSM overlap integrals. Second, by a direct matching with the Whittaker function. The corrected overlap integrals are then used for the {sup 7}Be(p,{gamma}){sup 8}B S-factor calculation. We study the convergence of the S-factor with respect to the NCSM HO frequency and the model space size. Our S-factor results are in agreement with recent direct measurement data. We also test the spectroscopic factors and the corrected overlap integrals from the NCSM in describing the momentum distributions in knockout reactions with {sup 8}B projectiles. A good agreement with the available experimental data is also found, attesting the overall consistency of the calculations.
Ab initio-aided CALPHAD thermodynamic modeling of the Sn-Pb binary system under current stressing.
Lin, Shih-kang; Yeh, Chao-kuei; Xie, Wei; Liu, Yu-chen; Yoshimura, Masahiro
2013-01-01
Soldering is an ancient process, having been developed 5000 years ago. It remains a crucial process with many modern applications. In electronic devices, electric currents pass through solder joints. A new physical phenomenon--the supersaturation of solders under high electric currents--has recently been observed. It involves (1) un-expected supersaturation of the solder matrix phase, and (2) the formation of unusual "ring-shaped" grains. However, the origin of these phenomena is not yet understood. Here we provide a plausible explanation of these phenomena based on the changes in the phase stability of Pb-Sn solders. Ab initio-aided CALPHAD modeling is utilized to translate the electric current-induced effect into the excess Gibbs free energies of the phases. Hence, the phase equilibrium can be shifted by current stressing. The Pb-Sn phase diagrams with and without current stressing clearly demonstrate the change in the phase stabilities of Pb-Sn solders under current stressing. PMID:24060995
Ab initio-aided CALPHAD thermodynamic modeling of the Sn-Pb binary system under current stressing
Lin, Shih-kang; Yeh, Chao-kuei; Xie, Wei; Liu, Yu-chen; Yoshimura, Masahiro
2013-01-01
Soldering is an ancient process, having been developed 5000 years ago. It remains a crucial process with many modern applications. In electronic devices, electric currents pass through solder joints. A new physical phenomenon – the supersaturation of solders under high electric currents – has recently been observed. It involves (1) un-expected supersaturation of the solder matrix phase, and (2) the formation of unusual “ring-shaped” grains. However, the origin of these phenomena is not yet understood. Here we provide a plausible explanation of these phenomena based on the changes in the phase stability of Pb-Sn solders. Ab initio-aided CALPHAD modeling is utilized to translate the electric current-induced effect into the excess Gibbs free energies of the phases. Hence, the phase equilibrium can be shifted by current stressing. The Pb-Sn phase diagrams with and without current stressing clearly demonstrate the change in the phase stabilities of Pb-Sn solders under current stressing. PMID:24060995
NASA Astrophysics Data System (ADS)
Motegi, Kyosuke; Nakajima, Takahito; Hirao, Kimihiko; Seijo, Luis
2001-04-01
A relativistic ab initio model potential (AIMP) for Pt, Au, and Hg atoms has been developed using a relativistic scheme by eliminating small components (RESC) in which the 5p, 5d, and 6s electrons are treated explicitly. The quality of new RESC-AIMP has been tested by calculating the spectroscopic properties of the hydrides of these elements using the Hartree-Fock and coupled cluster with singles and doubles (CCSD) methods. The agreement with reference all-electron RESC calculations is excellent. The RESC-AIMP method is applied successfully in the investigation of the spectroscopic constants of Au2 and Hg2 using the CCSD method with a perturbative estimate of the contributions of triples. The ground state of Pt2 is also determined by RESC-AIMP with the second-order complete active space perturbation method. The results show that scalar relativistic effects on the valence properties are well described by the RESC-AIMP method. The effect on the basis set superposition error on the spectroscopic constants is also examined.
NASA Astrophysics Data System (ADS)
Buchachenko, A. A.; Stolyarov, A. V.; Szczȩśniak, M. M.; Chałasiński, G.
2012-09-01
The coefficients at the lowest-order electrostatic, induction, and dispersion terms of the anisotropic long-range potential between the two KRb(1Σ+) molecules are evaluated through the static and dynamic molecular properties using the ab initio coupled cluster techniques. Adiabatic channel potentials for the ground-state molecules are obtained and used for the numerical quantum capture probability calculations in the spirit of the statistical adiabatic channel models. Capture rate coefficients for indistinguishable (polarized) and distinguishable (unpolarized) molecules at temperatures below 10 μK agree well with those computed with the simple isotropic dispersion R-6 potential, but underestimate the measured ones [Ospelkaus et al., Science 327, 853 (2010), 10.1126/science.1184121] up to a factor of 3. Preliminary assessment of the effects of higher-order long-range terms, retardation of dispersion forces, and magnetic dipole-dipole interaction does not offer any clear perspectives for drastic improvement of the capture approximation for the reactions studied.
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.
SurfKin: an ab initio kinetic code for modeling surface reactions.
Le, Thong Nguyen-Minh; Liu, Bin; Huynh, Lam K
2014-10-01
In this article, we describe a C/C++ program called SurfKin (Surface Kinetics) to construct microkinetic mechanisms for modeling gas-surface reactions. Thermodynamic properties of reaction species are estimated based on density functional theory calculations and statistical mechanics. Rate constants for elementary steps (including adsorption, desorption, and chemical reactions on surfaces) are calculated using the classical collision theory and transition state theory. Methane decomposition and water-gas shift reaction on Ni(111) surface were chosen as test cases to validate the code implementations. The good agreement with literature data suggests this is a powerful tool to facilitate the analysis of complex reactions on surfaces, and thus it helps to effectively construct detailed microkinetic mechanisms for such surface reactions. SurfKin also opens a possibility for designing nanoscale model catalysts. PMID:25111729
Influence of ab initio chemistry models on simulations of the Ionian atmosphere
NASA Astrophysics Data System (ADS)
Parsons, Neal; Levin, Deborah A.; Walker, Andrew C.; Moore, Chris H.; Goldstein, David B.; Varghese, Philip L.; Trafton, Laurence
2014-09-01
There is significant scientific interest in simulating the unique atmospheric conditions on the Jovian moon Io that range from cold surface temperatures to hyperthermal interactions which possibly supply the Jovian plasma torus. The Direct Simulation Monte Carlo (DSMC) method is well suited to model the rarefied, predominantly SO2, Ionian atmosphere. High speed collisions between SO2 and the hypervelocity O atoms and ions that compose the plasma torus are a significant mechanism in determining the composition of the atmosphere; therefore, high-fidelity modeling of their interactions is crucial to the accuracy of such simulations. Typically, the Total Collision Energy (TCE) model is used to determine molecular dissociation probabilities and the Variable Hard Sphere (VHS) model is used to determine collision cross sections. However, the parameters for each of these baseline models are based on low-temperature experimental data and thus have unknown reliability for the hyperthermal conditions in the Ionian atmosphere. Recently, Molecular Dynamics/Quasi-Classical Trajectory (MD/QCT) studies have been conducted to generate accurate collision and chemistry models for the SO2-O collision pair in order to replace the baseline models. However, the influence of MD/QCT models on Ionian simulations compared to the previously used models is not well understood. In this work, 1D simulations are conducted using both the MD/QCT-based and baseline models in order to determine the effect of MD/QCT models on Ionian simulations. It is found that atmospheric structure predictions are highly sensitive to the chemistry and collision models. Specifically, the MD/QCT model predicts approximately half the SO2 atmospheric dissociation due to O and O+ bombardment compared to TCE models, and also predicts a temperature rise due to plasma heating further from the Ionian surface than the existing baseline methodologies. These findings indicate that the accurate MD/QCT chemistry and collision
NASA Astrophysics Data System (ADS)
Yung, Y. Y.; Choi, B. H.; Tang, K. T.
1980-01-01
Three dimensional quantum mechanical calculations are carried out for the reactive scattering of D+H2→DH+H on the ab initio potential energy surface calculated by Liu and Siegbahn and fitted by Truhlar and Horowitz. The differential and total cross sections as well as the S matrix elements are obtained from the adiabatic distorted wave method. Threshold energy, cross sections and product distributions over final states are all in good agreement with experimental measurements. Results are also compared with the corresponding ones obtained on the Porter-Karplus and the Yates-Lester semi-empirical surfaces.
Ab initio cluster model study of the chemisorptions of CO on low-index platinum surfaces
Curulla, D.; Clotet, A.; Ricart, J.M.; Illas, F.
1999-06-24
A systematic theoretical study of the adsorption of CO on the Pt{l_brace}100{r_brace}, Pt{l_brace}110{r_brace}, and Pt{l_brace}111{r_brace} surfaces is presented. The calculated equilibrium geometries and vibrational frequencies have been found to be rather independent of the cluster model chosen to represent the surface. However, calculated interaction energies are found to be very sensitive to the surface cluster model. The analysis of the chemisorption bond has been carried out by means of the constrained space orbital variation, CSOV, and of projection operator techniques. These analysis reveal that the bonding interactions are dominated by the {pi}-back-donation although {sigma}-donation plays a significant role. It is also clearly shown that all bonding mechanisms, other than Pauli repulsion, but specially {pi}-back-donation, contribute to the observed red shift. However, the {pi}-back-donation contribution to the red shift is very similar for CO on different sites. Hence, {pi}-back-donation cannot be the mechanism responsible for the observed difference for the CO vibrational frequency on on-top and bridge sites. The CSOV decomposition reveals that the leading term contributing to this difference in vibrational frequency of chemisorbed CO is the initial Pauli repulsion or wall effect; this is a new, important and unexpected conclusion.
Alexandrov, Vitali Y.; Rosso, Kevin M.
2015-01-01
Goethite (α-FeOOH) surfaces represent one of the most ubiquitous redox-active interfaces in the environment, playing an important role in biogeochemical metal cycling and contaminant residence in the subsurface. Fe(II)-catalyzed recrystallization of goethite is a fundamental process in this context, but the proposed Fe(II)aq-Fe(III)goethite electron and iron atom exchange mechanism of recrystallization remains poorly understood at the atomic level. We examine the adsorption of aqueous Fe(II) and subsequent interfacial electron transfer (ET) between adsorbed Fe(II) and structural Fe(III) at the (110) and (021) goethite surfaces using density functional theory calculations including Hubbard U corrections (DFT+U) aided by ab initio molecular dynamics simulations. We investigate various surface sites for the adsorption of Fe2+(H2O)6 in different coordination environments. Calculated energies for adsorbed complexes at both surfaces favor monodentate complexes with reduced 4- and 5-fold coordination over higher-dentate structures and 6- fold coordination. The hydrolysis of H2O ligands is observed for some pre-ET adsorbed Fe(II) configurations. ET from the adsorbed Fe(II) into the goethite lattice is calculated to be energetically uphill always, but simultaneous proton transfer from H2O ligands of the adsorbed complexes to the surface oxygen species stabilizes post-ET states. We find that surface defects such as oxygen vacancies near the adsorption site also can stabilize post-ET states, enabling the Fe(II)aq-Fe(III)goethite interfacial electron transfer reaction implied from experiments to proceed.
Ab Initio Modeling of Bulk and Intragranular Diffusion in Ni Alloys
Alexandrov, Vitali Y.; Sushko, Maria L.; Schreiber, Daniel K.; Bruemmer, Stephen M.; Rosso, Kevin M.
2015-05-07
importance for understanding mechanisms of grain boundary (GB) oxidation causing environmental degradation and cracking of Ni-base structural alloys. In this study, first-principles calculations of vacancy-mediated diffusion are performed across a wide series of alloying elements commonly used in Ni-based superalloys, as well as interstitial diffusion of atomic oxygen and sulfur in the bulk, at the (111) surface, <110> symmetric tilt GBs of Ni corresponding to model low- (Σ=3/(111)) and high-energy (Σ=9/(221)) GBs. A substantial enhancement of diffusion is found for all species at the high-energy GB as compared to the bulk and the low-energy GB, with Cr, Mn and Ti exhibiting remarkably small activation barriers (<0.1 eV; ~10 times lower than in the bulk). Calculations also show that the bulk diffusion mechanism and kinetics differ for oxygen and sulfur, with oxygen having a faster mobility and preferentially diffusing through the tetrahedral interstitial sites in Ni matrix where it can be trapped in a local minimum.
Ab Initio Modeling of Bulk and Intragranular Diffusion in Ni Alloys.
Alexandrov, Vitaly; Sushko, Maria L; Schreiber, Daniel K; Bruemmer, Stephen M; Rosso, Kevin M
2015-05-01
Knowledge of solid-state and interfacial species diffusion kinetics is of paramount importance for understanding mechanisms of grain boundary (GB) oxidation causing environmental degradation and cracking of Ni-base structural alloys. In this study, first-principles calculations of vacancy-mediated diffusion are performed across a wide series of alloying elements commonly used in Ni-based superalloys, as well as interstitial diffusion of atomic oxygen and sulfur in the bulk, at the (111) surface, ⟨110⟩ symmetric tilt GBs of Ni corresponding to model low- (Σ = 3/(111)) and high-energy (Σ = 9/(221)) GBs. A substantial enhancement of diffusion is found for all species at the high-energy GB as compared with the bulk and the low-energy GB, with Cr, Mn, and Ti exhibiting remarkably small activation barriers (<0.1 eV; ~10 times lower than in the bulk). Calculations also show that the bulk diffusion mechanism and kinetics differ for oxygen and sulfur, with oxygen having a faster mobility and preferentially diffusing through the tetrahedral interstitial sites in Ni matrix, where it can be trapped in a local minimum. PMID:26263324
Mo;ecular Modeling of Biogenic Manganese Oxides Using ab Initio Density Functional Theory
NASA Astrophysics Data System (ADS)
Oconnor, M.; Sposito, G.; Refson, K.
2003-12-01
Layer type manganese oxides with short-range crystalline order (birnessites) are produced by many species of bacteria.Deposits of these oxides form a highly reactive catalytic surface that plays a major role in the destruction and sequestration of organic compounds and metals.Biogenic oxides also contain vacant Mn(IV)sites;these sites,with their associated negative charge, are the probable main cause of the high sorptive reactivity of the oxide surfaces. In order to acquire a deeper understanding of the molecular mechanisms involved in these processes, a model of a biogenic oxide was built and its structure was optimized using the CASTEP three-dimensional periodic system computational package. The resulting crystal structure shows good agreement with EXAFS data from crystals formed by a strain of the common soil and freshwater bacterium, Pseudomonas putida. The greatest challenge in modeling Mn oxides (like other transition metal oxides)comes in dealing with the electronic factors that lead to their magnetic and catalytic properties: they are highly correlated systems where the spin must be taken into account in order to obtain accurate predictions of their properties.
Obtaining model parameters for real materials from ab-initio calculations: Heisenberg exchange
NASA Astrophysics Data System (ADS)
Korotin, Dmitry; Mazurenko, Vladimir; Anisimov, Vladimir; Streltsov, Sergey
An approach to compute exchange parameters of the Heisenberg model in plane-wave based methods is presented. This calculation scheme is based on the Green's function method and Wannier function projection technique. It was implemented in the framework of the pseudopotential method and tested on such materials as NiO, FeO, Li2MnO3, and KCuF3. The obtained exchange constants are in a good agreement with both the total energy calculations and experimental estimations for NiO and KCuF3. In the case of FeO our calculations explain the pressure dependence of the Néel temperature. Li2MnO3 turns out to be a Slater insulator with antiferromagnetic nearest neighbor exchange defined by the spin splitting. The proposed approach provides a unique way to analyze magnetic interactions, since it allows one to calculate orbital contributions to the total exchange coupling and study the mechanism of the exchange coupling. The work was supported by a grant from the Russian Scientific Foundation (Project No. 14-22-00004).
Magnetism in Sr2CrMoO6 : A combined ab initio and model study
NASA Astrophysics Data System (ADS)
Sanyal, Prabuddha; Halder, Anita; Si, Liang; Wallerberger, Markus; Held, Karsten; Saha-Dasgupta, Tanusri
2016-07-01
Using a combination of first-principles density functional theory (DFT) calculations and exact diagonalization studies of a first-principles derived model, we carry out a microscopic analysis of the magnetic properties of the half-metallic double perovskite compound Sr2CrMoO6 , a sister compound of the much discussed material Sr2FeMoO6 . The electronic structure of Sr2CrMoO6 , though appearing similar to Sr2FeMoO6 at first glance, shows nontrivial differences with that of Sr2FeMoO6 on closer examination. In this context, our study highlights the importance of charge transfer energy between the two transition metal sites. The change in charge transfer energy due to a shift of Cr d states in Sr2CrMoO6 compared to Fe d in Sr2FeMoO6 suppresses the hybridization between Cr t2 g and Mo t2 g. This strongly weakens the hybridization-driven mechanism of magnetism discussed for Sr2FeMoO6 . Our study reveals that, nonetheless, the magnetic transition temperature of Sr2CrMoO6 remains high since an additional superexchange contribution to magnetism arises with a finite intrinsic moment developed at the Mo site. We further discuss the situation in comparison to another related double perovskite compound, Sr2CrWO6 . We also examine the effect of correlation beyond DFT, using dynamical mean field theory.
Ab initio modeling of carbohydrates: on the proper selection of computational methods and basis sets
Technology Transfer Automated Retrieval System (TEKTRAN)
With the development of faster computer hardware and quantum mechanical software it has become more feasible to study large carbohydrate molecules via quantum mechanical methods. In the past, studies of carbohydrates were restricted to empirical/semiempirical methods and Hartree Fock. In the last ...
Phenylalanine ab initio models for the simulation of skin natural moisturizing factor
NASA Astrophysics Data System (ADS)
Carvalho, B. G.; Raniero, L. J.; Martin, A. A.; Favero, P. P.
2013-04-01
In this study, we evaluated models that can be used to simulate amino acids in biological environments via density functional theory (DFT). The goal was to obtain realistic representations that combine computational economy and result quality when compared to experimental data. We increased the complexity of the models by using a model of an amino acid in a vacuum, followed by a water-solvated amino acid model. To consider pH variation, we simulated zwitterionic and nonionic amino acid configurations. The amino acid chosen for testing was phenylalanine, an aromatic amino acid present in high concentrations in the natural moisturizing factor of skin that plays a fundamental role in ultraviolet protection and vitiligo disease. To validate the models, vibrational modes and electronic properties were calculated and compared to experimental results.
Symplectic Symmetry and the Ab Initio No-Core Shell Model
Draayer, Jerry P.; Dytrych, Tomas; Sviratcheva, Kristina D.; Bahri, Chairul; Vary, James P.; /Iowa State U. /LLNL, Livermore /SLAC
2007-03-14
The symplectic symmetry of eigenstates for the 0{sub gs}{sup +} in {sup 16}O and the 0{sub gs}{sup +} and lowest 2{sup +} and 4{sup +} configurations of {sup 12}C that are well-converged within the framework of the no-core shell model with the JISP16 realistic interaction is examined. These states are found to project at the 85-90% level onto very few symplectic representations including the most deformed configuration, which confirms the importance of a symplectic no-core shell model and reaffirms the relevance of the Elliott SU(3) model upon which the symplectic scheme is built.
Evidence for Symplectic Symmetry in Ab Initio No-Core Shell Model Results for Light Nuclei
Dytrych, Tomas; Sviratcheva, Kristina D.; Bahri, Chairul; Draayer, Jerry P.; Vary, James P.; /Iowa State U. /LLNL, Livermore /SLAC
2007-04-24
Clear evidence for symplectic symmetry in low-lying states of {sup 12}C and {sup 16}O is reported. Eigenstates of {sup 12}C and {sup 16}O, determined within the framework of the no-core shell model using the JISP16 NN realistic interaction, typically project at the 85-90% level onto a few of the most deformed symplectic basis states that span only a small fraction of the full model space. The results are nearly independent of whether the bare or renormalized effective interactions are used in the analysis. The outcome confirms Elliott's SU(3) model which underpins the symplectic scheme, and above all, points to the relevance of a symplectic no-core shell model that can reproduce experimental B(E2) values without effective charges as well as deformed spatial modes associated with clustering phenomena in nuclei.
Ab Initio Study of 40Ca with an Importance Truncated No-Core Shell Model
Roth, R; Navratil, P
2007-05-22
We propose an importance truncation scheme for the no-core shell model, which enables converged calculations for nuclei well beyond the p-shell. It is based on an a priori measure for the importance of individual basis states constructed by means of many-body perturbation theory. Only the physically relevant states of the no-core model space are considered, which leads to a dramatic reduction of the basis dimension. We analyze the validity and efficiency of this truncation scheme using different realistic nucleon-nucleon interactions and compare to conventional no-core shell model calculations for {sup 4}He and {sup 16}O. Then, we present the first converged calculations for the ground state of {sup 40}Ca within no-core model spaces including up to 16{h_bar}{Omega}-excitations using realistic low-momentum interactions. The scheme is universal and can be easily applied to other quantum many-body problems.
High throughput ab initio modeling of charge transport for bio-molecular-electronics
NASA Astrophysics Data System (ADS)
Bruque, Nicolas Alexander
2009-12-01
Self-assembled nanostructures, composed of inorganic and organic materials, have multiple applications in the fields of engineering and nanotechnology. Experimental research using nanoscaled materials, such as semiconductor/metallic nanocrystals, nanowires (NW), and carbon nanotube (CNT)-molecular systems have potential applications in next generation nano electronic devices. Many of these molecular systems exhibit electronic device functionality. However, experimental analytical techniques to determine how the chemistry and geometry affects electron transport through these devices does not yet exist. Using theory and modeling, one can approximate the chemistry and geometry at the atomic level and also determine how the chemistry and geometry governs electron current. Nanoelectronic devices however, contain several thousand atoms which makes quantum modeling difficult. Popular atomistic modeling approaches are capable of handling small molecular systems, which are of scientific interest, but have little engineering value. The lack of large scale modeling tools has left the scientific and engineering community with a limited ability to understand, explore, and design complex systems of engineering interest. To address these issues, I have developed a high performance general quantum charge transport model based on the non-equilibrium Green function (NEGF) formalism using density functional theory (DFT) as implemented in the FIREBALL software. FIREBALL is a quantum molecular dynamics code which has demonstrated the ability to model large molecular systems. This dissertation project of integrating NEGF into FIREBALL provides researchers with a modeling tool capable of simulating charge current in large inorganic/organic systems. To provide theoretical support for experimental efforts, this project focused on CNT-molecular systems, which includes the discovery of a CNT-molecular resonant tunneling diode (RTD) for electronic circuit applications. This research also
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
Effective Operators Within the Ab Initio No-Core Shell Model
Stetcu, I; Barrett, B R; Navratil, P; Vary, J P
2004-11-30
We implement an effective operator formalism for general one- and two-body operators, obtaining results consistent with the no-core shell model (NCSM) wave functions. The Argonne V8' nucleon-nucleon potential was used in order to obtain realistic wave functions for {sup 4}He, {sup 6}Li and {sup 12}C. In the NCSM formalism, we compute electromagnetic properties using the two-body cluster approximation for the effective operators and obtain results which are sensitive to the range of the bare operator. To illuminate the dependence on the range, we employ a Gaussian two-body operator of variable range, finding weak renormalization of long range operators (e.g., quadrupole) in a fixed model space. This is understood in terms of the two-body cluster approximation which accounts mainly for short-range correlations. Consequently, short range operators, such as the relative kinetic energy, will be well renormalized in the two-body cluster approximation.
The ab initio model potential method. Second series transition metal elements
Barandiaran, Z.; Seijo, L. ); Huzinaga, S. )
1990-10-15
The {ital ab} {ital initio} core method potential model (AIMP) has already been presented in its nonrelativistic version and applied to the main group and first series transition metal elements (J. Chem. Phys. {bold 86}, 2132 (1987); {bold 91}, 7011 (1989)). In this paper we extend the AIMP method to include relativistic effects within the Cowan--Griffin approximation and we present relativistic Zn-like core model potentials and valence basis sets, as well as their nonrelativistic Zn-like core and Kr-like core counterparts. The pilot molecular calculations on YO, TcO, AgO, and AgH reveal that the 4{ital p} orbital is indeed a core orbital only at the end part of the series, whereas the 4{ital s} orbital can be safely frozen from Y to Cd. The all-electron and model potential results agree in 0.01--0.02 A in {ital R}{sub {ital e}} and 25--50 cm{sup {minus}1} in {bar {nu}}{sub {ital e}} if the same type of valence part of the basis set is used. The comparison of the relativistic results on AgH with those of the all-electron Dirac--Fock calculations by Lee and McLean is satisfactory: the absolute value of {ital R}{sub {ital e}} is reproduced within the 0.01 A margin and the relativistic contraction of 0.077 A is also very well reproduced (0.075 A). Finally, the relative magnitude of the effects of the core orbital change, mass--velocity potential, and Darwin potential on the net relativistic effects are analyzed in the four molecules studied.
Ab initio transport coefficients of Ar+ ions in Ar for cold plasma jet modeling
NASA Astrophysics Data System (ADS)
Chicheportiche, A.; Lepetit, B.; Gadéa, F. X.; Benhenni, M.; Yousfi, M.; Kalus, R.
2014-06-01
Collision cross sections and transport coefficients are calculated for Ar+ ions, in the ground state 2P3/2 and in the metastable state 2P1/2, colliding with their parent gas. Differential and integral collision cross sections are obtained using a numerical integration of the nuclear Schrödinger equation for several published interaction potentials. The Cohen-Schneider semi-empirical model is used for the inclusion of the spin-orbit interaction. The corresponding differential collision cross sections are then used in an optimized Monte Carlo code to calculate the ion transport coefficients for each initial ion state over a wide range of reduced electric field. Ion swarm data results are then compared with available experimental data for different proportions of ions in each state. This allows us to identify the most reliable interaction potential which reproduces ion transport coefficients falling within the experimental error bars. Such ion transport data will be used in electrohydrodynamic and chemical kinetic models of the low temperature plasma jet to quantify and to tune the active species production for a better use in biomedical applications.
Ab initio transport coefficients of Ar⁺ ions in Ar for cold plasma jet modeling.
Chicheportiche, A; Lepetit, B; Gadéa, F X; Benhenni, M; Yousfi, M; Kalus, R
2014-06-01
Collision cross sections and transport coefficients are calculated for Ar{+} ions, in the ground state {2}P_{3/2} and in the metastable state {2}P_{1/2}, colliding with their parent gas. Differential and integral collision cross sections are obtained using a numerical integration of the nuclear Schrödinger equation for several published interaction potentials. The Cohen-Schneider semi-empirical model is used for the inclusion of the spin-orbit interaction. The corresponding differential collision cross sections are then used in an optimized Monte Carlo code to calculate the ion transport coefficients for each initial ion state over a wide range of reduced electric field. Ion swarm data results are then compared with available experimental data for different proportions of ions in each state. This allows us to identify the most reliable interaction potential which reproduces ion transport coefficients falling within the experimental error bars. Such ion transport data will be used in electrohydrodynamic and chemical kinetic models of the low temperature plasma jet to quantify and to tune the active species production for a better use in biomedical applications. PMID:25019899
Ab initio approach to the non-perturbative scalar Yukawa model
NASA Astrophysics Data System (ADS)
Li, Yang; Karmanov, V. A.; Maris, P.; Vary, J. P.
2015-09-01
We report on the first non-perturbative calculation of the scalar Yukawa model in the single-nucleon sector up to four-body Fock sector truncation (one "scalar nucleon" and three "scalar pions"). The light-front Hamiltonian approach with a systematic non-perturbative renormalization is applied. We study the n-body norms and the electromagnetic form factor. We find that the one- and two-body contributions dominate up to coupling α ≈ 1.7. As we approach the coupling α ≈ 2.2, we discover that the four-body contribution rises rapidly and overtakes the two- and three-body contributions. By comparing with lower sector truncations, we show that the form factor converges with respect to the Fock sector expansion.
Ab initio electronic structure study of a model water splitting dimer complex.
Fernando, Amendra; Aikens, Christine M
2015-12-28
A model manganese dimer electrocatalyst bridged by μ-OH ligands is used to investigate changes in spin states that may occur during water oxidation. We have employed restricted open-shell Hartree-Fock (ROHF), second-order Møller-Plesset perturbation theory (MP2), complete active space self-consistent field (CASSCF), and multireference second-order Møller-Plesset perturbation theory (MRMP2) calculations to investigate this system. Multiconfigurational methods like CASSCF and MRMP2 are appropriate methods to study these systems with antiferromagnetically-coupled electrons. Orbital occupations and distributions have been closely analyzed to understand the electronic details and contributions to the water splitting from manganese and oxygen atoms. The presence of Mn(IV)O˙ radical moieties has been observed in this catalytic pathway. Multiple nearly degenerate excited states were found close to the ground state in all structures. This suggests competing potential energy landscapes near the ground state may influence the reactivity of manganese complexes such as the dimers studied in this work. PMID:26593689
NASA Astrophysics Data System (ADS)
Ramírez, Benjamín R.; Ghoniem, Nasr; Po, Giacomo
2012-09-01
We develop a model of cross-slip in face-centered cubic (fcc) metals based on an extension of the Peierls-Nabarro representation of the dislocation core. The dissociated core is described by a group of parametric fractional Volterra dislocations, subject to their mutual elastic interaction and a lattice-restoring force. The elastic interaction between them is computed from a nonsingular expression, while the lattice force is derived from the γ surface obtained directly from ab initio calculations. Using a network-based formulation of dislocation dynamics, the dislocation core structure is not restricted to be planar, and the activation energy is determined for a path where the core has three-dimensional equilibrium configurations. We show that the activation energy for cross-slip in Cu is 1.9eV when the core is represented by only two Shockley partials, while this value converges to 1.43eV when the core is distributed over a bundle of 20 Volterra partial fractional dislocations. The results of the model compare favorably with the experimental value of 1.15±0.37eV [J. Bonneville and B. Escaig, Acta Metall.AMETAR0001-616010.1016/0001-6160(79)90170-6 27, 1477 (1979)]. We also show that the cross-slip activation energy decreases significantly when the core is in a particular local stress field. Results are given for a representative uniform “Escaig” stress and for the nonuniform stress field at the head of a dislocation pileup. A local homogeneous stress field is found to result in a significant reduction of the cross-slip energy. Additionally, for a nonhomogeneous stress field at the head of a five-dislocation pileup compressed against a Lomer-Cottrell junction, the cross-slip energy is found to decrease to 0.62eV. The relatively low values of the activation energy in local stress fields predicted by the proposed model suggest that cross-slip events are energetically more favorable in strained fcc crystals.
NASA Astrophysics Data System (ADS)
Mei, Yuan; Sherman, David M.; Liu, Weihua; Etschmann, Barbara; Testemale, Denis; Brugger, Joël
2015-02-01
The solubility of zinc minerals in hydrothermal fluids is enhanced by chloride complexation of Zn2+. Thermodynamic models of these complexation reactions are central to models of Zn transport and ore formation. However, existing thermodynamic models, derived from solubility measurements, are inconsistent with spectroscopic measurements of Zn speciation. Here, we used ab initio molecular dynamics simulations (with the PBE exchange-correlation functional) to predict the speciation of Zn-Cl complexes from 25 to 600 °C. We also obtained in situ XAS measurements of Zn-Cl solutions at 30-600 °C. Qualitatively, the simulations reproduced the main features derived from in situ XANES and EXAFS measurements: octahedral to tetrahedral transition with increasing temperature and salinity, stability of ZnCl42- at high chloride concentration up to ⩾500 °C, and increasing stability of the trigonal planar [ZnCl3]- complex at high temperature. Having confirmed the dominant species, we directly determined the stability constants for the Zn-Cl complexes using thermodynamic integration along constrained Zn-Cl distances in a series of MD simulations. We corrected our stability constants to infinite dilution using the b-dot model for the activity coefficients of the solute species. In order to compare the ab initio results with experiments, we need to re-model the existing solubility data using the species we identified in our MD simulations. The stability constants derived from refitting published experimental data are in reasonable agreement with those we obtained using ab initio MD simulations. Our new thermodynamic model accurately predicts the experimentally observed changes in ZnO(s) and ZnCO3(s) solubility as a function of chloride concentration from 200 (Psat) to 600 °C (2000 bar). This study demonstrates that metal speciation and geologically useful stability constants can be derived for species in hydrothermal fluids from ab initio MD simulations even at the generalized
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.
Clerc, Daryl G
2016-07-21
An ab initio approach was used to study the molecular-level interactions that connect gene-mutation to changes in an organism׳s phenotype. The study provides new insights into the evolutionary process and presents a simplification whereby changes in phenotypic properties may be studied in terms of the binding affinities of the chemical interactions affected by mutation, rather than by correlation to the genes. The study also reports the role that nonlinear effects play in the progression of organs, and how those effects relate to the classical theory of evolution. Results indicate that the classical theory of evolution occurs as a special case within the ab initio model - a case having two attributes. The first attribute: proteins and promoter regions are not shared among organs. The second attribute: continuous limiting behavior exists in the physical properties of organs as well as in the binding affinity of the associated chemical interactions, with respect to displacements in the chemical properties of proteins and promoter regions induced by mutation. Outside of the special case, second-order coupling contributions are significant and nonlinear effects play an important role, a result corroborated by analyses of published activity levels in binding and transactivation assays. Further, gradations in the state of perfection of an organ may be small or large depending on the type of mutation, and not necessarily closely-separated as maintained by the classical theory. Results also indicate that organs progress with varying degrees of interdependence, the likelihood of successful mutation decreases with increasing complexity of the affected chemical system, and differences between the ab initio model and the classical theory increase with increasing complexity of the organism. PMID:27029513
NASA Astrophysics Data System (ADS)
Farahani, Pooria; Lundberg, Marcus; Karlsson, Hans O.
2013-11-01
The SN2 substitution reactions at phosphorus play a key role in organic and biological processes. Quantum molecular dynamics simulations have been performed to study the prototype reaction Cl-+PH2Cl→ClPH2+Cl-, using one and two-dimensional models. A potential energy surface, showing an energy well for a transition complex, was generated using ab initio electronic structure calculations. The one-dimensional model is essentially reflection free, whereas the more realistic two-dimensional model displays involved resonance structures in the reaction probability. The reaction rate is almost two orders of magnitude smaller for the two-dimensional compared to the one-dimensional model. Energetic errors in the potential energy surface is estimated to affect the rate by only a factor of two. This shows that for these types of reactions it is more important to increase the dimensionality of the modeling than to increase the accuracy of the electronic structure calculation.
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
Reply to Comment on "Ab Initio Study of 40Ca with an Importance Truncated No-Core Shell Model"
Roth, R; Navratil, P
2008-01-04
In their comment on our recent Letter [1] Dean et al. [2] criticize the calculations for the ground-state energy of {sup 40}Ca within the importance truncated no-core shell model (NCSM). In particular they address the role of configurations beyond the 3p3h level, which have not been included in the {sup 40}Ca calculations for large N{sub max} {h_bar}{Omega} model spaces. Before responding to this point, the following general statements are in order. For the atomic nucleus as a self-bound system, translational invariance is an important symmetry. The only possibility to preserve translational invariance when working with a Slater determinant basis is to use the harmonic oscillator (HO) basis in conjunction with a basis truncation according to the total HO excitation energy, i.e. N{sub max} {h_bar}{Omega}, as done in the ab initio NCSM. This is important not only for obtaining proper binding or excitation energies, but also for a correct extraction of physical wavefunctions. The spurious center-of-mass components can be exactly removed only if the HO basis and the N{sub max} {h_bar}{Omega} truncation are employed. The minimal violation of the translational invariance was one of the main motivations for developing the importance-truncation scheme introduced in the Letter. In this scheme, we start with the complete N{sub max} {h_bar}{Omega} HO basis space and select important configurations via perturbation theory. All symmetries are under control and our importance-truncated NCSM calculations are completely variational and provide an upper bound of the ground-state energy of the system. The restriction to the 3p3h level, made for computational reasons in the N{sub max} > 8 calculations for {sup 40}Ca, is not inherent to the importance truncation scheme. The explicit inclusion of 4p4h configurations--though computationally more demanding--is straight-forward, even for the largest N{sub max} {h_bar}{Omega} model spaces discussed. To demonstrate this fact we have
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
NASA Astrophysics Data System (ADS)
Wang, Yaocen; Takeuchi, Akira; Makino, Akihiro; Liang, Yunye; Kawazoe, Yoshiyuki
2014-05-01
Iron-based amorphous and nano-crystalline alloys have attracted a growing interest due to their potential in the application of magnetic coil production. However, fundamental understanding of the nano-crystallization mechanisms and magnetic features in the amorphous structure are still lack of knowledge. In the present work, we performed ab initio molecular dynamics simulation to clarify the ionic and electronic structure in atomic scale, and to derive the origin of the good magnetic property of Fe85Si2B8P4Cu1 amorphous alloy. The simulation gave a direct evidence of the Cu-P bonding preference in the amorphous alloy, which may promote nucleation in nano-crystallization process. On the other hand, the electron transfer and the band/orbital features in the amorphous alloy suggests that alloying elements with large electronegativity and the potential to expand Fe disordered matrix are preferred for enhancing the magnetization.
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
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
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.
Matsubara, Hiroshi; Horvat, Sonia M; Schiesser, Carl H
2003-04-01
Ab initio calculations using 6-311G**, cc-pVDZ, aug-cc-pVDZ, and a (valence) double-zeta pseudopotential (DZP) basis sets, with (MP2, QCISD, CCSD(T)) and without (UHF) the inclusion of electron correlation, and density functional (B3LYO) calculations predict that homolytic substitution reactions of the methyl radical at the silicon atom in disilane can proceed via both backside and frontside attack mechanisms. At the highest level of theory (CCSD(T)/aug-cc-pVDZ//MP2/aug-cc-pVDZ), energy barriers (delta E) of 47.4 and 48.6 kJ mol-1 are calculated for the backside and frontside reactions respectively. Similar results are obtained for reactions involving germanium and tin with energy barriers (delta E) of between 46.5 and 67.3, and 41.0 and 73.3 kJ mol-1 for the backside and frontside mechanisms, respectively. These data suggest that homolytic substitution reactions of methyl radical at silicon, germanium, and tin can proceed via either homolytic substitution mechanism. PMID:12926395
NASA Astrophysics Data System (ADS)
Zeng, Xiancheng; Hu, Hao; Hu, Xiangqian; Yang, Weitao
2009-04-01
A quantum mechanical/molecular mechanical minimum free energy path (QM/MM-MFEP) method was developed to calculate the redox free energies of large systems in solution with greatly enhanced efficiency for conformation sampling. The QM/MM-MFEP method describes the thermodynamics of a system on the potential of mean force surface of the solute degrees of freedom. The molecular dynamics (MD) sampling is only carried out with the QM subsystem fixed. It thus avoids "on-the-fly" QM calculations and thus overcomes the high computational cost in the direct QM/MM MD sampling. In the applications to two metal complexes in aqueous solution, the new QM/MM-MFEP method yielded redox free energies in good agreement with those calculated from the direct QM/MM MD method. Two larger biologically important redox molecules, lumichrome and riboflavin, were further investigated to demonstrate the efficiency of the method. The enhanced efficiency and uncompromised accuracy are especially significant for biochemical systems. The QM/MM-MFEP method thus provides an efficient approach to free energy simulation of complex electron transfer reactions.
Zeng Xiancheng; Hu Hao; Hu Xiangqian; Yang Weitao
2009-04-28
A quantum mechanical/molecular mechanical minimum free energy path (QM/MM-MFEP) method was developed to calculate the redox free energies of large systems in solution with greatly enhanced efficiency for conformation sampling. The QM/MM-MFEP method describes the thermodynamics of a system on the potential of mean force surface of the solute degrees of freedom. The molecular dynamics (MD) sampling is only carried out with the QM subsystem fixed. It thus avoids 'on-the-fly' QM calculations and thus overcomes the high computational cost in the direct QM/MM MD sampling. In the applications to two metal complexes in aqueous solution, the new QM/MM-MFEP method yielded redox free energies in good agreement with those calculated from the direct QM/MM MD method. Two larger biologically important redox molecules, lumichrome and riboflavin, were further investigated to demonstrate the efficiency of the method. The enhanced efficiency and uncompromised accuracy are especially significant for biochemical systems. The QM/MM-MFEP method thus provides an efficient approach to free energy simulation of complex electron transfer reactions.
NASA Astrophysics Data System (ADS)
Cottom, J.; Gruber, G.; Hadley, P.; Koch, M.; Pobegen, G.; Aichinger, T.; Shluger, A.
2016-05-01
Electrically detected magnetic resonance (EDMR) is a powerful technique for the observation and categorization of paramagnetic defects within semiconductors. The interpretation of the recorded EDMR spectra has long proved to be challenging. Here, defect spectra are identified by comparing EDMR measurements with extensive ab initio calculations. The defect identification is based upon the defect symmetry and the form of the hyperfine (HF) structure. A full description is given of how an accurate spectrum can be generated from the theoretical data by considering some thousand individual HF contributions out of some billion possibilities. This approach is illustrated with a defect observed in nitrogen implanted silicon carbide (SiC). Nitrogen implantation is a high energy process that gives rise to a high defect concentration. The majority of these defects are removed during the dopant activation anneal, shifting the interstitial nitrogen to the desired substitutional lattice sites, where they act as shallow donors. EDMR shows that a deep-level defect persists after the dopant activation anneal. This defect is characterized as having a g c ∥ B = 2.0054 ( 4 ) and g c ⊥ B = 2.0006 ( 4 ) , with pronounced hyperfine shoulder peaks with a 13 G peak to peak separation. The nitrogen at a carbon site next to a silicon vacancy ( N C V Si ) center is identified as the persistent deep-level defect responsible for the observed EDMR signal and the associated dopant deactivation.
Atomistic and Ab initio modeling of CaAl2O4 high-pressure polymorphs under Earth's mantle conditions
NASA Astrophysics Data System (ADS)
Eremin, N. N.; Grechanovsky, A. E.; Marchenko, E. I.
2016-05-01
Semi-empirical and ab initio theoretical investigation of crystal structure geometry, interatomic distances, phase densities and elastic properties for some CaAl2O4 phases under pressures up to 200 GPa was performed. Two independent simulation methods predicted the appearance of a still unknown super-dense CaAl2O4 modification. In this structure, the Al coordination polyhedron might be described as distorted one with seven vertices. Ca atoms were situated inside polyhedra with ten vertices and Ca-O distances from 1.96 to 2.49 Å. It became the densest modification under pressures of 170 GPa (density functional theory prediction) or 150 GPa (semi-empirical prediction). Both approaches indicated that this super-dense CaAl2O4 modification with a "stuffed α-PbO2" type structure could be a probable candidate for mutual accumulation of Ca and Al in the lower mantle. The existence of this phase can be verified experimentally using high pressure techniques.
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.
Qian, X.; Nimlos, M. R.; Davis, M.; Johnson, D. K.; Himmel, M. E.
2005-01-01
Ab initio molecular dynamics simulations were employed to investigate, with explicit solvent water molecules, {beta}-d-glucose and {beta}-d-xylose degradation mechanisms in acidic media. The rate-limiting step in sugar degradation was found to be protonation of the hydroxyl groups on the sugar ring. We found that the structure of water molecules plays a significant role in the acidic sugar degradation pathways. Firstly, a water molecule competes with the hydroxyl group on the sugar ring for protons. Secondly, water forms hydrogen bonds with the hydroxyl groups on the sugar rings, thus weakening the C-C and C-O bonds (each to a different degree). Note that the reaction pathways could be altered due to the change of relative stability of the C-C and C-O bonds. Thirdly, water molecules that are hydrogen-bonded to sugar hydroxyls could easily extract a proton from the reaction intermediate, terminating the reaction. Indeed, the sugar degradation pathway is complex due to multiple protonation probabilities and the surrounding water structure. Our experimental data support multiple sugar acidic degradation pathways.
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
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.
NASA Astrophysics Data System (ADS)
Nattino, Francesco; Galparsoro, Oihana; Costanzo, Francesca; Díez Muiño, Ricardo; Alducin, Maite; Kroes, Geert-Jan
2016-06-01
Accurately modeling surface temperature and surface motion effects is necessary to study molecule-surface reactions in which the energy dissipation to surface phonons can largely affect the observables of interest. We present here a critical comparison of two methods that allow to model such effects, namely, the ab initio molecular dynamics (AIMD) method and the generalized Langevin oscillator (GLO) model, using the dissociation of N2 on W(110) as a benchmark. AIMD is highly accurate as the surface atoms are explicitly part of the dynamics, but this advantage comes with a large computational cost. The GLO model is much more computationally convenient, but accounts for lattice motion effects in a very approximate way. Results show that, despite its simplicity, the GLO model is able to capture the physics of the system to a large extent, returning dissociation probabilities which are in better agreement with AIMD than static-surface results. Furthermore, the GLO model and the AIMD method predict very similar energy transfer to the lattice degrees of freedom in the non-reactive events, and similar dissociation dynamics.
NASA Astrophysics Data System (ADS)
Cragnolini, Tristan; Derreumaux, Philippe; Pasquali, Samuela
2015-06-01
RNA molecules are essential cellular machines performing a wide variety of functions for which a specific three-dimensional structure is required. Over the last several years, the experimental determination of RNA structures through x-ray crystallography and NMR seems to have reached a plateau in the number of structures resolved each year, but as more and more RNA sequences are being discovered, the need for structure prediction tools to complement experimental data is strong. Theoretical approaches to RNA folding have been developed since the late nineties, when the first algorithms for secondary structure prediction appeared. Over the last 10 years a number of prediction methods for 3D structures have been developed, first based on bioinformatics and data-mining, and more recently based on a coarse-grained physical representation of the systems. In this review we are going to present the challenges of RNA structure prediction and the main ideas behind bioinformatic approaches and physics-based approaches. We will focus on the description of the more recent physics-based phenomenological models and on how they are built to include the specificity of the interactions of RNA bases, whose role is critical in folding. Through examples from different models, we will point out the strengths of physics-based approaches, which are able not only to predict equilibrium structures, but also to investigate dynamical and thermodynamical behavior, and the open challenges to include more key interactions ruling RNA folding.
Cragnolini, Tristan; Derreumaux, Philippe; Pasquali, Samuela
2015-06-17
RNA molecules are essential cellular machines performing a wide variety of functions for which a specific three-dimensional structure is required. Over the last several years, the experimental determination of RNA structures through x-ray crystallography and NMR seems to have reached a plateau in the number of structures resolved each year, but as more and more RNA sequences are being discovered, the need for structure prediction tools to complement experimental data is strong. Theoretical approaches to RNA folding have been developed since the late nineties, when the first algorithms for secondary structure prediction appeared. Over the last 10 years a number of prediction methods for 3D structures have been developed, first based on bioinformatics and data-mining, and more recently based on a coarse-grained physical representation of the systems. In this review we are going to present the challenges of RNA structure prediction and the main ideas behind bioinformatic approaches and physics-based approaches. We will focus on the description of the more recent physics-based phenomenological models and on how they are built to include the specificity of the interactions of RNA bases, whose role is critical in folding. Through examples from different models, we will point out the strengths of physics-based approaches, which are able not only to predict equilibrium structures, but also to investigate dynamical and thermodynamical behavior, and the open challenges to include more key interactions ruling RNA folding. PMID:25993396
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.
Li, Zhe; Wu, Yinuo; Feng, Ling-Jun; Wu, Ruibo; Luo, Hai-Bin
2014-12-01
Phosphodiesterases (PDEs) are the sole enzymes hydrolyzing the important second messengers cGMP and cAMP and have been identified as therapeutic targets for several diseases. The most successful examples are PDE5 inhibitors (i.e., sildenafil and tadalafil), which have been approved for the treatment of male erectile dysfunction and pulmonary hypertension. However, the side effects mostly due to nonselective inhibition toward other PDE isoforms, set back the clinical usage of PDE5 inhibitors. Until now, the exact catalytic mechanism of the substrate cGMP by PDE5 is still unclear. Herein, the first computational study on the catalytic hydrolysis mechanism of cGMP for PDE5 (catalytic domain) is performed by employing the state-of-the-art ab initio quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations. Our simulations show a SN2 type reaction procedure via a highly dissociated transition state with a reaction barrier of 8.88 kcal/mol, which is quite different from the previously suggested hydrolysis mechanism of cAMP for PDE4. Furthermore, the subsequent ligand exchange and the release of the product GMP have also been investigated by binding energy analysis and MD simulations. It is deduced that ligand exchange would be the rate-determining step of the whole reaction, which is consistent with many previous experimental results. The obtained mechanistic insights should be valuable for not only the rational design of more specific inhibitors toward PDE5 but also understanding the general hydrolysis mechanism of cGMP-specific PDEs. PMID:26583228
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
NASA Astrophysics Data System (ADS)
Abadias, G.; Kanoun, M. B.; Goumri-Said, S.; Koutsokeras, L.; Dub, S. N.; Djemia, Ph.
2014-10-01
The structure, phase stability, and mechanical properties of ternary alloys of the Zr-Ta-N system are investigated by combining thin-film growth and ab initio calculations. Zr1-xTaxN films with 0≤x≤1 were deposited by reactive magnetron cosputtering in Ar +N2 plasma discharge and their structural properties characterized by x-ray diffraction. We considered both ordered and disordered alloys, using supercells and special quasirandom structure approaches, to account for different possible metal atom distributions on the cation sublattice. Density functional theory within the generalized gradient approximation was employed to calculate the electronic structure as well as predict the evolution of the lattice parameter and key mechanical properties, including single-crystal elastic constants and polycrystalline elastic moduli, of ternary Zr1-xTaxN compounds with cubic rocksalt structure. These calculated values are compared with experimental data from thin-film measurements using Brillouin light scattering and nanoindentation tests. We also study the validity of Vegard's empirical rule and the effect of growth-dependent stresses on the lattice parameter. The thermal stability of these Zr1-xTaxN films is also studied, based on their structural and mechanical response upon vacuum annealing at 850 °C for 3 h. Our findings demonstrate that Zr1-xTaxN alloys with Ta fraction 0.51⩽x⩽0.78 exhibit enhanced toughness, while retaining high hardness ˜30 GPa, as a result of increased valence electron concentration and phase stability tuning. Calculations performed for disordered or ordered structures both lead to the same conclusion regarding the mechanical behavior of these nitride alloys, in agreement with recent literature findings [H. Kindlund, D. G. Sangiovanni, L. Martinez-de-Olcoz, J. Lu, J. Jensen, J. Birch, I. Petrov, J. E. Greene, V. Chirita, and L. Hultman, APL Materials 1, 042104 (2013), 10.1063/1.4822440].
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.
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.
NASA Astrophysics Data System (ADS)
Tu, Qing; Lange, Björn; Lopes, J. Marcelo J.; Zauscher, Stefan; Blum, Volker
Contact resonance AFM is demonstrated as a powerful tool for mapping differences in the mechanical properties of 2D materials and heterostructures, permitting to resolve surface and subsurface structural differences of different domains. Measured contact resonance frequencies are related to the contact stiffness of the combined tip-sample system. Based on first principles predicted elastic properties and a continuum approach to model the mechanical impedance, we find contact stiffness ratios between different domains of few-layer graphene on 3C-SiC(111) in excellent agreement with experiment. We next demonstrate that the approach is able to quantitatively resolve differences between other 2D materials domains, e.g., for h-BN, MoS2 and MoO3 on graphene on SiC. We show that the combined effect of several materials parameters, especially the in-plane elastic properties and the layer thickness, determines the contact stiffness, therefore boosting the sensitivity even if the out-of-plane elastic properties are similar.
Tachikawa, Hiroto
2014-06-01
The mechanism of dissolution of the Li(+) ion in an electrolytic solvent is investigated by the direct ab initio molecular dynamics (AIMD) method. Lithium fluoroborate (Li(+)BF4(-)) and ethylene carbonate (EC) are examined as the origin of the Li(+) ion and the solvent molecule, respectively. This salt is widely utilized as the electrolyte in the lithium ion secondary battery. The binding of EC to the Li(+) moiety of the Li(+)BF4(-) salt is exothermic, and the binding energies at the CAM-B3LYP/6-311++G(d,p) level for n=1, 2, 3, and 4, where n is the number of EC molecules binding to the Li(+) ion, (EC)n(Li(+)BF4(-)), are calculated to be 91.5, 89.8, 87.2, and 84.0 kcal mol(-1) (per EC molecule), respectively. The intermolecular distances between Li(+) and the F atom of BF4(-) are elongated: 1.773 Å (n=0), 1.820 Å (n=1), 1.974 Å (n=2), 1.942 Å (n=3), and 4.156 Å (n=4). The atomic bond populations between Li(+) and the F atom for n=0, 1, 2, 3, and 4 are 0.202, 0.186, 0.150, 0.038, and 0.0, respectively. These results indicate that the interaction of Li(+) with BF4(-) becomes weaker as the number of EC molecules is increased. The direct AIMD calculation for n=4 shows that EC reacts spontaneously with (EC)3(Li(+)BF4(-)) and the Li(+) ion is stripped from the salt. The following substitution reaction takes place: EC+(EC)3(Li(+)BF4(-))→(EC)4Li(+)-(BF4(-)). The reaction mechanism is discussed on the basis of the theoretical results. PMID:24616076
NASA Astrophysics Data System (ADS)
Arab, Fahima; Sahraoui, F. Ali; Haddadi, Khelifa; Bouhemadou, Abdelmadjid; Louail, Layachi
2016-05-01
Structural stability and mechanical and thermodynamic properties of the orthorhombic and trigonal MgSiN2 polymorphs (or-MgSiN2 and tr-MgSiN2) were investigated through density functional theory and quasi-harmonic Debye model (QHDM). Our calculations show that or-MgSiN2 is energetically the stable polymorph at low pressure, in agreement with previous experimental and theoretical study. Under pressure, a crystallographic transition from the orthorhombic structure to the trigonal one occurs around 25, 17.45 and 19.05 GPa as obtained from the generalized gradient approximation of Perdew-Wang (GGA-PW91), the generalized gradient approximation parameterized recently by Perdew et al (GGA-PBEsol) and the local density approximation developed by Ceperley and Alder and parameterized by Perdew and Zunger (LDA-CAPZ), respectively. Single-crystalline and polycrystalline elastic constants and related properties, namely Vickers hardness, acoustic Grüneisen parameter, minimum thermal conductivity, isotropic sound velocities and Debye temperature, were numerically estimated for both or-MgSiN2 and tr-MgSiN2. We have showed that the hardness of tr-MgSiN2 is comparable to that of the harder materials like c-BN and B6O. Temperature and pressure dependencies of volume, bulk modulus, thermal expansion, Grüneisen parameter, heat capacities and Debye temperature were investigated using QHDM.
Comas-Vives, Aleix
2016-03-14
In this contribution, realistic amorphous SiO2 models of 2.1 × 2.1 nm with silanol densities ranging 1.1-7.2 OH per nm(2) are obtained by means of ab initio calculations via the dehydroxylation of a fully hydroxylated silica surface. The dehydroxyation process is considered to take place via direct condensation of adjacent silanol groups and silica migration steps. The latter reconstructions are needed in order to obtain highly dehydroxylated silica surfaces with favorable energetics and without the formation of defects. The obtained surface phase diagram of different silica models as a function of temperature and PH2O is able to correctly describe the silanol density under different conditions, and the IR spectroscopic signatures of the silanols are in qualitative agreement with the experiment. The amorphous silica models presented here have a high degree of heterogeneity as found from the big variability obtained in the energetics of the dehydroxylation steps. It was also found that the resulting average Si-O distance of the newly formed siloxane bridges serves as a descriptor of the strain introduced in the silica surface. All these factors can be crucial in order to simulate the activity of catalysts grafted onto silica with different silanol densities, especially the one containing ca. 1 OH per nm(2), which can serve as a model for the SiO2 surface pretreated under high vacuum and at 700 °C. PMID:26898649
Ab initio investigation of helium in Y2Ti2O7: Mobility and effects on mechanical properties
NASA Astrophysics Data System (ADS)
Danielson, T.; Tea, E.; Hin, C.
2016-08-01
Oxide nanoclusters (NCs) in nanostructured ferritic alloys (NFAs) are known to be efficient trapping sites for the transmutation product helium. In this study, the migration barriers and potential energy surfaces of helium in Y2Ti2O7 are presented to explain the mobility of helium through oxide NCs and shed light on the accumulation of helium and the trapping mechanisms of the oxides. A complex tunnel-shaped potential energy surface is identified and gives rise to relatively large migration barriers. Subsequently, the effect of helium accumulation on the mechanical properties of Y2Ti2O7 oxide nanoclusters is investigated and it is shown that the mechanical properties of the oxide do not significantly degrade as helium accumulates.
NASA Astrophysics Data System (ADS)
Passler, Peter P.; Rode, Bernd M.
2015-10-01
Hydration properties of trivalent thulium and ytterbium ions in aqueous solution are investigated via quantum mechanical charge field molecular dynamics (QMCF-MD) simulations. The QMCF-MD formalism is a special type of QM/MM simulation, where the chemically most relevant part of the system - in this case the ion with its first and second hydration shells - is treated by quantum mechanics. The mean ionsbnd O distances and the average coordination numbers of the first hydration shells are compared with experimental EXAFS data. Mean ligand residence times, vibrational frequencies and force constants were evaluated to characterise the dynamics of the systems.
Baroni, Simone; Navratil, Petr; Quaglioni, Sofia
2013-03-26
In this study, we introduce a unified approach to nuclear bound and continuum states based on the coupling of the no-core shell model (NCSM), a bound-state technique, with the no-core shell model/resonating group method (NCSM/RGM), a nuclear scattering technique. This new ab initio method, no-core shell model with continuum (NCSMC), leads to convergence properties superior to either NCSM or NCSM/RGM while providing a balanced approach to different classes of states. In the NCSMC, the ansatz for the many-nucleon wave function includes (i) a square-integrable A-nucleon component expanded in a complete harmonic oscillator basis and (ii) a binary-cluster component with asymptoticmore » boundary conditions that can properly describe weakly bound states, resonances, and scattering. The Schrödinger equation is transformed into a system of coupled-channel integral-differential equations that we solve using a modified microscopic R-matrix formalism within a Lagrange mesh basis. We demonstrate the usefulness of the approach by investigating the unbound 7He nucleus.« less
Ab initio study on the electronic and mechanical properties of ReB and ReC
Zhao Erjun; Wang Jinping; Meng Jian; Wu Zhijian
2009-04-15
The structural, electronic, and mechanical properties of ReB and ReC have been studied by use of the density functional theory. For each compound, six structures are considered, i.e., hexagonal WC, NiAs, wurtzite, cubic NaCl, CsCl, and zinc-blende type structures. The results indicate that for ReB and ReC, WC type structure is energetically the most stable among the considered structures, followed by NiAs type structure. ReB-WC (i.e., ReB in WC type structure) and ReB-NiAs are both thermodynamically and mechanically stable. ReC-WC and ReC-NiAs are mechanically stable and becomes thermodynamically stable above 35 and 55 GPa, respectively. The estimated hardness from shear modulus is 34 GPa for ReB-WC, 28 GPa for ReB-NiAs, 35 GPa for ReC-WC and 37 GPa for ReC-NiAs, indicating that they are potential candidates to be ultra-incompressible and hard materials. - Graphical Abstract: The enthalpy-pressure diagram for ReC from the selected structures. ReC-WC and ReC-NiAs become thermodynamically stable above 35 and 55 GPa, respectively. ReC-WC is the most stable phase.
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.
{sup 7}Be(p,{gamma}){sup 8}B S factor from ab initio no-core shell model wave functions
Navratil, P.; Bertulani, C.A.; Caurier, E.
2006-06-15
Nuclear structure of {sup 7}Be, {sup 8}B, and {sup 7,8}Li is studied within the ab initio no-core shell model (NCSM). Starting from high-precision nucleon-nucleon (NN) interactions, wave functions of {sup 7}Be and {sup 8}B bound states are obtained in basis spaces up to 10({Dirac_h}/2{pi}){omega} and used to calculate channel cluster form factors (overlap integrals) of the {sup 8}B ground state with {sup 7}Be+p. Due to the use of the harmonic oscillator (HO) basis, the overlap integrals have incorrect asymptotic properties. We fix this problem in two alternative ways. First, by a Woods-Saxon potential solution fit to the interior of the NCSM overlap integrals. Second, by a direct matching with the Whittaker function. The corrected overlap integrals are then used for the {sup 7}Be(p,{gamma}){sup 8}B S-factor calculation. We study the convergence of the S factor with respect to the NCSM HO frequency and the model space size. Our S factor results agree with recent direct measurement data. We also test the spectroscopic factors and the corrected overlap integrals from the NCSM in describing the momentum distributions in knockout reactions with {sup 8}B projectiles. A good agreement with the available experimental data is also found, attesting to the overall consistency of the calculations.
Westphal, Eduard; Pliego, Josefredo R. Jr.
2005-08-15
The solvation of the lithium and sodium ions in dimethyl sulfoxide solution was theoretically investigated using ab initio calculations coupled with the hybrid cluster-continuum model, a quasichemical theory of solvation. We have investigated clusters of ions with up to five dimethyl sulfoxide (DMSO) molecules, and the bulk solvent was described by a dielectric continuum model. Our results show that the lithium and sodium ions have four and five DMSO molecules into the first coordination shell, and the calculated solvation free energies are -135.5 and -108.6 kcal mol{sup -1}, respectively. These data suggest a solvation free energy value of -273.2 kcal mol{sup -1} for the proton in dimethyl sulfoxide solution, a value that is more negative than the present uncertain experimental value. This and previous studies on the solvation of ions in water solution indicate that the tetraphenylarsonium tetraphenylborate assumption is flawed and the absolute value of the free energy of transfer of ions from water to DMSO solution is higher than the present experimental values.
NASA Astrophysics Data System (ADS)
Zhang, Yi; Zhao, Yusheng; Chen, Changfeng
2013-04-01
Recently, a family of halogen-based Li-rich antiperovskites were synthesized [J. Am. Chem. Soc.JACSAT0002-786310.1021/ja305709z 134, 15042 (2012)], and the measured superionic conductivity makes these materials promising candidates as solid electrolytes for applications in Li-ion and Li-air batteries. This discovery raises several pressing issues on the fundamental physics concerning the thermodynamic and electrochemical stability of the synthesized materials and the mechanism of the observed superionic Li+ transport. Here, we study the reported antiperovskites Li3OCl, Li3OBr, and their mixed compounds using first-principles density functional theory and molecular dynamics simulations. Our calculations show that these materials are thermodynamically metastable. Their large electronic band gaps and chemical stability against electrodes suggest the excellent electrochemical performance, which bodes well for the use in potentially harsh working conditions in practical battery applications. The calculated low activation enthalpy for Li-ion migration well below the crystal melting temperature and superionic transport near the Li sublattice melting state explain the experimentally observed phenomena. Our study identifies mobile Li vacancies and anion disorder as the primary driving mechanisms for superionic Li+ conductivity in the antiperovskites. This work unveils essential working principles of the Li-rich antiperovskites, which are crucial to further exploration, development, and application of these and other charge-inverted materials with tailored properties.
NASA Astrophysics Data System (ADS)
Prasetyo, Niko; Armunanto, Ria
2016-05-01
Structures and dynamics of Ag+ in 18.6% aqueous ammonia have been studied using Quantum Mechanical Charge Field Molecular Dynamics (QMCF-MD) simulation at the Hartree-Fock (HF) level theory employing LANL2DZ ECP basis set for Ag+ and Dunning DZP for solvent molecules. Structural properties are in excellent agreement with previous QM/MM and experiments studies. [Ag(NH3)2(H2O)3]+ was found as dominant species during simulation time. For 20 ps of simulation time, a labile first solvation shell was observed with both fast ammonia and water ligands exchanges. QMCF-MD framework describes first solvation shell more labile than conventional QM/MM MD simulation.
NASA Astrophysics Data System (ADS)
Lu, Xiu-hui; Han, Jun-feng; Yu, Hai-bin; Lian, Zhen-xia
2010-10-01
The mechanism of the cycloaddition reaction of forming a silicic bis-heterocyclic compound between singlet dimethylmethylenesilylene (Me2C=Si:) and ethene has been investigated with the CCSD(T)//MP2/6-31G* method. From the potential energy profile, it can be predicted that, this reaction has one dominant channel. The presented rule of this dominant channel: the 3p unoccupied orbital of Si in dimethylmethylenesilylene and the π orbital of ethene forming the π→p donor-acceptor bond, resulting in the formation of three-membered ring intermediate (INT1); INT1 then isomerizes to a four-membered ring silylene (P2), which is driven by ring-enlargement effect; due to sp3 hybridization of Si atom in P2, P2 further combines with ethene to form a silicic bis-heterocyclic compound.
Dawson, Matthew L.; Varner, Mychel E.; Perraud, Véronique; Ezell, Michael J.; Gerber, R. Benny; Finlayson-Pitts, Barbara J.
2012-01-01
Airborne particles affect human health and significantly influence visibility and climate. A major fraction of these particles result from the reactions of gaseous precursors to generate low-volatility products such as sulfuric acid and high-molecular weight organics that nucleate to form new particles. Ammonia and, more recently, amines, both of which are ubiquitous in the environment, have also been recognized as important contributors. However, accurately predicting new particle formation in both laboratory systems and in air has been problematic. During the oxidation of organosulfur compounds, gas-phase methanesulfonic acid is formed simultaneously with sulfuric acid, and both are found in particles in coastal regions as well as inland. We show here that: (i) Amines form particles on reaction with methanesulfonic acid, (ii) water vapor is required, and (iii) particle formation can be quantitatively reproduced by a semiempirical kinetics model supported by insights from quantum chemical calculations of likely intermediate clusters. Such an approach may be more broadly applicable in models of outdoor, indoor, and industrial settings where particles are formed, and where accurate modeling is essential for predicting their impact on health, visibility, and climate. PMID:23090988
Dawson, Matthew L; Varner, Mychel E; Perraud, Véronique; Ezell, Michael J; Gerber, R Benny; Finlayson-Pitts, Barbara J
2012-11-13
Airborne particles affect human health and significantly influence visibility and climate. A major fraction of these particles result from the reactions of gaseous precursors to generate low-volatility products such as sulfuric acid and high-molecular weight organics that nucleate to form new particles. Ammonia and, more recently, amines, both of which are ubiquitous in the environment, have also been recognized as important contributors. However, accurately predicting new particle formation in both laboratory systems and in air has been problematic. During the oxidation of organosulfur compounds, gas-phase methanesulfonic acid is formed simultaneously with sulfuric acid, and both are found in particles in coastal regions as well as inland. We show here that: (i) Amines form particles on reaction with methanesulfonic acid, (ii) water vapor is required, and (iii) particle formation can be quantitatively reproduced by a semiempirical kinetics model supported by insights from quantum chemical calculations of likely intermediate clusters. Such an approach may be more broadly applicable in models of outdoor, indoor, and industrial settings where particles are formed, and where accurate modeling is essential for predicting their impact on health, visibility, and climate. PMID:23090988
Arismendi-Arrieta, Daniel J; Riera, Marc; Bajaj, Pushp; Prosmiti, Rita; Paesani, Francesco
2016-03-01
New potential energy functions (i-TTM) describing the interactions between halide ions and water molecules are reported. The i-TTM potentials are derived from fits to electronic structure data and include an explicit treatment of two-body repulsion, electrostatics, and dispersion energy. Many-body effects are represented through classical polarization within an extended Thole-type model. By construction, the i-TTM potentials are compatible with the flexible and fully ab initio MB-pol potential, which has recently been shown to accurately predict the properties of water from the gas to the condensed phase. The accuracy of the i-TTM potentials is assessed through extensive comparisons with CCSD(T)-F12, DF-MP2, and DFT data as well as with results obtained with common polarizable force fields for X(-)(H2O)n clusters with X(-) = F(-), Cl(-), Br(-), and I(-), and n = 1-8. By construction, the new i-TTM potentials will enable direct simulations of vibrational spectra of halide-water systems from clusters to bulk and interfaces. PMID:26560189
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
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.
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.
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
NASA Astrophysics Data System (ADS)
Benlamari, S.; Amara Korba, S.; Lakel, S.; Meradji, H.; Ghemid, S.; El Haj Hassan, F.
2016-01-01
The structural, elastic, thermal and electronic properties of perovskite hydrides SrLiH3 and SrPdH3 have been investigated using the all-electron full-potential linear augmented plane wave (FP-LAPW) method based on the density functional theory (DFT). For the exchange-correlation potential, local-density approximation (LDA) and generalized gradient approximation (GGA) have been used to calculate theoretical lattice parameters, bulk modulus, and its pressure derivative. The present results are in good agreement with available theoretical and experimental data. The three independent elastic constants (C11, C12 and C44) are also reported. From electronic band structure and density of states (DOSs), it is found that SrLiH3 is an insulator characterized by an indirect gap of 3.48 eV, while SrPdH3 is metallic with a calculated DOSs at Fermi energy of 0.745 states/eV-unit cell. Poisson’s ratio (σ), Young’s modulus (E), shear modulus (G), anisotropy factor (A), average sound velocities (vm) and density (ρ) of these compounds are also estimated for the first time. The Debye temperature is deduced from the average sound velocity. Variation of elastic constants and bulk modulus of these compounds as a function of pressure is also reported. Pressure and thermal effects on some macroscopic properties are predicted using the quasi-harmonic Debye model.
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
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
NASA Astrophysics Data System (ADS)
Li, Xiaojie; Schönecker, Stephan; Li, Ruihuan; Li, Xiaoqing; Wang, Yuanyuan; Zhao, Jijun; Johansson, Börje; Vitos, Levente
2016-07-01
To examine the effect of neutron transmutation on tungsten as the first wall material of fusion reactors, the elastic properties of W1‑x‑y Re x Os y (0 ⩽ x, y ⩽ 6%) random alloys in body centered cubic (bcc) structure are investigated systematically using the all-electron exact muffin-tin orbitals (EMTO) method in combination with the coherent-potential approximation (CPA). The calculated lattice constant and elastic properties of pure W are consistent with available experiments. Both Os and Re additions reduce the lattice constant and increase the bulk modulus of W, with Os having the stronger effect. The polycrystalline shear modulus, Young’s modulus and the Debye temperature increase (decrease) with the addition of Re (Os). Except for C 11, the other elastic parameters including C 12, C 44, Cauchy pressure, Poisson ratio, B/G, increase as a function of Re and Os concentration. The variations of the latter three parameters and the trend in the ratio of cleavage energy to shear modulus for the most dominant slip system indicate that the ductility of the alloy enhances with increasing Re and Os content. The calculated elastic anisotropy of bcc W slightly increases with the concentration of both alloying elements. The estimated melting temperatures of the W–Re–Os alloy suggest that Re or Os addition will reduce the melting temperature of pure W solid. The classical Labusch–Nabarro model for solid-solution hardening predicts larger strengthening effects in W1‑y Os y than in W1‑x Re x . A strong correlation between C‧ and the fcc–bcc structural energy difference for W1‑x‑y Re x Os y is revealed demonstrating that canonical band structure dictates the alloying effect on C‧. The structural energy difference is exploited to estimate the alloying effect on the ideal tensile strength in the [0 0 1] direction.
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.
Plotnikov, Nikolay V.; Prasad, B. Ram; Chakrabarty, Suman; Chu, Zhen T.; Warshel, Arieh
2013-01-01
Understanding the nature of the free energy surfaces for phosphate hydrolysis is a prerequisite for understanding the corresponding key chemical reactions in biology. Here the challenge has been to move to careful ab initio QM/MM (QM(ai)/MM) free energy calculations, where obtaining converging results is very demanding and computationally expensive. This work describes such calculations, focusing on the free energy surface for the hydrolysis of phosphate monoesters, paying a special attention to the comparison between the one water (1W) and two water (2W) paths for the proton transfer (PT) step. This issue has been explored before by energy minimization with implicit solvent models and by non-systematic QM/MM energy minimization, as well as by non-systematic free energy mapping. However, no study has provided the needed reliable 2D (3D) surfaces which are necessary for reaching concrete conclusions. Our study generated in a systematic way the 2D (3D) free energy maps for several relevant systems, comparing the results of QM(ai)/MM and QM(ai)/implicit solvent surfaces, and provides an advanced description of the relevant energetics. It is found that the 1W path for the hydrolysis of methyl diphosphate (MDP) trianion is 6–9 kcal/mol higher than the 2W path. This difference becomes slightly larger in the presence of Mg2+ ion, since this ion reduces the pKa of the conjugated acid form of the phosphate oxygen that accepts the proton. Interestingly, the BLYP approach (which has been used extensively in some studies) gives much smaller difference between the 1W and 2W activation barriers. At any rate, it is worth to point out that the 2W transition state for the PT is not much higher that the common plateau that serves as the starting point of both the 1W and 2W PT paths. Thus, the calculated catalytic effects of proteins based on the 2W PT mechanistic models are not expected to be different from the catalytic effects predicted using the 1W PT mechanistic models
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
Kemege, Kyle E.; Hickey, John M.; Lovell, Scott; Battaile, Kevin P.; Zhang, Yang; Hefty, P. Scott
2012-02-13
Chlamydia trachomatis is a medically important pathogen that encodes a relatively high percentage of proteins with unknown function. The three-dimensional structure of a protein can be very informative regarding the protein's functional characteristics; however, determining protein structures experimentally can be very challenging. Computational methods that model protein structures with sufficient accuracy to facilitate functional studies have had notable successes. To evaluate the accuracy and potential impact of computational protein structure modeling of hypothetical proteins encoded by Chlamydia, a successful computational method termed I-TASSER was utilized to model the three-dimensional structure of a hypothetical protein encoded by open reading frame (ORF) CT296. CT296 has been reported to exhibit functional properties of a divalent cation transcription repressor (DcrA), with similarity to the Escherichia coli iron-responsive transcriptional repressor, Fur. Unexpectedly, the I-TASSER model of CT296 exhibited no structural similarity to any DNA-interacting proteins or motifs. To validate the I-TASSER-generated model, the structure of CT296 was solved experimentally using X-ray crystallography. Impressively, the ab initio I-TASSER-generated model closely matched (2.72-{angstrom} C{alpha} root mean square deviation [RMSD]) the high-resolution (1.8-{angstrom}) crystal structure of CT296. Modeled and experimentally determined structures of CT296 share structural characteristics of non-heme Fe(II) 2-oxoglutarate-dependent enzymes, although key enzymatic residues are not conserved, suggesting a unique biochemical process is likely associated with CT296 function. Additionally, functional analyses did not support prior reports that CT296 has properties shared with divalent cation repressors such as Fur.
NASA Astrophysics Data System (ADS)
Ulian, Gianfranco; Tosoni, Sergio; Valdrè, Giovanni
2013-11-01
The quantum chemical characterization of solid state systems is conducted with many different approaches, among which the adoption of periodic boundary conditions to deal with three-dimensional infinite condensed systems. This method, coupled to the Density Functional Theory (DFT), has been proved successful in simulating a huge variety of solids. Only in relatively recent years this ab initio quantum-mechanic approach has been used for the investigation of layer silicate structures and minerals. In the present work, a systematic comparison of different DFT functionals (GGA-PBEsol and hybrid B3LYP) and basis sets (plane waves and all-electron Gaussian-type orbitals) on the geometry, energy, and phonon properties of a model layer silicate, talc [Mg3Si4O10(OH)2], is presented. Long range dispersion is taken into account by DFT+D method. Results are in agreement with experimental data reported in literature, with minimal deviation given by the GTO/B3LYP-D* method regarding both axial lattice parameters and interaction energy and by PW/PBE-D for the unit-cell volume and angular values. All the considered methods adequately describe the experimental talc infrared spectrum.
Ulian, Gianfranco; Tosoni, Sergio; Valdrè, Giovanni
2013-11-28
The quantum chemical characterization of solid state systems is conducted with many different approaches, among which the adoption of periodic boundary conditions to deal with three-dimensional infinite condensed systems. This method, coupled to the Density Functional Theory (DFT), has been proved successful in simulating a huge variety of solids. Only in relatively recent years this ab initio quantum-mechanic approach has been used for the investigation of layer silicate structures and minerals. In the present work, a systematic comparison of different DFT functionals (GGA-PBEsol and hybrid B3LYP) and basis sets (plane waves and all-electron Gaussian-type orbitals) on the geometry, energy, and phonon properties of a model layer silicate, talc [Mg3Si4O10(OH)2], is presented. Long range dispersion is taken into account by DFT+D method. Results are in agreement with experimental data reported in literature, with minimal deviation given by the GTO∕B3LYP-D* method regarding both axial lattice parameters and interaction energy and by PW/PBE-D for the unit-cell volume and angular values. All the considered methods adequately describe the experimental talc infrared spectrum. PMID:24289338
Ulian, Gianfranco; Valdrè, Giovanni; Tosoni, Sergio
2013-11-28
The quantum chemical characterization of solid state systems is conducted with many different approaches, among which the adoption of periodic boundary conditions to deal with three-dimensional infinite condensed systems. This method, coupled to the Density Functional Theory (DFT), has been proved successful in simulating a huge variety of solids. Only in relatively recent years this ab initio quantum-mechanic approach has been used for the investigation of layer silicate structures and minerals. In the present work, a systematic comparison of different DFT functionals (GGA-PBEsol and hybrid B3LYP) and basis sets (plane waves and all-electron Gaussian-type orbitals) on the geometry, energy, and phonon properties of a model layer silicate, talc [Mg{sub 3}Si{sub 4}O{sub 10}(OH){sub 2}], is presented. Long range dispersion is taken into account by DFT+D method. Results are in agreement with experimental data reported in literature, with minimal deviation given by the GTO/B3LYP-D* method regarding both axial lattice parameters and interaction energy and by PW/PBE-D for the unit-cell volume and angular values. All the considered methods adequately describe the experimental talc infrared spectrum.