Ab Initio Many-Body Calculations Of Nucleon-Nucleus Scattering
Quaglioni, S; Navratil, P
2008-12-17
We develop a new ab initio many-body approach capable of describing simultaneously both bound and scattering states in light nuclei, by combining the resonating-group method with the use of realistic interactions, and a microscopic and consistent description of the nucleon clusters. This approach preserves translational symmetry and Pauli principle. We outline technical details and present phase shift results for neutron scattering on {sup 3}H, {sup 4}He and {sup 10}Be and proton scattering on {sup 3,4}He, using realistic nucleon-nucleon (NN) potentials. Our A = 4 scattering results are compared to earlier ab initio calculations. We find that the CD-Bonn NN potential in particular provides an excellent description of nucleon-{sup 4}He S-wave phase shifts. We demonstrate that a proper treatment of the coupling to the n-{sup 10}Be continuum is successful in explaining the parity-inverted ground state in {sup 11}Be.
Ab initio many-body calculations of light nuclei neutron and proton scattering
NASA Astrophysics Data System (ADS)
Quaglioni, Sofia
2008-10-01
One of the greatest challenges of nuclear physics today is the development of a quantitative microscopic theory of low-energy reactions on light nuclei. At the same time, technical progress on the theoretical front is urgent to match the major experimental advances in the study of exotic nuclei at the radioactive beam facilities. We build a new ab initio many-body approachootnotetextS. Quaglioni and P. Navratil, arXiv:0804.1560. capable of describing simultaneously both bound and scattering states in light nuclei, by combining the resonating-group methodootnotetextY. C. Tang et al., Phys. Rep. 47, 167 (1978); K. Langanke and H. Friedrich, Advances in Nuclear Physics, chapter 4., Plenum, New York, 1987. with the ab initio no-core shell model.ootnotetextP. Navratil, J. P. Vary, and B. R. Barrett, Phys. Rev. Lett. 84, 5728 (2000); Phys. Rev. C 62, 054311 (2000).. In this way, we complement a microscopic-cluster technique with the use of realistic interactions, and a microscopic and consistent description of the nucleon clusters, while preserving Pauli principle and translational symmetry. I will present results for neutron and proton scattering on light nuclei, including n- and p-^4He phase shifts, and low-lying states of one-neutron halo p-shell nuclei, obtained using realistic nucleon-nucleon potentials. In particular, I will address the parity inversion of the ^11Be ground state.
Ab initio many-body calculations of nucleon scattering on ^16O
NASA Astrophysics Data System (ADS)
Navratil, Petr; Quaglioni, Sofia; Roth, Robert
2008-10-01
We develop a new ab initio many-body approachootnotetextS. Quaglioni and P. Navratil, arXiv:0804.1560. capable of describing simultaneously both bound and scattering states in light nuclei, by combining the resonating-group methodootnotetextY. C. Tang et al., Phys. Rep. 47, 167 (1978); K. Langanke and H. Friedrich, Advances in Nuclear Physics, Plenum, New York, 1987. with the ab initio no-core shell model (NCSM).ootnotetextP. Navratil, J. P. Vary, and B. R. Barrett, Phys. Rev. Lett. 84, 5728 (2000); Phys. Rev. C 62, 054311 (2000). In this way, we complement a microscopic-cluster technique with the use of realistic interactions, and a microscopic and consistent description of the nucleon clusters, while preserving Pauli principle and translational symmetry. We will present results for low-energy nucleon scattering on ^16O and for A=17 bound states obtained using realistic nucleon-nucleon potentials. The ^16O wave functions are calculated within the importance-truncated NCSMootnotetextR. Roth and P. Navratil, Phys. Rev. Lett. 99, 092501 (2007). that allows the use of model spaces up to 18φ and ultimately enables to reach convergence of phase-shifts and other observables. Prepared by LLNL under Contract DE-AC52-07NA27344. Support from the U.S. DOE/SC/NP (Work Proposal No. SCW0498), and from the U. S. Department of Energy Grant DE-FC02-07ER41457 is acknowledged.
An AB Initio Relativistic Many-Body Calculation of Mercury 6
NASA Astrophysics Data System (ADS)
Cai, Ziyong
1990-01-01
In this work, the energy levels and widths of mercury 6p^2 resonance states are investigated, completing the second branch of our whole ab initio relativistic many-body calculation procedure. ^1 Relativistic zeroth-order wavefunctions are obtained by solving the multi-configurational Dirac -Fock (MCDF) equations.^2 The Breit effects are included. Relativistic configuration interaction (RCI) calculations are then carried out for the localized portion of the wavefunction^3 for which we included single and double excitations to account for both valence-shell correlations and the most important parts of the core-valence correlations, 5d-valence correlations. A total of 4 sets of virtual basis functions are used in the RCI calculations. Two virtual sets are used for valence shell correlations, both from vs_{1/2 } to vg_{9/2}. Another two sets are used for core-valence correlations, each ranging from vs_{1/2} to vf _{7/2}. To evaluate the continuum effects, relativistic continuum wavefunctions are calculated by using Perger's program^4 in the presence of a frozen core generated from Grant's program. ^5 Wavefunctions of the open channel discrete states are generated using Desclaux's program. ^2 The interactions between the resonances and the open channels are then calculated by using the method of the configuration interaction in continuum (CIC), ^6 for which a program, PVINT, ^7 is implemented to give both energy shifts and widths for the resonance states caused by the open channels. References. (1) D. R. Beck and Z. Cai, Phys. Rev. A 41, 301 (1990). (2) J. P. Desclaux, Comput. Phys. Commun 9, 31 (1975). (3) D. R. Beck, Program RELCOR, 1987, unpublished. (4) W. F. Perger and V. S. Karighattam, Program to generate continua and R^{k } integrals, 1990, Submitted to Comp. Phys. Commun. (5) I. P. Grant, B. J. McKenzie, P. H. Norrington, D. F. Mayers and N. C. Pyper, Comput. Phys. Commun. 21, 207 (1980). (6) U. Fano, Phys. Rev. 124, 1866 (1961). (7) Z. Cai, Program PVINT, 1990
The many-body Wigner Monte Carlo method for time-dependent ab-initio quantum simulations
Sellier, J.M. Dimov, I.
2014-09-15
The aim of ab-initio approaches is the simulation of many-body quantum systems from the first principles of quantum mechanics. These methods are traditionally based on the many-body Schrödinger equation which represents an incredible mathematical challenge. In this paper, we introduce the many-body Wigner Monte Carlo method in the context of distinguishable particles and in the absence of spin-dependent effects. Despite these restrictions, the method has several advantages. First of all, the Wigner formalism is intuitive, as it is based on the concept of a quasi-distribution function. Secondly, the Monte Carlo numerical approach allows scalability on parallel machines that is practically unachievable by means of other techniques based on finite difference or finite element methods. Finally, this method allows time-dependent ab-initio simulations of strongly correlated quantum systems. In order to validate our many-body Wigner Monte Carlo method, as a case study we simulate a relatively simple system consisting of two particles in several different situations. We first start from two non-interacting free Gaussian wave packets. We, then, proceed with the inclusion of an external potential barrier, and we conclude by simulating two entangled (i.e. correlated) particles. The results show how, in the case of negligible spin-dependent effects, the many-body Wigner Monte Carlo method provides an efficient and reliable tool to study the time-dependent evolution of quantum systems composed of distinguishable particles.
Predictive Nuclear Many-Body Theory with Ab Initio Methods: A Brief Survey and A Look Ahead
NASA Astrophysics Data System (ADS)
Hergert, Heiko
2015-10-01
The reach of ab initio many-body techniques has increased tremendously in recent years, owing to new developments in many-body theory as well as advances in their numerical implementation. Coupled Cluster, Self-Consistent Green's Function, and In-Medium Similarity Renormalization Group (IM-SRG) calculations are routinely performed for isotopes in the A ~ 100 region. Moreover, these techniques have been extended to tackle open-shell nuclei, either directly or through the auxiliary step of deriving valence-space interactions for use with existing Shell Model technology. One of the most powerful aspects of ab initio methods is their capability to provide results for energies and other observables with systematic uncertainties. Together with new accurate nuclear forces (and operators) derived from Chiral Effective Field Theory, they provide a consistent framework--and a road map--for a predictive description of nuclei. This will have a critical impact on the search for the limits of nuclear existence, tests of fundamental symmetries (e.g., the search for neutrinoless double beta decay), our understanding of quenching and effective charges in phenomenological Shell Model calculations etc. Using the Multi-Reference IM-SRG as a representative example, I will survey the current capabilities of ab initio methods with an emphasis on uncertainty quantification, highlight successes in the description of ground-state properties and spectra, and preview upcoming developments like the construction of consistent transition operators.
Ab initio many-body calculations of nucleon-4He scattering with three-nucleon forces
Hupin, Guillaume; Langhammer, Joachim; Navratil, Petr; Quaglioni, Sofia; Calci, Angelo; Roth, Robert
2013-11-27
We extend the ab initio no-core shell model/resonating-group method to include three-nucleon (3N) interactions for the description of nucleon-nucleus collisions. We outline the formalism, give algebraic expressions for the 3N-force integration kernels, and discuss computational aspects of two alternative implementations. The extended theoretical framework is then applied to nucleon-4He elastic scattering using similarity-renormalization-group (SRG)-evolved nucleon-nucleon plus 3N potentials derived from chiral effective field theory. We analyze the convergence properties of the calculated phase shifts and explore their dependence upon the SRG evolution parameter. We include up to six excited states of the 4He target and find significant effects from themore » inclusion of the chiral 3N force, e.g., it enhances the spin-orbit splitting between the 3/2– and 1/2– resonances and leads to an improved agreement with the phase shifts obtained from an accurate R-matrix analysis of the five-nucleon experimental data. As a result, we find remarkably good agreement with measured differential cross sections at various energies below the d+3H threshold, while analyzing powers manifest larger deviations from experiment for certain energies and angles.« less
Ab initio nuclear many-body perturbation calculations in the Hartree-Fock basis
NASA Astrophysics Data System (ADS)
Hu, B. S.; Xu, F. R.; Sun, Z. H.; Vary, J. P.; Li, T.
2016-07-01
Starting from realistic nuclear forces, the chiral N3LO and JISP16, we have applied many-body perturbation theory (MBPT) to the structure of closed-shell nuclei, 4He and 16O. The two-body N3LO interaction is softened by a similarity renormalization group transformation while JISP16 is adopted without renormalization. The MBPT calculations are performed within the Hartree-Fock (HF) bases. The angular momentum coupled scheme is used, which can reduce the computational task. Corrections up to the third order in energy and up to the second order in radius are evaluated. Higher-order corrections in the HF basis are small relative to the leading-order perturbative result. Using the antisymmetrized Goldstone diagram expansions of the wave function, we directly correct the one-body density for the calculation of the radius, rather than calculate corrections to the occupation probabilities of single-particle orbits as found in other treatments. We compare our results with other methods where available and find good agreement. This supports the conclusion that our methods produce reasonably converged results with these interactions. We also compare our results with experimental data.
NASA Astrophysics Data System (ADS)
Kong, Bo; Zhang, Yachao
2016-07-01
The electronic structures of the cubic GdH3 are extensively investigated using the ab initio many-body GW calculations treating the Gd 4f electrons either in the core (4f-core) or in the valence states (4f-val). Different degrees of quasiparticle (QP) self-consistent calculations with the different starting points are used to correct the failures of the GGA/GGA + U/HSE03 calculations. In the 4f-core case, GGA + G0W0 calculations give a fundamental band gap of 1.72 eV, while GGA+ GW0 or GGA + GW calculations present a larger band gap. In the 4f-val case, the nonlocal exchange-correlation (xc) functional HSE03 can account much better for the strong localization of the 4f states than the semilocal or Hubbard U corrected xc functional in the Kohn-Sham equation. We show that the fundamental gap of the antiferromagnetic (AFM) or ferromagnetic (FM) GdH3 can be opened up by solving the QP equation with improved starting point of eigenvalues and wave functions given by HSE03. The HSE03 + G0W0 calculations present a fundamental band gap of 2.73 eV in the AFM configuration, and the results of the corresponding GW0 and GW calculations are 2.89 and 3.03 eV, respectively. In general, for the cubic structure, the fundamental gap from G0W0 calculations in the 4f-core case is the closest to the real result. By G0W0 calculations in the 4f-core case, we find that H or Gd defects can strongly affect the band structure, especially the H defects. We explain the mechanism in terms of the possible electron correlation on the hydrogen site. Under compression, the insulator-to-metal transition in the cubic GdH3 occurs around 40 GPa, which might be a satisfied prediction.
NASA Astrophysics Data System (ADS)
Jorissen, Kevin; Story, Shauna; Rehr, John
2014-03-01
Accurate calculations of x-ray absorption spectra (XAS) often require linking several materials science codes. To reduce the complexity and support the hardware requirements of such calculations, we have virtualized XAS modeling workflows using a Cloud-based approach, with interfacing and configuration of codes handled by developers, and virtual HPC resources allocated on demand. When coupled to user-friendly GUIs this puts powerful multi-code simulations in the hands of general users. For instance, FEFF users can improve XAS interpretation and analysis using accurate ab initio Debye-Waller factors and self energy from the ABINIT DFT code, rather than semi-empirical models. Additionally, such workflows allow robust automation of large-scale calculation sets such as the Materials Project where our approach could enable a theoretical spectroscopy database of many thousands of structures for systematic study of materials. Supported by NSF-1216716.
NASA Astrophysics Data System (ADS)
France-Lanord, Arthur; Soukiassian, Patrick; Glattli, Christian; Wimmer, Erich
2016-03-01
In an effort to extend the reach of current ab initio calculations to simulations requiring millions of configurations for complex systems such as heterostructures, we have parameterized the third-generation Charge Optimized Many-Body (COMB3) potential using solely ab initio total energies, forces, and stress tensors as input. The quality and the predictive power of the new forcefield are assessed by computing properties including the cohesive energy and density of SiO2 polymorphs, surface energies of alpha-quartz, and phonon densities of states of crystalline and amorphous phases of SiO2. Comparison with data from experiments, ab initio calculations, and molecular dynamics simulations using published forcefields including BKS (van Beest, Kramer, and van Santen), ReaxFF, and COMB2 demonstrates an overall improvement of the new parameterization. The computed temperature dependence of the thermal conductivity of crystalline alpha-quartz and the Kapitza resistance of the interface between crystalline Si(001) and amorphous silica is in excellent agreement with experiment, setting the stage for simulations of complex nanoscale heterostructures.
France-Lanord, Arthur; Soukiassian, Patrick; Glattli, Christian; Wimmer, Erich
2016-03-14
In an effort to extend the reach of current ab initio calculations to simulations requiring millions of configurations for complex systems such as heterostructures, we have parameterized the third-generation Charge Optimized Many-Body (COMB3) potential using solely ab initio total energies, forces, and stress tensors as input. The quality and the predictive power of the new forcefield are assessed by computing properties including the cohesive energy and density of SiO2 polymorphs, surface energies of alpha-quartz, and phonon densities of states of crystalline and amorphous phases of SiO2. Comparison with data from experiments, ab initio calculations, and molecular dynamics simulations using published forcefields including BKS (van Beest, Kramer, and van Santen), ReaxFF, and COMB2 demonstrates an overall improvement of the new parameterization. The computed temperature dependence of the thermal conductivity of crystalline alpha-quartz and the Kapitza resistance of the interface between crystalline Si(001) and amorphous silica is in excellent agreement with experiment, setting the stage for simulations of complex nanoscale heterostructures. PMID:26979702
Many-body ab initio diffusion quantum Monte Carlo applied to the strongly correlated oxide NiO
Mitra, Chandrima; Krogel, Jaron T.; Santana, Juan A.; Reboredo, Fernando A.
2015-10-28
We present a many-body diffusion quantum Monte Carlo (DMC) study of the bulk and defect properties of NiO. We find excellent agreement with experimental values, within 0.3%, 0.6%, and 3.5% for the lattice constant, cohesive energy, and bulk modulus, respectively. The quasiparticle bandgap was also computed, and the DMC result of 4.72 (0.17) eV compares well with the experimental value of 4.3 eV. Furthermore, DMC calculations of excited states at the L, Z, and the gamma point of the Brillouin zone reveal a flat upper valence band for NiO, in good agreement with Angle Resolved Photoemission Spectroscopy results. To study defect properties, we evaluated the formation energies of the neutral and charged vacancies of oxygen and nickel in NiO. A formation energy of 7.2 (0.15) eV was found for the oxygen vacancy under oxygen rich conditions. For the Ni vacancy, we obtained a formation energy of 3.2 (0.15) eV under Ni rich conditions. These results confirm that NiO occurs as a p-type material with the dominant intrinsic vacancy defect being Ni vacancy.
Ab initio many-body calculations of the (3)H(d,n)(4)He and (3)He(d,p)(4)He fusion reactions.
Navrátil, Petr; Quaglioni, Sofia
2012-01-27
We apply the ab initio no-core shell model combined with the resonating-group method approach to calculate the cross sections of the (3)H(d,n)(4)He and (3)He(d,p)(4)He fusion reactions. These are important reactions for the big bang nucleosynthesis and the future of energy generation on Earth. Starting from a selected similarity-transformed chiral nucleon-nucleon interaction that accurately describes two-nucleon data, we performed many-body calculations that predict the S factor of both reactions. Virtual three-body breakup effects are obtained by including excited pseudostates of the deuteron in the calculation. Our results are in satisfactory agreement with experimental data and pave the way for microscopic investigations of polarization and electron-screening effects, of the (3)H(d,γn)(4)He bremsstrahlung and other reactions relevant to fusion research. PMID:22400830
Willow, Soohaeng Yoo; Salim, Michael A.; Kim, Kwang S.; Hirata, So
2015-01-01
A direct, simultaneous calculation of properties of a liquid using an ab initio electron-correlated theory has long been unthinkable. Here we present structural, dynamical, and response properties of liquid water calculated by ab initio molecular dynamics using the embedded-fragment spin-component-scaled second-order many-body perturbation method with the aug-cc-pVDZ basis set. This level of theory is chosen as it accurately and inexpensively reproduces the water dimer potential energy surface from the coupled-cluster singles, doubles, and noniterative triples with the aug-cc-pVQZ basis set, which is nearly exact. The calculated radial distribution function, self-diffusion coefficient, coordinate number, and dipole moment, as well as the infrared and Raman spectra are in excellent agreement with experimental results. The shapes and widths of the OH stretching bands in the infrared and Raman spectra and their isotropic-anisotropic Raman noncoincidence, which reflect the diverse local hydrogen-bond environment, are also reproduced computationally. The simulation also reveals intriguing dynamic features of the environment, which are difficult to probe experimentally, such as a surprisingly large fluctuation in the coordination number and the detailed mechanism by which the hydrogen donating water molecules move across the first and second shells, thereby causing this fluctuation. PMID:26400690
Willow, Soohaeng Yoo; Salim, Michael A; Kim, Kwang S; Hirata, So
2015-01-01
A direct, simultaneous calculation of properties of a liquid using an ab initio electron-correlated theory has long been unthinkable. Here we present structural, dynamical, and response properties of liquid water calculated by ab initio molecular dynamics using the embedded-fragment spin-component-scaled second-order many-body perturbation method with the aug-cc-pVDZ basis set. This level of theory is chosen as it accurately and inexpensively reproduces the water dimer potential energy surface from the coupled-cluster singles, doubles, and noniterative triples with the aug-cc-pVQZ basis set, which is nearly exact. The calculated radial distribution function, self-diffusion coefficient, coordinate number, and dipole moment, as well as the infrared and Raman spectra are in excellent agreement with experimental results. The shapes and widths of the OH stretching bands in the infrared and Raman spectra and their isotropic-anisotropic Raman noncoincidence, which reflect the diverse local hydrogen-bond environment, are also reproduced computationally. The simulation also reveals intriguing dynamic features of the environment, which are difficult to probe experimentally, such as a surprisingly large fluctuation in the coordination number and the detailed mechanism by which the hydrogen donating water molecules move across the first and second shells, thereby causing this fluctuation. PMID:26400690
NASA Astrophysics Data System (ADS)
Reyes-Lillo, Sebastian E.; Rangel, Tonatiuh; Bruneval, Fabien; Neaton, Jeffrey B.
2016-07-01
The Ruddlesden-Popper (RP) homologous series Srn +1TinO3 n +1 provides a useful template for the study and control of the effects of dimensionality and quantum confinement on the excited state properties of the complex oxide SrTiO3. We use ab initio many-body perturbation theory within the G W approximation and the Bethe-Salpeter equation approach to calculate quasiparticle energies and absorption spectra of Srn +1TinO3 n +1 for n =1 -5 and ∞ . Our computed direct and indirect optical gaps are in excellent agreement with spectroscopic measurements. The calculated optical spectra reproduce the main experimental features and reveal excitonic structure near the gap edge. We find that electron-hole interactions are important across the series, leading to significant exciton binding energies that increase for small n and reach a value of 330 meV for n =1 , a trend attributed to increased quantum confinement. We find that the lowest-energy singlet exciton of Sr2TiO4 (n =1 ) localizes in the two-dimensional plane defined by the TiO2 layer, and we explain the origin of its localization.
Ab initio many-body calculations of nucleon-^{4}He scattering with three-nucleon forces
Hupin, Guillaume; Langhammer, Joachim; Navratil, Petr; Quaglioni, Sofia; Calci, Angelo; Roth, Robert
2013-11-27
We extend the ab initio no-core shell model/resonating-group method to include three-nucleon (3N) interactions for the description of nucleon-nucleus collisions. We outline the formalism, give algebraic expressions for the 3N-force integration kernels, and discuss computational aspects of two alternative implementations. The extended theoretical framework is then applied to nucleon-^{4}He elastic scattering using similarity-renormalization-group (SRG)-evolved nucleon-nucleon plus 3N potentials derived from chiral effective field theory. We analyze the convergence properties of the calculated phase shifts and explore their dependence upon the SRG evolution parameter. We include up to six excited states of the ^{4}He target and find significant effects from the inclusion of the chiral 3N force, e.g., it enhances the spin-orbit splitting between the 3/2^{–} and 1/2^{–} resonances and leads to an improved agreement with the phase shifts obtained from an accurate R-matrix analysis of the five-nucleon experimental data. As a result, we find remarkably good agreement with measured differential cross sections at various energies below the d+^{3}H threshold, while analyzing powers manifest larger deviations from experiment for certain energies and angles.
Coccia, Emanuele; Varsano, Daniele; Guidoni, Leonardo
2016-01-01
In this letter, we report the singlet ground state structure of the full carotenoid peridinin by means of variational Monte Carlo (VMC) calculations. The VMC relaxed geometry has an average bond length alternation of 0.1165(10) Å, larger than the values obtained by DFT (PBE, B3LYP, and CAM-B3LYP) and shorter than that calculated at the Hartree–Fock (HF) level. TDDFT and EOM-CCSD calculations on a reduced peridinin model confirm the HOMO–LUMO major contribution of the Bu+-like (S2) bright excited state. Many Body Green’s Function Theory (MBGFT) calculations of the vertical excitation energy of the Bu+-like state for the VMC structure (VMC/MBGFT) provide an excitation energy of 2.62 eV, in agreement with experimental results in n-hexane (2.72 eV). The dependence of the excitation energy on the bond length alternation in the MBGFT and TDDFT calculations with different functionals is discussed. PMID:26580027
Misquitta, Alston J; Stone, Anthony J
2016-09-13
Creating accurate, analytic atom-atom potentials for small organic molecules from first principles can be a time-consuming and computationally intensive task, particularly if we also require them to include explicit polarization terms, which are essential in many systems. We describe how the CamCASP suite of programs can be used to generate such potentials using some of the most accurate electronic structure methods currently applicable. We derive the long-range terms from monomer properties and determine the short-range anisotropy parameters by a novel and robust method based on the iterated stockholder atom approach. Using these techniques, we develop distributed multipole models for the electrostatic, polarization, and dispersion interactions in the pyridine dimer and develop a series of many-body potentials for the pyridine system. Even the simplest of these potentials exhibits root mean square errors of only about 0.6 kJ mol(-1) for the low-energy pyridine dimers, significantly surpassing the best empirical potentials. Our best model is shown to support eight stable minima, four of which have not been reported before in the literature. Further, the functional form can be made systematically more elaborate so as to improve the accuracy without a significant increase in the human-time spent in their generation. We investigate the effects of anisotropy, rank of multipoles, and choice of polarizability and dispersion models. PMID:27467814
Liu, Jie; Herbert, John M
2016-02-01
We introduce a "pair-pair" approximation to the generalized many-body expansion (pp-GMBE) as an approximation to a traditional four-body expansion, the latter of which is accurate but quickly becomes numerically unstable and ultimately intractable as the number of "bodies" (fragments) increases. The pp-GMBE method achieves a good balance between accuracy and efficiency by defining significant fragment pairs and then fragment quartets. An efficient fragmentation scheme is introduced for proteins such that the largest subsystems contain about 60 atoms. Application of the pp-GMBE method to proteins with as many as 70 amino acids (1142 atoms) reveals that pp-GMBE energies are quite faithful to those obtained when the same level of density functional theory is applied to the entire macromolecule. When combined with embedding charges obtained from natural population analysis, the pp-GMBE approach affords absolute energies that differ by 1-3 kcal/mol from full supersystem results, but it yields conformational energy profiles that are practically indistinguishable from the supersystem calculation at the same level of theory. PMID:26730608
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
Collective rotation from ab initio theory
NASA Astrophysics Data System (ADS)
Caprio, Mark A.; Maris, Pieter; Vary, James P.
2015-10-01
The challenge of ab initio nuclear theory is to quantitatively predict the complex and highly-correlated behavior of the nuclear many-body system, starting from the underlying internucleon interactions. We may now seek to understand the wealth of nuclear collective phenomena through ab initio approaches. No-core configuration interaction (NCCI) calculations for p-shell nuclei give rise to rotational bands, as evidenced by rotational patterns for excitation energies, electromagnetic moments, and electromagnetic transitions. In this talk, the intrinsic structure of these bands is discussed, and the predicted rotational bands are compared to experiment. Supported by the US DOE under Award Nos. DE-FG02-95ER-40934, DESC0008485 (SciDAC/NUCLEI), and DE-FG02-87ER40371 and the US NSF under Award No. 0904782. Computational resources provided by NERSC (US DOE Contract No. DE-AC02-05CH11231).
AB INITIO 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.
Collective rotation from ab initio theory
NASA Astrophysics Data System (ADS)
Caprio, M. A.; Maris, P.; Vary, J. P.; Smith, R.
2015-08-01
Through ab initio approaches in nuclear theory, we may now seek to quantitatively understand the wealth of nuclear collective phenomena starting from the underlying internucleon interactions. No-core configuration interaction (NCCI) calculations for p-shell nuclei give rise to rotational bands, as evidenced by rotational patterns for excitation energies, electromagnetic moments and electromagnetic transitions. In this review, NCCI calculations of 7-9Be are used to illustrate and explore ab initio rotational structure, and the resulting predictions for rotational band properties are compared with experiment. We highlight the robustness of ab initio rotational predictions across different choices for the internucleon interaction.
Ab initio phonon limited transport
NASA Astrophysics Data System (ADS)
Verstraete, Matthieu
We revisit the thermoelectric (TE) transport properties of two champion materials, PbTe and SnSe, using fully first principles methods. In both cases the performance of the material is due to subtle combinations of structural effects, scattering, and phase space reduction. In PbTe anharmonic effects are completely opposite to the predicted quasiharmonic evolution of phonon frequencies and to frequently (and incorrectly) cited extrapolations of experiments. This stabilizes the material at high T, but also tends to enhance its thermal conductivity, in a non linear manner, above 600 Kelvin. This explains why PbTe is in practice limited to room temperature applications. SnSe has recently been shown to be the most efficient TE material in bulk form. This is mainly due to a strongly enhanced carrier concentration and electrical conductivity, after going through a phase transition from 600 to 800 K. We calculate the transport coefficients as well as the defect concentrations ab initio, showing excellent agreement with experiment, and elucidating the origin of the double phase transition as well as the new charge carriers. AH Romero, EKU Gross, MJ Verstraete, and O Hellman PRB 91, 214310 (2015) O. Hellman, IA Abrikosov, and SI Simak, PRB 84 180301 (2011)
Navráail, Petr; Quaglioni, Sofia
2012-01-24
We apply the ab initio no-core shell model combined with the resonating-group method approach to calculate the cross sections of the ^{3}H(d,n)^{4}He and ^{3}He(d,p)^{4}He fusion reactions. These are important reactions for the big bang nucleosynthesis and the future of energy generation on Earth. Starting from a selected similarity-transformed chiral nucleon-nucleon interaction that accurately describes two-nucleon data, we performed many-body calculations that predict the S factor of both reactions. Virtual three-body breakup effects are obtained by including excited pseudostates of the deuteron in the calculation. Lastly, our results are in satisfactory agreement with experimental data and pave the way for microscopic investigations of polarization and electron-screening effects, of the ^{3}H(d,γn)^{4}He bremsstrahlung and other reactions relevant to fusion research.
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.
Zabidi, Noriza Ahmad; Kassim, Hasan Abu; Shrivastava, Keshav N.
2008-05-20
Polonium is the only element with a simple cubic (sc) crystal structure. Atoms in solid polonium sit at the corners of a simple cubic unit cell and no where else. Polonium has a valence electron configuration 6s{sup 2}6p{sup 4} (Z = 84). The low temperature {alpha}-phase transforms into the rhombohedral (trigonal) {beta} structure at {approx}348 K. The sc {alpha}-Po unit cell constant is a = 3.345 A. The beta form of polonium ({beta}-Po) has the lattice parameters, a{sub R} = 3.359 A and a rhombohedral angle 98 deg. 13'. We have performed an ab initio electronic structure calculation by using the density functional theory. We have performed the calculation with and without spin-orbit (SO) coupling by using both the LDA and the GGA for the exchange-correlations. The k-points in a simple cubic BZ are determined by R (0.5, 0.5, 0.5), {gamma} (0, 0, 0), X (0.5, 0, 0), M (0.5, 0.5, 0) and {gamma} (0, 0, 0). Other directions of k-points are {gamma} (0, 0, 0), X (0.5, 0, 0), R (0.5, 0.5, 0.5) and {gamma} (0, 0, 0). The SO splittings of p states at the {gamma} point in the GGA+SO scheme for {alpha}-Po are 0.04 eV and 0.02 eV while for the {beta}-Po these are 0.03 eV and 0.97 eV. We have also calculated the vibrational spectra for the unit cells in both the structures. We find that exchanging of a Po atom by Pb atom produces several more bands and destabilizes the {beta} phase.
Ab initio calculations of reactions with light nuclei
NASA Astrophysics Data System (ADS)
Quaglioni, Sofia; Hupin, Guillaume; Calci, Angelo; Navrátil, Petr; Roth, Robert
2016-03-01
An ab initio (i.e., from first principles) theoretical framework capable of providing a unified description of the structure and low-energy reaction properties of light nuclei is desirable to further our understanding of the fundamental interactions among nucleons, and provide accurate predictions of crucial reaction rates for nuclear astrophysics, fusion-energy research, and other applications. In this contribution we review ab initio calculations for nucleon and deuterium scattering on light nuclei starting from chiral two- and three-body Hamiltonians, obtained within the framework of the ab initio no-core shell model with continuum. This is a unified approach to nuclear bound and scattering states, in which square-integrable energy eigenstates of the A-nucleon system are coupled to (A-a)+a target-plus-projectile wave functions in the spirit of the resonating group method to obtain an efficient description of the many-body nuclear dynamics both at short and medium distances and at long ranges.
Many-Body Models for Molecular Nanomagnets
NASA Astrophysics Data System (ADS)
Chiesa, A.; Carretta, S.; Santini, P.; Amoretti, G.; Pavarini, E.
2013-04-01
We present a flexible and effective ab initio scheme to build many-body models for molecular nanomagnets, and to calculate magnetic exchange couplings and zero-field splittings. It is based on using localized Foster-Boys orbitals as a one-electron basis. We apply this scheme to three paradigmatic systems, the antiferromagnetic rings Cr8 and Cr7Ni, and the single-molecule magnet Fe4. In all cases we identify the essential magnetic interactions and find excellent agreement with experiments.
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
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.
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
Ab initio infrared and Raman spectra
NASA Technical Reports Server (NTRS)
Fredkin, D. R.; White, S. R.; Wilson, K. R.; Komornicki, A.
1983-01-01
It is pointed out that with increased computer power and improved computational techniques, such as the gradients developed in recent years, it is becoming practical to compute spectra ab initio, from the fundamental constants of nature, for systems of increasing complexity. The present investigation has the objective to explore several possible ab initio approaches to spectra, giving particular attention to infrared and nonresonance Raman. Two approaches are discussed. The sequential approach, in which first the electronic part and then later the nuclear part of the Born-Oppenheimer approximation is solved, is appropriate for small systems. The simultaneous approach, in which the electronic and nuclear parts are solved at the same time, is more appropriate for many-atom systems. A review of the newer quantum gradient techniques is provided, and the infrared and Raman spectral band contours for the water molecule are computed.
NASA Astrophysics Data System (ADS)
Thiessen, P. A.; Treder, H.-J.
Jedes initium wird durch experimenta crucis zum eventus. Jedes theoretisch interpretierbare ex-eventu-Resultat führt auf ein neues Initium. Gerade dies ist die gemeinsame Aussage von Atomistik, Quantenmechanik und Relativitätstheorie.Translated AbstractAb initio vel ex eventu. IIEvery initium becomes an eventus by experimenta crucis. Every theoretically interpretable ex-eventu result leads to a new initium. Right this is the joint assertion of atomism, quantum mechanics, and relativity.
Ab initio Bogoliubov coupled cluster theory
NASA Astrophysics Data System (ADS)
Signoracci, Angelo; Hagen, Gaute; Duguet, Thomas
2014-09-01
Coupled cluster (CC) theory has become a standard method in nuclear theory for realistic ab initio calculations of medium mass nuclei, but remains limited by its requirement of a Slater determinant reference state which reasonably approximates the nuclear system of interest. Extensions of the method, such as equation-of-motion CC, permit the calculation of nuclei with one or two nucleons added or removed from a doubly magic core, yet still only a few dozen nuclei are accessible with modern computational restrictions. In order to extend the applicability of ab initio methods to open-shell systems, the superfluid nature of nuclei must be taken into account. By utilizing Bogoliubov algebra and employing spontaneous symmetry breaking with respect to particle number conservation, superfluid systems can be treated by a single reference state. An ab initio theory to include correlations on top of a Bogoliubov reference state has been developed in the guise of standard CC theory. The formalism and first results of this Bogoliubov coupled cluster theory will be presented to demonstrate the applicability of the method.
Ab initio tight-binding Hamiltonian for transition metal dichalcogenides
NASA Astrophysics Data System (ADS)
Fang, Shiang; Kuate Defo, Rodrick; Shirodkar, Sharmila N.; Lieu, Simon; Tritsaris, Georgios A.; Kaxiras, Efthimios
2015-11-01
We present an accurate ab initio tight-binding Hamiltonian for the transition metal dichalcogenides, MoS2, MoSe2, WS2, WSe2, with a minimal basis (the d orbitals for the metal atoms and p orbitals for the chalcogen atoms) based on a transformation of the Kohn-Sham density functional theory Hamiltonian to a basis of maximally localized Wannier functions. The truncated tight-binding Hamiltonian, with only on-site, first, and partial second neighbor interactions, including spin-orbit coupling, provides a simple physical picture and the symmetry of the main band-structure features. Interlayer interactions between adjacent layers are modeled by transferable hopping terms between the chalcogen p orbitals. The full-range tight-binding Hamiltonian can be reduced to hybrid-orbital k .p effective Hamiltonians near the band extrema that capture important low-energy excitations. These ab initio Hamiltonians can serve as the starting point for applications to interacting many-body physics including optical transitions and Berry curvature of bands, of which we give some examples.
Three-cluster dynamics within an ab initio framework
Quaglioni, Sofia; Romero-Redondo, Carolina; Navratil, Petr
2013-09-26
In this study, we introduce a fully antisymmetrized treatment of three-cluster dynamics within the ab initio framework of the no-core shell model/resonating-group method. Energy-independent nonlocal interactions among the three nuclear fragments are obtained from realistic nucleon-nucleon interactions and consistent ab initio many-body wave functions of the clusters. The three-cluster Schrödinger equation is solved with bound-state boundary conditions by means of the hyperspherical-harmonic method on a Lagrange mesh. We discuss the formalism in detail and give algebraic expressions for systems of two single nucleons plus a nucleus. Using a soft similarity-renormalization-group evolved chiral nucleon-nucleon potential, we apply the method to amore » 4He+n+n description of 6He and compare the results to experiment and to a six-body diagonalization of the Hamiltonian performed within the harmonic-oscillator expansions of the no-core shell model. Differences between the two calculations provide a measure of core (4He) polarization effects.« less
Ab initio study of hot electrons in GaAs.
Bernardi, Marco; Vigil-Fowler, Derek; Ong, Chin Shen; Neaton, Jeffrey B; Louie, Steven G
2015-04-28
Hot carrier dynamics critically impacts the performance of electronic, optoelectronic, photovoltaic, and plasmonic devices. Hot carriers lose energy over nanometer lengths and picosecond timescales and thus are challenging to study experimentally, whereas calculations of hot carrier dynamics are cumbersome and dominated by empirical approaches. In this work, we present ab initio calculations of hot electrons in gallium arsenide (GaAs) using density functional theory and many-body perturbation theory. Our computed electron-phonon relaxation times at the onset of the Γ, L, and X valleys are in excellent agreement with ultrafast optical experiments and show that the ultrafast (tens of femtoseconds) hot electron decay times observed experimentally arise from electron-phonon scattering. This result is an important advance to resolve a controversy on hot electron cooling in GaAs. We further find that, contrary to common notions, all optical and acoustic modes contribute substantially to electron-phonon scattering, with a dominant contribution from transverse acoustic modes. This work provides definitive microscopic insight into hot electrons in GaAs and enables accurate ab initio computation of hot carriers in advanced materials. PMID:25870287
Unified ab initio approaches to nuclear structure and reactions
NASA Astrophysics Data System (ADS)
Navrátil, Petr; Quaglioni, Sofia; Hupin, Guillaume; Romero-Redondo, Carolina; Calci, Angelo
2016-05-01
The description of nuclei starting from the constituent nucleons and the realistic interactions among them has been a long-standing goal in nuclear physics. In addition to the complex nature of the nuclear forces, with two-, three- and possibly higher many-nucleon components, one faces the quantum-mechanical many-nucleon problem governed by an interplay between bound and continuum states. In recent years, significant progress has been made in ab initio nuclear structure and reaction calculations based on input from QCD-employing Hamiltonians constructed within chiral effective field theory. After a brief overview of the field, we focus on ab initio many-body approaches—built upon the no-core shell model—that are capable of simultaneously describing both bound and scattering nuclear states, and present results for resonances in light nuclei, reactions important for astrophysics and fusion research. In particular, we review recent calculations of resonances in the 6He halo nucleus, of five- and six-nucleon scattering, and an investigation of the role of chiral three-nucleon interactions in the structure of 9Be. Further, we discuss applications to the 7Be {({{p}},γ )}8{{B}} radiative capture. Finally, we highlight our efforts to describe transfer reactions including the 3H{({{d}},{{n}})}4He fusion.
Ab initio no core full configuration approach for light nuclei
NASA Astrophysics Data System (ADS)
Kim, Youngman; Shin, Ik Jae; Maris, Pieter; Vary, James P.; Forssén, Christian; Rotureau, Jimmy
2014-07-01
Comprehensive understanding of the structure and reactions of light nuclei poses theoretical and computational challenges. Still, a number of ab initio approaches have been developed to calculate the properties of atomic nuclei using fundamental interactions among nucleons. Among them, we work with the ab initio no core full configuration (NCFC) method and ab initio no core Gamow Shell Model (GSM). We first review these approaches and present some recent results.
Ab initio no core full configuration approach for light nuclei
NASA Astrophysics Data System (ADS)
Kim, Youngman; Shin, Ik Jae; Maris, Pieter; Vary, James P.; Forssén, Christian; Rotureau, Jimmy
2015-10-01
Comprehensive understanding of the structure and reactions of light nuclei poses theoretical and computational challenges. Still, a number of ab initio approaches have been developed to calculate the properties of atomic nuclei using fundamental interactions among nucleons. Among them, we work with the ab initio no core full configuration (NCFC) method and ab initio no core Gamow Shell Model (GSM). We first review these approaches and present some recent results.
Ab initio melting curve of osmium
NASA Astrophysics Data System (ADS)
Burakovsky, L.; Burakovsky, N.; Preston, D. L.
2015-11-01
The melting curve of osmium up to a pressure P of 500 GPa is obtained from an extensive suite of ab initio quantum molecular dynamics (QMD) simulations using the Z method. The ab initio P =0 melting point of Os is 3370 ±75 K; this range encompasses all of the available data in the literature and corroborates the conclusion of J. W. Arblaster [Platinum Metals Rev. 49, 166 (2005)], 10.1595/147106705X70264 that the melting temperature of pure Os is 3400 ±50 K and that the 3300 K typically quoted in the literature is the melting point of impure Os. The T =0 equation of state (EOS) of Os and the P dependence of the optimized c /a ratio for the hexagonal unit cell, both to pressures ˜900 GPa, are obtained in the ab initio approach as validation of its use. Although excellent agreement with the available experimental data (P ≲80 GPa) is found, it is the third-order Birch-Murnaghan EOS with B0'=5 rather than the more widely accepted B0'=4 that describes the QMD data to higher pressures, in agreement with the more recent experimental EOS by Godwal et al. The theoretical melting curve of Os obtained earlier by Joshi et al. is shown to be inconsistent with our QMD results, and the possible reason for this discrepancy is suggested. Regularities in the melting curves of Os and five other third-row transition metals (Ta, W, Re, Pt, Au) could be used to estimate the currently unknown melting curves of Hf and Ir.
Ab initio quantum chemistry: Methodology and applications
Friesner, Richard A.
2005-01-01
This Perspective provides an overview of state-of-the-art ab initio quantum chemical methodology and applications. The methods that are discussed include coupled cluster theory, localized second-order Moller–Plesset perturbation theory, multireference perturbation approaches, and density functional theory. The accuracy of each approach for key chemical properties is summarized, and the computational performance is analyzed, emphasizing significant advances in algorithms and implementation over the past decade. Incorporation of a condensed-phase environment by means of mixed quantum mechanical/molecular mechanics or self-consistent reaction field techniques, is presented. A wide range of illustrative applications, focusing on materials science and biology, are discussed briefly. PMID:15870212
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.
Towards Accurate Ab Initio Predictions of the Spectrum of Methane
NASA Technical Reports Server (NTRS)
Schwenke, David W.; Kwak, Dochan (Technical Monitor)
2001-01-01
We have carried out extensive ab initio calculations of the electronic structure of methane, and these results are used to compute vibrational energy levels. We include basis set extrapolations, core-valence correlation, relativistic effects, and Born- Oppenheimer breakdown terms in our calculations. Our ab initio predictions of the lowest lying levels are superb.
Discovering chemistry with an ab initio nanoreactor
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S.; Martínez, Todd J.
2014-11-02
Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provides detailed physical insight. While theoretical and computational studies have generally focused on specific processes or mechanistic hypotheses, recent methodological and computational advances harken the advent of their principal role in discovery. Here we report the development and application of the ab initio nanoreactor – a highly accelerated, first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor we show new pathways for glycine synthesis frommore » primitive compounds proposed to exist on the early Earth, providing new insight into the classic Urey-Miller experiment. Ultimately, these results highlight the emergence of theoretical and computational chemistry as a tool for discovery in addition to its traditional role of interpreting experimental findings.« less
Ab initio molar volumes and Gaussian radii.
Parsons, Drew F; Ninham, Barry W
2009-02-12
Ab initio molar volumes are calculated and used to derive radii for ions and neutral molecules using a spatially diffuse model of the electron distribution with Gaussian spread. The Gaussian radii obtained can be used for computation of nonelectrostatic ion-ion dispersion forces that underlie Hofmeister specific ion effects. Equivalent hard-sphere radii are also derived, and these are in reasonable agreement with crystalline ionic radii. The Born electrostatic self-energy is derived for a Gaussian model of the electronic charge distribution. It is shown that the ionic volumes used in electrostatic calculations of strongly hydrated cosmotropic ions ought best to include the first hydration shell. Ionic volumes for weakly hydrated chaotropic metal cations should exclude electron overlap (in electrostatic calculations). Spherical radii are calculated as well as nonisotropic ellipsoidal radii for nonspherical ions, via their nonisotropic static polarizability tensors. PMID:19140766
Guiding ab initio calculations by alchemical derivatives.
to Baben, M; Achenbach, J O; von Lilienfeld, O A
2016-03-14
We assess the concept of alchemical transformations for predicting how a further and not-tested change in composition would change materials properties. This might help to guide ab initio calculations through multidimensional property-composition spaces. Equilibrium volumes, bulk moduli, and relative lattice stability of fcc and bcc 4d transition metals Zr, Nb, Mo, Tc, Ru, Rh, Pd, and Ag are calculated using density functional theory. Alchemical derivatives predict qualitative trends in lattice stability while equilibrium volumes and bulk moduli are predicted with less than 9% and 28% deviation, respectively. Predicted changes in equilibrium volume and bulk moduli for binary and ternary mixtures of Rh-Pd-Ag are in qualitative agreement even for predicted bulk modulus changes as large as +100% or -50%. Based on these results, it is suggested that alchemical transformations could be meaningful for enhanced sampling in the context of virtual high-throughput materials screening projects. PMID:26979677
Ab initio determination of light hadron masses.
Dürr, S; Fodor, Z; Frison, J; Hoelbling, C; Hoffmann, R; Katz, S D; Krieg, S; Kurth, T; Lellouch, L; Lippert, T; Szabo, K K; Vulvert, G
2008-11-21
More than 99% of the mass of the visible universe is made up of protons and neutrons. Both particles are much heavier than their quark and gluon constituents, and the Standard Model of particle physics should explain this difference. We present a full ab initio calculation of the masses of protons, neutrons, and other light hadrons, using lattice quantum chromodynamics. Pion masses down to 190 mega-electron volts are used to extrapolate to the physical point, with lattice sizes of approximately four times the inverse pion mass. Three lattice spacings are used for a continuum extrapolation. Our results completely agree with experimental observations and represent a quantitative confirmation of this aspect of the Standard Model with fully controlled uncertainties. PMID:19023076
Guiding ab initio calculations by alchemical derivatives
NASA Astrophysics Data System (ADS)
to Baben, M.; Achenbach, J. O.; von Lilienfeld, O. A.
2016-03-01
We assess the concept of alchemical transformations for predicting how a further and not-tested change in composition would change materials properties. This might help to guide ab initio calculations through multidimensional property-composition spaces. Equilibrium volumes, bulk moduli, and relative lattice stability of fcc and bcc 4d transition metals Zr, Nb, Mo, Tc, Ru, Rh, Pd, and Ag are calculated using density functional theory. Alchemical derivatives predict qualitative trends in lattice stability while equilibrium volumes and bulk moduli are predicted with less than 9% and 28% deviation, respectively. Predicted changes in equilibrium volume and bulk moduli for binary and ternary mixtures of Rh-Pd-Ag are in qualitative agreement even for predicted bulk modulus changes as large as +100% or -50%. Based on these results, it is suggested that alchemical transformations could be meaningful for enhanced sampling in the context of virtual high-throughput materials screening projects.
Discovering chemistry with an ab initio nanoreactor.
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S; Martínez, Todd J
2014-12-01
Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provide detailed physical insight. Although theoretical and computational studies have generally focused on specific processes or mechanistic hypotheses, recent methodological and computational advances harken the advent of their principal role in discovery. Here we report the development and application of the ab initio nanoreactor--a highly accelerated first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor, we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, which provide new insight into the classic Urey-Miller experiment. These results highlight the emergence of theoretical and computational chemistry as a tool for discovery, in addition to its traditional role of interpreting experimental findings. PMID:25411881
Discovering chemistry with an ab initio nanoreactor
NASA Astrophysics Data System (ADS)
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S.; Martínez, Todd J.
2014-12-01
Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provide detailed physical insight. Although theoretical and computational studies have generally focused on specific processes or mechanistic hypotheses, recent methodological and computational advances harken the advent of their principal role in discovery. Here we report the development and application of the ab initio nanoreactor—a highly accelerated first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor, we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, which provide new insight into the classic Urey-Miller experiment. These results highlight the emergence of theoretical and computational chemistry as a tool for discovery, in addition to its traditional role of interpreting experimental findings.
Ab Initio Modeling of Molecular Radiation
NASA Technical Reports Server (NTRS)
Jaffe, Richard; Schwenke, David
2014-01-01
Radiative emission from excited states of atoms and molecules can comprise a significant fraction of the total heat flux experienced by spacecraft during atmospheric entry at hypersonic speeds. For spacecraft with ablating heat shields, some of this radiative flux can be absorbed by molecular constituents in the boundary layer that are formed by the ablation process. Ab initio quantum mechanical calculations are carried out to predict the strengths of these emission and absorption processes. This talk will describe the methods used in these calculations using, as examples, the 4th positive emission bands of CO and the 1g+ 1u+ absorption in C3. The results of these calculations are being used as input to NASA radiation modeling codes like NeqAir, HARA and HyperRad.
Ab Initio Calculation of the Hoyle State
Epelbaum, Evgeny; Krebs, Hermann; Lee, Dean; Meissner, Ulf-G.
2011-05-13
The Hoyle state plays a crucial role in the helium burning of stars heavier than our Sun and in the production of carbon and other elements necessary for life. This excited state of the carbon-12 nucleus was postulated by Hoyle as a necessary ingredient for the fusion of three alpha particles to produce carbon at stellar temperatures. Although the Hoyle state was seen experimentally more than a half century ago nuclear theorists have not yet uncovered the nature of this state from first principles. In this Letter we report the first ab initio calculation of the low-lying states of carbon-12 using supercomputer lattice simulations and a theoretical framework known as effective field theory. In addition to the ground state and excited spin-2 state, we find a resonance at -85(3) MeV with all of the properties of the Hoyle state and in agreement with the experimentally observed energy.
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
NASA Astrophysics Data System (ADS)
Martinez, Todd
Traditional approaches for modeling chemical reaction networks such as those involved in combustion have focused on identifying individual reactions and using theoretical approaches to explore the underlying mechanisms. Recent advances involving graphical processing units (GPUs), commodity products developed for the videogaming industry, have made it possible to consider a distinct approach wherein one attempts to discover chemical reactions and mechanisms. We provide a brief summary of these developments and then discuss the concept behind the ``ab initio nanoreactor'' which explores the space of possible chemical reactions and molecular species for a given stoichiometry. The nanoreactor concept is exemplified with an example to the Urey-Miller reaction network which has been previously advanced as a potential model for prebiotic chemistry. We briefly discuss some of the future directions envisioned for the development of this nanoreactor concept.
Discovering chemistry with an ab initio nanoreactor
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S.; Martínez, Todd J.
2014-01-01
Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provides detailed physical insight. While theoretical and computational studies have generally focused on specific processes or mechanistic hypotheses, recent methodological and computational advances harken the advent of their principal role in discovery. Here we report the development and application of the ab initio nanoreactor – a highly accelerated, first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, providing new insight into the classic Urey-Miller experiment. These results highlight the emergence of theoretical and computational chemistry as a tool for discovery in addition to its traditional role of interpreting experimental findings. PMID:25411881
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
Ab initio derivation of model energy density functionals
NASA Astrophysics Data System (ADS)
Dobaczewski, Jacek
2016-08-01
I propose a simple and manageable method that allows for deriving coupling constants of model energy density functionals (EDFs) directly from ab initio calculations performed for finite fermion systems. A proof-of-principle application allows for linking properties of finite nuclei, determined by using the nuclear nonlocal Gogny functional, to the coupling constants of the quasilocal Skyrme functional. The method does not rely on properties of infinite fermion systems but on the ab initio calculations in finite systems. It also allows for quantifying merits of different model EDFs in describing the ab initio results.
Ab initio effective interactions for s d -shell valence nucleons
NASA Astrophysics Data System (ADS)
Dikmen, E.; Lisetskiy, A. F.; Barrett, B. R.; Maris, P.; Shirokov, A. M.; Vary, J. P.
2015-06-01
We perform ab initio no-core shell-model calculations for A =18 and 19 nuclei in a 4 ℏ Ω , or Nmax=4 , model space by using the effective JISP16 and chiral N3LO nucleon-nucleon potentials and transform the many-body effective Hamiltonians into the 0 ℏ Ω model space to construct the A -body effective Hamiltonians in the s d shell. We separate the A -body effective Hamiltonians with A =18 and A =19 into inert core, one-, and two-body components. Then we use these core, one-, and two-body components to perform standard shell-model calculations for the A =18 and A =19 systems with valence nucleons restricted to the s d shell. Finally, we compare the standard shell-model results in the 0 ℏ Ω model space with the exact no-core shell-model results in the 4 ℏ Ω model space for the A =18 and A =19 systems and find good agreement.
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.
Microsolvation of methyl hydrogen peroxide: Ab initio quantum chemical approach
NASA Astrophysics Data System (ADS)
Kulkarni, Anant D.; Rai, Dhurba; Bartolotti, Libero J.; Pathak, Rajeev K.
2009-08-01
Methyl hydrogen peroxide (MHP), one of the simplest organic hydroperoxides, is a strong oxidant, with enhanced activity in aqueous ambience. The present study investigates, at the molecular level, the role of hydrogen bonding that is conducive to cluster formation of MHP with water molecules from its peroxide end, with the methyl group remaining hydrophobic for up to five water molecules. Ab initio quantum chemical computations on MHP⋯(H2O)n, [n =1-5] are performed at second order Møller-Plesset (MP2) perturbation theory employing the basis sets 6-31G(d,p) and 6-311++G(2d,2p) to study the cluster formation of MHP with water molecules from its peroxide end and hydrophobic hydration due to the methyl group. Successive addition of water molecules alters the hydrogen bonding pattern, which leads to changes in overall cluster geometry and in turn to IR vibrational frequency shifts. Molecular co-operativity in these clusters is gauged directly through a detailed many-body interaction energy analysis. Molecular electrostatic potential maps are shown to have a bearing on predicting further growth of these clusters, which is duly corroborated through sample calculations for MHP⋯(H2O)8. Further, a continuum solvation model calculation for energetically stable clusters suggests that this study should serve as a precursor for pathways to aqueous solvation of MHP.
Ab initio molecular dynamics simulation of liquid water by quantum Monte Carlo
Zen, Andrea; Luo, Ye Mazzola, Guglielmo Sorella, Sandro; Guidoni, Leonardo
2015-04-14
Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article, we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous density functional theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab initio simulations of complex chemical systems.
Large-scale ab initio configuration interaction calculations for light nuclei
NASA Astrophysics Data System (ADS)
Maris, Pieter; Metin Aktulga, H.; Caprio, Mark A.; Çatalyürek, Ümit V.; Ng, Esmond G.; Oryspayev, Dossay; Potter, Hugh; Saule, Erik; Sosonkina, Masha; Vary, James P.; Yang, Chao; Zhou, Zheng
2012-12-01
In ab-initio Configuration Interaction calculations, the nuclear wavefunction is expanded in Slater determinants of single-nucleon wavefunctions and the many-body Schrodinger equation becomes a large sparse matrix problem. The challenge is to reach numerical convergence to within quantified numerical uncertainties for physical observables using finite truncations of the infinite-dimensional basis space. We discuss strategies for constructing and solving the resulting large sparse matrix eigenvalue problems on current multicore computer architectures. Several of these strategies have been implemented in the code MFDn, a hybrid MPI/OpenMP Fortran code for ab-initio nuclear structure calculations that can scale to 100,000 cores and more. Finally, we will conclude with some recent results for 12C including emerging collective phenomena such as rotational band structures using SRG evolved chiral N3LO interactions.
Ab Initio: And a New Era of Airline Pilot Training.
ERIC Educational Resources Information Center
Gesell, Laurence E.
1995-01-01
Expansion of air transportation and decreasing numbers seeking pilot training point to a shortage of qualified pilots. Ab initio training, in which candidates with no flight time are trained to air transport proficiency, could resolve the problem. (SK)
Phonocatalysis. An ab initio simulation experiment
NASA Astrophysics Data System (ADS)
Kim, Kwangnam; Kaviany, Massoud
2016-06-01
Using simulations, we postulate and show that heterocatalysis on large-bandgap semiconductors can be controlled by substrate phonons, i.e., phonocatalysis. With ab initio calculations, including molecular dynamic simulations, the chemisorbed dissociation of XeF6 on h-BN surface leads to formation of XeF4 and two surface F/h-BN bonds. The reaction pathway and energies are evaluated, and the sorption and reaction emitted/absorbed phonons are identified through spectral analysis of the surface atomic motion. Due to large bandgap, the atomic vibration (phonon) energy transfer channels dominate and among them is the match between the F/h-BN covalent bond stretching and the optical phonons. We show that the chemisorbed dissociation (the pathway activation ascent) requires absorption of large-energy optical phonons. Then using progressively heavier isotopes of B and N atoms, we show that limiting these high-energy optical phonons inhibits the chemisorbed dissociation, i.e., controllable phonocatalysis.
Ab initio two-component Ehrenfest dynamics
NASA Astrophysics Data System (ADS)
Ding, Feizhi; Goings, Joshua J.; Liu, Hongbin; Lingerfelt, David B.; Li, Xiaosong
2015-09-01
We present an ab initio two-component Ehrenfest-based mixed quantum/classical molecular dynamics method to describe the effect of nuclear motion on the electron spin dynamics (and vice versa) in molecular systems. The two-component time-dependent non-collinear density functional theory is used for the propagation of spin-polarized electrons while the nuclei are treated classically. We use a three-time-step algorithm for the numerical integration of the coupled equations of motion, namely, the velocity Verlet for nuclear motion, the nuclear-position-dependent midpoint Fock update, and the modified midpoint and unitary transformation method for electronic propagation. As a test case, the method is applied to the dissociation of H2 and O2. In contrast to conventional Ehrenfest dynamics, this two-component approach provides a first principles description of the dynamics of non-collinear (e.g., spin-frustrated) magnetic materials, as well as the proper description of spin-state crossover, spin-rotation, and spin-flip dynamics by relaxing the constraint on spin configuration. This method also holds potential for applications to spin transport in molecular or even nanoscale magnetic devices.
Ab Initio 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 two-component Ehrenfest dynamics
Ding, Feizhi; Goings, Joshua J.; Liu, Hongbin; Lingerfelt, David B.; Li, Xiaosong
2015-09-21
We present an ab initio two-component Ehrenfest-based mixed quantum/classical molecular dynamics method to describe the effect of nuclear motion on the electron spin dynamics (and vice versa) in molecular systems. The two-component time-dependent non-collinear density functional theory is used for the propagation of spin-polarized electrons while the nuclei are treated classically. We use a three-time-step algorithm for the numerical integration of the coupled equations of motion, namely, the velocity Verlet for nuclear motion, the nuclear-position-dependent midpoint Fock update, and the modified midpoint and unitary transformation method for electronic propagation. As a test case, the method is applied to the dissociation of H{sub 2} and O{sub 2}. In contrast to conventional Ehrenfest dynamics, this two-component approach provides a first principles description of the dynamics of non-collinear (e.g., spin-frustrated) magnetic materials, as well as the proper description of spin-state crossover, spin-rotation, and spin-flip dynamics by relaxing the constraint on spin configuration. This method also holds potential for applications to spin transport in molecular or even nanoscale magnetic devices.
Ab initio molecular crystal structures, spectra, and phase diagrams.
Hirata, So; Gilliard, Kandis; He, Xiao; Li, Jinjin; Sode, Olaseni
2014-09-16
Conspectus Molecular crystals are chemists' solids in the sense that their structures and properties can be understood in terms of those of the constituent molecules merely perturbed by a crystalline environment. They form a large and important class of solids including ices of atmospheric species, drugs, explosives, and even some organic optoelectronic materials and supramolecular assemblies. Recently, surprisingly simple yet extremely efficient, versatile, easily implemented, and systematically accurate electronic structure methods for molecular crystals have been developed. The methods, collectively referred to as the embedded-fragment scheme, divide a crystal into monomers and overlapping dimers and apply modern molecular electronic structure methods and software to these fragments of the crystal that are embedded in a self-consistently determined crystalline electrostatic field. They enable facile applications of accurate but otherwise prohibitively expensive ab initio molecular orbital theories such as Møller-Plesset perturbation and coupled-cluster theories to a broad range of properties of solids such as internal energies, enthalpies, structures, equation of state, phonon dispersion curves and density of states, infrared and Raman spectra (including band intensities and sometimes anharmonic effects), inelastic neutron scattering spectra, heat capacities, Gibbs energies, and phase diagrams, while accounting for many-body electrostatic (namely, induction or polarization) effects as well as two-body exchange and dispersion interactions from first principles. They can fundamentally alter the role of computing in the studies of molecular crystals in the same way ab initio molecular orbital theories have transformed research practices in gas-phase physical chemistry and synthetic chemistry in the last half century. In this Account, after a brief summary of formalisms and algorithms, we discuss applications of these methods performed in our group as compelling
Ab initio no core calculations of light nuclei and preludes to Hamiltonian quantum field theory
Vary, J. P.; Maris, P.; Honkanen, H.; Li, J.; Shirokov, A. M.; Brodsky, S. J.; Harindranath, A.
2009-12-17
Recent advances in ab initio quantum many-body methods and growth in computer power now enable highly precise calculations of nuclear structure. The precision has attained a level sufficient to make clear statements on the nature of 3-body forces in nuclear physics. Total binding energies, spin-dependent structure effects, and electroweak properties of light nuclei play major roles in pinpointing properties of the underlying strong interaction. Eventually, we anticipate a theory bridge with immense predictive power from QCD through nuclear forces to nuclear structure and nuclear reactions. Light front Hamiltonian quantum field theory offers an attractive pathway and we outline key elements.
Ab initio no core calculations of light nuclei and preludes to Hamiltonian quantum field theory
Vary, J.P.; Maris, P.; Shirokov, A.M.; Honkanen, H.; li, J.; Brodsky, S.J.; Harindranath, A.; Teramond, G.F.de; /Costa Rica U.
2009-08-03
Recent advances in ab initio quantum many-body methods and growth in computer power now enable highly precise calculations of nuclear structure. The precision has attained a level sufficient to make clear statements on the nature of 3-body forces in nuclear physics. Total binding energies, spin-dependent structure effects, and electroweak properties of light nuclei play major roles in pinpointing properties of the underlying strong interaction. Eventually,we anticipate a theory bridge with immense predictive power from QCD through nuclear forces to nuclear structure and nuclear reactions. Light front Hamiltonian quantum field theory offers an attractive pathway and we outline key elements.
Skutterudites under pressure: An ab initio study
Ram, Swetarekha; Kanchana, V.; Valsakumar, M. C.
2014-03-07
Ab initio results on the band structure, density of states, and Fermi surface (FS) properties of LaRu{sub 4}X{sub 12} (X = P, As, Sb) are presented at ambient pressure as well as under compression. The analysis of density of states reveals the major contribution at the Fermi level to be mainly from the Ru-d and X-p states. We have a complicated Fermi surface with both electron and hole characters for all the three compounds which is derived mainly from the Ru-d and X-p states. There is also a simpler FS with hole character derived from the P-p{sub z} orbital for LaRu{sub 4}P{sub 12} and Ru-d{sub z{sup 2}} orbital in the case of As and Sb containing compounds. More interestingly, Fermi surface nesting feature is observed only in the case of the LaRu{sub 4}P{sub 12}. Under compression, we observe the topology of the complicated FS sheet of LaRu{sub 4}As{sub 12} to change around V/V{sub 0} = 0.85, leading to a behaviour similar to that of a multiband superconductor, and in addition, we have two more hole pockets centered around Γ at V/V{sub 0} = 0.8 for the same compound. Apart from this, we find the hole pocket to vanish at V/V{sub 0} = 0.8 in the case of LaRu{sub 4}Sb{sub 12} and the opening of the complicated FS sheet gets reduced. The de Haas van Alphen calculation shows the number of extremal orbits in the complicated sheet to change in As and Sb containing compounds under compression, where we also observe the FS topology to change.
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.
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.
Ab Initio Studies of Stratospheric Ozone Depletion Chemistry
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Head-Gordon, Martin; Langhoff, Stephen R. (Technical Monitor)
1995-01-01
An overview of the current understanding of ozone depletion chemistry, particularly with regards the formation of the so-called Antarctic ozone hole, will be presented together with an outline as to how ab initio quantum chemistry can be used to further our understanding of stratospheric chemistry. The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of fluorine, chlorine, bromine and nitrogen oxide species will be demonstrated by presentation of some example studies. The ab initio results will be shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the theoretical results are shown to fill in the gaps and to resolve experimental controversies. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of halogen oxide species will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of experimental studies.
NASA Astrophysics Data System (ADS)
Song, Yu-Zhi; Li, Yong-Qing; Gao, Shou-Bao; Meng, Qing-Tian
2014-01-01
A globally accurate potential energy surface is reported for the electronic ground-state HLi2 by fitting ab initio energies to double many-body expansion formalism. The total 3726 ab initio energies used to map the HLi2 potential energy surface are calculated using the multi-reference configuration interaction method, with their dynamical correlation being semiempirically corrected by the double many-body expansion-scaled external correlation method. The current potential energy surface generates an excellent fit of the ab initio energies, showing a small root-mean squared derivation of 0.636 kcal mol-1. The topographical features of the HLi2 potential energy surface are examined in detail, which concludes that the H + Li2(X 1Σg) → Li + LiH(X 1Σ) reaction is essentially barrierless and the exothermicity is calculated to be 33.668 kcal mol-1, thus corroborates the available experimental and theoretical results.
Dominant Modes in Light Nuclei - Ab Initio View of Emergent Symmetries
NASA Astrophysics Data System (ADS)
Draayer, J. P.; Dytrych, T.; Launey, K. D.; Dreyfuss, A. C.; Langr, D.
2015-01-01
An innovative symmetry-guided concept is discussed with a focus on emergent symmetry patterns in complex nuclei. In particular, the ab initio symmetry-adapted no-core shell model (SA-NCSM), which capitalizes on exact as well as partial symmetries that underpin the structure of nuclei, provides remarkable insight into how simple symmetry patterns emerge in the many-body nuclear dynamics from first principles. This ab initio view is complemented by a fully microscopic no-core symplectic shell-model framework (NCSpM), which, in turn, informs key features of the primary physics responsible for the emergent phenomena of large deformation and alpha-cluster substructures in studies of the challenging Hoyle state in Carbon-12 and enhanced collectivity in intermediate-mass nuclei. Furthermore, by recognizing that deformed configurations often dominate the low-energy regime, the SA-NCSM provides a strategy for determining the nature of bound states of nuclei in terms of a relatively small subspace of the symmetry-reorganized complete model space, which opens new domains of nuclei for ab initio investigations, namely, the intermediate-mass region, including isotopes of Ne, Mg, and Si.
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
Ab initio computations of photodissociation products of CFC alternatives
Tai, S.; Illinger, K.H.; Kenny, J.E.
1995-12-31
Ab initio computations, have already been used to examine the energetics of the photodissociation of stratospheric chlorofluorocarbons. Our awn research has investigated the ab initio computation of vibrational frequencies and infrared intensities of CF{sub 3}CH{sub 2}F, CF{sub 3}CF{sub 2}H, and CF{sub 3}CH{sub 3}; continuing research will attempt to expand these computations to the energetics of the photodissociation of these molecules, since sane of the most common types of chlorofluorocarbon substitutes are hydrofluoroethanes.
Ab Initio Structure Analysis Using Laboratory Powder Diffraction Data
NASA Astrophysics Data System (ADS)
Sasaki, Akito
Today, laboratory X-ray diffractometers are seeing increasingly wide use in the ab initio crystal structure analysis of organic powder samples. This is because optics and optical devices have been improved, making it possible to obtain precise integrated intensities of reflections in high 2-theta ranges. Another reason is that one can use direct-space methods, which do not require “high-resolution diffraction data”, much more easily than before. Described here are some key points to remember when performig ab initio crystal structure analysis using powder diffraction data from organic compounds.
Takahasi, S. ); Curtiss, L.A.; Gosztola, D.; Koura, N. ); Loong, C.K.; Saboungi, M.L. . Materials Science Div.)
1993-04-01
The Raman and neutron scattering spectra of 46 mol% AlCl[sub 3] -54 mol% 1-ethyl-3-methyl imidazolium chloride (EMIC) and 67 mol% AlCl[sub 3] - 33 mol% EMIC melts are presented. Ab initio molecular orbital calculations have been carried out on structures of chloroaluminate anion and EMI cation and the interaction between anion and cation.
Towards SiC Surface Functionalization: An Ab Initio Study
Cicero, G; Catellani, A
2005-01-28
We present a microscopic model of the interaction and adsorption mechanism of simple organic molecules on SiC surfaces as obtained from ab initio molecular dynamics simulations. Our results open the way to functionalization of silicon carbide, a leading candidate material for bio-compatible devices.
Ab initio pseudopotential band calculation of organic conductors
Ishibashi, Shoji; Kohyama, Masanori
1999-12-01
The authors have calculated the band structures of organic conductors TTF-TCNQ and {beta}-(BEDT-TTF){sub 2}I{sub 3} using the ab initio plane-wave pseudopotential method within the local-density approximation (LDA). The Fermi-surface shape and the origin of bands near the Fermi level are investigated for each compound.
Multiple time step integrators in ab initio molecular dynamics
Luehr, Nathan; Martínez, Todd J.; Markland, Thomas E.
2014-02-28
Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties associated with splitting the ab initio potential into fast and slowly varying components. Here we present two schemes that enable efficient time-scale separation in ab initio calculations: one based on fragment decomposition and the other on range separation of the Coulomb operator in the electronic Hamiltonian. We demonstrate for both water clusters and a solvated hydroxide ion that multiple time-scale molecular dynamics allows for outer time steps of 2.5 fs, which are as large as those obtained when such schemes are applied to empirical potentials, while still allowing for bonds to be broken and reformed throughout the dynamics. This permits computational speedups of up to 4.4x, compared to standard Born-Oppenheimer ab initio molecular dynamics with a 0.5 fs time step, while maintaining the same energy conservation and accuracy.
Motif based Hessian matrixfor ab initio geometry optimization ofnanostructures
Zhao, Zhengji; Wang, Lin-Wang; Meza, Juan
2006-04-05
A simple method to estimate the atomic degree Hessian matrixof a nanosystem is presented. The estimated Hessian matrix, based on themotif decomposition of the nanosystem, can be used to accelerate abinitio atomic relaxations with speedups of 2 to 4 depending on the sizeof the system. In addition, the programing implementation for using thismethod in a standard ab initio package is trivial.
The In-Medium Similarity Renormalization Group: A novel ab initio method for nuclei
NASA Astrophysics Data System (ADS)
Hergert, H.; Bogner, S. K.; Morris, T. D.; Schwenk, A.; Tsukiyama, K.
2016-03-01
We present a comprehensive review of the In-Medium Similarity Renormalization Group (IM-SRG), a novel ab initio method for nuclei. The IM-SRG employs a continuous unitary transformation of the many-body Hamiltonian to decouple the ground state from all excitations, thereby solving the many-body problem. Starting from a pedagogical introduction of the underlying concepts, the IM-SRG flow equations are developed for systems with and without explicit spherical symmetry. We study different IM-SRG generators that achieve the desired decoupling, and how they affect the details of the IM-SRG flow. Based on calculations of closed-shell nuclei, we assess possible truncations for closing the system of flow equations in practical applications, as well as choices of the reference state. We discuss the issue of center-of-mass factorization and demonstrate that the IM-SRG ground-state wave function exhibits an approximate decoupling of intrinsic and center-of-mass degrees of freedom, similar to Coupled Cluster (CC) wave functions. To put the IM-SRG in context with other many-body methods, in particular many-body perturbation theory and non-perturbative approaches like CC, a detailed perturbative analysis of the IM-SRG flow equations is carried out. We conclude with a discussion of ongoing developments, including IM-SRG calculations with three-nucleon forces, the multi-reference IM-SRG for open-shell nuclei, first non-perturbative derivations of shell-model interactions, and the consistent evolution of operators in the IM-SRG. We dedicate this review to the memory of Gerry Brown, one of the pioneers of many-body calculations of nuclei.
Gravitational Many-Body Problem
Makino, J.
2008-04-29
In this paper, we briefly review some aspects of the gravitational many-body problem, which is one of the oldest problems in the modern mathematical science. Then we review our GRAPE project to design computers specialized to this problem.
Ab initio Monte Carlo investigation of small lithium clusters.
Srinivas, S.
1999-06-16
Structural and thermal properties of small lithium clusters are studied using ab initio-based Monte Carlo simulations. The ab initio scheme uses a Hartree-Fock/density functional treatment of the electronic structure combined with a jump-walking Monte Carlo sampling of nuclear configurations. Structural forms of Li{sub 8} and Li{sub 9}{sup +} clusters are obtained and their thermal properties analyzed in terms of probability distributions of the cluster potential energy, average potential energy and configurational heat capacity all considered as a function of the cluster temperature. Details of the gradual evolution with temperature of the structural forms sampled are examined. Temperatures characterizing the onset of structural changes and isomer coexistence are identified for both clusters.
Recent progress in ab initio density matrix renormalization group methodology
NASA Astrophysics Data System (ADS)
Hachmann, Johannes; Dorando, Jonathan J.; Kin-Lic Chan, Garnet
2008-03-01
We present some recent developments in the ab initio density matrix renormalization group (DMRG) method for quantum chemical problems, in particular our local, quadratic scaling algorithm [1] for low dimensional systems. This method is particularly suited for the description of strong nondynamic correlation, and allows us to compute numerically exact (FCI) correlated energies for large active spaces, up to one order of magnitude larger then can be done by conventional CASCI techniques. Other features of this method are its inherent multireference nature, compactness, variational results, size-consistency and size-extensivity. In addition we will review the problems (predominantly organic electronic materials) on which we applied the ab initio DMRG: 1) metal-insulator transition in hydrogen chains [1] 2) all-trans polyacetylene [1] 3) acenes [2] 4) polydiacetylenes [3]. References [1] Hachmann, Cardoen, Chan, JCP 125 (2006), 144101. [2] Hachmann, Dorando, Avil'es, Chan, JCP 127 (2007), 134309. [3] unpublished.
Understanding phonon transport in thermoelectric materials using ab initio approaches
NASA Astrophysics Data System (ADS)
Broido, David
Good thermoelectric materials have low phonon thermal conductivity, kph. Accurate theories to describe kph are important components in developing predictive models of thermoelectric efficiency that can help guide synthesis and measurement efforts. We have developed ab initio approaches to calculate kph, in which phonon modes and phonon scattering rates are computed using interatomic force constants determined from density functional theory, and a full solution of the Boltzmann transport equation for phonons is implemented. A recent approach to calculate interatomic force constants using ab initio molecular dynamics has yielded a good description of the thermal properties of Bi2Te3. But, the complexity of new promising candidate thermoelectric materials introduces computational challenges in assessing their thermal properties. An example is germanane, a germanium based hydrogen-terminated layered semiconductor, which we will discuss in this talk.
Towards AB Initio Calculation of the Circular Dichroism of Peptides
NASA Astrophysics Data System (ADS)
Molteni, E.; Onida, G.; Tiana, G.
2012-08-01
In this work we plan to use ab initio spectroscopy calculations to compute circular dichroism (CD) spectra of peptides. CD provides information on protein secondary structure content; peptides, instead, remain difficult to address, due to their tendency to adopt multiple conformations in equilibrium. Therefore peptides are an interesting test-case for ab initio calculation of CD spectra. As a first application, we focus on the (83-92) fragment of HIV-1 protease, which is known to be involved in the folding and dimerization of this protein. As a preliminary step, we performed classical molecular dynamics (MD) simulations, in order to obtain a set of representative conformers of the peptide. Then, on some of the obtained conformations, we calculated absorption spectra at the independent particle, RPA and TDLDA levels, showing the presence of charge transfer excitations, and their influence on spectral features.
Ab-initio calculations on melting of thorium
NASA Astrophysics Data System (ADS)
Mukherjee, D.; Sahoo, B. D.; Joshi, K. D.; Kaushik, T. C.; Gupta, Satish C.
2016-05-01
Ab-initio molecular dynamics study has been performed on face centered cubic structured thorium to determine its melting temperature at room pressure. The ion-electron interaction potential energy calculated as a function of temperature for three volumes (a0)3 and (1.02a0)3 and (1.04a0)3 increases gradually with temperature and undergoes a sharp jump at ~2200 K, ~2100 K and ~1800 K, respectively. Here, a0 = 5.043 Å is the equilibrium lattice parameter at 0 K obtained from ab-initio calculations. These jumps in interaction energy are treated as due to the onset of melting and corresponding temperatures as melting point. The melting point of 2100 K is close to the experimental value of 2023K. Further, the same has been verified by plotting the atomic arrangement evolved at various temperatures and corresponding pair correlation functions.
Ab initio molecular dynamics: concepts, recent developments, and future trends.
Iftimie, Radu; Minary, Peter; Tuckerman, Mark E
2005-05-10
The methodology of ab initio molecular dynamics, wherein finite-temperature dynamical trajectories are generated by using forces computed "on the fly" from electronic structure calculations, has had a profound influence in modern theoretical research. Ab initio molecular dynamics allows chemical processes in condensed phases to be studied in an accurate and unbiased manner, leading to new paradigms in the elucidation of microscopic mechanisms, rationalization of experimental data, and testable predictions of new phenomena. The purpose of this work is to give a brief introduction to the technique and to review several important recent developments in the field. Several illustrative examples showing the power of the technique have been chosen. Perspectives on future directions in the field also will be given. PMID:15870204
Exploring complex chemical reactions by ab-initio simulation
NASA Astrophysics Data System (ADS)
Parrinello, Michele
1998-03-01
Recent progress in the ab-initio molecular dynamics method and the power of parallel computing, allow the detailed study of complex chemical reaction of great industrial relevance. We illustrate this unprecedented capability by investigating the second generation Ziegler-Natta catalytic process. In this inhomogeneous catalyst, a polymerization reaction is induced by TiCl4 molecules deposited on an MgCl2 solid support. A density functional based ab-initio molecular dynamics calculation conducted with a minimum of initial assumption allows to understand the nature of the catalytic center and to determine the reaction path with the associated free energy barrier. Furthermore our calculation can explain in a nontrivial way the stereo-selectivity of the process.
Ab Initio Electronic Relaxation Times and Transport in Noble Metals
NASA Astrophysics Data System (ADS)
Mustafa, Jamal I.; Bernardi, Marco; Neaton, Jeffrey B.; Louie, Steven G.
Relaxation times employed to study electron transport in metals are typically assumed to be constants and obtained empirically using the Drude model. Here, we employ ab initio calculations to compute the electron-phonon relaxation times of Cu, Ag, and Au, and find that they vary significantly on the Fermi surface, spanning ~15 -45 fs. We compute room temperature resistivities in excellent agreement with experiment by combining GW bandstructures, Wannier-interpolated band velocities, and ab initio relaxation times. Our calculations are compared to other approximations used for the relaxation times. Additionally, an importance sampling scheme is introduced to speed up the convergence of resistivity and transport calculations by sampling directly points on the Fermi surface. This work was supported by NSF Grant No. DMR15-1508412 and U.S. DOE under Contract No. DE-AC02-05CH11231. Computational resources have been provided by DOE at LBNL's NERSC facility.
Ab Initio Calculations Of Light-Ion Reactions
Navratil, P; Quaglioni, S; Roth, R; Horiuchi, W
2012-03-12
The exact treatment of nuclei starting from the constituent nucleons and the fundamental interactions among them has been a long-standing goal in nuclear physics. In addition to the complex nature of nuclear forces, one faces the quantum-mechanical many-nucleon problem governed by an interplay between bound and continuum states. In recent years, significant progress has been made in ab initio nuclear structure and reaction calculations based on input from QCD employing Hamiltonians constructed within chiral effective field theory. In this contribution, we present one of such promising techniques capable of describing simultaneously both bound and scattering states in light nuclei. By combining the resonating-group method (RGM) with the ab initio no-core shell model (NCSM), we complement a microscopic cluster approach with the use of realistic interactions and a microscopic and consistent description of the clusters. We discuss applications to light nuclei scattering, radiative capture and fusion reactions.
Spin-orbit decomposition of ab initio nuclear wave functions
NASA Astrophysics Data System (ADS)
Johnson, Calvin W.
2015-03-01
Although the modern shell-model picture of atomic nuclei is built from single-particle orbits with good total angular momentum j , leading to j -j coupling, decades ago phenomenological models suggested that a simpler picture for 0 p -shell nuclides can be realized via coupling of the total spin S and total orbital angular momentum L . I revisit this idea with large-basis, no-core shell-model calculations using modern ab initio two-body interactions and dissect the resulting wave functions into their component L - and S -components. Remarkably, there is broad agreement with calculations using the phenomenological Cohen-Kurath forces, despite a gap of nearly 50 years and six orders of magnitude in basis dimensions. I suggest that L -S decomposition may be a useful tool for analyzing ab initio wave functions of light nuclei, for example, in the case of rotational bands.
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.
GAUSSIAN 76: An ab initio Molecular Orbital Program
DOE R&D Accomplishments Database
Binkley, J. S.; Whiteside, R.; Hariharan, P. C.; Seeger, R.; Hehre, W. J.; Lathan, W. A.; Newton, M. D.; Ditchfield, R.; Pople, J. A.
1978-01-01
Gaussian 76 is a general-purpose computer program for ab initio Hartree-Fock molecular orbital calculations. It can handle basis sets involving s, p and d-type Gaussian functions. Certain standard sets (STO-3G, 4-31G, 6-31G*, etc.) are stored internally for easy use. Closed shell (RHF) or unrestricted open shell (UHF) wave functions can be obtained. Facilities are provided for geometry optimization to potential minima and for limited potential surface scans.
The study of molecular spectroscopy by ab initio methods
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.
1991-01-01
This review illustrates the potential of theory for solving spectroscopic problems. The accuracy of approximate techniques for including electron correlation have been calibrated by comparison with full configuration-interaction calculations. Examples of the application of ab initio calculations to vibrational, rotational, and electronic spectroscopy are given. It is shown that the state-averaged, complete active space self-consistent field, multireference configuration-interaction procedure provides a good approach for treating several electronic states accurately in a common molecular orbital basis.
NASA Astrophysics Data System (ADS)
Kang, Youngho; Jeon, Sang Ho; Cho, Youngmi; Han, Seungwu
2016-01-01
We investigate the vertical ionization potential (IP) and electron affinity (EA) of organic semiconductors in the solid state that govern the optoelectrical property of organic devices using a fully ab initio way. The present method combines the density functional theory and many-body perturbation theory based on G W approximations. To demonstrate the accuracy of this approach, we carry out calculations on several prototypical organic molecules. Since IP and EA depend on the molecular orientation at the surface, the molecular geometry of the surface is explicitly considered through the slab model. The computed IP and EA are in reasonable and consistent agreements with spectroscopic data on organic surfaces with various molecular arrangements. However, the transport gaps are slightly underestimated in calculations, which can be explained by different screening effects between surface and bulk regions.
Steinmann, Casper; Fedorov, Dmitri G.; Jensen, Jan H.
2013-01-01
We extend the Effective Fragment Molecular Orbital (EFMO) method to the frozen domain approach where only the geometry of an active part is optimized, while the many-body polarization effects are considered for the whole system. The new approach efficiently mapped out the entire reaction path of chorismate mutase in less than four days using 80 cores on 20 nodes, where the whole system containing 2398 atoms is treated in the ab initio fashion without using any force fields. The reaction path is constructed automatically with the only assumption of defining the reaction coordinate a priori. We determine the reaction barrier of chorismate mutase to be kcal mol−1 for MP2/cc-pVDZ and for MP2/cc-pVTZ in an ONIOM approach using EFMO-RHF/6-31G(d) for the high and low layers, respectively. PMID:23593259
Palummo, Maurizia; Hogan, Conor; Sottile, Francesco; Bagalá, Paolo; Rubio, Angel
2009-08-28
We present a theoretical investigation of electronic and optical properties of free-base porphyrins based on density functional theory and many-body perturbation theory. The electronic levels of free-base porphine (H(2)P) and its phenyl derivative, free-base tetraphenylporphyrin (H(2)TPP) are calculated using the ab initio GW approximation for the self-energy. The approach is found to yield results that compare favorably with the available photoemission spectra. The excitonic nature of the optical peaks is revealed by solving the Bethe-Salpeter equation, which provides an accurate description of the experimental absorption spectra. The lowest triplet transition energies are in good agreement with the measured values. PMID:19725603
Global ab initio ground-state potential energy surface of N4
NASA Astrophysics Data System (ADS)
Paukku, Yuliya; Yang, Ke R.; Varga, Zoltan; Truhlar, Donald G.
2013-07-01
We present a global ground-state potential energy surface for N4 suitable for treating high-energy vibrational-rotational energy transfer and collision-induced dissociation in N2-N2 collisions. To obtain the surface, complete active space second-order perturbation theory calculations were performed for the ground singlet state with an active space of 12 electrons in 12 orbitals and the maug-cc-pVTZ triple zeta basis set. About 17 000 ab initio data points have been calculated for the N4 system, distributed along nine series of N2 + N2 geometries and three series of N3 + N geometries. The six-dimensional ground-state potential energy surface is fitted using least-squares fits to the many-body component of the electronic energies based on permutationally invariant polynomials in bond order variables.
Velaga, Srinath C; Anderson, Brian J
2014-01-16
Gas hydrate deposits are receiving increased attention as potential locations for CO2 sequestration, with CO2 replacing the methane that is recovered as an energy source. In this scenario, it is very important to correctly characterize the cage occupancies of CO2 to correctly assess the sequestration potential as well as the methane recoverability. In order to predict accurate cage occupancies, the guest–host interaction potential must be represented properly. Earlier, these potential parameters were obtained by fitting to experimental equilibrium data and these fitted parameters do not match with those obtained by second virial coefficient or gas viscosity data. Ab initio quantum mechanical calculations provide an independent means to directly obtain accurate intermolecular potentials. A potential energy surface (PES) between H2O and CO2 was computed at the MP2/aug-cc-pVTZ level and corrected for basis set superposition error (BSSE), an error caused due to the lower basis set, by using the half counterpoise method. Intermolecular potentials were obtained by fitting Exponential-6 and Lennard-Jones 6-12 models to the ab initio PES, correcting for many-body interactions. We denoted this model as the “VAS” model. Reference parameters for structure I carbon dioxide hydrate were calculated using the VAS model (site–site ab initio intermolecular potentials) as Δμ(w)(0) = 1206 ± 2 J/mol and ΔH(w)(0) = 1260 ± 12 J/mol. With these reference parameters and the VAS model, pure CO2 hydrate equilibrium pressure was predicted with an average absolute deviation of less than 3.2% from the experimental data. Predictions of the small cage occupancy ranged from 32 to 51%, and the large cage is more than 98% occupied. The intermolecular potentials were also tested by calculating the pure CO2 density and diffusion of CO2 in water using molecular dynamics simulations. PMID:24328234
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics
Makhov, Dmitry V.; Shalashilin, Dmitrii V.; Glover, William J.; Martinez, Todd J.
2014-08-07
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as “cloning,” in analogy to the “spawning” procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, “trains,” as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions.
Ab Initio Nuclear Structure and Reaction Calculations for Rare Isotopes
Draayer, Jerry P.
2014-09-28
We have developed a novel ab initio symmetry-adapted no-core shell model (SA-NCSM), which has opened the intermediate-mass region for ab initio investigations, thereby providing an opportunity for first-principle symmetry-guided applications to nuclear structure and reactions for nuclear isotopes from the lightest p-shell systems to intermediate-mass nuclei. This includes short-lived proton-rich nuclei on the path of X-ray burst nucleosynthesis and rare neutron-rich isotopes to be produced by the Facility for Rare Isotope Beams (FRIB). We have provided ab initio descriptions of high accuracy for low-lying (including collectivity-driven) states of isotopes of Li, He, Be, C, O, Ne, Mg, Al, and Si, and studied related strong- and weak-interaction driven reactions that are important, in astrophysics, for further understanding stellar evolution, X-ray bursts and triggering of s, p, and rp processes, and in applied physics, for electron and neutrino-nucleus scattering experiments as well as for fusion ignition at the National Ignition Facility (NIF).
A Complete and Accurate Ab Initio Repeat Finding Algorithm.
Lian, Shuaibin; Chen, Xinwu; Wang, Peng; Zhang, Xiaoli; Dai, Xianhua
2016-03-01
It has become clear that repetitive sequences have played multiple roles in eukaryotic genome evolution including increasing genetic diversity through mutation, changes in gene expression and facilitating generation of novel genes. However, identification of repetitive elements can be difficult in the ab initio manner. Currently, some classical ab initio tools of finding repeats have already presented and compared. The completeness and accuracy of detecting repeats of them are little pool. To this end, we proposed a new ab initio repeat finding tool, named HashRepeatFinder, which is based on hash index and word counting. Furthermore, we assessed the performances of HashRepeatFinder with other two famous tools, such as RepeatScout and Repeatfinder, in human genome data hg19. The results indicated the following three conclusions: (1) The completeness of HashRepeatFinder is the best one among these three compared tools in almost all chromosomes, especially in chr9 (8 times of RepeatScout, 10 times of Repeatfinder); (2) in terms of detecting large repeats, HashRepeatFinder also performed best in all chromosomes, especially in chr3 (24 times of RepeatScout and 250 times of Repeatfinder) and chr19 (12 times of RepeatScout and 60 times of Repeatfinder); (3) in terms of accuracy, HashRepeatFinder can merge the abundant repeats with high accuracy. PMID:26272474
Ab Initio Study of Defect Properties in YPO4
Gao, Fei; Xiao, Haiyan Y.; Zhou, Yungang; Devanathan, Ramaswami; Hu, Shenyang Y.; Li, Yulan; Sun, Xin; Khaleel, Mohammad A.
2012-03-01
Ab initio methods based on density functional theory have been used to calculate the formation energies of intrinsic defects, including vacancies, interstitials, antisites and Frenkel pairs in YPO4 under the O-rich and Y2O3-rich, and the O-rich and Y-rich conditions. The larger size of the yttrium atom may give rise to higher formation energy of the phosphorus antisite defect. In general, the formation energies of anion interstitials are much smaller than those of cation interstitials for both conditions considered. It is of greatly interest to find that the relative stabilities among the same types of interstitials are independent of the reference states. The most stable configuration for oxygen interstitials is an O-O split interstitial near the Ta site, while the most stable configuration for cation interstitials is a tetrahedral interstitial near the Ta site. The cation split interstitials are unfavorable in YPO4, with much higher formation energies. Furthermore, the properties of Frenkel pairs are compared with those calculated using empirical potentials. The results reveal that both ab initio and empirical potential calculations show a similar trend in the formation energies of Frenkel pairs, but the formation energies obtained by empirical potentials are much larger than those calculated by ab initio method.
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics
NASA Astrophysics Data System (ADS)
Makhov, Dmitry V.; Glover, William J.; Martinez, Todd J.; Shalashilin, Dmitrii V.
2014-08-01
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as "cloning," in analogy to the "spawning" procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions.
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics.
Makhov, Dmitry V; Glover, William J; Martinez, Todd J; Shalashilin, Dmitrii V
2014-08-01
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as "cloning," in analogy to the "spawning" procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions. PMID:25106573
Can an ab initio three-body virial equation describe the mercury gas phase?
Wiebke, J; Wormit, M; Hellmann, R; Pahl, E; Schwerdtfeger, P
2014-03-27
We report a sixth-order ab initio virial equation of state (EOS) for mercury. The virial coefficients were determined in the temperature range from 500 to 7750 K using a three-body approximation to the N-body interaction potential. The underlying two-body and three-body potentials were fitted to highly accurate Coupled-Cluster interaction energies of Hg2 (Pahl, E.; Figgen, D.; Thierfelder, C.; Peterson, K. A.; Calvo, F.; Schwerdtfeger, P. J. Chem. Phys. 2010, 132, 114301-1) and equilateral-triangular configurations of Hg3. We find the virial coefficients of order four and higher to be negative and to have large absolute values over the entire temperature range considered. The validity of our three-body, sixth-order EOS seems to be limited to small densities of about 1.5 g cm(-3) and somewhat higher densities at higher temperatures. Termwise analysis and comparison to experimental gas-phase data suggest a small convergence radius of the virial EOS itself as well as a failure of the three-body interaction model (i.e., poor convergence of the many-body expansion for mercury). We conjecture that the nth-order term of the virial EOS is to be evaluated from the full n-body interaction potential for a quantitative picture. Consequently, an ab initio three-body virial equation cannot describe the mercury gas phase. PMID:24547987
4He+n+n continuum within an ab initio framework
Romero-Redondo, Carolina; Quaglioni, Sofia; Navratil, Petr; Hupin, Guillaume
2014-07-16
In this study, the low-lying continuum spectrum of the 6He nucleus is investigated for the first time within an ab initio framework that encompasses the 4He+n+n three-cluster dynamics characterizing its lowest decay channel. This is achieved through an extension of the no-core shell model combined with the resonating-group method, in which energy-independent nonlocal interactions among three nuclear fragments can be calculated microscopically, starting from realistic nucleon-nucleon interactions and consistent ab initio many-body wave functions of the clusters. The three-cluster Schrödinger equation is solved with three-body scattering boundary conditions by means of the hyperspherical-harmonics method on a Lagrange mesh. Using amore » soft similarity-renormalization-group evolved chiral nucleon-nucleon potential, we find the known Jπ = 2+ resonance as well as a result consistent with a new low-lying second 2+ resonance recently observed at GANIL at ~2.6 MeV above the He6 ground state. We also find resonances in the 2–, 1+, and 0– channels, while no low-lying resonances are present in the 0+ and 1– channels.« less
Ab Initio Study of Hot Carriers in the First Picosecond after Sunlight Absorption in Silicon
NASA Astrophysics Data System (ADS)
Bernardi, Marco; Vigil-Fowler, Derek; Lischner, Johannes; Neaton, Jeffrey B.; Louie, Steven G.
2014-06-01
Hot carrier thermalization is a major source of efficiency loss in solar cells. Because of the subpicosecond time scale and complex physics involved, a microscopic characterization of hot carriers is challenging even for the simplest materials. We develop and apply an ab initio approach based on density functional theory and many-body perturbation theory to investigate hot carriers in semiconductors. Our calculations include electron-electron and electron-phonon interactions, and require no experimental input other than the structure of the material. We apply our approach to study the relaxation time and mean free path of hot carriers in Si, and map the band and k dependence of these quantities. We demonstrate that a hot carrier distribution characteristic of Si under solar illumination thermalizes within 350 fs, in excellent agreement with pump-probe experiments. Our work sheds light on the subpicosecond time scale after sunlight absorption in Si, and constitutes a first step towards ab initio quantification of hot carrier dynamics in materials.
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.
Ab initio quantum Monte Carlo calculations of ground-state properties of manganese's oxides
NASA Astrophysics Data System (ADS)
Sharma, Vinit; Krogel, Jaron T.; Kent, P. R. C.; Reboredo, Fernando A.
One of the critical scientific challenges of contemporary research is to obtain an accurate theoretical description of the electronic properties of strongly correlated systems such as transition metal oxides and rare-earth compounds, since state-of-art ab-initio methods based on approximate density functionals are not always sufficiently accurate. Quantum Monte Carlo (QMC) methods, which use statistical sampling to evaluate many-body wave functions, have the potential to answer this challenge. Owing to the few fundamental approximations made and the direct treatment of electron correlation, QMC methods are among the most accurate electronic structure methods available to date. We assess the accuracy of the diffusion Monte Carlo method in the case of rocksalt manganese oxide (MnO). We study the electronic properties of this strongly-correlated oxide, which has been identified as a suitable candidate for many applications ranging from catalysts to electronic devices. ``This work was supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division.'' Ab initio quantum Monte Carlo calculations of ground-state properties of manganese's oxides.
Three-cluster dynamics within an ab initio framework
Quaglioni, Sofia; Romero-Redondo, Carolina; Navratil, Petr
2013-09-26
In this study, we introduce a fully antisymmetrized treatment of three-cluster dynamics within the ab initio framework of the no-core shell model/resonating-group method. Energy-independent nonlocal interactions among the three nuclear fragments are obtained from realistic nucleon-nucleon interactions and consistent ab initio many-body wave functions of the clusters. The three-cluster Schrödinger equation is solved with bound-state boundary conditions by means of the hyperspherical-harmonic method on a Lagrange mesh. We discuss the formalism in detail and give algebraic expressions for systems of two single nucleons plus a nucleus. Using a soft similarity-renormalization-group evolved chiral nucleon-nucleon potential, we apply the method to a ^{4}He+n+n description of ^{6}He and compare the results to experiment and to a six-body diagonalization of the Hamiltonian performed within the harmonic-oscillator expansions of the no-core shell model. Differences between the two calculations provide a measure of core (^{4}He) polarization effects.
Lefkidis, Georgios; Zhang, G P; Hübner, W
2009-11-20
In an ultrafast laser-induced magnetization-dynamics scenario we demonstrate for the first time an exact microscopic spin-switch mechanism. Combining ab initio electronic many-body theory and quantum optics analysis we show in detail how the coherently induced material polarization for every elementary process leads to angular-momentum exchange between the light and the irradiated antiferromagnetic NiO. Thus we answer the long-standing question where the angular momentum goes. The calculation also predicts a dynamic Kerr effect, which provides a signature for monitoring spin dynamics, by simply measuring the transient rotation and ellipticity of the reflected light. PMID:20366066
NASA Astrophysics Data System (ADS)
Duguet, T.; Bender, M.; Ebran, J.-P.; Lesinski, T.; Somà, V.
2015-12-01
This programmatic paper lays down the possibility to reconcile the necessity to resum many-body correlations into the energy kernel with the fact that safe multi-reference energy density functional (EDF) calculations cannot be achieved whenever the Pauli principle is not enforced, as is for example the case when many-body correlations are parametrized under the form of empirical density dependencies. Our proposal is to exploit a newly developed ab initio many-body formalism to guide the construction of safe, explicitly correlated and systematically improvable parametrizations of the off-diagonal energy and norm kernels that lie at the heart of the nuclear EDF method. The many-body formalism of interest relies on the concepts of symmetry breaking and restoration that have made the fortune of the nuclear EDF method and is, as such, amenable to this guidance. After elaborating on our proposal, we briefly outline the project we plan to execute in the years to come.
Chasman, R.R.
1995-08-01
In the past few years, we developed many-body variational wave functions that allow one to treat pairing and particle-hole two-body interactions on an equal footing. The complexity of these wave functions depends on the number of levels included in the valence space, but does not depend on the number of nucleons in the system. By using residual interaction strengths (e.g. the quadrupole interaction strength or pairing interaction strength) as generator coordinates, one gets many different wave functions, each having a different expectation value for the relevant interaction mode. These wave functions are particularly useful when one is dealing with a situation in which the mean-field approximation is inadequate. Because the same basis states are used in the construction of the many-body wave functions, it is possible to calculate overlaps and interaction matrix elements for the many-body wave functions (which are not in general orthogonal) easily. The valence space can contain a large number of single-particle basis states, when there are constants of motion that can be used to break the levels up into groups. We added a cranking term to the many-body Hamiltonian and modified the projection procedure to get states of good signature before variation. In our present implementation, each group is limited to eight pairs of single-particle levels. We are working on ways of increasing the number of levels that can be included in each group. We are also working on including particle-particle residual interaction modes, in addition to pairing, in our Hamiltonian.
Macromolecular ab initio phasing enforcing secondary and tertiary structure
Millán, Claudia; Sammito, Massimo; Usón, Isabel
2015-01-01
Ab initio phasing of macromolecular structures, from the native intensities alone with no experimental phase information or previous particular structural knowledge, has been the object of a long quest, limited by two main barriers: structure size and resolution of the data. Current approaches to extend the scope of ab initio phasing include use of the Patterson function, density modification and data extrapolation. The authors’ approach relies on the combination of locating model fragments such as polyalanine α-helices with the program PHASER and density modification with the program SHELXE. Given the difficulties in discriminating correct small substructures, many putative groups of fragments have to be tested in parallel; thus calculations are performed in a grid or supercomputer. The method has been named after the Italian painter Arcimboldo, who used to compose portraits out of fruit and vegetables. With ARCIMBOLDO, most collections of fragments remain a ‘still-life’, but some are correct enough for density modification and main-chain tracing to reveal the protein’s true portrait. Beyond α-helices, other fragments can be exploited in an analogous way: libraries of helices with modelled side chains, β-strands, predictable fragments such as DNA-binding folds or fragments selected from distant homologues up to libraries of small local folds that are used to enforce nonspecific tertiary structure; thus restoring the ab initio nature of the method. Using these methods, a number of unknown macromolecules with a few thousand atoms and resolutions around 2 Å have been solved. In the 2014 release, use of the program has been simplified. The software mediates the use of massive computing to automate the grid access required in difficult cases but may also run on a single multicore workstation (http://chango.ibmb.csic.es/ARCIMBOLDO_LITE) to solve straightforward cases. PMID:25610631
Ab initio study of neutron drops with chiral Hamiltonians
NASA Astrophysics Data System (ADS)
Potter, H. D.; Fischer, S.; Maris, P.; Vary, J. P.; Binder, S.; Calci, A.; Langhammer, J.; Roth, R.
2014-12-01
We report ab initio calculations for neutron drops in a 10 MeV external harmonic-oscillator trap using chiral nucleon-nucleon plus three-nucleon interactions. We present total binding energies, internal energies, radii and odd-even energy differences for neutron numbers N = 2- 18 using the no-core shell model with and without importance truncation. Furthermore, we present total binding energies for N = 8 , 16 , 20 , 28 , 40 , 50 obtained in a coupled-cluster approach. Comparisons with quantum Monte Carlo results, where available, using Argonne v8‧ with three-nucleon interactions reveal important dependences on the chosen Hamiltonian.
Ab initio vibrational and dielectric properties of Y V O
NASA Astrophysics Data System (ADS)
Vali, R.
2009-10-01
For the yttrium orthovanadate Y V O with a tetragonal zircon-type structure, the first complete set of Raman-active and IR-active phonon modes has been calculated using ab initio density functional perturbation theory. The calculated IR reflectivity spectra are in good agreement with available experimental data. We report the calculated frequencies of three Raman-active modes that could not be detected experimentally and a new assignment of the experimental Raman data. The contributions of each IR-active phonon modes to static dielectric tensor have been determined.
Ab Initio Computation of the Energies of Circular Quantum Dots
Lohne, M. Pedersen; Hagen, Gaute; Hjorth-Jensen, M.; Kvaal, S.; Pederiva, F.
2011-01-01
We perform coupled-cluster and diffusion Monte Carlo calculations of the energies of circular quantum dots up to 20 electrons. The coupled-cluster calculations include triples corrections and a renormalized Coulomb interaction defined for a given number of low-lying oscillator shells. Using such a renormalized Coulomb interaction brings the coupled-cluster calculations with triples correlations in excellent agreement with the diffusion Monte Carlo calculations. This opens up perspectives for doing ab initio calculations for much larger systems of electrons.
The implementation of ab initio quantum chemistry calculations on transporters.
Cooper, M D; Hillier, I H
1991-06-01
The RHF and geometry optimization sections of the ab initio quantum chemistry code, GAMESS, have been optimized for a network of parallel microprocessors, Inmos T800-20 transputers, using both indirect and direct SCF techniques. The results indicate great scope for implementation of such codes on small parallel computer systems, very high efficiencies having been achieved, particularly in the cases of direct SCF and geometry optimization with large basis sets. The work, although performed upon one particular parallel system, the Meiko Computing Surface, is applicable to a wide range of parallel systems with both shared and distributed memory. PMID:1919615
Ab initio electronic properties of dual phosphorus monolayers in silicon
2014-01-01
In the midst of the epitaxial circuitry revolution in silicon technology, we look ahead to the next paradigm shift: effective use of the third dimension - in particular, its combination with epitaxial technology. We perform ab initio calculations of atomically thin epitaxial bilayers in silicon, investigating the fundamental electronic properties of monolayer pairs. Quantitative band splittings and the electronic density are presented, along with effects of the layers’ relative alignment and comments on disordered systems, and for the first time, the effective electronic widths of such device components are calculated. PMID:25246862
Ab Initio Calculations Applied to Problems in Metal Ion Chemistry
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.; Partridge, Harry; Arnold, James O. (Technical Monitor)
1994-01-01
Electronic structure calculations can provide accurate spectroscopic data (such as molecular structures) vibrational frequencies, binding energies, etc.) that have been very useful in explaining trends in experimental data and in identifying incorrect experimental measurements. In addition, ab initio calculations. have given considerable insight into the many interactions that make the chemistry of transition metal systems so diverse. In this review we focus on cases where calculations and experiment have been used to solve interesting chemical problems involving metal ions. The examples include cases where theory was used to differentiate between disparate experimental values and cases where theory was used to explain unexpected experimental results.
Ab initio study of hydrogen on beryllium surfaces
NASA Astrophysics Data System (ADS)
Bachurin, D. V.; Vladimirov, P. V.
2015-11-01
Static ab initio calculations were performed for five principal hexagonal close-packed beryllium surfaces: basal, prismatic (type I and II) and pyramidal (type I and II). The basal plane was found to be the most energetically favorable, while the energies of the prismatic (type I) and pyramidal (type I) planes were slightly higher followed by the type II planes. Beryllium is known to show extreme interlayer distance relaxation near the surface. Up to five outermost atomic layers were involved in surface relaxation. The presence of hydrogen on the beryllium surfaces led to a noticeable reduction of the surface energy.
An improved ab initio structure for fluorine peroxide (FOOF)
NASA Astrophysics Data System (ADS)
Mack, Hans-Georg; Oberhammer, Heinz
1988-03-01
Ab initio calculations with the 6-31G* and Dunning (9s5p/4s2p) basis sets augmented with p and d functions at various levels of theory (RHF, MP2, MP3, and MP4) were carried out on F 2O 2. The best result was obtained at the MP2 level with the Dunning basis plus one set of d functions on fluorine and two sets of d functions on oxygen. These calculations reproduce the experimental bond lengths to within 0.01 Å and the angles to within the experimental uncertainties.
Ab-initio study of transition metal hydrides
Sharma, Ramesh; Shukla, Seema Dwivedi, Shalini Sharma, Yamini
2014-04-24
We have performed ab initio self consistent calculations based on Full potential linearized augmented plane wave (FP-LAPW) method to investigate the optical and thermal properties of yttrium hydrides. From the band structure and density of states, the optical absorption spectra and specific heats have been calculated. The band structure of Yttrium metal changes dramatically due to hybridization of Y sp orbitals with H s orbitals and there is a net charge transfer from metal to hydrogen site. The electrical resistivity and specific heats of yttrium hydrides are lowered but the thermal conductivity is slightly enhanced due to increase in scattering from hydrogen sites.
Potential Dependence of Electrochemical Barriers from ab Initio Calculations.
Chan, Karen; Nørskov, Jens K
2016-05-01
We present a simple and computationally efficient method to determine the potential dependence of the activation energies for proton-electron transfer from a single ab initio barrier calculation. We show that the potential dependence of the activation energy is given by the partial charge transferred at the transition state. The method is evaluated against the potential dependence determined explicitly through multiple calculations at varying potential. We show that the transfer coefficient is given by the charge transferred from the initial to transition state, which has significant implications for electrochemical kinetics. PMID:27088442
Ab initio Study of He Stability in hcp-Ti
Dai, Yunya; Yang, Li; Peng, SM; Long, XG; Gao, Fei; Zu, Xiaotao T.
2010-12-20
The stability of He in hcp-Ti was studied using ab initio method based on density functional theory. The results indicate that a single He atom prefers to occupy the tetrahedral site rather than the octahedral site. The interaction of He defects with Ti atoms has been used to explain the relative stabilities of He point defects in hcp-Ti. The relative stability of He defects in hcp-Ti is useful for He clustering and bubble nucleation in metal tritides, which provides the basis for development of improved atomistic models.
NASA Astrophysics Data System (ADS)
Boyle, J. J.; Pindzola, M. S.
2005-11-01
Preface; Contributors; Introduction; Part I. Atomic Structure: 1. Development of atomic many-body theory Ingvar Lindgren; 2. Relativistic MBPT for highly charged ions W. R. Johnson; 3. Parity nonconservation in atoms S. A. Blundell, W. R. Johnson, and J. Sapirstein; Part II. Photoionization of Atoms: 4. Single photoionization processes J. J. Boyle, and M. D. Kutzner; 5. Photoionization dominated by double excitation T. N. Chang; 6. Direct double photoionization in atoms Z. W. Liu; 7. Photoelectron angular distributions Steven T. Manson; Part III. A. Atomic Scattering - General Considerations: 8. The many-body approach to electron-atom collisions M. Ya Amusia; 9. Theoretical aspects of electron impact ionization P. L. Altick; Part III. B. Atomic Scattering - Low-Order Applications: 10. Perturbation series methods D. H. Madison; 11. Target dependence of the triply differential cross section Cheng Pan and Anthony F. Starace; 12. Overview of Thomas processes for fast mass transfer J. H. McGuire, Jack C. Straton and T. Ishihara; Part III. C. Atomic Scattering - All-Order Applications: 13. R-matrix Theory: Some Recent Applications Philip G. Burke: 14. Electron scattering: application of Dirac R-matrix theory Wasantha Wijesundera, Ian Grant and Patrick Norrington; 15. Close coupling and distorted-wave theory D. C. Griffin and M. S. Pindzola; Appendix: Units and notation; References; Index.
NASA Astrophysics Data System (ADS)
Boyle, J. J.; Pindzola, M. S.
1998-09-01
Preface; Contributors; Introduction; Part I. Atomic Structure: 1. Development of atomic many-body theory Ingvar Lindgren; 2. Relativistic MBPT for highly charged ions W. R. Johnson; 3. Parity nonconservation in atoms S. A. Blundell, W. R. Johnson, and J. Sapirstein; Part II. Photoionization of Atoms: 4. Single photoionization processes J. J. Boyle, and M. D. Kutzner; 5. Photoionization dominated by double excitation T. N. Chang; 6. Direct double photoionization in atoms Z. W. Liu; 7. Photoelectron angular distributions Steven T. Manson; Part III. A. Atomic Scattering - General Considerations: 8. The many-body approach to electron-atom collisions M. Ya Amusia; 9. Theoretical aspects of electron impact ionization P. L. Altick; Part III. B. Atomic Scattering - Low-Order Applications: 10. Perturbation series methods D. H. Madison; 11. Target dependence of the triply differential cross section Cheng Pan and Anthony F. Starace; 12. Overview of Thomas processes for fast mass transfer J. H. McGuire, Jack C. Straton and T. Ishihara; Part III. C. Atomic Scattering - All-Order Applications: 13. R-matrix Theory: Some Recent Applications Philip G. Burke: 14. Electron scattering: application of Dirac R-matrix theory Wasantha Wijesundera, Ian Grant and Patrick Norrington; 15. Close coupling and distorted-wave theory D. C. Griffin and M. S. Pindzola; Appendix: Units and notation; References; Index.
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
Oxidation of GaN: An ab initio thermodynamic approach
NASA Astrophysics Data System (ADS)
Jackson, Adam J.; Walsh, Aron
2013-10-01
GaN is a wide-band-gap semiconductor used in high-efficiency light-emitting diodes and solar cells. The solid is produced industrially at high chemical purities by deposition from a vapor phase, and oxygen may be included at this stage. Oxidation represents a potential path for tuning its properties without introducing more exotic elements or extreme processing conditions. In this work, ab initio computational methods are used to examine the energy potentials and electronic properties of different extents of oxidation in GaN. Solid-state vibrational properties of Ga, GaN, Ga2O3, and a single substitutional oxygen defect have been studied using the harmonic approximation with supercells. A thermodynamic model is outlined which combines the results of ab initio calculations with data from experimental literature. This model allows free energies to be predicted for arbitrary reaction conditions within a wide process envelope. It is shown that complete oxidation is favorable for all industrially relevant conditions, while the formation of defects can be opposed by the use of high temperatures and a high N2:O2 ratio.
Ab initio thermodynamic model for magnesium carbonates and hydrates.
Chaka, Anne M; Felmy, Andrew R
2014-09-01
An ab initio thermodynamic framework for predicting properties of hydrated magnesium carbonate minerals has been developed using density-functional theory linked to macroscopic thermodynamics through the experimental chemical potentials for MgO, water, and CO2. Including semiempirical dispersion via the Grimme method and small corrections to the generalized gradient approximation of Perdew, Burke, and Ernzerhof for the heat of formation yields a model with quantitative agreement for the benchmark minerals brucite, magnesite, nesquehonite, and hydromagnesite. The model shows how small differences in experimental conditions determine whether nesquehonite, hydromagnesite, or magnesite is the result of laboratory synthesis from carbonation of brucite, and what transformations are expected to occur on geological time scales. Because of the reliance on parameter-free first-principles methods, the model is reliably extensible to experimental conditions not readily accessible to experiment and to any mineral composition for which the structure is known or can be hypothesized, including structures containing defects, substitutions, or transitional structures during solid state transformations induced by temperature changes or processes such as water, CO2, or O2 diffusion. Demonstrated applications of the ab initio thermodynamic framework include an independent means to evaluate differences in thermodynamic data for lansfordite, predicting the properties of Mg analogues of Ca-based hydrated carbonates monohydrocalcite and ikaite, which have not been observed in nature, and an estimation of the thermodynamics of barringtonite from the stoichiometry and a single experimental observation. PMID:24679248
Ab Initio Thermodynamic Model for Magnesium Carbonates and Hydrates
Chaka, Anne M.; Felmy, Andrew R.
2014-03-28
An ab initio thermodynamic framework for predicting properties of hydrated magnesium carbonate minerals has been developed using density-functional theory linked to macroscopic thermodynamics through the experimental chemical potentials for MgO, water, and CO2. Including semiempirical dispersion via the Grimme method and small corrections to the generalized gradient approximation of Perdew, Burke, and Ernzerhof for the heat of formation yields a model with quantitative agreement for the benchmark minerals brucite, magnesite, nesquehonite, and hydromagnesite. The model shows how small differences in experimental conditions determine whether nesquehonite, hydromagnesite, or magnesite is the result of laboratory synthesis from carbonation of brucite, and what transformations are expected to occur on geological time scales. Because of the reliance on parameter-free first principles methods, the model is reliably extensible to experimental conditions not readily accessible to experiment and to any mineral composition for which the structure is known or can be hypothesized, including structures containing defects, substitutions, or transitional structures during solid state transformations induced by temperature changes or processes such as water, CO2, or O2 diffusion. Demonstrated applications of the ab initio thermodynamic framework include an independent means to evaluate differences in thermodynamic data for lansfordite, predicting the properties of Mg analogs of Ca-based hydrated carbonates monohydrocalcite and ikaite which have not been observed in nature, and an estimation of the thermodynamics of barringtonite from the stoichiometry and a single experimental observation.
Entropy of Liquid Water from Ab Initio Molecular Dynamics
NASA Astrophysics Data System (ADS)
Spanu, Leonardo; Zhang, Cui; Galli, Giulia
2012-02-01
The debate on the structural properties of water has been mostly based on the calculation of pair correlation functions. However, the simulation of thermodynamic and spectroscopic quantities may be of great relevance for the characterization of liquid water properties. We have computed the entropy of liquid water using a two-phase thermodynamic model and trajectories generated by ab initio molecular dynamics simulations [1]. In an attempt to better understand the performance of several density functionals in simulating liquid water, we have performed ab initio molecular dynamics using semilocal, hybrid [2] and van der Waals density functionals [3]. We show that in all cases, at the experimental equilibrium density and at temperatures in the vicinity of 300 K, the computed entropies are underestimated, with respect to experiment, and the liquid exhibits a degree of tetrahedral order higher than in experiments. We also discuss computational strategies to simulate spectroscopic properties of water, including infrared and Raman spectra.[4pt] [1] C.Zhang, L.Spanu and G.Galli, J.Phys.Chem. B 2011 (in press)[0pt] [2] C.Zhang, D.Donadio, F.Gygi and G.Galli, J. Chem. Theory Comput. 7, 1443 (2011)[0pt] [3] C.Zhang, J.Wu, G.Galli and F.Gygi, J. Chem. Theory Comput. 7, 3061 (2011)
AB INITIO SIMULATIONS FOR MATERIAL PROPERTIES ALONG THE JUPITER ADIABAT
French, Martin; Becker, Andreas; Lorenzen, Winfried; Nettelmann, Nadine; Bethkenhagen, Mandy; Redmer, Ronald; Wicht, Johannes
2012-09-15
We determine basic thermodynamic and transport properties of hydrogen-helium-water mixtures for the extreme conditions along Jupiter's adiabat via ab initio simulations, which are compiled in an accurate and consistent data set. In particular, we calculate the electrical and thermal conductivity, the shear and longitudinal viscosity, and diffusion coefficients of the nuclei. We present results for associated quantities like the magnetic and thermal diffusivity and the kinematic shear viscosity along an adiabat that is taken from a state-of-the-art interior structure model. Furthermore, the heat capacities, the thermal expansion coefficient, the isothermal compressibility, the Grueneisen parameter, and the speed of sound are calculated. We find that the onset of dissociation and ionization of hydrogen at about 0.9 Jupiter radii marks a region where the material properties change drastically. In the deep interior, where the electrons are degenerate, many of the material properties remain relatively constant. Our ab initio data will serve as a robust foundation for applications that require accurate knowledge of the material properties in Jupiter's interior, e.g., models for the dynamo generation.
Ab initio prediction of the critical thickness of a precipitate
NASA Astrophysics Data System (ADS)
Sampath, S.; Janisch, R.
2013-09-01
Segregation and precipitation of second phases in metals and metallic alloys is an important phenomenon that has a strong influence on the mechanical properties of the material. Models exist that describe the growth of coherent, semi-coherent and incoherent precipitates. One important parameter of these models is the energy of the interface between matrix and precipitate. In this work we apply ab initio density functional theory calculations to obtain this parameter and to understand how it depends on chemical composition and mechanical strain at the interface. Our example is a metastable Mo-C phase, the body-centred tetragonal structure, which exists as a semi-coherent precipitate in body-centred cubic molybdenum. The interface of this precipitate is supposed to change from coherent to semi-coherent during the growth of the precipitate. We predict the critical thickness of the precipitate by calculating the different contributions to a semi-coherent interface energy by means of ab initio density functional theory calculations. The parameters in our model include the elastic strain energy stored in the precipitate, as well as a misfit dislocation energy that depends on the dislocation core width and the dislocation spacing. Our predicted critical thickness agrees well with experimental observations.
Ab initio dynamics of the cytochrome P450 hydroxylation reaction
NASA Astrophysics Data System (ADS)
Elenewski, Justin E.; Hackett, John C.
2015-02-01
The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis.
Ab initio dynamics of the cytochrome P450 hydroxylation reaction
Elenewski, Justin E.; Hackett, John C
2015-02-14
The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis.
Development of Novel Analytical Method for Ab Initio Powder Structural Analysis
NASA Astrophysics Data System (ADS)
Sakata, Makoto; Nishibori, Eiji; Sawa, Hiroshi
Genetic Algorithm (GA) applied to ab initio structure determination from synchrotron powder diffraction is described. It seems to have an advantage over other real space methods for ab initio structure determination because of the existence of schema theorem. As an example, the case of Prednisolone Succinate is shown in some detail. Future development of GA in crystallography is briefly described.
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.
Song, Y Z; Varandas, A J C
2009-04-01
A single-sheeted potential energy surface is reported for the electronic ground-state of H(2)S by fitting accurate multireference configuration interaction energies calculated using aug-cc-pVTZ and aug-cc-pVQZ basis sets with extrapolation of the electron correlation energy to the complete basis set limit, plus extrapolation to the complete basis set limit of the complete-active-space self-consistent field energy. A switching function formalism has been used to warrant the correct behavior at the H(2)(X (1)Sigma(g) (+))+S((1)D) and SH(X (2)Pi)+H((2)S) dissociation limits. The topographical features of the novel global potential energy surface are examined in detail, with the former being used for exploratory quasiclassical trajectory calculations of the thermal rate constant for the S((1)D)+H(2), S((1)D)+D(2), and S((1)D)+HD reactions at room temperature. A comparison with other available potential energy surfaces as well as kinetics data is also provided. PMID:19355742
Ab initio electronic stopping power of protons in bulk materials
NASA Astrophysics Data System (ADS)
Shukri, Abdullah Atef; Bruneval, Fabien; Reining, Lucia
2016-01-01
The electronic stopping power is a crucial quantity for ion irradiation: it governs the deposited heat, the damage profile, and the implantation depth. Whereas experimental data are readily available for elemental solids, the data are much more scarce for compounds. Here we develop a fully ab initio computational scheme based on linear response time-dependent density-functional theory to predict the random electronic stopping power (RESP) of materials without any empirical fitting. We show that the calculated RESP compares well with experimental data, when at full convergence, with the inclusion of the core states and of the exchange correlation. We evaluate the unexpectedly limited magnitude of the nonlinear terms in the RESP by comparing with other approaches based on the time propagation of time-dependent density-functional theory. Finally, we check the validity of a few empirical rules of thumbs that are commonly used to estimate the electronic stopping power.
Isofulminic acid, HONC: Ab initio theory and microwave spectroscopy.
Mladenović, Mirjana; Lewerenz, Marius; McCarthy, Michael C; Thaddeus, Patrick
2009-11-01
Isofulminic acid, HONC, the most energetic stable isomer of isocyanic acid HNCO, higher in energy by 84 kcal/mol, has been detected spectroscopically by rotational spectroscopy supported by coupled cluster electronic structure calculations. The fundamental rotational transitions of the normal, carbon-13, oxygen-18, and deuterium isotopic species have been detected in the centimeter band in a molecular beam by Fourier transform microwave spectroscopy, and rotational constants and nitrogen and deuterium quadrupole coupling constants have been derived. The measured constants agree well with those predicted by ab initio calculations. A number of other electronic and spectroscopic parameters of isofulminic acid, including the dipole moment, vibrational frequencies, infrared intensities, and centrifugal distortion constants have been calculated at a high level of theory. Isofulminic acid is a good candidate for astronomical detection with radio telescopes because it is highly polar and its more stable isomers (HNCO, HOCN, and HCNO) have all been identified in space. PMID:19895013
High-throughput ab-initio dilute solute diffusion database
Wu, Henry; Mayeshiba, Tam; Morgan, Dane
2016-01-01
We demonstrate automated generation of diffusion databases from high-throughput density functional theory (DFT) calculations. A total of more than 230 dilute solute diffusion systems in Mg, Al, Cu, Ni, Pd, and Pt host lattices have been determined using multi-frequency diffusion models. We apply a correction method for solute diffusion in alloys using experimental and simulated values of host self-diffusivity. We find good agreement with experimental solute diffusion data, obtaining a weighted activation barrier RMS error of 0.176 eV when excluding magnetic solutes in non-magnetic alloys. The compiled database is the largest collection of consistently calculated ab-initio solute diffusion data in the world. PMID:27434308
Ab initio engineering of materials with stacked hexagonal tin frameworks
Shao, Junping; Beaufils, Clément; Kolmogorov, Aleksey N.
2016-01-01
The group-IV tin has been hypothesized to possess intriguing electronic properties in an atom-thick hexagonal form. An attractive pathway of producing sizable 2D crystallites of tin is based on deintercalation of bulk compounds with suitable tin frameworks. Here, we have identified a new synthesizable metal distannide, NaSn2, with a 3D stacking of flat hexagonal layers and examined a known compound, BaSn2, with buckled hexagonal layers. Our ab initio results illustrate that despite being an exception to the 8-electron rule, NaSn2 should form under pressures easily achievable in multi-anvil cells and remain (meta)stable under ambient conditions. Based on calculated Z2 invariants, the predicted NaSn2 may display topologically non-trivial behavior and the known BaSn2 could be a strong topological insulator. PMID:27387140
Ab initio calculation of the shock Hugoniot of bulk silicon
NASA Astrophysics Data System (ADS)
Strickson, Oliver; Artacho, Emilio
2016-03-01
We describe how ab initio molecular dynamics can be used to determine the Hugoniot locus (states accessible by a shock wave) for materials with a number of stable phases, and with an approximate treatment of plasticity and yield, without having to simulate these phenomena directly. We consider the case of bulk silicon, with forces from density-functional theory, up to 70 GPa. The fact that shock waves can split into multiple waves due to phase transitions or yielding is taken into account here by specifying the strength of any preceding waves explicitly based on their yield strain. Points corresponding to uniaxial elastic compression along three crystal axes and a number of postshock phases are given, including a plastically yielded state, approximated by an isotropic stress configuration following an elastic wave of predetermined strength. The results compare well to existing experimental data for shocked silicon.
Interatomic Coulombic decay widths of helium trimer: Ab initio calculations
Kolorenč, Přemysl; Sisourat, Nicolas
2015-12-14
We report on an extensive study of interatomic Coulombic decay (ICD) widths in helium trimer computed using a fully ab initio method based on the Fano theory of resonances. Algebraic diagrammatic construction for one-particle Green’s function is utilized for the solution of the many-electron problem. An advanced and universal approach to partitioning of the configuration space into discrete states and continuum subspaces is described and employed. Total decay widths are presented for all ICD-active states of the trimer characterized by one-site ionization and additional excitation of an electron into the second shell. Selected partial decay widths are analyzed in detail, showing how three-body effects can qualitatively change the character of certain relaxation transitions. Previously unreported type of three-electron decay processes is identified in one class of the metastable states.
Interatomic Coulombic decay widths of helium trimer: Ab initio calculations.
Kolorenč, Přemysl; Sisourat, Nicolas
2015-12-14
We report on an extensive study of interatomic Coulombic decay (ICD) widths in helium trimer computed using a fully ab initio method based on the Fano theory of resonances. Algebraic diagrammatic construction for one-particle Green's function is utilized for the solution of the many-electron problem. An advanced and universal approach to partitioning of the configuration space into discrete states and continuum subspaces is described and employed. Total decay widths are presented for all ICD-active states of the trimer characterized by one-site ionization and additional excitation of an electron into the second shell. Selected partial decay widths are analyzed in detail, showing how three-body effects can qualitatively change the character of certain relaxation transitions. Previously unreported type of three-electron decay processes is identified in one class of the metastable states. PMID:26671378
Ab initio study of helium behavior in titanium tritides
Liang, J. H.; Dai, Yunya; Yang, Li; Peng, SM; Fan, K. M.; Long, XG; Zhou, X. S.; Zu, Xiaotao; Gao, Fei
2013-03-01
Ab initio calculations based on density functional theory have been performed to investigate the relative stability of titanium tritides and the helium behavior in stable titanium tritides. The results show that the β-phase TiT1.5 without two tritium along the [100] direction (TiT1.5[100]) is more stable than other possible structures. The stability of titanium tritides decrease with the increased generation of helium in TiT1.5[100]. In addition, helium generated by tritium decay prefers locating at a tetrahedral site, and favorably migrates between two neighbor vacant tetrahedral sites through an intermediate octahedral site in titanium tritides, with a migration energy of 0.23 eV. Furthermore, helium is easily accumulated on a (100) plane in β-phase TiT1.5[100].
The ab-initio density matrix renormalization group in practice
Olivares-Amaya, Roberto; Hu, Weifeng; Sharma, Sandeep; Yang, Jun; Chan, Garnet Kin-Lic; Nakatani, Naoki
2015-01-21
The ab-initio density matrix renormalization group (DMRG) is a tool that can be applied to a wide variety of interesting problems in quantum chemistry. Here, we examine the density matrix renormalization group from the vantage point of the quantum chemistry user. What kinds of problems is the DMRG well-suited to? What are the largest systems that can be treated at practical cost? What sort of accuracies can be obtained, and how do we reason about the computational difficulty in different molecules? By examining a diverse benchmark set of molecules: π-electron systems, benchmark main-group and transition metal dimers, and the Mn-oxo-salen and Fe-porphine organometallic compounds, we provide some answers to these questions, and show how the density matrix renormalization group is used in practice.
Ab initio engineering of materials with stacked hexagonal tin frameworks.
Shao, Junping; Beaufils, Clément; Kolmogorov, Aleksey N
2016-01-01
The group-IV tin has been hypothesized to possess intriguing electronic properties in an atom-thick hexagonal form. An attractive pathway of producing sizable 2D crystallites of tin is based on deintercalation of bulk compounds with suitable tin frameworks. Here, we have identified a new synthesizable metal distannide, NaSn2, with a 3D stacking of flat hexagonal layers and examined a known compound, BaSn2, with buckled hexagonal layers. Our ab initio results illustrate that despite being an exception to the 8-electron rule, NaSn2 should form under pressures easily achievable in multi-anvil cells and remain (meta)stable under ambient conditions. Based on calculated Z2 invariants, the predicted NaSn2 may display topologically non-trivial behavior and the known BaSn2 could be a strong topological insulator. PMID:27387140
Ab Initio Calculations of Excited Carrier Dynamics in Gallium Nitride
NASA Astrophysics Data System (ADS)
Jhalani, Vatsal; Bernardi, Marco
Bulk wurtzite GaN is the primary material for blue light-emission technology. The radiative processes in GaN are regulated by the dynamics of excited (or so-called ``hot'') carriers, through microscopic processes not yet completely understood. We present ab initio calculations of electron-phonon (e-ph) scattering rates for hot carriers in GaN. Our work combines density functional theory to compute the electronic states, and density functional perturbation theory to obtain the phonon dispersions and e-ph coupling matrix elements. These quantities are interpolated on fine Brillouin zone grids with maximally localized Wannier functions, to converge the e-ph scattering rates within 5 eV of the band edges. We resolve the contribution of the different phonon modes to the total scattering rate, and study the impact on the relaxation times of the long-range Fröhlich interaction due to the longitudinal-optical phonon modes.
Highly anisotropic thermal conductivity of arsenene: An ab initio study
NASA Astrophysics Data System (ADS)
Zeraati, Majid; Vaez Allaei, S. Mehdi; Abdolhosseini Sarsari, I.; Pourfath, Mahdi; Donadio, Davide
2016-02-01
Elemental two-dimensional (2D) materials exhibit intriguing heat transport and phononic properties. Here we have investigated the lattice thermal conductivity of newly proposed arsenene, the 2D honeycomb structure of arsenic, using ab initio calculations. Solving the Boltzmann transport equation for phonons, we predict a highly anisotropic thermal conductivity of 30.4 and 7.8 W/mK along the zigzag and armchair directions, respectively, at room temperature. Our calculations reveal that phonons with mean free paths between 20 nm and 1 μ m provide the main contribution to the large thermal conductivity in the zigzag direction; mean free paths of phonons contributing to heat transport in the armchair directions range between 20 and 100 nm. The obtained anisotropic thermal conductivity and feasibility of synthesis, in addition to high electron mobility reported elsewhere, make arsenene a promising material for nanoelectronic applications and thermal management.
Ab initio study of II-(VI)2 dichalcogenides.
Olsson, P; Vidal, J; Lincot, D
2011-10-12
The structural stabilities of the (Zn,Cd)(S,Se,Te)(2) dichalcogenides have been determined ab initio. These compounds are shown to be stable in the pyrite phase, in agreement with available experiments. Structural parameters for the ZnTe(2) pyrite semiconductor compound proposed here are presented. The opto-electronic properties of these dichalcogenide compounds have been calculated using quasiparticle GW theory. Bandgaps, band structures and effective masses are proposed as well as absorption coefficients and refraction indices. The compounds are all indirect semiconductors with very flat conduction band dispersion and high absorption coefficients. The work functions and surface properties are predicted. The Te and Se based compounds could be of interest as absorber materials in photovoltaic applications. PMID:21937783
Ab initio study of guanine damage by hydroxyl radical.
Chaban, Galina M; Wang, Dunyou; Huo, Winifred M
2015-01-15
Multiconfigurational ab initio methods are used in this study to examine two initial reactions that take place during the OH radical attack of the DNA base guanine: a ring opening reaction and a hydrogen transfer reaction. The same reactions are also studied in the presence of a single water molecule. The ring opening reaction has a moderate barrier height of ∼20-25 kcal/mol that is relatively insensitive to the presence of water. The barrier of the H-transfer reaction, on the other hand, is lowered from ∼50 to ∼22 kcal/mol when one water molecule is added, thus becoming comparable to the barrier height of the ring opening reaction. PMID:25517252
Vibrational and ab initio molecular dynamics studies of bradykinin
NASA Astrophysics Data System (ADS)
Święch, Dominika; Kubisiak, Piotr; Andrzejak, Marcin; Borowski, Piotr; Proniewicz, Edyta
2016-07-01
In this study, the comprehensive theoretical and experimental investigations of Raman (RS) and infrared absorption (IR) spectra of bradykinin (BK) are presented. The ab initio Born-Oppenheimer molecular dynamics (BOMD) calculations, in the presence of water molecules that form the first coordination sphere, were used for conformational analysis of the BK structure. Based on the Density Functional Theory (DFT) calculations at the B3LYP/6-311G(d) level the vibrational spectra were interpreted. The calculated frequencies were scaled by means of the effective scaling frequency factor (ESFF) method. The theoretical data, which confirm the compact structure of BK in the presence of the water molecules revealed the remarkable effect of the intermolecular hydrogen bonding on the BK structural properties.
Ab initio water pair potential with flexible monomers.
Jankowski, Piotr; Murdachaew, Garold; Bukowski, Robert; Akin-Ojo, Omololu; Leforestier, Claude; Szalewicz, Krzysztof
2015-03-26
A potential energy surface for the water dimer with explicit dependence on monomer coordinates is presented. The surface was fitted to a set of previously published interaction energies computed on a grid of over a quarter million points in the 12-dimensional configurational space using symmetry-adapted perturbation theory and coupled-cluster methods. The present fit removes small errors in published fits, and its accuracy is critically evaluated. The minimum and saddle-point structures of the potential surface were found to be very close to predictions from direct ab initio optimizations. The computed second virial coefficients agreed well with experimental values. At low temperatures, the effects of monomer flexibility in the virial coefficients were found to be much smaller than the quantum effects. PMID:25687650
Reactive Monte Carlo sampling with an ab initio potential
NASA Astrophysics Data System (ADS)
Leiding, Jeff; Coe, Joshua D.
2016-05-01
We present the first application of reactive Monte Carlo in a first-principles context. The algorithm samples in a modified NVT ensemble in which the volume, temperature, and total number of atoms of a given type are held fixed, but molecular composition is allowed to evolve through stochastic variation of chemical connectivity. We discuss general features of the method, as well as techniques needed to enhance the efficiency of Boltzmann sampling. Finally, we compare the results of simulation of NH3 to those of ab initio molecular dynamics (AIMD). We find that there are regions of state space for which RxMC sampling is much more efficient than AIMD due to the "rare-event" character of chemical reactions.
Efficient Ab initio Modeling of Random Multicomponent Alloys
NASA Astrophysics Data System (ADS)
Jiang, Chao; Uberuaga, Blas P.
2016-03-01
We present in this Letter a novel small set of ordered structures (SSOS) method that allows extremely efficient ab initio modeling of random multicomponent alloys. Using inverse II-III spinel oxides and equiatomic quinary bcc (so-called high entropy) alloys as examples, we demonstrate that a SSOS can achieve the same accuracy as a large supercell or a well-converged cluster expansion, but with significantly reduced computational cost. In particular, because of this efficiency, a large number of quinary alloy compositions can be quickly screened, leading to the identification of several new possible high-entropy alloy chemistries. The SSOS method developed here can be broadly useful for the rapid computational design of multicomponent materials, especially those with a large number of alloying elements, a challenging problem for other approaches.
Ab initio H2O in realistic hydrophilic confinement.
Allolio, Christoph; Klameth, Felix; Vogel, Michael; Sebastiani, Daniel
2014-12-15
A protocol for the ab initio construction of a realistic cylindrical pore in amorphous silica, serving as a geometric nanoscale confinement for liquids and solutions, is presented. Upon filling the pore with liquid water at different densities, the structure and dynamics of the liquid inside the confinement can be characterized. At high density, the pore introduces long-range oscillations into the water density profile, which makes the water structure unlike that of the bulk across the entire pore. The tetrahedral structure of water is also affected up to the second solvation shell of the pore wall. Furthermore, the effects of the confinement on hydrogen bonding and diffusion, resulting in a weakening and distortion of the water structure at the pore walls and a slowdown in diffusion, are characterized. PMID:25208765
Ab initio quantum dynamics using coupled-cluster.
Kvaal, Simen
2012-05-21
The curse of dimensionality (COD) limits the current state-of-the-art ab initio propagation methods for non-relativistic quantum mechanics to relatively few particles. For stationary structure calculations, the coupled-cluster (CC) method overcomes the COD in the sense that the method scales polynomially with the number of particles while still being size-consistent and extensive. We generalize the CC method to the time domain while allowing the single-particle functions to vary in an adaptive fashion as well, thereby creating a highly flexible, polynomially scaling approximation to the time-dependent Schrödinger equation. The method inherits size-consistency and extensivity from the CC method. The method is dubbed orbital-adaptive time-dependent coupled-cluster, and is a hierarchy of approximations to the now standard multi-configurational time-dependent Hartree method for fermions. A numerical experiment is also given. PMID:22612082
Efficient Ab initio Modeling of Random Multicomponent Alloys.
Jiang, Chao; Uberuaga, Blas P
2016-03-11
We present in this Letter a novel small set of ordered structures (SSOS) method that allows extremely efficient ab initio modeling of random multicomponent alloys. Using inverse II-III spinel oxides and equiatomic quinary bcc (so-called high entropy) alloys as examples, we demonstrate that a SSOS can achieve the same accuracy as a large supercell or a well-converged cluster expansion, but with significantly reduced computational cost. In particular, because of this efficiency, a large number of quinary alloy compositions can be quickly screened, leading to the identification of several new possible high-entropy alloy chemistries. The SSOS method developed here can be broadly useful for the rapid computational design of multicomponent materials, especially those with a large number of alloying elements, a challenging problem for other approaches. PMID:27015491
Ab initio engineering of materials with stacked hexagonal tin frameworks
NASA Astrophysics Data System (ADS)
Shao, Junping; Beaufils, Clément; Kolmogorov, Aleksey N.
2016-07-01
The group-IV tin has been hypothesized to possess intriguing electronic properties in an atom-thick hexagonal form. An attractive pathway of producing sizable 2D crystallites of tin is based on deintercalation of bulk compounds with suitable tin frameworks. Here, we have identified a new synthesizable metal distannide, NaSn2, with a 3D stacking of flat hexagonal layers and examined a known compound, BaSn2, with buckled hexagonal layers. Our ab initio results illustrate that despite being an exception to the 8-electron rule, NaSn2 should form under pressures easily achievable in multi-anvil cells and remain (meta)stable under ambient conditions. Based on calculated Z2 invariants, the predicted NaSn2 may display topologically non-trivial behavior and the known BaSn2 could be a strong topological insulator.
High-throughput ab-initio dilute solute diffusion database.
Wu, Henry; Mayeshiba, Tam; Morgan, Dane
2016-01-01
We demonstrate automated generation of diffusion databases from high-throughput density functional theory (DFT) calculations. A total of more than 230 dilute solute diffusion systems in Mg, Al, Cu, Ni, Pd, and Pt host lattices have been determined using multi-frequency diffusion models. We apply a correction method for solute diffusion in alloys using experimental and simulated values of host self-diffusivity. We find good agreement with experimental solute diffusion data, obtaining a weighted activation barrier RMS error of 0.176 eV when excluding magnetic solutes in non-magnetic alloys. The compiled database is the largest collection of consistently calculated ab-initio solute diffusion data in the world. PMID:27434308
Ab initio calculations of grain boundaries in bcc metals
NASA Astrophysics Data System (ADS)
Scheiber, Daniel; Pippan, Reinhard; Puschnig, Peter; Romaner, Lorenz
2016-03-01
In this study, we compute grain boundary (GB) properties for a large set of GBs in bcc transition metals with a special focus on W, Mo and Fe using ab initio density functional theory (DFT) and semi-empirical second nearest neighbour modified embedded atom method (2NN-MEAM) potentials. The GB properties include GB energies, surface energies, GB excess volume and work of separation, which we analyse and then compare to experimental data. We find that the used 2NN-MEAM potentials can predict general trends of GB properties, but do not always reproduce the GB ground state structure and energy found with DFT. In particular, our results explain the experimental finding that W and Mo prefer intergranular fracture, while other bcc metals prefer transgranular cleavage.
Transport coefficients in diamond from ab-initio calculations
NASA Astrophysics Data System (ADS)
Löfâs, Henrik; Grigoriev, Anton; Isberg, Jan; Ahuja, Rajeev
2013-03-01
By combining the Boltzmann transport equation with ab-initio electronic structure calculations, we obtain transport coefficients for boron-doped diamond. We find the temperature dependence of the resistivity and the hall coefficients in good agreement with experimental measurements. Doping in the samples is treated via the rigid band approximation and scattering is treated in the relaxation time approximation. In contrast to previous results, the acoustic phonon scattering is the dominating scattering mechanism for the considered doping range. At room temperature, we find the thermopower, S, in the range 1-1.6 mV/K and the power factor, S2σ, in the range 0.004-0.16 μW /cm K2.
Ab initio potential energy surface and rovibrational states of HBO
NASA Astrophysics Data System (ADS)
Ha, Tae-Kyu; Makarewicz, Jan
1999-01-01
The potential energy surface describing the large-amplitude motion of H around the BO core in the HBO molecule has been determined from ab initio calculations. This surface has been sampled by a set of 170 grid points from a two-dimensional space defined by the stretching and the bending coordinates of the H nucleus. At each grid point, the BO bond length has been optimized using the second-order Møller-Plesset perturbation theory with the basis set aug-cc-pVTZ. The surface has a local minimum for the linear as well as the bent configuration of HBO. A low energy barrier to the linear configuration BOH causes a large-amplitude motion and a strong rovibrational interaction in the molecule. Its rovibrational dynamics is different from the dynamics in bent or quasilinear triatomic molecules.
Ab Initio Screening of CO2-philic Groups.
Tian, Ziqi; Saito, Tomonori; Jiang, De-En
2015-04-23
Ab initio calculations were used to identify CO2-philic groups. Over 55 neutral molecules were screened for CO2 affinity via binding energetics. It is found that poly(ethylene oxide)s (PEO) oligomers with more than three repeating units are good CO2-binding groups, consistent with the high-performance of PEO-based materials for CO2/N2 separation. More interestingly, two triazole groups linked with a methylene chain are also excellent for CO2 binding with a favorable interaction of more than 28 kJ/mol, indicating that polymers or covalent-organic frameworks (COFs) with triazoles may be utilized for CO2 capture. This work provides a useful guide to introduce promising organic groups into polymeric membranes and COFs for CO2/N2 separation media. PMID:25825811
Ab initio calculation of thermodynamic properties of silicon
NASA Astrophysics Data System (ADS)
Wei, Siqing; Li, Changlin; Chou, M. Y.
1994-11-01
We present a fully ab initio calculation of the thermodynamic properties for silicon within the quasiharmonic approximation, making use of volume-dependent phonon frequencies obtained from pseudopotential local-density calculations. The temperature dependence of the thermal-expansion coefficient, specific heat (at constant volume), and other related quantities are studied. We confirm that the thermal-expansion coefficient behaves differently in three temperature regions: positive for temperature below 15 K, negative between 15 and 125 K, and positive again above 125 K. This finding agrees with experiment. The abnormal (negative) thermal-expansion coefficient at low temperatures is explained through a detailed study of mode Grüneisen parameters. Both specific-heat and thermal-expansion-coefficient values calculated are in excellent agreement with experiment up to a few hundred kelvin.
XMVB: a program for ab initio nonorthogonal valence bond computations.
Song, Lingchun; Mo, Yirong; Zhang, Qianer; Wu, Wei
2005-04-15
An ab initio nonorthogonal valence bond program, called XMVB, is described in this article. The XMVB package uses Heitler-London-Slater-Pauling (HLSP) functions as state functions, and calculations can be performed with either all independent state functions for a molecule or preferably a few selected important state functions. Both our proposed paired-permanent-determinant approach and conventional Slater determinant expansion algorithm are implemented for the evaluation of the Hamiltonian and overlap matrix elements among VB functions. XMVB contains the capabilities of valence bond self-consistent field (VBSCF), breathing orbital valence bond (BOVB), and valence bond configuration interaction (VBCI) computations. The VB orbitals, used to construct VB functions, can be defined flexibly in the calculations depending on particular applications and focused problems, and they may be strictly localized, delocalized, or bonded-distorted (semidelocalized). The parallel version of XMVB based on MPI (Message Passing Interface) is also available. PMID:15704237
Ab initio theory of NMR chemical shifts in solids
Louie, S.G. |
1997-12-31
A new formalism for ab initio calculation of the orbital magnetic susceptibility and the NMR chemical shifts in solids and liquids is presented. The approach can be applied to periodic systems such as crystals, surfaces or polymers, and with a supercell technique, to nonperiodic systems such as amorphous materials, liquids, or solids with defects. The formalism is based on the density functional theory in the local density approximation and makes use of a generalized f-sum rule to eliminate the divergent terms that plagued previous theories. Calculations have been successfully carried out for the diamagnetic susceptibility of a number of insulators and for the NMR chemical shifts of a variety of systems including free molecules, ionic crystals, hydrogen-bonded materials and amorphous carbon.
Ab initio Potential Energy Surface for H-H2
NASA Technical Reports Server (NTRS)
Partridge, Harry; Bauschlicher, Charles W., Jr.; Stallcop, James R.; Levin, Eugene
1993-01-01
Ab initio calculations employing large basis sets are performed to determine an accurate potential energy surface for H-H2 interactions for a broad range of separation distances. At large distances, the spherically averaged potential determined from the calculated energies agrees well with the corresponding results determined from dispersion coefficients; the van der Waals well depth is predicted to be 75 +/- (mu)E(sub h). Large basis sets have also been applied to reexamine the accuracy of theoretical repulsive potential energy surfaces. Multipolar expansions of the computed H-H2 potential energy surface are reported for four internuclear separation distances (1.2, 1.401, 1.449, and 1.7a(sub 0) of the hydrogen molecule. The differential elastic scattering cross section calculated from the present results is compared with the measurements from a crossed beam experiment.
Ab initio correlated calculations of rare-gas dimer quadrupoles
NASA Astrophysics Data System (ADS)
Donchev, Alexander G.
2007-10-01
This paper reports ab initio calculations of rare gas ( RG=Kr , Ar, Ne, and He) dimer quadrupoles at the second order of Møller-Plesset perturbation theory (MP2). The study reveals the crucial role of the dispersion contribution to the RG2 quadrupole in the neighborhood of the equilibrium dimer separation. The magnitude of the dispersion quadrupole is found to be much larger than that predicted by the approximate model of Hunt. As a result, the total MP2 quadrupole moment is significantly smaller than was assumed in virtually all previous related studies. An analytical model for the distance dependence of the RG2 quadrupole is proposed. The model is based on the effective-electron approach of Jansen, but replaces the original Gaussian approximation to the electron density in an RG atom by an exponential one. The role of the nonadditive contribution in RG3 quadrupoles is discussed.
Ab initio methods for nuclear properties - a computational physics approach
NASA Astrophysics Data System (ADS)
Maris, Pieter
2011-04-01
A microscopic theory for the structure and reactions of light nuclei poses formidable challenges for high-performance computing. Several ab-initio methods have now emerged that provide nearly exact solutions for some nuclear properties. The ab-initio no-core full configuration (NCFC) approach is based on basis space expansion methods and uses Slater determinants of single-nucleon basis functions to express the nuclear wave function. In this approach, the quantum many-particle problem becomes a large sparse matrix eigenvalue problem. The eigenvalues of this matrix give us the binding energies, and the corresponding eigenvectors the nuclear wave functions. These wave functions can be employed to evaluate experimental quantities. In order to reach numerical convergence for fundamental problems of interest, the matrix dimension often exceeds 1 billion, and the number of nonzero matrix elements may saturate available storage on present-day leadership class facilities. I discuss different strategies for distributing and solving this large sparse matrix on current multicore computer architectures, including methods to deal with with memory bottleneck. Several of these strategies have been implemented in the code MFDn, which is a parallel fortran code for nuclear structure calculations. I will show scaling behavior and compare the performance of the pure MPI version with the hybrid MPI/OpenMP code on Cray XT4 and XT5 platforms. For large core counts (typically 5,000 and above), the hybrid version is more efficient than pure MPI. With this code, we have been able to predict properties of the unstable nucleus 14F, which have since been confirmed by experiments. I will also give an overview of other recent results for nuclei in the A = 6 to 16 range with 2- and 3-body interactions. Supported in part by US DOE Grant DE-FC02-09ER41582.
4He Thermophysical Properties: New Ab Initio Calculations
Hurly, John J.; Mehl, James B.
2007-01-01
Since 2000, atomic physicists have reduced the uncertainty of the helium-helium “ab initio” potential; for example, from approximately 0.6 % to 0.1 % at 4 bohr, and from 0.8 % to 0.1 % at 5.6 bohr. These results led us to: (1) construct a new inter-atomic potential ϕ07, (2) recalculate values of the second virial coefficient, the viscosity, and the thermal conductivity of 4He from 1 K to 10,000 K, and (3), analyze the uncertainties of the thermophysical properties that propagate from the uncertainty of ϕ07 and from the Born-Oppenheimer approximation of the electron-nucleon quantum mechanical system. We correct minor errors in a previous publication [J. J. Hurly and M. R. Moldover, J. Res. Nat. Inst. Standards Technol. 105, 667 (2000)] and compare our results with selected data published after 2000. The ab initio results tabulated here can serve as standards for the measurement of thermophysical properties. PMID:27110456
Tuning Many-Body Interactions in Graphene: The Effects of Doping on Excitons and Carrier Lifetimes
NASA Astrophysics Data System (ADS)
Mak, Kin Fai; da Jornada, Felipe H.; He, Keliang; Deslippe, Jack; Petrone, Nicholas; Hone, James; Shan, Jie; Louie, Steven G.; Heinz, Tony F.
2014-05-01
The optical properties of graphene are strongly affected by electron-electron (e-e) and electron-hole (e-h) interactions. Here we tune these many-body interactions through varying the density of free charge carriers. Measurements from the infrared to the ultraviolet reveal significant changes in the optical conductivity of graphene for both electron and hole doping. The shift, broadening, and modification in shape of the saddle-point exciton resonance reflect strong screening of the many-body interactions by the carriers, as well as changes in quasiparticle lifetimes. Ab initio calculations by the GW Bethe-Salpeter equation method, which take into account the modification of both the repulsive e-e and the attractive e-h interactions, provide excellent agreement with experiment. Understanding the optical properties and high-energy carrier dynamics of graphene over a wide range of doping is crucial for both fundamental graphene physics and for emerging applications of graphene in photonics.
Operator evolution for ab initio theory of light nuclei
NASA Astrophysics Data System (ADS)
Schuster, Micah; Quaglioni, Sofia; Johnson, Calvin; Jurgenson, Eric; Navrátil, Petr
2014-09-01
The past two decades have seen a revolution in ab initio calculations of nuclear properties. One key element has been the development of a rigorous effective interaction theory, applying unitary transformations to soften the nuclear Hamiltonian and hence accelerate the convergence as a function of the model space size. For consistency, however, one ought to apply the same transformation to other operators when calculating transitions and mean values from the eigenstates of the renormalized Hamiltonian. Working in a translationally invariant harmonic oscillator basis for the two- and three-nucleon systems, we evolve the Hamiltonian, square radius, and total dipole strength operators by the similarity renormalization group (SRG). The inclusion of up to three-body matrix elements in the 4He nucleus all but completely restores the invariance of the expectation values under the transformation. We also consider a Gaussian operator with adjustable range; short ranges have the largest absolute renormalization when including two- and three-body induced terms, while at long ranges the induced three-body contribution takes on increased relative importance. The past two decades have seen a revolution in ab initio calculations of nuclear properties. One key element has been the development of a rigorous effective interaction theory, applying unitary transformations to soften the nuclear Hamiltonian and hence accelerate the convergence as a function of the model space size. For consistency, however, one ought to apply the same transformation to other operators when calculating transitions and mean values from the eigenstates of the renormalized Hamiltonian. Working in a translationally invariant harmonic oscillator basis for the two- and three-nucleon systems, we evolve the Hamiltonian, square radius, and total dipole strength operators by the similarity renormalization group (SRG). The inclusion of up to three-body matrix elements in the 4He nucleus all but completely restores
Heats of Segregation of BCC Binaries from ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2004-01-01
We compare dilute-limit heats of segregation for selected BCC transition metal binaries computed using ab initio and quantum approximate energy methods. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent LMTO-based parameters. Quantum approximate segregation energies are computed with and without atomistic relaxation, while the ab initio calculations are performed without relaxation. Results are discussed within the context of a segregation model driven by strain and bond-breaking effects. We compare our results with full-potential quantum calculations and with available experimental results.
Ab initio calculation and anharmonic force field of hypochlorous acid, HOCl
NASA Astrophysics Data System (ADS)
Halonen, L.; Ha, T.-K.
1988-03-01
Ab initio calculations on HOCl have been performed at the third-order Møller-Plesset perturbation theory level to determine the equilibrium structure and the anharmonic force field. An empirical anharmonic force field based on the ab initio results is obtained using available experimental vibration-rotation data. Four of the six harmonic and six of the ten cubic force constants have been determined experimentally, the remaining values being fixed at the ab initio values. A good fit to the experimental vibration-rotation data of four isotopic species is obtained.
Ab initio quasiparticle band structure of ABA and ABC-stacked graphene trilayers
NASA Astrophysics Data System (ADS)
Menezes, Marcos G.; Capaz, Rodrigo B.; Louie, Steven G.
2014-01-01
We obtain the quasiparticle band structure of ABA and ABC-stacked graphene trilayers through ab initio density-functional theory (DFT) and many-body quasiparticle calculations within the GW approximation. To interpret our results, we fit the DFT and GW π bands to a low-energy tight-binding model, which is found to reproduce very well the observed features near the K point. The values of the extracted hopping parameters are reported and compared with available theoretical and experimental data. For both stackings, the self-energy corrections lead to a renormalization of the Fermi velocity, an effect also observed in previous calculations on monolayer graphene. They also increase the separation between the higher-energy bands, which is proportional to the nearest-neighbor interlayer hopping parameter γ1. Both features are brought to closer agreement with experiment through the self-energy corrections. Finally, other effects, such as trigonal warping, electron-hole asymmetry, and energy gaps, are discussed in terms of the associated parameters.
Ab initio calculation of the electronic absorption spectrum of liquid water
NASA Astrophysics Data System (ADS)
Martiniano, Hugo F. M. C.; Galamba, Nuno; Cabral, Benedito J. Costa
2014-04-01
The electronic absorption spectrum of liquid water was investigated by coupling a one-body energy decomposition scheme to configurations generated by classical and Born-Oppenheimer Molecular Dynamics (BOMD). A Frenkel exciton Hamiltonian formalism was adopted and the excitation energies in the liquid phase were calculated with the equation of motion coupled cluster with single and double excitations method. Molecular dynamics configurations were generated by different approaches. Classical MD were carried out with the TIP4P-Ew and AMOEBA force fields. The BLYP and BLYP-D3 exchange-correlation functionals were used in BOMD. Theoretical and experimental results for the electronic absorption spectrum of liquid water are in good agreement. Emphasis is placed on the relationship between the structure of liquid water predicted by the different models and the electronic absorption spectrum. The theoretical gas to liquid phase blue-shift of the peak positions of the electronic absorption spectrum is in good agreement with experiment. The overall shift is determined by a competition between the O-H stretching of the water monomer in liquid water that leads to a red-shift and polarization effects that induce a blue-shift. The results illustrate the importance of coupling many-body energy decomposition schemes to molecular dynamics configurations to carry out ab initio calculations of the electronic properties in liquid phase.
Ab initio calculation of the electronic absorption spectrum of liquid water
Martiniano, Hugo F. M. C.; Galamba, Nuno; Cabral, Benedito J. Costa
2014-04-28
The electronic absorption spectrum of liquid water was investigated by coupling a one-body energy decomposition scheme to configurations generated by classical and Born-Oppenheimer Molecular Dynamics (BOMD). A Frenkel exciton Hamiltonian formalism was adopted and the excitation energies in the liquid phase were calculated with the equation of motion coupled cluster with single and double excitations method. Molecular dynamics configurations were generated by different approaches. Classical MD were carried out with the TIP4P-Ew and AMOEBA force fields. The BLYP and BLYP-D3 exchange-correlation functionals were used in BOMD. Theoretical and experimental results for the electronic absorption spectrum of liquid water are in good agreement. Emphasis is placed on the relationship between the structure of liquid water predicted by the different models and the electronic absorption spectrum. The theoretical gas to liquid phase blue-shift of the peak positions of the electronic absorption spectrum is in good agreement with experiment. The overall shift is determined by a competition between the O–H stretching of the water monomer in liquid water that leads to a red-shift and polarization effects that induce a blue-shift. The results illustrate the importance of coupling many-body energy decomposition schemes to molecular dynamics configurations to carry out ab initio calculations of the electronic properties in liquid phase.
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.
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.
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.
Ab initio quasiparticle bandstructure of ABA and ABC-stacked graphene trilayers
NASA Astrophysics Data System (ADS)
Menezes, Marcos; Capaz, Rodrigo; Louie, Steven
2013-03-01
We obtain the quasiparticle band structure of ABA and ABC-stacked graphene trilayers through ab initio density functional theory (DFT) and many-body quasiparticle calculations within the GW approximation. To interpret our results, we fit the DFT and GW π bands to a low energy tight-binding model, which is found to reproduce very well the observed features near the K point. The values of the extracted hopping parameters are reported and compared with available theoretical and experimental data. For both stackings, the quasiparticle corrections lead to a renormalization of the Fermi velocity, an effect also observed in previous calculations on monolayer graphene. They also increase the separation between the higher energy bands, which is proportional to the nearest neighbor interlayer hopping parameter γ1. Both features are brought to closer agreement with experiment through the quasiparticle corrections. Finally, other effects, such as trigonal warping, electron-hole assymetry and energy gaps are discussed in terms of the associated parameters. This work was supported by the Brazilian funding agencies: CAPES, CNPq, FAPERJ and INCT-Nanomateriais de Carbono. It was also supported by NSF grant No. DMR10-1006184 and U.S. DOE under Contract No. DE-AC02-05CH11231.
Towards an ab initio description of correlated materials
NASA Astrophysics Data System (ADS)
Yee, Chuck-Hou
Strongly-correlated materials are a rich playground for physical phenomena, exhibiting complex phase diagrams with many competing orders. Ab initio insights into materials combined with physical ideas provide the ability to identify the organizing principles driving the correlated electronic behavior and pursue first-principles design of new compounds. Realistic modeling of correlated materials is an active area of research, especially with the recent merger of density functional theory (DFT) with dynamical mean-field theory (DMFT). This thesis is structured in two parts. The first describes the methods and algorithmic developments which drive advances in DFT+DMFT. In Ch. 2 and 3, we provide an overview of the two foundational theories, DMFT and DFT. In the second half of Ch. 3, we describe some of the principles guiding the combination of the two theories to form DFT+DMFT. In Ch. 4, we describe the algorithm lying at the heart of modern DFT+DMFT implementations, the hybridization expansion formulation of continuous-time quantum monte carlo (CTQMC) for the general Anderson impurity problem, as well as a fast rejection algorithm for speeding-up the local trace evaluation. The final chapter in the methods section describes an algorithm for direct sampling of the partition function, and thus the free energy and entropy, of simple Anderson impurity models within CTQMC. The second part of the thesis is a collection of applications of our ab initio approach to key correlated materials. We first apply our method to plutonium binary alloys (Ch. 6), which when supplemented with slave-boson mean-field theory, allows us to understand the observed photoemission spectra. Ch. 7 describes the computation of spectra and optical conductivity for rare-earth nickelates grown as epitaxial thin films. In the final two chapters, we turn our attention to the high-temperature superconductors. In the first, we show that the charge-transfer energy is a key chemical variable which controls
Ab-initio calculations of low-lying excited states of water clusters (H2O)n, n = 2-6
NASA Astrophysics Data System (ADS)
Zvereva, Natalja A.; Ippolitov, Ivan I.
1997-03-01
In recent years there has been a great deal of interest in excited states of the water molecule with a number of experimental and theoretical studies reported. The umber of spectroscopic studies of complexes involving hydrogen bonding has greatly increased. However, most of this research has involved studies of the IR and Raman vibrational spectra of the molecule hydrogen bond formation strongly perturbs the vibrational transitions, but the changes in the electronic spectrum can also be observed. The interaction between water molecules is dominated by hydrogen bonding and the level of the ab initio theory is used because it accurately reproduces the major components of the interaction energy. Ab initio studies of cyclic water clusters (H2O)n, n equals 2 minus 6 and analysis of many-body interactions for ground states of these clusters has been done. The magnitudes of the two-through six-body energy terms and their contribution to the interaction energy of small ring water clusters has been computed at the Hartree-Fock (HF) and second-through fourth-order many-body perturbation (MP2, MP4) levels of theory. In this paper, we investigated the lowest S1 excited electronic states of the cyclic water clusters (H2O)n, n equals 2 minus 6 from the point of changes in the absorption spectra for many-body interactions and examined their trends as the cluster grew.
Aiming for benchmark accuracy with the many-body expansion.
Richard, Ryan M; Lao, Ka Un; Herbert, John M
2014-09-16
Conspectus The past 15 years have witnessed an explosion of activity in the field of fragment-based quantum chemistry, whereby ab initio electronic structure calculations are performed on very large systems by decomposing them into a large number of relatively small subsystem calculations and then reassembling the subsystem data in order to approximate supersystem properties. Most of these methods are based, at some level, on the so-called many-body (or "n-body") expansion, which ultimately requires calculations on monomers, dimers, ..., n-mers of fragments. To the extent that a low-order n-body expansion can reproduce supersystem properties, such methods replace an intractable supersystem calculation with a large number of easily distributable subsystem calculations. This holds great promise for performing, for example, "gold standard" CCSD(T) calculations on large molecules, clusters, and condensed-phase systems. The literature is awash in a litany of fragment-based methods, each with their own working equations and terminology, which presents a formidable language barrier to the uninitiated reader. We have sought to unify these methods under a common formalism, by means of a generalized many-body expansion that provides a universal energy formula encompassing not only traditional n-body cluster expansions but also methods designed for macromolecules, in which the supersystem is decomposed into overlapping fragments. This formalism allows various fragment-based methods to be systematically classified, primarily according to how the fragments are constructed and how higher-order n-body interactions are approximated. This classification furthermore suggests systematic ways to improve the accuracy. Whereas n-body approaches have been thoroughly tested at low levels of theory in small noncovalent clusters, we have begun to explore the efficacy of these methods for large systems, with the goal of reproducing benchmark-quality calculations, ideally meaning complete
Ab Initio Studies of Halogen and Nitrogen Oxide Species of Interest in Stratospheric Chemistry
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Langhoff, Stephen R. (Technical Monitor)
1995-01-01
The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of fluorine, chlorine, bromine and nitrogen oxide species will be demonstrated by presentation of some example studies. The ab initio results are shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the theoretical results are shown to fill in the gaps and to resolve experimental controversies. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of halogen oxide species will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of experimental studies.
Ab initio Bogoliubov coupled cluster theory for open-shell nuclei
Signoracci, Angelo J.; Duguet, Thomas; Hagen, Gaute; Jansen, G. R.
2015-06-29
Background: Ab initio many-body methods have been developed over the past 10 yr to address closed-shell nuclei up to mass A≈130 on the basis of realistic two- and three-nucleon interactions. A current frontier relates to the extension of those many-body methods to the description of open-shell nuclei. Several routes to address open-shell nuclei are currently under investigation, including ideas that exploit spontaneous symmetry breaking. Purpose: Singly open-shell nuclei can be efficiently described via the sole breaking of U(1) gauge symmetry associated with particle-number conservation as a way to account for their superfluid character. While this route was recently followed within the framework of self-consistent Green's function theory, the goal of the present work is to formulate a similar extension within the framework of coupled cluster theory. Methods: We formulate and apply Bogoliubov coupled cluster (BCC) theory, which consists of representing the exact ground-state wave function of the system as the exponential of a quasiparticle excitation cluster operator acting on a Bogoliubov reference state. Equations for the ground-state energy and the cluster amplitudes are derived at the singles and doubles level (BCCSD) both algebraically and diagrammatically. The formalism includes three-nucleon forces at the normal-ordered two-body level. The first BCC code is implemented in m scheme, which will permit the treatment of doubly open-shell nuclei via the further breaking of SU(2) symmetry associated with angular momentum conservation. Results: Proof-of-principle calculations in an N_{max}=6 spherical harmonic oscillator basis for ^{16,18}O and ^{18}Ne in the BCCD approximation are in good agreement with standard coupled cluster results with the same chiral two-nucleon interaction, while ^{20}O and ^{20}Mg display underbinding relative to experiment. The breaking of U(1) symmetry, monitored by computing the variance
Ab initio Bogoliubov coupled cluster theory for open-shell nuclei
Signoracci, Angelo J.; Duguet, Thomas; Hagen, Gaute; Jansen, G. R.
2015-06-29
Background: Ab initio many-body methods have been developed over the past 10 yr to address closed-shell nuclei up to mass A≈130 on the basis of realistic two- and three-nucleon interactions. A current frontier relates to the extension of those many-body methods to the description of open-shell nuclei. Several routes to address open-shell nuclei are currently under investigation, including ideas that exploit spontaneous symmetry breaking. Purpose: Singly open-shell nuclei can be efficiently described via the sole breaking of U(1) gauge symmetry associated with particle-number conservation as a way to account for their superfluid character. While this route was recently followed withinmore » the framework of self-consistent Green's function theory, the goal of the present work is to formulate a similar extension within the framework of coupled cluster theory. Methods: We formulate and apply Bogoliubov coupled cluster (BCC) theory, which consists of representing the exact ground-state wave function of the system as the exponential of a quasiparticle excitation cluster operator acting on a Bogoliubov reference state. Equations for the ground-state energy and the cluster amplitudes are derived at the singles and doubles level (BCCSD) both algebraically and diagrammatically. The formalism includes three-nucleon forces at the normal-ordered two-body level. The first BCC code is implemented in m scheme, which will permit the treatment of doubly open-shell nuclei via the further breaking of SU(2) symmetry associated with angular momentum conservation. Results: Proof-of-principle calculations in an Nmax=6 spherical harmonic oscillator basis for 16,18O and 18Ne in the BCCD approximation are in good agreement with standard coupled cluster results with the same chiral two-nucleon interaction, while 20O and 20Mg display underbinding relative to experiment. The breaking of U(1) symmetry, monitored by computing the variance associated with the particle-number operator, is
Zhang, Chunfang; Fu, Mingkai; Shen, Zhitao; Ma, Haitao; Bian, Wensheng
2014-06-21
A new global ab initio potential energy surface (called ZMB-a) for the 1(1)A' state of the C((1)D)+H2 reactive system has been constructed. This is based upon ab initio calculations using the internally contracted multireference configuration interaction approach with the aug-cc-pVQZ basis set, performed at about 6300 symmetry unique geometries. Accurate analytical fits are generated using many-body expansions with the permutationally invariant polynomials, except that the fit of the deep well region is taken from our previous fit. The ZMB-a surface is unique in the accurate description of the regions around conical intersections (CIs) and of van der Waals (vdW) interactions. The CIs between the 1(1)A' and 2(1)A' states cause two kinds of barriers on the ZMB-a surface: one is in the linear H-CH dissociation direction with a barrier height of 9.07 kcal/mol, which is much higher than those on the surfaces reported before; the other is in the C((1)D) collinearly attacking H2 direction with a barrier height of 12.39 kcal/mol. The ZMB-a surface basically reproduces our ab initio calculations in the vdW interaction regions, and supports a linear C-HH vdW complex in the entrance channel, and two vdW complexes in the exit channel, at linear CH-H and HC-H geometries, respectively. PMID:24952535
Lead-Chalcogenides Under Pressure: Ab-Initio Study
NASA Astrophysics Data System (ADS)
Gupta, Dinesh C.; Hamid, Idris
ab-initio calculations using fully relativistic pseudo-potential have been performed to investigate the high pressure phase transition, elastic and electronic properties of lead-chalcogenides including the less known lead polonium. The calculated ground state parameters, for the rock-salt structure show good agreement with the experimental data. The enthalpy calculations show that these materials undergo a first-order phase transition from rock-salt to CsCl structure at 19.4, 15.5, 11.5 and 7.3 GPa for PbS, PbSe, PbTe and PbPo, respectively. Present calculations successfully predicted the location of the band gap at L-point of Brillouin zone as well as the value of the band gap in every case at ambient pressure. It is observed that unlike other lead-chalcogenides, PbPo is semi-metal at ambient pressure. The pressure variation of the energy gap indicates that these materials metalized under high pressures. For this purpose, the electronic structure of these materials has also been computed in parent as well as in high pressure phase.
Exploring the free energy surface using ab initio molecular dynamics.
Samanta, Amit; Morales, Miguel A; Schwegler, Eric
2016-04-28
Efficient exploration of configuration space and identification of metastable structures in condensed phase systems are challenging from both computational and algorithmic perspectives. In this regard, schemes that utilize a set of pre-defined order parameters to sample the relevant parts of the configuration space [L. Maragliano and E. Vanden-Eijnden, Chem. Phys. Lett. 426, 168 (2006); J. B. Abrams and M. E. Tuckerman, J. Phys. Chem. B 112, 15742 (2008)] have proved useful. Here, we demonstrate how these order-parameter aided temperature accelerated sampling schemes can be used within the Born-Oppenheimer and the Car-Parrinello frameworks of ab initio molecular dynamics to efficiently and systematically explore free energy surfaces, and search for metastable states and reaction pathways. We have used these methods to identify the metastable structures and reaction pathways in SiO2 and Ti. In addition, we have used the string method [W. E, W. Ren, and E. Vanden-Eijnden, Phys. Rev. B 66, 052301 (2002); L. Maragliano et al., J. Chem. Phys. 125, 024106 (2006)] within the density functional theory to study the melting pathways in the high pressure cotunnite phase of SiO2 and the hexagonal closed packed to face centered cubic phase transition in Ti. PMID:27131525
Ab initio study of MoS2 nanotube bundles
NASA Astrophysics Data System (ADS)
Verstraete, Matthieu; Charlier, Jean-Christophe
2003-07-01
Recently, the synthesis of a new phase of MoS2I1/3 stoichiometry was reported [M. Remskar, A. Mrzel, Z. Skraba, A. Jesih, M. Ceh, J. Demšar, P. Stadelmann, F. Lévy, and D. Mihailovic, Science 292, 479 (2001)]. Electron microscope images and diffraction data were interpreted to indicate bundles of sub-nanometer-diameter single-wall MoS2 nanotubes. After experimental characterization, the structure was attributed to an assembly of “armchair” nanotubes with interstitial iodine. Using first-principles total-energy calculations, bundles of MoS2 nanotubes with different topologies and stoichiometries are investigated. All of the systems are strongly metallic. Configurations with “zigzag” structures are found to be more stable energetically than the “armchair” ones, though all of the structures have similar stabilities. After relaxation, there remain several candidates which give a lattice parameter in relative agreement with experiment. Further, spin-polarized calculations indicate that a structure with armchair tubes iodine atoms in their center acquires a very large spontaneous magnetic moment of 12μB, while the other structures are nonmagnetic. Our ab initio calculations show that in most of the other structures, the tubes are very strongly bound together, and that the compounds should be considered as a crystal, rather than as a bundle of tubes in the habitual sense.
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
Engineering Room-temperature Superconductors Via ab-initio Calculations
NASA Astrophysics Data System (ADS)
Gulian, Mamikon; Melkonyan, Gurgen; Gulian, Armen
The BCS, or bosonic model of superconductivity, as Little and Ginzburg have first argued, can bring in superconductivity at room temperatures in the case of high-enough frequency of bosonic mode. It was further elucidated by Kirzhnitset al., that the condition for existence of high-temperature superconductivity is closely related to negative values of the real part of the dielectric function at finite values of the reciprocal lattice vectors. In view of these findings, the task is to calculate the dielectric function for real materials. Then the poles of this function will indicate the existence of bosonic excitations which can serve as a "glue" for Cooper pairing, and if the frequency is high enough, and the dielectric matrix is simultaneously negative, this material is a good candidate for very high-Tc superconductivity. Thus, our approach is to elaborate a methodology of ab-initio calculation of the dielectric function of various materials, and then point out appropriate candidates. We used the powerful codes (TDDF with the DP package in conjunction with ABINIT) for computing dielectric responses at finite values of the wave vectors in the reciprocal lattice space. Though our report is concerned with the particular problem of superconductivity, the application range of the data processing methodology is much wider. The ability to compute the dielectric function of existing and still non-existing (though being predicted!) materials will have many more repercussions not only in fundamental sciences but also in technology and industry.
Ab initio simulations of pseudomorphic silicene and germanene bidimensional heterostructures
NASA Astrophysics Data System (ADS)
Debernardi, Alberto; Marchetti, Luigi
2016-06-01
Among the novel two-dimensional (2D) materials, silicene and germanene, which are two honeycomb crystal structures composed of a monolayer of Si and Ge, respectively, have attracted the attention of material scientists because they combine the advantages of the new 2D ultimate-scaled electronics with their compatibility with industrial processes presently based on Si and Ge. We envisage pseudomorphic lateral heterostructures based on ribbons of silicene and germanene, which are the 2D analogs of conventional 3D Si/Ge superlattices and quantum wells. In spite of the considerable lattice mismatch (˜4 % ) between free-standing silicene and germanene, our ab initio simulations predict that, considering striped 2D lateral heterostructures made by alternating silicene and germanene ribbons of constant width, the silicene/germanene junction remains pseudomorphic—i.e., it maintains lattice-matched edges—up to critical ribbon widths that can reach some tens of nanometers. Such critical widths are one order of magnitude larger than the critical thickness measured in 3D pseudomorphic Si/Ge heterostructures and the resolution of state-of-the-art lithography, thus enabling the possibility of lithography patterned silicene/germanene junctions. We computed how the strain produced by the pseudomorphic growth modifies the crystal structure and electronic bands of the ribbons, providing a mechanism for band-structure engineering. Our results pave the way for lithography patterned lateral heterostructures that can serve as the building blocks of novel 2D electronics.
Ab Initio Investigation of NH_3-O_2 Exciplex
NASA Astrophysics Data System (ADS)
Haupert, L. M.; Simpson, G.; Slipchenko, L. V.
2010-06-01
In their recent investigation of fluorescence from poly(amido amine) (PAMAM) dendrimers, Chu and Imae suggested an exciplex composed of tertiary amine and oxygen molecules might be responsible for fluorescence in PAMAM dendrimers. In this work, we present an ab initio investigation of the electronic structure of a possible ammonia-oxygen exciplex model system using equation-of-motion coupled cluster techniques. Geometry optimization of the triplet ground state produced a weakly bound state with an equilibrium separation of ˜ 3.5 Å, and an excited state geometry scan revealed a bound, excited triplet state with an equilibrium separation of 2.02 Å, consistent with results of earlier PM3 work by Juranic et al. The energy gap between the triplet ground state and first triplet excited state of the exciplex at 2.02 Å is 412.8 nm, lending support to the exciplex hypothesis. C.-C. Chu, and T. Imae, Macromol. Rapid. Commun., 30, 89-93 (2009). I. Juranic, H. S. Rzepa, and Y. MinYan, J. Chem. Soc. Perkin Trans., 2 (1990)
Amide N-oxides: an ab initio molecular orbital study
NASA Astrophysics Data System (ADS)
Greenberg, Arthur; DuBois, Thomas D.
2001-06-01
There are no known examples of amide N-oxides. The present study employs ab initio molecular orbital calculations at the 6-3G ∗ level to explore potential target molecules in this class. Bridgehead bicyclic lactams appear to be attractive targets for oxidation to form the corresponding N-oxides because they have reduced (or zero) amide resonance energy. The amide N-oxide linkage is predicted to have a ca. 9-10 kcal/mol rotational barrier due to eclipsing of nonbonded oxygen atoms in the transition state. The linkage has a nearly flat conformational ( ΦON-CO) profile in the range 120-240° and this suggests that a very sterically hindered acyclic amide N-oxide may be a practical synthetic target. The oxidation of strained amides is calculated to be highly exothermic if dimethyldioxirane is employed. This reagent is predicted to react appreciably exothermically with normal, stable amides such as N, N-dimethylacetamide, thus offering the potential for generating and studying such relatively unstable amide N-oxides at low temperatures.
FTIR, Raman spectra and ab initio calculations of 2-mercaptobenzothiazole.
Rai, Amareshwar K; Singh, Rachana; Singh, K N; Singh, V B
2006-02-01
FTIR and Raman spectra of a rubber vulcanization accelerator, 2-mercaptobenzothiazole (MBT), were recorded in the solid phase. The harmonic vibrational wavenumbers, for both the toutomeric forms of MBT, as well as for its dimeric complex, have been calculated, using ab initio RHF and density functional B3LYP methods invoking different basis sets upto RHF/6-31G** and B3LYP/6-31G** and the results were compared with the experimental values. Conformational studies have been also carried out regarding its toutomeric monomer forms and its dimer form. With all the basis sets the thione form of MBT (II) is predicted to be more stable than thiol form (I) and dimeric conformation (III) is predicted to be more stable with monomeric conformations (I) and (II). Vibrational assignments have been made, and it has been found that the calculated normal mode frequencies of dimeric conformation (III) are required for the analysis of IR and Raman bands of the MBT. The predicted shift in NH- stretching vibration towards the lower wave number side with the B3LYP/6-31G** calculations for the most stable dimer form (III), is in better agreement with experimental results. The intermolecular sulfur-nitrogen distance in N-H...S hydrogen bond was found to be 3.35 angstroms from these calculations, is also in agreement to the experimental value. PMID:16098806
Ab initio simulations on rutile-based titania nanowires
NASA Astrophysics Data System (ADS)
Zhukovskii, Yu F.; Evarestov, R. A.
2012-08-01
The rod symmetry groups for monoperiodic (1D) nanostructures have been applied for construction of models for bulk-like TiO2 nanowires (NWs) cut from a rutile-based 3D crystal along the chosen [001] and [110] directions of crystallographic axes. In this study, we have considered nanowires described by both the Ti-atom centered rotation axes as well as the hollow site centered axes passing through the interstitial positions between the Ti and O atoms closest to the axes. The most stable [001]-oriented TiO2 NWs with rhombic cross sections are found to display the energetically preferable {110} facets only while the nanowires with quasi-square sections across the [110] axis are formed by the alternating { 1bar 10 } and {001} facets. For simulations on rutile-based nanowires possessing different diameters for each NW type, we have performed large-scale ab initio Density Functional Theory (DFT) and hybrid DFT-Hartree Fock (DFT-HF) calculations with total geometry optimization within the Generalized Gradient Approximation (GGA) in the form of the Perdew-Becke-Ernzenhof (PBE) exchange-correlation functionals (PBE and PBE0, respectively), using the formalism of linear combination of localized atomic functions (LCAO). We have simulated both structural and electronic properties of TiO2 NWs depending both on orientation and position of symmetry axes as well as on diameter and morphology of nanowires.
Ab initio studies of niobium defects in uranium
Xiang, S; Huang, H; Hsiung, L
2007-06-01
Uranium (U), with the addition of small amount of niobium (Nb), is stainless. The Nb is fully miscible with the high temperature phase of U and tends to segregate upon cooling below 647 C. The starting point of segregation is the configuration of Nb substitutional or interstitial defects. Using density-functional-theory based ab initio calculations, the authors find that the formation energy of a single vacancy is 1.08 eV, that of Nb substitution is 0.59 eV, that of Nb interstitial at octahedral site is 1.58 eV, and that of Nb interstitial at tetrahedral site is 2.35 eV; all with reference to a reservoir of {gamma} phase U and pure Nb. The formation energy of Nb defects correlates with the local perturbation of electron distribution; higher formation energy to larger perturbation. Based on this study, Nb atoms thermodynamically prefer to occupy substitutional sites in {gamma} phase U, and they prefer to be in individual substitutional defects than clusters.
Melting curves of metals by ab initio calculations
NASA Astrophysics Data System (ADS)
Minakov, Dmitry; Levashov, Pavel
2015-06-01
In this work we used several ab initio approaches to reproduce melting curves and discussed their abilities, advantages and drawbacks. We used quasiharmonic appoximation and Lindemann criterion to build melting curves in wide region of pressures. This approach allows to calculate the total free energy of electrons and phonons, so it is possible to obtain all thermodynamic properties in the crystalline state. We also used quantum molecular dynamics simulations to investigate melting at various pressures. We explored the size-effect of the heat until it melts (HUM) method in detail. Special attention was paid to resolve the boundaries of the melting region on density. All calculations were performed for aluminum, copper and gold. Results were in good agreement with available experimental data. Also we studied the influence of electronic temperature on melting curves. It turned out that the melting temperature increased with the rise of electron temperature at normal density and had non-monotonic behavior at higher densities. This work is supported by the Ministry of Education and Science of the Russian Federation (Project No. 3.522.2014/K).
Ab initio simulation of gap discrete breathers in strained graphene
NASA Astrophysics Data System (ADS)
Lobzenko, I. P.; Chechin, G. M.; Bezuglova, G. S.; Baimova, Yu. A.; Korznikova, E. A.; Dmitriev, S. V.
2016-03-01
The methods of the density functional theory were used for the first time for the simulation of discrete breathers in graphene. It is demonstrated that breathers can exist with frequencies lying in the gap of the phonon spectrum, induced by uniaxial tension of a monolayer graphene sheet in the "zigzag" direction (axis X), polarized in the "armchair" direction (axis Y). The found gap breathers are highly localized dynamic objects, the core of which is formed by two adjacent carbon atoms located on the Y axis. The atoms surrounding the core vibrate at much lower amplitudes along both the axes ( X and Y). The dependence of the frequency of these breathers on amplitude is found, which shows a soft type of nonlinearity. No breathers of this type were detected in the gap induced by stretching along the Y axis. It is shown that the breather vibrations may be approximated by the Morse oscillators, the parameters of which are determined from ab initio calculations. The results are of fundamental importance, as molecular dynamics calculations based on empirical potentials cannot serve as a reliable proof of the existence of breathers in crystals.
Volumic omit maps in ab initio dual-space phasing.
Oszlányi, Gábor; Sütő, András
2016-07-01
Alternating-projection-type dual-space algorithms have a clear construction, but are susceptible to stagnation and, thus, inefficient for solving the phase problem ab initio. To improve this behaviour new omit maps are introduced, which are real-space perturbations applied periodically during the iteration process. The omit maps are called volumic, because they delete some predetermined subvolume of the unit cell without searching for atomic regions or analysing the electron density in any other way. The basic algorithms of positivity, histogram matching and low-density elimination are tested by their solution statistics. It is concluded that, while all these algorithms based on weak constraints are practically useless in their pure forms, appropriate volumic omit maps can transform them to practically useful methods. In addition, the efficiency of the already useful reflector-type charge-flipping algorithm can be further improved. It is important that these results are obtained by using non-sharpened structure factors and without any weighting scheme or reciprocal-space perturbation. The mathematical background of volumic omit maps and their expected applications are also discussed. PMID:27357850
Ab initio study of optical excitations in VO2
NASA Astrophysics Data System (ADS)
Coulter, John; Gali, Adam; Manousakis, Efstratios
2014-03-01
Motivated by recent experimental efforts to fabricate p-n junctions from transition metal oxides (TMOs) and a recent theoretical study claiming TMOs to be good absorbers and promising materials for efficient carrier multiplication, we study the optical properties of a prototypical TMO, the insulator M1 phase of vanadium dioxide (VO2), by ab initio methods. We applied the Bethe-Salpeter equations (BSE) to calculate the optical properties, starting from self-consistent GW quasi-particle energy levels and states. In contrast to expectations, the exciton binding energy obtained by BSE is in good agreement with the experiment. We find that the electron-electron interaction is very strong which makes this material promising for efficient carrier multiplication that might lead to an enhanced efficiency in photo-voltaics applications. To illustrate this more quantitatively, we calculated the impact ionization rate within the independent quasiparticle approximation, and find that the rate is significantly higher than silicon in the region of highest solar intensity, due to the strong multiple carrier excitations.
Ab initio studies of phoshorene island single electron transistor.
Ray, S J; Venkata Kamalakar, M; Chowdhury, R
2016-05-18
Phosphorene is a newly unveiled two-dimensional crystal with immense potential for nanoelectronic and optoelectronic applications. Its unique electronic structure and two dimensionality also present opportunities for single electron devices. Here we report the behaviour of a single electron transistor (SET) made of a phosphorene island, explored for the first time using ab initio calculations. We find that the band gap and the charging energy decrease monotonically with increasing layer numbers due to weak quantum confinement. When compared to two other novel 2D crystals such as graphene and MoS2, our investigation reveals larger adsorption energies of gas molecules on phosphorene, which indicates better a sensing ability. The calculated charge stability diagrams show distinct changes in the presence of an individual molecule which can be applied to detect the presence of different molecules with sensitivity at a single molecular level. The higher charging energies of the molecules within the SET display operational viability at room temperature, which is promising for possible ultra sensitive detection applications. PMID:27093536
Ab initio studies of phosphorene island single electron transistor
NASA Astrophysics Data System (ADS)
Ray, S. J.; Venkata Kamalakar, M.; Chowdhury, R.
2016-05-01
Phosphorene is a newly unveiled two-dimensional crystal with immense potential for nanoelectronic and optoelectronic applications. Its unique electronic structure and two dimensionality also present opportunities for single electron devices. Here we report the behaviour of a single electron transistor (SET) made of a phosphorene island, explored for the first time using ab initio calculations. We find that the band gap and the charging energy decrease monotonically with increasing layer numbers due to weak quantum confinement. When compared to two other novel 2D crystals such as graphene and MoS2, our investigation reveals larger adsorption energies of gas molecules on phosphorene, which indicates better a sensing ability. The calculated charge stability diagrams show distinct changes in the presence of an individual molecule which can be applied to detect the presence of different molecules with sensitivity at a single molecular level. The higher charging energies of the molecules within the SET display operational viability at room temperature, which is promising for possible ultra sensitive detection applications.
An Ab Initio Based Potential Energy Surface for Water
NASA Technical Reports Server (NTRS)
Partridge, Harry; Schwenke, David W.; Langhoff, Stephen R. (Technical Monitor)
1996-01-01
We report a new determination of the water potential energy surface. A high quality ab initio potential energy surface (PES) and dipole moment function of water have been computed. This PES is empirically adjusted to improve the agreement between the computed line positions and those from the HITRAN 92 data base. The adjustment is small, nonetheless including an estimate of core (oxygen 1s) electron correlation greatly improves the agreement with experiment. Of the 27,245 assigned transitions in the HITRAN 92 data base for H2(O-16), the overall root mean square (rms) deviation between the computed and observed line positions is 0.125/cm. However the deviations do not correspond to a normal distribution: 69% of the lines have errors less than 0.05/cm. Overall, the agreement between the line intensities computed in the present work and those contained in the data base is quite good, however there are a significant number of line strengths which differ greatly.
Ab initio description of the exotic unbound 7He nucleus
Baroni, Simone; Navratil, Petr; Quaglioni, Sofia
2013-01-11
In this study, the neutron-rich unbound 7He nucleus has been the subject of many experimental investigations. While the ground-state 3/2– resonance is well established, there is a controversy concerning the excited 1/2– resonance reported in some experiments as low lying and narrow (ER~1 MeV, Γ≤1 MeV) while in others as very broad and located at a higher energy. This issue cannot be addressed by ab initio theoretical calculations based on traditional bound-state methods. We introduce a new unified approach to nuclear bound and continuum states based on the coupling of the no-core shell model, a bound-state technique, with the no-coremore » shell model combined with the resonating-group method, a nuclear scattering technique. Our calculations describe the ground-state resonance in agreement with experiment and, at the same time, predict a broad 1/2– resonance above 2 MeV.« less
Ab initio investigation of grain boundary cohesion in Al alloys
NASA Astrophysics Data System (ADS)
Zhang, Shengjun; Kontsevoi, Oleg Y.; Freeman, A. J.; Olson, G. B.
2010-03-01
Strength and hardness of aluminum alloys can be substantially increased by alloying with Mg, Zn, Cu, Si, and other elements. The main drawback of Al alloys is their susceptibility to stress corrosion cracking, which is caused by alloying impurities segregated at grain boundaries. We investigated the embrittling and cohesion-enhancing effects of impurities on a σ5(012)[100] grain boundary in Al by means of the full-potential linearized augmented plane-wave (FLAPW) method within the framework of the Rice-Wang thermodynamic model and within the ab initio tensile test approach. We calculated segregation energies, analyzed local atomic configurations, electronic structures and spatial charge density distributions around segregated impurities, and identified the roles of atomic size and the bonding behavior of the impurity with the surrounding Al atoms. The results show that He, H and Na are strong embrittlers, Zn is a weak embrittler, while Sc, B, Cu and Mg are cohesion enhancers. We further evaluated the effect of co-alloying with two or more elements on grain boundary strength. This work provides a fundamental basis for the design of high strength Al alloys.
Ab Initio Simulation of the Photoelectron Spectrum for Methoxy Radical
NASA Astrophysics Data System (ADS)
Cheng, Lan; Weichman, Marissa L.; Kim, Jongjin B.; Ichino, Takatoshi; Neumark, Daniel; Stanton, John F.
2015-06-01
A theoretical simulation of the photoelectron spectrum for the ground state of methoxy radical is reported based on the quasidiabatic model Hamiltonian originally proposed by Köppel, Domcke, and Cederbaum. The parameters in the model Hamiltonian have been obtained from ab initio coupled-cluster calculations. The linear and quadratic force constants have been calculated using equation-of-motion coupled-cluster ionization potential method with the singles, doubles, and triples (EOMIP-CCSDT) truncation scheme together with atomic natural orbital basis sets of triple-zeta quality (ANO1). The cubic and quartic force constants have been obtained from EOMIP-CCSD calculations with ANO basis sets of double-zeta quality (ANO0), and the spin-orbit coupling constant has been computed at the EOMIP-CCSD/pCVTZ level. The nuclear Schroedinger equation has been solved using the Lanzcos algorithm to obtain vibronic energy levels as well as the corresponding intensities. The simulated spectrum compares favorably with the recent high-resolution slow electron velocity-map imaging experiment for vibronic levels up to 2000 cm-1.
Ab initio simulations of peptide-mineral interactions
NASA Astrophysics Data System (ADS)
Hug, Susanna; Hunter, Graeme K.; Goldberg, Harvey; Karttunen, Mikko
We performed Car-Parrinello Molecular Dynamics (CPMD) simulations of two amino acids, aspartic acid (Asp) and phophoserine (pSer), on a calcium oxalate monohydrate (COM) surface as a model of the interactions of phosphoproteins with biominerals. In our earlier work using in vitro experiments and classical Molecular Dynamics (MD) simulations we have demonstrated the importance of phosphorylation of serine on the interactions of osteopontin (OPN) with COM. We used configurations from our previous classical MD simulations as a starting point for the ab initio simulations. In the case of Asp we found that the α-carboxyl and amine groups form temporary close contacts with the surface. For the dipeptide Asp-pSer the carboxyl groups form permanent close contacts with the surface and the distances of its other functional groups do not vary much. We show how the interaction of carboxyl groups with COM crystal is established and confirm the importance of phosphorylation in mediating the interactions between COM surfaces and OPN.
Ab initio study of the phenol-water cation radical
NASA Astrophysics Data System (ADS)
Hobza, Pavel; Burcl, Rudolf; Špirko, Vladimír; Dopfer, Otto; Müller-Dethlefs, Klaus; Schlag, Edward W.
1994-07-01
The phenol-water cation radical has been investigated by ab initio theory using the spin-restricted open-shell Hartree-Fock and spin-restricted open-shell second-order Møller-Plesset theories with 3-21G*(O) and 6-31G* basis sets. The full geometrical optimization was performed for several hydrogen-bonded structures and one hemibonded structure. Clearly, the most stable structure has been found for Cs symmetry with the linear hydrogen bond between the proton of the OH group of the phenol cation radical and the oxygen of the water, and the water hydrogens pointing away from the phenyl ring. For this structure harmonic (and for some intermolecular modes anharmonic) vibrational frequencies have been computed for various isotopic complexes. The computed shifts of phenol-localized intramolecular modes on complexation and on deuteration as well as the calculated intermolecular frequencies of the different isotopic complexes allow for an assignment of vibrational frequencies observed in the experimental zero-kinetic-energy (ZEKE) photoelectron spectra. Five out of a possible six intermolecular vibrations and several intramolecular modes have been assigned, including the 18b vibration which shows a strong blue shift in frequency upon complexation. Structure and properties of the phenol-water cation radical are compared with those of the corresponding neutral complex.
Ab initio study of the bonding in diatomic nickel
Noell, J.O.; Newton, M.D.; Hay, P.J.; Martin, R.L.; Bobrowicz, F.W.
1980-09-01
Hartree--Fock, GVB, and configuration interaction calculations were performed for diatomic nickel using an ab initio effective core potential. A basis set specifically optimized for the /sup 3/D state of atomic nickel is found to be far superior to the more common basis obtained from the /sup 3/F atomic state. Correlation effects are found to be significant in determining the bond energy. In particular, the two electrons of the s--s bond must be appropriately correlated. In addition, correlation effects which one would interpret as being principally intra-atomic in character are found to have a marked effect on the molecular properties. The theoretically predicted bond dissociation energy (D/sub e/) of 43.4 kcal/mol is significantly lower than the experimental estimate of 55 +- 5 kcal/mol. However, molecular partition functions calculated using the present results indicate that the experimental value should be revised downward to a value of approx.46 +- 5 kcal/mol, in good agreement with our calculations. An interatomic distance of 4.27 bohr is computed and compared with experimental estimates. Spectroscopic parameters for dipole-allowed transitions from the ground state were determined from SCF and GVB calculations and discussed in relation to the experimentally observed visible and ultraviolet spectra attributed to Ni/sub 2/.
An efficient approach to ab initio Monte Carlo simulation
Leiding, Jeff; Coe, Joshua D.
2014-01-21
We present a Nested Markov chain Monte Carlo (NMC) scheme for building equilibrium averages based on accurate potentials such as density functional theory. Metropolis sampling of a reference system, defined by an inexpensive but approximate potential, was used to substantially decorrelate configurations at which the potential of interest was evaluated, thereby dramatically reducing the number needed to build ensemble averages at a given level of precision. The efficiency of this procedure was maximized on-the-fly through variation of the reference system thermodynamic state (characterized here by its inverse temperature β{sup 0}), which was otherwise unconstrained. Local density approximation results are presented for shocked states of argon at pressures from 4 to 60 GPa, where—depending on the quality of the reference system potential—acceptance probabilities were enhanced by factors of 1.2–28 relative to unoptimized NMC. The optimization procedure compensated strongly for reference potential shortcomings, as evidenced by significantly higher speedups when using a reference potential of lower quality. The efficiency of optimized NMC is shown to be competitive with that of standard ab initio molecular dynamics in the canonical ensemble.
Ab initio Raman spectroscopy of water under extreme conditions
NASA Astrophysics Data System (ADS)
Rozsa, Viktor; Pan, Ding; Wan, Quan; Galli, Giulia
Water exhibits one of the most complex phase diagrams of any binary compound. Despite extensive studies, the melting lines of high-pressure ice phases remain very controversial, with reports differing by hundreds of Kelvin. The boundary between ice VII and liquid phase is particularly disputed, with recent work exploring plasticity and amorphization mediating the transition. Raman measurements are often used to fingerprint melting, yet their interpretation is difficult without atomistic modeling. Here, we report a study of high P/T water where we computed Raman spectra using a method combining ab initio molecular dynamics and density functional perturbation theory, as implemented in the Qbox code. Spectra were computed for the liquid at 10 and 20 GPa, both at 1000 K, and for solid ice VII (20 GPa, 500 K). Decomposing the spectra into inter and intra molecular contributions provided insight into the dynamics of the hydrogen-bonded network at extreme conditions. The relevance of our simulation results for models of water in Earth, Uranus, and Neptune will be discussed, and an interpretation of existing experiments at high pressure will be presented.
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.
Accurate ab initio vibrational energies of methyl chloride.
Owens, Alec; Yurchenko, Sergei N; Yachmenev, Andrey; Tennyson, Jonathan; Thiel, Walter
2015-06-28
Two new nine-dimensional potential energy surfaces (PESs) have been generated using high-level ab initio theory for the two main isotopologues of methyl chloride, CH3 (35)Cl and CH3 (37)Cl. The respective PESs, CBS-35( HL), and CBS-37( HL), are based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set (CBS) limit, and incorporate a range of higher-level (HL) additive energy corrections to account for core-valence electron correlation, higher-order coupled cluster terms, scalar relativistic effects, and diagonal Born-Oppenheimer corrections. Variational calculations of the vibrational energy levels were performed using the computer program TROVE, whose functionality has been extended to handle molecules of the form XY 3Z. Fully converged energies were obtained by means of a complete vibrational basis set extrapolation. The CBS-35( HL) and CBS-37( HL) PESs reproduce the fundamental term values with root-mean-square errors of 0.75 and 1.00 cm(-1), respectively. An analysis of the combined effect of the HL corrections and CBS extrapolation on the vibrational wavenumbers indicates that both are needed to compute accurate theoretical results for methyl chloride. We believe that it would be extremely challenging to go beyond the accuracy currently achieved for CH3Cl without empirical refinement of the respective PESs. PMID:26133427
Accurate ab initio vibrational energies of methyl chloride
NASA Astrophysics Data System (ADS)
Owens, Alec; Yurchenko, Sergei N.; Yachmenev, Andrey; Tennyson, Jonathan; Thiel, Walter
2015-06-01
Two new nine-dimensional potential energy surfaces (PESs) have been generated using high-level ab initio theory for the two main isotopologues of methyl chloride, CH335Cl and CH337Cl. The respective PESs, CBS-35 HL, and CBS-37 HL, are based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set (CBS) limit, and incorporate a range of higher-level (HL) additive energy corrections to account for core-valence electron correlation, higher-order coupled cluster terms, scalar relativistic effects, and diagonal Born-Oppenheimer corrections. Variational calculations of the vibrational energy levels were performed using the computer program TROVE, whose functionality has been extended to handle molecules of the form XY 3Z. Fully converged energies were obtained by means of a complete vibrational basis set extrapolation. The CBS-35 HL and CBS-37 HL PESs reproduce the fundamental term values with root-mean-square errors of 0.75 and 1.00 cm-1, respectively. An analysis of the combined effect of the HL corrections and CBS extrapolation on the vibrational wavenumbers indicates that both are needed to compute accurate theoretical results for methyl chloride. We believe that it would be extremely challenging to go beyond the accuracy currently achieved for CH3Cl without empirical refinement of the respective PESs.
Ab initio description of p-shell hypernuclei.
Wirth, Roland; Gazda, Daniel; Navrátil, Petr; Calci, Angelo; Langhammer, Joachim; Roth, Robert
2014-11-01
We present the first ab initio calculations for p-shell single-Λ hypernuclei. For the solution of the many-baryon problem, we develop two variants of the no-core shell model with explicit Λ and Σ(+),Σ(0),Σ(-) hyperons including Λ-Σ conversion, optionally supplemented by a similarity renormalization group transformation to accelerate model-space convergence. In addition to state-of-the-art chiral two- and three-nucleon interactions, we use leading-order chiral hyperon-nucleon interactions and a recent meson-exchange hyperon-nucleon interaction. We validate the approach for s-shell hypernuclei and apply it to p-shell hypernuclei, in particular to (Λ)(7)Li, (Λ)(9)Be, and (Λ)(13)C. We show that the chiral hyperon-nucleon interactions provide ground-state and excitation energies that generally agree with experiment within the cutoff dependence. At the same time we demonstrate that hypernuclear spectroscopy provides tight constraints on the hyperon-nucleon interactions. PMID:25415901
Ab initio SCF calculations on hydrogen bonded cresol isomers
NASA Astrophysics Data System (ADS)
Pohl, M.; Kleinermanns, K.
1988-12-01
Ab initio GAUSSIAN 80 calculations with two different basis sets (STO-3G and 4 31 G*) were performed on hydrogen bonded cresol isomers for comparison with experimental data from free jet fluorescence excitation spectroscopy. For m-cresol, the calculated barriers for hindered internal rotation of the OH-group and the CH3-group are in good agreement with experiment. The calculations show the trans-linear configuration of p-cresol· B-clusters ( B = H2O, CH3OH) to be more stable than the all-planar configuration. This agrees with CI calculations and microwave spectroscopic investigations of the water dimer. Calculations of both the intermolecular stretch and bend frequencies of p-cresol· B-clusters show little dependence on the all-planar or trans-linear configuration but a strong dependence on the choice of the basis set. With the minimal basis set STO-3G, the vibrational energies are generally too high. The agreement between the calculated vibrational frequencies from the 4 31 G* basis set and the experimental values is fair.
Predicting lattice thermal conductivity with help from ab initio methods
NASA Astrophysics Data System (ADS)
Broido, David
2015-03-01
The lattice thermal conductivity is a fundamental transport parameter that determines the utility a material for specific thermal management applications. Materials with low thermal conductivity find applicability in thermoelectric cooling and energy harvesting. High thermal conductivity materials are urgently needed to help address the ever-growing heat dissipation problem in microelectronic devices. Predictive computational approaches can provide critical guidance in the search and development of new materials for such applications. Ab initio methods for calculating lattice thermal conductivity have demonstrated predictive capability, but while they are becoming increasingly efficient, they are still computationally expensive particularly for complex crystals with large unit cells . In this talk, I will review our work on first principles phonon transport for which the intrinsic lattice thermal conductivity is limited only by phonon-phonon scattering arising from anharmonicity. I will examine use of the phase space for anharmonic phonon scattering and the Grüneisen parameters as measures of the thermal conductivities for a range of materials and compare these to the widely used guidelines stemming from the theory of Liebfried and Schölmann. This research was supported primarily by the NSF under Grant CBET-1402949, and by the S3TEC, an Energy Frontier Research Center funded by the US DOE, office of Basic Energy Sciences under Award No. DE-SC0001299.
Ab initio calculation of double ionization of atoms
Serov, V. V.
2013-02-15
The Solov'ev-Vinitsky method was used to perform an ab initio calculation of the triple-differential cross section for the double single-photon photoionization of helium for the case of equal emitted-electron energies. A Gaussian width {gamma} describing angular electron-electron correlations at the total electron energy E taking values in range between 0.1 and 100 eV was obtained for this cross section. The results agree with available experimental data, but they raise a doubt as to whether the well-known Wannier law {gamma} {proportional_to} E{sup 1/4} is applicable at experimentally accessible energies. The Gaussian width {gamma} was investigated as a function of the total emitted-electron energy for targets that have a strongly asymmetric configuration of the initial state-specifically, a negative atomic-hydrogen ion H{sup -} and heliumin the 1s2s{sup 1}S and 1s3s{sup 1}S excited states. It was found that this function, {gamma}(E), had a maximum at low energies. It was also shown that, at low energies, the dependence of the double-differential cross section on the angle between the emitted-electron momenta for the targets indicated above differed substantially from the Gaussian dependence, featuring maxima whose number was equal to the number of radial nodes in the initial state. This opens new possibilities for a qualitative analysis of the electron structure of targets.
Ab initio calculation of infrared intensities for hydrogen peroxide
NASA Technical Reports Server (NTRS)
Rogers, J. D.; Hillman, J. J.
1982-01-01
Results of an ab initio SCF quantum mechanical study are used to derive estimates for the infrared intensities of the fundamental vibrations of hydrogen peroxide. Atomic polar tensors (APTs) were calculated on the basis of a 4-31G basis set, and used to derive absolute intensities for the vibrational transitions. Comparison of the APTs calculated for H2O2 with those previously obtained for H2O and CH3OH, and of the absolute intensities derived from the H2O2 APTs with those derived from APTs transferred from H2O and CH3OH, reveals the sets of values to differ by no more than a factor of two, supporting the validity of the theoretical calculation. Values of the infrared intensities obtained correspond to A1 = 14.5 km/mol, A2 = 0.91 km/mol, A3 = 0.058 km/mol, A4 = 123 km/mol, A5 = 46.2 km/mol, and A6 = 101 km/mol. Charge, charge flux and overlap contributions to the dipole moment derivatives are also computed.
Ab initio molecular dynamics calculations of ion hydration free energies
Leung, Kevin; Rempe, Susan B.; Lilienfeld, O. Anatole von
2009-05-28
We apply ab initio molecular dynamics (AIMD) methods in conjunction with the thermodynamic integration or '{lambda}-path' technique to compute the intrinsic hydration free energies of Li{sup +}, Cl{sup -}, and Ag{sup +} ions. Using the Perdew-Burke-Ernzerhof functional, adapting methods developed for classical force field applications, and with consistent assumptions about surface potential ({phi}) contributions, we obtain absolute AIMD hydration free energies ({Delta}G{sub hyd}) within a few kcal/mol, or better than 4%, of Tissandier et al.'s [J. Phys. Chem. A 102, 7787 (1998)] experimental values augmented with the SPC/E water model {phi} predictions. The sums of Li{sup +}/Cl{sup -} and Ag{sup +}/Cl{sup -} AIMD {Delta}G{sub hyd}, which are not affected by surface potentials, are within 2.6% and 1.2 % of experimental values, respectively. We also report the free energy changes associated with the transition metal ion redox reaction Ag{sup +}+Ni{sup +}{yields}Ag+Ni{sup 2+} in water. The predictions for this reaction suggest that existing estimates of {Delta}G{sub hyd} for unstable radiolysis intermediates such as Ni{sup +} may need to be extensively revised.
Ab initio study of ice catalyzation of HOCl + HCl reaction
Zhou, Y.F.; Liu, C.B.
2000-06-15
The observations by Farman et al. revealed remarkable depletions in the total atmospheric ozone content in Antarctica. The observed total ozone decreased smoothing during the spring season from about 1975. Satellite observations have proved Antarctic ozone depletions over a very extended region, in general agreement with the local ground-based data of Farman et al. It was suggested that heterogeneous reactions occurring on particles in polar stratospheric clouds (PSCs) play a central role in the depletion of stratospheric ozone. Experiments proved that the reaction of HOCl + HCl was very slow in the gas phase, but on ice surface it was rapid. In this work the ice catalysis of HOCl + HCl reaction was investigated by using ab initio molecular orbital theory. The authors applied the Hartree-Fock self-consistent field and the second-order Moeller-Plesset perturbation theory with the basis sets of 6-31G* to the model system. The complexes and transition state were obtained along the reaction with and without the presence of ice surface. By comparing the results, a possible catalyzation mechanism of ice on the reaction is proposed.
Exploring the free energy surface using ab initio molecular dynamics
NASA Astrophysics Data System (ADS)
Samanta, Amit; Morales, Miguel A.; Schwegler, Eric
2016-04-01
Efficient exploration of configuration space and identification of metastable structures in condensed phase systems are challenging from both computational and algorithmic perspectives. In this regard, schemes that utilize a set of pre-defined order parameters to sample the relevant parts of the configuration space [L. Maragliano and E. Vanden-Eijnden, Chem. Phys. Lett. 426, 168 (2006); J. B. Abrams and M. E. Tuckerman, J. Phys. Chem. B 112, 15742 (2008)] have proved useful. Here, we demonstrate how these order-parameter aided temperature accelerated sampling schemes can be used within the Born-Oppenheimer and the Car-Parrinello frameworks of ab initio molecular dynamics to efficiently and systematically explore free energy surfaces, and search for metastable states and reaction pathways. We have used these methods to identify the metastable structures and reaction pathways in SiO2 and Ti. In addition, we have used the string method [W. E, W. Ren, and E. Vanden-Eijnden, Phys. Rev. B 66, 052301 (2002); L. Maragliano et al., J. Chem. Phys. 125, 024106 (2006)] within the density functional theory to study the melting pathways in the high pressure cotunnite phase of SiO2 and the hexagonal closed packed to face centered cubic phase transition in Ti.
Accurate ab initio vibrational energies of methyl chloride
Owens, Alec; Yurchenko, Sergei N.; Yachmenev, Andrey; Tennyson, Jonathan; Thiel, Walter
2015-06-28
Two new nine-dimensional potential energy surfaces (PESs) have been generated using high-level ab initio theory for the two main isotopologues of methyl chloride, CH{sub 3}{sup 35}Cl and CH{sub 3}{sup 37}Cl. The respective PESs, CBS-35{sup HL}, and CBS-37{sup HL}, are based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set (CBS) limit, and incorporate a range of higher-level (HL) additive energy corrections to account for core-valence electron correlation, higher-order coupled cluster terms, scalar relativistic effects, and diagonal Born-Oppenheimer corrections. Variational calculations of the vibrational energy levels were performed using the computer program TROVE, whose functionality has been extended to handle molecules of the form XY {sub 3}Z. Fully converged energies were obtained by means of a complete vibrational basis set extrapolation. The CBS-35{sup HL} and CBS-37{sup HL} PESs reproduce the fundamental term values with root-mean-square errors of 0.75 and 1.00 cm{sup −1}, respectively. An analysis of the combined effect of the HL corrections and CBS extrapolation on the vibrational wavenumbers indicates that both are needed to compute accurate theoretical results for methyl chloride. We believe that it would be extremely challenging to go beyond the accuracy currently achieved for CH{sub 3}Cl without empirical refinement of the respective PESs.
Ab Initio Potential Energy Surface for H-H2
NASA Technical Reports Server (NTRS)
Patridge, Harry; Bauschlicher, Charles W., Jr.; Stallcop, James R.; Levin, Eugene
1993-01-01
Ab initio calculations employing large basis sets are performed to determine an accurate potential energy surface for H-H2 interactions for a broad range of separation distances. At large distances, the spherically averaged potential determined from the calculated energies agrees well with the corresponding results determined from dispersion coefficients; the van der Waals well depth is predicted to be 75 +/- 3 micro E(h). Large basis sets have also been applied to reexamine the accuracy of theoretical repulsive potential energy surfaces (25-70 kcal/mol above the H-H2 asymptote) at small interatomic separations; the Boothroyd, Keogh, Martin, and Peterson (BKMP) potential energy surface is found to agree with results of the present calculations within the expected uncertainty (+/- 1 kcal/mol) of the fit. Multipolar expansions of the computed H-H2 potential energy surface are reported for four internuclear separation distances (1.2, 1.401, 1.449, and 1.7a(0)) of the hydrogen molecule. The differential elastic scattering cross section calculated from the present results is compared with the measurements from a crossed beam experiment.
^{4}He+n+n continuum within an ab initio framework
Romero-Redondo, Carolina; Quaglioni, Sofia; Navratil, Petr; Hupin, Guillaume
2014-07-16
In this study, the low-lying continuum spectrum of the ^{6}He nucleus is investigated for the first time within an ab initio framework that encompasses the ^{4}He+n+n three-cluster dynamics characterizing its lowest decay channel. This is achieved through an extension of the no-core shell model combined with the resonating-group method, in which energy-independent nonlocal interactions among three nuclear fragments can be calculated microscopically, starting from realistic nucleon-nucleon interactions and consistent ab initio many-body wave functions of the clusters. The three-cluster Schrödinger equation is solved with three-body scattering boundary conditions by means of the hyperspherical-harmonics method on a Lagrange mesh. Using a soft similarity-renormalization-group evolved chiral nucleon-nucleon potential, we find the known J^{π} = 2^{+} resonance as well as a result consistent with a new low-lying second 2^{+} resonance recently observed at GANIL at ~2.6 MeV above the He6 ground state. We also find resonances in the 2^{–}, 1^{+}, and 0^{–} channels, while no low-lying resonances are present in the 0^{+} and 1^{–} channels.
NASA Astrophysics Data System (ADS)
Kritayakornupong, Chinapong; Hannongbua, Supot
2007-01-01
The structural and dynamical properties of high-spin Ru 2+ in aqueous solution have been theoretically studied using molecular dynamics (MD) simulations. The conventional MD simulation based on pair potentials gives the overestimated average first shell coordination number of 9, whereas the value of 5.9 was observed when the three-body corrected function was included. A combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation has been performed to take into account the many-body effects on the hydration shell structure of Ru 2+. The most important region, the first hydration shell, was treated by ab initio quantum mechanics at UHF level using the SBKJC VDZ ECP basis set for Ru 2+ and the 6-31G ∗ basis sets for water. An exact coordination number of 6 for the first hydration shell was obtained from the QM/MM simulation. The QM/MM simulation predicts the average Ru 2+-O distance of 2.42 Å for the first hydration shell, whereas the values of 2.34 and 2.46 Å are resulted from the pair potentials without and with the three-body corrected simulations, respectively. Several other structural properties representing position and orientation of the solvate molecules were evaluated for describing the hydration shell structure of the Ru 2+ ion in dilute aqueous solution. A mean residence time of 7.1 ps was obtained for water ligands residing in the second hydration shell.
Efficient conformational space exploration in ab initio protein folding simulation
Ullah, Ahammed; Ahmed, Nasif; Pappu, Subrata Dey; Shatabda, Swakkhar; Ullah, A. Z. M. Dayem; Rahman, M. Sohel
2015-01-01
Ab initio protein folding simulation largely depends on knowledge-based energy functions that are derived from known protein structures using statistical methods. These knowledge-based energy functions provide us with a good approximation of real protein energetics. However, these energy functions are not very informative for search algorithms and fail to distinguish the types of amino acid interactions that contribute largely to the energy function from those that do not. As a result, search algorithms frequently get trapped into the local minima. On the other hand, the hydrophobic–polar (HP) model considers hydrophobic interactions only. The simplified nature of HP energy function makes it limited only to a low-resolution model. In this paper, we present a strategy to derive a non-uniform scaled version of the real 20×20 pairwise energy function. The non-uniform scaling helps tackle the difficulty faced by a real energy function, whereas the integration of 20×20 pairwise information overcomes the limitations faced by the HP energy function. Here, we have applied a derived energy function with a genetic algorithm on discrete lattices. On a standard set of benchmark protein sequences, our approach significantly outperforms the state-of-the-art methods for similar models. Our approach has been able to explore regions of the conformational space which all the previous methods have failed to explore. Effectiveness of the derived energy function is presented by showing qualitative differences and similarities of the sampled structures to the native structures. Number of objective function evaluation in a single run of the algorithm is used as a comparison metric to demonstrate efficiency. PMID:26361554
Cosmic-Ray Modulation: an Ab Initio Approach
NASA Astrophysics Data System (ADS)
Engelbrecht, N. E.; Burger, R. A.
2014-10-01
A better understanding of cosmic-ray modulation in the heliosphere can only be gained through a proper understanding of the effects of turbulence on the diffusion and drift of cosmic rays. We present an ab initio model for cosmic-ray modulation, incorporating for the first time the results yielded by a two-component turbulence transport model. This model is solved for periods of minimum solar activity, utilizing boundary values chosen so that model results are in fair to good agreement with spacecraft observations of turbulence quantities, not only in the solar ecliptic plane but also along the out-of-ecliptic trajectory of the Ulysses spacecraft. These results are employed as inputs for modelled slab and 2D turbulence energy spectra. The latter spectrum is chosen based on physical considerations, with a drop-off at the very lowest wavenumbers commencing at the 2D outerscale. There currently exist no models or observations for this quantity, and it is the only free parameter in this study. The modelled turbulence 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 cosmic-ray drifts are modelled in a self-consistent way, employing a recently developed model for drift along the wavy current sheet. The resulting diffusion coefficients and drift expressions are applied to the study of galactic cosmic-ray protons and antiprotons using a three-dimensional, steady-state cosmic-ray modulation code, and sample solutions in fair agreement with multiple spacecraft observations are presented.
Chiroptical properties of unsubstituted carbohydrates: Ab initio and semiempirical studies
NASA Astrophysics Data System (ADS)
Parra C., Alejandro
Ab initio calculations support assignment of the vacuum ultraviolet circular dichroism (CD) of simple saccharides to 11A 1 --> 21B1 and 11A 1 --> 11A2 transitions centered on the oxygen atoms of the acetal group treated as two weakly coupled ether chromophores. The calculations are consistent with assignments previously made on the basis of a deconvolution of CD spectra. Estimates of the oxygen centered contributions to magnetic transition dipole moments were made. Semiempirical calculations were performed to model the NaD molar optical rotation of 1,6- and 3,6- anhydrosugars. For 1,6-anhydrosugars, current parameters produce reasonable agreement with experimental values. For 3,6-anhydrosugars, modifications to the ether parameters had to be introduced. The most relevant included a reorientation of the bond-centered s-->s* transition dipole charges in the ether chromophore to a C2v orientation, and a shift from prolate polarizability ellipsoids to general ellipsoids. These changes result in good agreement with experimental Na D molar rotations for 3,6-anhydrosugars. A low energy CD band arises in 3,6- and 1,6-anhydrosugars when agreement with the experimental NaD molar rotations is achieved. It is proposed that this band is a real feature in the spectrum. The origin of the band is primarily the interaction between b1 symmetry- oriented transition dipoles in the COC groups with other transition dipoles in the molecule. Comparison with experimental spectra leads to an assignment of this band to 11A1 --> 21B1 transitions centered on the COC groups.
Efficient conformational space exploration in ab initio protein folding simulation.
Ullah, Ahammed; Ahmed, Nasif; Pappu, Subrata Dey; Shatabda, Swakkhar; Ullah, A Z M Dayem; Rahman, M Sohel
2015-08-01
Ab initio protein folding simulation largely depends on knowledge-based energy functions that are derived from known protein structures using statistical methods. These knowledge-based energy functions provide us with a good approximation of real protein energetics. However, these energy functions are not very informative for search algorithms and fail to distinguish the types of amino acid interactions that contribute largely to the energy function from those that do not. As a result, search algorithms frequently get trapped into the local minima. On the other hand, the hydrophobic-polar (HP) model considers hydrophobic interactions only. The simplified nature of HP energy function makes it limited only to a low-resolution model. In this paper, we present a strategy to derive a non-uniform scaled version of the real 20×20 pairwise energy function. The non-uniform scaling helps tackle the difficulty faced by a real energy function, whereas the integration of 20×20 pairwise information overcomes the limitations faced by the HP energy function. Here, we have applied a derived energy function with a genetic algorithm on discrete lattices. On a standard set of benchmark protein sequences, our approach significantly outperforms the state-of-the-art methods for similar models. Our approach has been able to explore regions of the conformational space which all the previous methods have failed to explore. Effectiveness of the derived energy function is presented by showing qualitative differences and similarities of the sampled structures to the native structures. Number of objective function evaluation in a single run of the algorithm is used as a comparison metric to demonstrate efficiency. PMID:26361554
Ab initio quantum transport calculations using plane waves
NASA Astrophysics Data System (ADS)
Garcia-Lekue, A.; Vergniory, M. G.; Jiang, X. W.; Wang, L. W.
2015-08-01
We present an ab initio method to calculate elastic quantum transport at the nanoscale. The method is based on a combination of density functional theory using plane wave nonlocal pseudopotentials and the use of auxiliary periodic boundary conditions to obtain the scattering states. The method can be applied to any applied bias voltage and the charge density and potential profile can either be calculated self-consistently, or using an approximated self-consistent field (SCF) approach. Based on the scattering states one can straightforwardly calculate the transmission coefficients and the corresponding electronic current. The overall scheme allows us to obtain accurate and numerically stable solutions for the elastic transport, with a computational time similar to that of a ground state calculation. This method is particularly suitable for calculations of tunneling currents through vacuum, that some of the nonequilibrium Greens function (NEGF) approaches based on atomic basis sets might have difficulty to deal with. Several examples are provided using this method from electron tunneling, to molecular electronics, to electronic devices: (i) On a Au nanojunction, the tunneling current dependence on the electrode-electrode distance is investigated. (ii) The tunneling through field emission resonances (FERs) is studied via an accurate description of the surface vacuum states. (iii) Based on quantum transport calculations, we have designed a molecular conformational switch, which can turn on and off a molecular junction by applying a perpendicular electric field. (iv) Finally, we have used the method to simulate tunnel field-effect transistors (TFETs) based on two-dimensional transition-metal dichalcogenides (TMDCs), where we have studied the performance and scaling limits of such nanodevices and proposed atomic doping to enhance the transistor performance.