1
|
Zhou X, Huang Z, He X. Diffusion Monte Carlo method for barrier heights of multiple proton exchanges and complexation energies in small water, ammonia, and hydrogen fluoride clusters. J Chem Phys 2024; 160:054103. [PMID: 38310472 DOI: 10.1063/5.0182164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2023] [Accepted: 01/09/2024] [Indexed: 02/05/2024] Open
Abstract
Proton exchange reactions are of key importance in many processes in water. However, it is nontrivial to achieve reliable barrier heights for multiple proton exchanges and complexation energies in hydrogen-bonded systems theoretically. Performance of the fixed-node diffusion quantum Monte Carlo (FN-DMC) with the single-Slater-Jastrow trial wavefunction on total energies, barrier heights of multiple proton exchanges, and complexation energies of small water, ammonia, and hydrogen fluoride clusters is investigated in this study. Effects of basis sets and those of locality approximation (LA), T-move approximation (T-move), and determinant localization approximation (DLA) schemes in dealing with the nonlocal part of pseudopotentials on FN-DMC results are evaluated. According to our results, diffuse basis functions are important in achieving reliable barrier heights and complexation energies with FN-DMC, although the cardinal number of the basis set is more important than diffuse basis functions on total energies of these systems. Our results also show that the time step bias with DLA and LA is smaller than T-move; however, the time step bias of DMC energies with respect to time steps using the T-move is roughly linear up to 0.06 a.u., while this is not the case with LA and DLA. Barrier heights and complexation energies with FN-DMC using these three schemes are always within chemical accuracy. Taking into account the fact that T-move and DLA are typically more stable than LA, FN-DMC calculations with the T-move or DLA scheme and basis sets containing diffuse basis functions are suggested for barrier heights of multiple proton exchanges and complexation energies of hydrogen-bonded clusters.
Collapse
Affiliation(s)
- Xiaojun Zhou
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, People's Republic of China
- School of Physics & Information Science, Shaanxi University of Science and Technology, Xi'an 710021, People's Republic of China
| | - Zhiru Huang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, People's Republic of China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, People's Republic of China
- Chongqing Key Laboratory of Precision Optics, Chongqing Institute of East China Normal University, Chongqing 401120, People's Republic of China
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, People's Republic of China
| |
Collapse
|
2
|
Novotný M, Dubecký M, Karlický F. Toward accurate modeling of structure and energetics of bulk hexagonal boron nitride. J Comput Chem 2024; 45:115-121. [PMID: 37737623 DOI: 10.1002/jcc.27222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Revised: 08/17/2023] [Accepted: 08/25/2023] [Indexed: 09/23/2023]
Abstract
Materials that exhibit both strong covalent and weak van der Waals interactions pose a considerable challenge to many computational methods, such as DFT. This makes assessing the accuracy of calculated properties, such as exfoliation energies in layered materials like hexagonal boron nitride (h-BN) problematic, when experimental data are not available. In this paper, we investigate the accuracy of equilibrium lattice constants and exfoliation energy calculation for various DFT-based computational approaches in bulk h-BN. We contrast these results with available experiments and reference fixed-node diffusion quantum Monte Carlo (QMC) results. From our reference QMC calculation, we obtained an exfoliation energy of - 33 ± 2 meV/atom (-0.38 ± 0.02 J/m2 ).
Collapse
Affiliation(s)
- Michal Novotný
- Department of Physics, Faculty of Science, University of Ostrava, Ostrava, Czech Republic
| | - Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, Ostrava, Czech Republic
- ATRI, Faculty of Materials Science and Technology in Trnava, Slovak University of Technology in Bratislava, Trnava, Slovakia
| | - František Karlický
- Department of Physics, Faculty of Science, University of Ostrava, Ostrava, Czech Republic
| |
Collapse
|
3
|
Ichibha T, Nikaido Y, Bennett MC, Krogel JT, Hongo K, Maezono R, Reboredo FA. Locality error free effective core potentials for 3d transition metal elements developed for the diffusion Monte Carlo method. J Chem Phys 2023; 159:164114. [PMID: 37888761 DOI: 10.1063/5.0175381] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Accepted: 10/05/2023] [Indexed: 10/28/2023] Open
Abstract
Pseudopotential locality errors have hampered the applications of the diffusion Monte Carlo (DMC) method in materials containing transition metals, in particular oxides. We have developed locality error free effective core potentials, pseudo-Hamiltonians, for transition metals ranging from Cr to Zn. We have modified a procedure published by some of us in Bennett et al. [J. Chem. Theory Comput. 18, 828 (2022)]. We carefully optimized our pseudo-Hamiltonians and achieved transferability errors comparable to the best semilocal pseudopotentials used with DMC but without incurring in locality errors. Our pseudo-Hamiltonian set (named OPH23) bears the potential to significantly improve the accuracy of many-body-first-principles calculations in fundamental science research of complex materials involving transition metals.
Collapse
Affiliation(s)
- Tom Ichibha
- School of Information Science, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Yutaka Nikaido
- School of Information Science, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| | - M Chandler Bennett
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Kenta Hongo
- Research Center for Advanced Computing Infrastructure, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| | - Ryo Maezono
- School of Information Science, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| | - Fernando A Reboredo
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| |
Collapse
|
4
|
Tyagi R, Zen A, Voora VK. Quantifying the Impact of Halogenation on Intermolecular Interactions and Binding Modes of Aromatic Molecules. J Phys Chem A 2023. [PMID: 37406194 DOI: 10.1021/acs.jpca.3c02291] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/07/2023]
Abstract
Halogenation of aromatic molecules is frequently used to modulate intermolecular interactions with ramifications for optoelectronic and mechanical properties. In this work, we accurately quantify and understand the nature of intermolecular interactions in perhalogenated benzene (PHB) clusters. Using benchmark binding energies from the fixed-node diffusion Monte Carlo (FN-DMC) method, we show that generalized Kohn-Sham semicanonical projected random phase approximation (GKS-spRPA) plus approximate exchange kernel (AKX) provides reliable interaction energies with mean absolute error (MAE) of 0.23 kcal/mol. Using the GKS-spRPA+AXK method, we quantify the interaction energies of several binding modes of PHB clusters ((C6X6)n; X = F, Cl, Br, I; n = 2, 3). For a given binding mode, the interaction energies increase 3-4 times from X = F to X = I; the X-X binding modes have energies in the range of 2-4 kcal/mol, while the π-π binding mode has interaction energies in the range of 4-12 kcal/mol. SAPT-DFT-based energy decomposition analysis is then used to show that the equilibrium geometries are dictated primarily by the dispersion and exchange interactions. Finally, we test the accuracy of several dispersion-corrected density functional approximations and show that only the r2SCAN-D4 method has a low MAE and correct long-range behavior, which makes it suitable for large-scale simulations and for developing structure-function relationships of halogenated aromatic systems.
Collapse
Affiliation(s)
- Ritaj Tyagi
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Andrea Zen
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
| | - Vamsee K Voora
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| |
Collapse
|
5
|
Huang B, von Lilienfeld OA, Krogel JT, Benali A. Toward DMC Accuracy Across Chemical Space with Scalable Δ-QML. J Chem Theory Comput 2023; 19:1711-1721. [PMID: 36857531 DOI: 10.1021/acs.jctc.2c01058] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2023]
Abstract
In the past decade, quantum diffusion Monte Carlo (DMC) has been demonstrated to successfully predict the energetics and properties of a wide range of molecules and solids by numerically solving the electronic many-body Schrödinger equation. With O(N3) scaling with the number of electrons N, DMC has the potential to be a reference method for larger systems that are not accessible to more traditional methods such as CCSD(T). Assessing the accuracy of DMC for smaller molecules becomes the stepping stone in making the method a reference for larger systems. We show that when coupled with quantum machine learning (QML)-based surrogate methods, the computational burden can be alleviated such that quantum Monte Carlo (QMC) shows clear potential to undergird the formation of high-quality descriptions across chemical space. We discuss three crucial approximations necessary to accomplish this: the fixed-node approximation, universal and accurate references for chemical bond dissociation energies, and scalable minimal amons-set-based QML (AQML) models. Numerical evidence presented includes converged DMC results for over 1000 small organic molecules with up to five heavy atoms used as amons and 50 medium-sized organic molecules with nine heavy atoms to validate the AQML predictions. Numerical evidence collected for Δ-AQML models suggests that already modestly sized QMC training data sets of amons suffice to predict total energies with near chemical accuracy throughout chemical space.
Collapse
Affiliation(s)
- Bing Huang
- Faculty of Physics, University of Vienna, Kolingasse 14-16, 1090 Vienna, Austria
| | - O Anatole von Lilienfeld
- Departments of Chemistry, Materials Science and Engineering, and Physics, University of Toronto, St. George Campus, Toronto, ON M5S 3H6, Canada.,Machine Learning Group, Technische Universität Berlin and Institute for the Foundations of Learning and Data, 10587 Berlin, Germany.,Vector Institute for Artificial Intelligence, Toronto, ON M5S 1M1, Canada
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Anouar Benali
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| |
Collapse
|
6
|
Kolesár V, Dubecký M. Accuracy of Noncovalent Interactions Involving d-Elements by the 1-Determinant Fixed-Node Diffusion Monte Carlo Method with Effective Core Potentials. J Chem Theory Comput 2023; 19:1170-1176. [PMID: 36751996 DOI: 10.1021/acs.jctc.2c00872] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
A critical assessment of effective core potential (ECP)-based single-determinant (SD) fixed-node diffusion quantum Monte Carlo (FNDMC) accuracy in prototypical noncovalent closed-shell systems involving d-elements is presented. Careful analysis of biases and elimination of possible bias sources leads to two findings of practical importance for SD FNDMC in these systems. First, in some systems (HCu:HCu, HCu:CuH), SD FNDMC reveals large biases of interaction energy differences (significantly exceeding the target 2% relative error) vs a reliable coupled-cluster CCSD(T)/CBS (complete basis set) reference. Second, the leading error of SD FNDMC with ECPs was attributed to a higher nuclear charge Z of d-group (pseudo) atoms, when compared to sp elements, in line with a previously reported finding that aggregate SD FNDMC bias tends to increase in systems with higher electronic densities. Therefore, SD FNDMC should only be used with caution in systems with a large Z.
Collapse
Affiliation(s)
- Vladimír Kolesár
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
| | - Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic.,ATRI, Faculty of Materials Science and Technology in Trnava, Slovak University of Technology in Bratislava, J. Bottu 25, 917 24 Trnava, Slovakia
| |
Collapse
|
7
|
Nagy PR, Gyevi-Nagy L, Lőrincz BD, Kállay M. Pursuing the basis set limit of CCSD(T) non-covalent interaction energies for medium-sized complexes: case study on the S66 compilation. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2109526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Péter R. Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - László Gyevi-Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Balázs D. Lőrincz
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Mihály Kállay
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| |
Collapse
|
8
|
Rao L, Wang F. Diffusion quantum Monte Carlo method on diradicals using single- and multi-determinant-Jastrow trial wavefunctions and different orbitals. J Chem Phys 2022; 156:124308. [DOI: 10.1063/5.0086606] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, the diffusion quantum Monte Carlo (DMC) method is employed to calculate the energies of singlet and triplet states for a series of organic diradicals and diatomic diradicals with π2 configuration. Single-determinant-Jastrow (SDJ) trial wavefunctions for triplet states, two-determinant-Jastrow (2DJ) trial wavefunctions for the singlet states, and multi-determinant-Jastrow (MDJ) trial wavefunctions are employed in DMC calculations using restricted open-shell B3LYP (ROB3LYP) orbitals, complete-active-space self-consistent field (CASSCF) orbitals, state-average CASSCF orbitals, or frozen-CASSCF orbitals. Our results show that DMC energies using either SDJ/2DJ or MDJ with ROB3LYP orbitals are close to or lower than those with the other orbitals for organic diradicals, while they are not very sensitive to the employed orbitals for diatomic diradicals. Furthermore, using MDJ can reduce DMC energies to some extent for most of the investigated organic diradicals and some diatomic diradicals. The importance of MDJ on DMC energies can be estimated based on the percentage of main determinants in the CASCI wavefunction. On the other hand, singlet–triplet gaps can be calculated reasonably with DMC using MDJ with a mean absolute error of less than 2 kcal/mol with all these orbitals. CASCI wavefunctions using density functional theory orbitals are preferred in constructing MDJ trial wavefunctions in practical DMC calculations since it is easier to obtain such wavefunctions than CASSCF methods.
Collapse
Affiliation(s)
- Lu Rao
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610065, People’s Republic of China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610065, People’s Republic of China
| |
Collapse
|
9
|
Zhou H, Scemama A, Wang G, Annaberdiyev A, Kincaid B, Caffarel M, Mitas L. A quantum Monte Carlo study of systems with effective core potentials and node nonlinearities. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2021.111402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
10
|
Kartashynska ES, Vollhardt D. Quantum chemical assessment of the molecular area corresponding to the onset of the LE-LC phase transition for amphiphilic 2D monolayers at the air/water interface. Phys Chem Chem Phys 2021; 23:25356-25364. [PMID: 34750603 DOI: 10.1039/d1cp03511h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
An approach for the assessment of the area per surfactant molecule in a monolayer at the onset of the LE-LC phase transition (Ac) is proposed based on the quantum chemical approach and a thermodynamic model for amphiphilic monolayers, which takes into account the nonideality of the mixing entropy. The values of the Gibbs' clusterization energy for small surfactant associates, as well as the geometric parameters of the monolayer unit cells, were used, previously calculated using the semiempirical PM3 method for eight classes of amphiphilic compounds: saturated and ethoxylated alcohols, saturated and unsaturated cis-carboxylic acids, α-hydroxylic and α-aminoacids, N-acyl-substituted alanine and dialkyl-substituted melamine. The obtained Ac values are in satisfactory agreement with the available experimental data. This allows using the proposed approach for prognostic purposes in the cases where there are no corresponding π-A isotherms for necessary surfactants, but there are calculated thermodynamic and structural parameters of its clusterization.
Collapse
Affiliation(s)
- Elena S Kartashynska
- L.M. Litvinenko Institute of Physical Organic and Coal Chemistry, 70 R. Luxemburg Str., Donetsk 83114, Ukraine.,Donetsk National Technical University, 58 Artema Str., Donetsk 83000, Ukraine
| | - Dieter Vollhardt
- Max Planck Institute of Polymer Research, Mainz D-55128, Germany.
| |
Collapse
|
11
|
Prayogo GI, Shin H, Benali A, Maezono R, Hongo K. Importance of Van der Waals Interactions in Hydrogen Adsorption on a Silicon-carbide Nanotube Revisited with vdW-DFT and Quantum Monte Carlo. ACS OMEGA 2021; 6:24630-24636. [PMID: 34604645 PMCID: PMC8482461 DOI: 10.1021/acsomega.1c03318] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/25/2021] [Accepted: 08/13/2021] [Indexed: 06/13/2023]
Abstract
Density functional theory (DFT) is a valuable tool for calculating adsorption energies toward designing materials for hydrogen storage. However, dispersion forces being absent from the local/semi-local theory, it remains unclear as to how the consideration of van der Waals (vdW) interactions affects such calculations. For the first time, we applied diffusion Monte Carlo (DMC) to evaluate the adsorption characteristics of a hydrogen molecule on a (5,5) armchair silicon-carbide nanotube (H2-SiCNT). Within the DFT framework, we benchmarked various exchange-correlation functionals, including those recently developed for treating dispersion or vdW interactions. We found that the vdW-corrected DFT methods agree well with DMC, whereas the local (semilocal) functional significantly over (under)-binds. Furthermore, we fully optimized the H2-SiCNT geometry within the DFT framework and investigated the correlation between the structure and charge density. The vdW contribution to the adsorption was found to be non-negligible at ∼1 kcal/mol per hydrogen molecule, which amounts to 9-29% of the ideal adsorption energy required for hydrogen storage applications.
Collapse
Affiliation(s)
- Genki I. Prayogo
- School
of Information Science, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| | - Hyeondeok Shin
- Computational
Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Anouar Benali
- Computational
Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Ryo Maezono
- School
of Information Science, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| | - Kenta Hongo
- Research
Center for Advanced Computing Infrastructure, JAIST, Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| |
Collapse
|
12
|
Doran AE, Qiu DL, Hirata S. Monte Carlo MP2-F12 for Noncovalent Interactions: The C 60 Dimer. J Phys Chem A 2021; 125:7344-7351. [PMID: 34433271 DOI: 10.1021/acs.jpca.1c05021] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green's functions in real space and imaginary time. These integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol-1 for the C60 dimer at the center distance of 9.748 Å.
Collapse
Affiliation(s)
- Alexander E Doran
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - David L Qiu
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
13
|
Ballesteros F, Dunivan S, Lao KU. Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA-ellipticine and buckycatcher-fullerene. J Chem Phys 2021; 154:154104. [PMID: 33887937 DOI: 10.1063/5.0042906] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA-ellipticine intercalation complex, and the buckycatcher-C60 complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller-Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA-ellipticine, and buckycatcher-C60, respectively. Such sizable discrepancy above "chemical accuracy" for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree-Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.
Collapse
Affiliation(s)
- Francisco Ballesteros
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, USA
| | - Shelbie Dunivan
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, USA
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, USA
| |
Collapse
|
14
|
Fanta R, Dubecký M. Noncovalent Interactions by the Quantum Monte Carlo Method: Strong Influence of Isotropic Jastrow Cutoff Radii. J Chem Theory Comput 2021; 17:4242-4249. [PMID: 34169721 DOI: 10.1021/acs.jctc.1c00467] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a paradigmatic example of a strong effect of Jastrow cutoff radii setup on the accuracy of noncovalent interaction energy differences within one-determinant Slater-Jastrow fixed-node diffusion Monte Carlo (1FNDMC) simulations using isotropic Jastrow terms and effective-core potentials. Analysis of total energies, absolute and relative errors, and local energy variance of energy differences vs the reference results suggests a simple procedure to marginalize the related biases. The presented data showcase improvements in dispersion-bounded systems within such a 1FNDMC method.
Collapse
Affiliation(s)
- Roman Fanta
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
| | - Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic.,ATRI, Faculty of Materials Science and Technology in Trnava, Slovak University of Technology in Bratislava, J. Bottu 25, 917 24 Trnava, Slovakia
| |
Collapse
|
15
|
Interactions between large molecules pose a puzzle for reference quantum mechanical methods. Nat Commun 2021; 12:3927. [PMID: 34168142 PMCID: PMC8225865 DOI: 10.1038/s41467-021-24119-3] [Citation(s) in RCA: 45] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2020] [Accepted: 06/02/2021] [Indexed: 02/05/2023] Open
Abstract
Quantum-mechanical methods are used for understanding molecular interactions throughout the natural sciences. Quantum diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] are state-of-the-art trusted wavefunction methods that have been shown to yield accurate interaction energies for small organic molecules. These methods provide valuable reference information for widely-used semi-empirical and machine learning potentials, especially where experimental information is scarce. However, agreement for systems beyond small molecules is a crucial remaining milestone for cementing the benchmark accuracy of these methods. We show that CCSD(T) and DMC interaction energies are not consistent for a set of polarizable supramolecules. Whilst there is agreement for some of the complexes, in a few key systems disagreements of up to 8 kcal mol-1 remain. These findings thus indicate that more caution is required when aiming at reproducible non-covalent interactions between extended molecules.
Collapse
|
16
|
Kartashynska E, Vysotsky YB, Vollhardt D, Fainerman VB. Relationship between the Bulk and Surface Basicity of Aliphatic Amines: A Quantum Chemical Approach. ACS OMEGA 2020; 5:32032-32039. [PMID: 33344857 PMCID: PMC7745424 DOI: 10.1021/acsomega.0c04939] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/09/2020] [Accepted: 11/19/2020] [Indexed: 05/23/2023]
Abstract
To assess the surface basicity constant (pK b) of aliphatic amine films, the use of a theoretical approach recently developed to evaluate the pK a of carboxylic acid monolayers on the water surface is tested. The present paper gives a new full picture of the change of acid-base properties of surfactants during their aggregation at the air/water interface. The exploited approach is simple because it does not involve the construction of thermodynamic cycles but uses the Gibbs energies of the formation and dimerization of surfactant monomers in neutral and ionized forms in the aqueous and gaseous phases. The quantum chemical semiempirical PM3 method is applied to perform calculations using a conductor-like screening model, which takes into account the aqueous phase. The calculation shows that aliphatic amines, as well as carboxylic acids, are characterized by a change of the value of the basicity/acidity constant during the film formation. The film formation of surfactants leads to a decrease in their acid-base properties, i.e., the surface pK a values of carboxylic acids and pK b values of amines increase. However, unlike carboxylic acids, there is practically no dependence of the surface pK b value on the alkyl chain length of the aliphatic amine, which is caused by almost identical contributions of one CH2 fragment to the solvation Gibbs energy of neutral and ionized monomers within the calculation error. The obtained results agree with existing experimental data.
Collapse
Affiliation(s)
- Elena
S. Kartashynska
- L.M.
Litvinenko Institute of Physical Organic and Coal Chemistry, 70 R. Luxemburg Str., 83114 Donetsk, Ukraine
| | - Yuri B. Vysotsky
- Donetsk
National Technical University, 58 Artema Str., 83000 Donetsk, Ukraine
| | - Dieter Vollhardt
- Max
Planck Institute of Polymer Research, D-55128 Mainz, Germany
| | | |
Collapse
|
17
|
Kitzmiller NL, Wolf ME, Turney JM, Schaefer HF. The HOX⋯SO 2 (X=F, Cl, Br, I) Binary Complexes: Implications for Atmospheric Chemistry. Chemphyschem 2020; 22:112-126. [PMID: 33090675 DOI: 10.1002/cphc.202000746] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2020] [Revised: 10/16/2020] [Indexed: 11/07/2022]
Abstract
Sulfur dioxide and hypohalous acids (HOX, X=F, Cl, Br, I) are ubiquitous molecules in the atmosphere that are central to important processes like seasonal ozone depletion, acid rain, and cloud nucleation. We present the first theoretical examination of the HOX⋯SO2 binary complexes and the associated trends due to halogen substitution. Reliable geometries were optimized at the CCSD(T)/aug-cc-pV(T+d)Z level of theory for HOF and HOCl complexes. The HOBr and HOI complexes were optimized at the CCSD(T)/aug-cc-pV(D+d)Z level of theory with the exception of the Br and I atoms which were modeled with an aug-cc-pwCVDZ-PP pseudopotential. 27 HOX⋯SO2 complexes were characterized and the focal point method was employed to produce CCSDT(Q)/CBS interaction energies. Natural Bond Orbital analysis and Symmetry Adapted Perturbation Theory were used to classify the nature of each principle interaction. The interaction energies of all HOX⋯SO2 complexes in this study ranged from 1.35 to 3.81 kcal mol-1 . The single-interaction hydrogen bonded complexes spanned a range of 2.62 to 3.07 kcal mol-1 , while the single-interaction halogen bonded complexes were far more sensitive to halogen substitution ranging from 1.35 to 3.06 kcal mol-1 , indicating that the two types of interactions are extremely competitive for heavier halogens. Our results provide insight into the interactions between HOX and SO2 which may guide further research of related systems.
Collapse
Affiliation(s)
- Nathaniel L Kitzmiller
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia, 30602
| | - Mark E Wolf
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia, 30602
| | - Justin M Turney
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia, 30602
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia, 30602
| |
Collapse
|
18
|
Benali A, Shin H, Heinonen O. Quantum Monte Carlo benchmarking of large noncovalent complexes in the L7 benchmark set. J Chem Phys 2020; 153:194113. [DOI: 10.1063/5.0026275] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Anouar Benali
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Hyeondeok Shin
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Olle Heinonen
- Material Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| |
Collapse
|
19
|
Per MC, Fletcher EK, Swann ET, Cleland DM. Reliable radical stabilization energies from diffusion Monte Carlo calculations. J Comput Chem 2020; 41:2378-2382. [DOI: 10.1002/jcc.26397] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2020] [Revised: 07/10/2020] [Accepted: 07/14/2020] [Indexed: 11/08/2022]
|
20
|
Červinka C, Štejfa V. Sublimation Properties of α,ω-Diamines Revisited from First-Principles Calculations. Chemphyschem 2020; 21:1184-1194. [PMID: 32243713 DOI: 10.1002/cphc.202000108] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2020] [Revised: 04/02/2020] [Indexed: 11/06/2022]
Abstract
Sublimation enthalpies of alkane-α,ω-diamines exhibit an odd-even pattern within their homologous series. First-principles calculations coupled with the quasi-harmonic approximation for crystals and with the conformation mixing model for the ideal gas are used to explain this phenomenon from the theoretical point of view. Crystals of the odd and even alkane-α,ω-diamines distinctly differ in their packing motifs. However, first-principles calculations indicate that it is a delicate interplay of the cohesive forces, phonons, molecular vibrations and conformational equilibrium which governs the odd-even pattern of the sublimation enthalpies within the homologous series. High molecular flexibility of the alkane-α,ω-diamines predetermines higher sensitivity of the computational model to the quality of the optimized geometries and relative conformational energies. Performance of high-throughput computational methods, such as the density functional tight binding (DFTB, GFN2-xTB) and the explicitly correlated dispersion-corrected Møller-Plesset perturbative method (MP2C-F12), are benchmarked against the consistent state-of-the-art calculations of conformational energies and interaction energies, respectively.
Collapse
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28, Prague 6, Czech Republic
| | - Vojtěch Štejfa
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28, Prague 6, Czech Republic
| |
Collapse
|
21
|
Sajid A, Ford MJ, Reimers JR. Single-photon emitters in hexagonal boron nitride: a review of progress. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2020; 83:044501. [PMID: 31846956 DOI: 10.1088/1361-6633/ab6310] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
This report summarizes progress made in understanding properties such as zero-phonon-line energies, emission and absorption polarizations, electron-phonon couplings, strain tuning and hyperfine coupling of single photon emitters in hexagonal boron nitride. The primary aims of this research are to discover the chemical nature of the emitting centres and to facilitate deployment in device applications. Critical analyses of the experimental literature and data interpretation, as well as theoretical approaches used to predict properties, are made. In particular, computational and theoretical limitations and challenges are discussed, with a range of suggestions made to overcome these limitations, striving to achieve realistic predictions concerning the nature of emitting centers. A symbiotic relationship is required in which calculations focus on properties that can easily be measured, whilst experiments deliver results in a form facilitating mass-produced calculations.
Collapse
Affiliation(s)
- A Sajid
- University of Technology Sydney, School of Mathematical and Physical Sciences, Ultimo, New South Wales 2007, Australia. CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark. Department of Physics, GC University Faisalabad, Allama Iqbal Road, 38000 Faisalabad, Pakistan. Author to whom any correspondence should be addressed
| | | | | |
Collapse
|
22
|
Qin KS, Ichibha T, Hongo K, Maezono R. Inconsistencies in ab initio evaluations of non-additive contributions of DNA stacking energies. Chem Phys 2020. [DOI: 10.1016/j.chemphys.2019.110554] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
23
|
Zhou X, Wang F. Singlet-triplet gaps in diradicals obtained with diffusion quantum Monte Carlo using a Slater-Jastrow trial wavefunction with a minimum number of determinants. Phys Chem Chem Phys 2019; 21:20422-20431. [PMID: 31501831 DOI: 10.1039/c9cp03045j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Diradicals are essential species in a wide range of chemical processes, whereas the computational study of their electronic structure often remains a challenge due to near-degeneracy of the frontier molecular orbitals. The fixed-node diffusion quantum Monte Carlo (FN-DMC) method is employed to calculate adiabatic energy gaps of some typical diradicals with the Slater-Jastrow trial wavefunction. The antisymmetrized part of the trial wavefunction is taken to be a linear combination of a minimum number of determinants using RB3LYP orbitals from the closed-shell singlet state or ROB3LYP orbitals from the triplet state. Our results show that using the two-determinant-Jastrow trial wavefunction is necessary to achieve reliable energy differences between closed-shell singlet states. The energy of the triplet state with MS = 1 is calculated to be lower than that with MS = 0 with FN-DMC even using trial wavefunctions with spin-pure states as their antisymmetrized parts and this difference is reduced with better orbitals. This indicates that the fixed-node error is smaller for the triplet state with MS = 1. Adiabatic energy gaps obtained from the present FN-DMC calculations are in reasonable agreement with available experimental values. Compared with results of the high level EOM-SF-CC method, energy gaps of FN-DMC with RB3LYP orbitals are slightly better than those using ROB3LYP orbitals and results of EOM-SF-CCSD. The present FN-DMC calculations using the simplest ansatz for the trial wavefunction can achieve reasonable results for these diradicals and they can readily be applied to large diradicals.
Collapse
Affiliation(s)
- Xiaojun Zhou
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, P. R. China.
| | | |
Collapse
|
24
|
Yang DC, Kim DY, Kim KS. Quantum Monte Carlo Study of the Water Dimer Binding Energy and Halogen-π Interactions. J Phys Chem A 2019; 123:7785-7791. [PMID: 31418568 DOI: 10.1021/acs.jpca.9b04072] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Halogen-π systems are involved with competition between halogen bonding and π-interaction. Using the diffusion quantum Monte Carlo (DMC) method, we compare the interaction of benzene with diatomic halogens (X2: Cl2/Br2) with the typical hydrogen bonding in the water dimer, taking into account explicit correlations of up to three bodies. The benzene-Cl2/Br2 binding energies (13.07 ± 0.42/16.62 ± 0.02 kJ/mol) attributed to both halogen bonding and dispersion are smaller than but comparable to the typical hydrogen bonding in the water dimer binding energy (20.88 ± 0.27 kJ/mol). All of the above values are in good agreement with those from the coupled-cluster with single, double, and noniterative triple excitations (CCSD(T)) results at the complete basis set limit (benzene-Cl2/Br2: 12.78/16.17 kJ/mol; water dimer: 21.0 kJ/mol).
Collapse
Affiliation(s)
- D ChangMo Yang
- Center for Superfunctional Materials, Department of Chemistry , Ulsan National Institute of Science and Technology , Ulsan 44919 , Republic of Korea
| | - Dong Yeon Kim
- Center for Superfunctional Materials, Department of Chemistry , Ulsan National Institute of Science and Technology , Ulsan 44919 , Republic of Korea
| | - Kwang S Kim
- Center for Superfunctional Materials, Department of Chemistry , Ulsan National Institute of Science and Technology , Ulsan 44919 , Republic of Korea
| |
Collapse
|
25
|
Dubecký M, Jurečka P, Mitas L, Ditte M, Fanta R. Toward Accurate Hydrogen Bonds by Scalable Quantum Monte Carlo. J Chem Theory Comput 2019; 15:3552-3557. [DOI: 10.1021/acs.jctc.9b00096] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Matúš Dubecký
- Department of Physics, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
- ATRI, Slovak University of Technology, Paulínska 16, 917 24 Trnava, Slovakia
| | - Petr Jurečka
- Department of Physical Chemistry, Palackỳ University Olomouc, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lubos Mitas
- Department of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Matej Ditte
- Department of Physics, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
| | - Roman Fanta
- Department of Physics, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic
| |
Collapse
|
26
|
Wineman-Fisher V, Al-Hamdani Y, Addou I, Tkatchenko A, Varma S. Ion-Hydroxyl Interactions: From High-Level Quantum Benchmarks to Transferable Polarizable Force Fields. J Chem Theory Comput 2019; 15:2444-2453. [PMID: 30830778 PMCID: PMC6598712 DOI: 10.1021/acs.jctc.8b01198] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Ion descriptors in molecular mechanics models are calibrated against reference data on ion-water interactions. It is then typically assumed that these descriptors will also satisfactorily describe interactions of ions with other functional groups, such as those present in biomolecules. However, several studies now demonstrate that this transferability assumption produces, in many different cases, large errors. Here we address this issue in a representative polarizable model and focus on transferability of cationic interactions from water to a series of alcohols. Both water and alcohols use hydroxyls for ion-coordination, and, therefore, this set of molecules constitutes the simplest possible case of transferability. We obtain gas phase reference data systematically from "gold-standard" quantum Monte Carlo and CCSD(T) methods, followed by benchmarked vdW-corrected DFT. We learn that the original polarizable model yields large gas phase water → alcohol transferability errors - the RMS and maximum errors are 2.3 and 5.1 kcal/mol, respectively. These errors are, nevertheless, systematic in that ion-alcohol interactions are overstabilized, and systematic errors typically imply that some essential physics is either missing or misrepresented. A comprehensive analysis shows that when both low- and high-field responses of ligand dipole polarization are described accurately, then transferability improves significantly - the RMS and maximum errors in the gas phase reduce, respectively, to 0.9 and 2.5 kcal/mol. Additionally, predictions of condensed phase transfer free energies also improve. Nevertheless, within the limits of the extrathermodynamic assumptions necessary to separate experimental estimates of salt dissolution into constituent cationic and anionic contributions, we note that the error in the condensed phase is systematic, which we attribute, at least, partially to the parametrization in long-range electrostatics. Overall, this work demonstrates a rational approach to boosting transferability of ionic interactions that will be applicable broadly to improving other polarizable and nonpolarizable models.
Collapse
Affiliation(s)
- Vered Wineman-Fisher
- Department of Cell Biology, Microbiology and Molecular Biology , University of South Florida , Tampa , Florida 33620 , United States
| | - Yasmine Al-Hamdani
- Physics and Materials Science Research Unit , University of Luxembourg , 162a avenue de la Fïancerie , Luxembourg City , L-1511 , Luxembourg
| | - Iqbal Addou
- Department of Cell Biology, Microbiology and Molecular Biology , University of South Florida , Tampa , Florida 33620 , United States
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit , University of Luxembourg , 162a avenue de la Fïancerie , Luxembourg City , L-1511 , Luxembourg
| | - Sameer Varma
- Department of Cell Biology, Microbiology and Molecular Biology , University of South Florida , Tampa , Florida 33620 , United States
| |
Collapse
|
27
|
Al-Hamdani YS, Tkatchenko A. Understanding non-covalent interactions in larger molecular complexes from first principles. J Chem Phys 2019; 150:010901. [PMID: 30621423 PMCID: PMC6910608 DOI: 10.1063/1.5075487] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2018] [Accepted: 12/05/2018] [Indexed: 01/02/2023] Open
Abstract
Non-covalent interactions pervade all matter and play a fundamental role in layered materials, biological systems, and large molecular complexes. Despite this, our accumulated understanding of non-covalent interactions to date has been mainly developed in the tens-of-atoms molecular regime. This falls considerably short of the scales at which we would like to understand energy trends, structural properties, and temperature dependencies in materials where non-covalent interactions have an appreciable role. However, as more reference information is obtained beyond moderately sized molecular systems, our understanding is improving and we stand to gain pertinent insights by tackling more complex systems, such as supramolecular complexes, molecular crystals, and other soft materials. In addition, accurate reference information is needed to provide the drive for extending the predictive power of more efficient workhorse methods, such as density functional approximations that also approximate van der Waals dispersion interactions. In this perspective, we discuss the first-principles approaches that have been used to obtain reference interaction energies for beyond modestly sized molecular complexes. The methods include quantum Monte Carlo, symmetry-adapted perturbation theory, non-canonical coupled cluster theory, and approaches based on the random-phase approximation. By considering the approximations that underpin each method, the most accurate theoretical references for supramolecular complexes and molecular crystals to date are ascertained. With these, we also assess a handful of widely used exchange-correlation functionals in density functional theory. The discussion culminates in a framework for putting into perspective the accuracy of high-level wavefunction-based methods and identifying future challenges.
Collapse
Affiliation(s)
- Yasmine S Al-Hamdani
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
28
|
Kartashynska E, Vysotsky Y, Fainerman V, Vollhardt D, Miller R. Quantum-chemical analysis of condensed monolayer phases of N-alkanoyl-substituted alanine at the air/water interface. Colloids Surf A Physicochem Eng Asp 2018. [DOI: 10.1016/j.colsurfa.2018.03.027] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
29
|
Song S, Kim MC, Sim E, Benali A, Heinonen O, Burke K. Benchmarks and Reliable DFT Results for Spin Gaps of Small Ligand Fe(II) Complexes. J Chem Theory Comput 2018; 14:2304-2311. [DOI: 10.1021/acs.jctc.7b01196] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Suhwan Song
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Min-Cheol Kim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | | | | | - Kieron Burke
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| |
Collapse
|
30
|
Reimers JR, Sajid A, Kobayashi R, Ford MJ. Understanding and Calibrating Density-Functional-Theory Calculations Describing the Energy and Spectroscopy of Defect Sites in Hexagonal Boron Nitride. J Chem Theory Comput 2018; 14:1602-1613. [DOI: 10.1021/acs.jctc.7b01072] [Citation(s) in RCA: 53] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jeffrey R. Reimers
- International Centre for Quantum and Molecular Structures and Department of Physics, Shanghai University, Shanghai 200444, China
- School of Mathematical and Physical Sciences, University of Technology Sydney, Ultimo, New South Wales 2007, Australia
| | - A. Sajid
- School of Mathematical and Physical Sciences, University of Technology Sydney, Ultimo, New South Wales 2007, Australia
- Department of Physics, GC University Faisalabad, Allama Iqbal Road, 38000 Faisalabad, Pakistan
| | - Rika Kobayashi
- National Computational Infrastructure, The Australian National University, Canberra, Austrailian Capital Territory 2600, Australia
| | - Michael J. Ford
- School of Mathematical and Physical Sciences, University of Technology Sydney, Ultimo, New South Wales 2007, Australia
| |
Collapse
|
31
|
Hongo K, Maezono R. A Computational Scheme To Evaluate Hamaker Constants of Molecules with Practical Size and Anisotropy. J Chem Theory Comput 2017; 13:5217-5230. [PMID: 28981266 DOI: 10.1021/acs.jctc.6b01159] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We propose a computational scheme to evaluate Hamaker constants, A, of molecules with practical sizes and anisotropies. Upon the increasing feasibility of diffusion Monte Carlo (DMC) methods to evaluate binding curves for such molecules to extract the constants, we discussed how to treat the averaging over anisotropy and how to correct the bias due to the nonadditivity. We have developed a computational procedure for dealing with the anisotropy and reducing statistical errors and biases in DMC evaluations, based on possible validations on predicted A. We applied the scheme to cyclohexasilane molecule, Si6H12, used in "printed electronics" fabrications, getting A ≈ 105 ± 2 zJ, being in plausible range supported even by other possible extrapolations. The scheme provided here would open a way to use handy ab initio evaluations to predict wettabilities as in the form of materials informatics over broader molecules.
Collapse
Affiliation(s)
- Kenta Hongo
- Research Center for Advanced Computing Infrastructure, JAIST , Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan.,Center for Materials Research by Information Integration, Research and Services Division of Materials Data and Integrated System, National Institute for Materials Science , Tsukuba 305-0047, Japan.,PRESTO, Japan Science and Technology Agency (JST) , Kawaguchi, Saitama 332-0012, Japan
| | - Ryo Maezono
- School of Information Science, JAIST , Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| |
Collapse
|
32
|
Dubecký M. Noncovalent Interactions by Fixed-Node Diffusion Monte Carlo: Convergence of Nodes and Energy Differences vs Gaussian Basis-Set Size. J Chem Theory Comput 2017; 13:3626-3635. [DOI: 10.1021/acs.jctc.7b00537] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701
03 Ostrava, Czech Republic
- ATRI, Faculty of Materials
Science and Technology, Slovak University of Technology, Paulínska
16, 917 24 Trnava, Slovakia
| |
Collapse
|
33
|
Krongchon K, Busemeyer B, Wagner LK. Accurate barrier heights using diffusion Monte Carlo. J Chem Phys 2017; 146:124129. [DOI: 10.1063/1.4979059] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Kittithat Krongchon
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Brian Busemeyer
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Lucas K. Wagner
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
34
|
Dubecký M. Bias cancellation in one-determinant fixed-node diffusion Monte Carlo: Insights from fermionic occupation numbers. Phys Rev E 2017; 95:033308. [PMID: 28415179 DOI: 10.1103/physreve.95.033308] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2016] [Indexed: 06/07/2023]
Abstract
The accuracy of the fixed-node diffusion Monte Carlo (FNDMC) depends on the node location of the supplied trial state Ψ_{T}. The practical FNDMC approaches available for large systems rely on compact yet effective Ψ_{T}, most often containing an explicitly correlated single Slater determinant (SD). However, SD nodes may be better suited to one system than to another, which may possibly lead to inaccurate FNDMC energy differences. It remains a challenge how to estimate nonequivalence or appropriateness of SDs. Here we use the differences of a measure based on the Euclidean distance between the natural orbital occupation number (NOON) vector of the SD and the exact solution in the NOON vector space, which can be viewed as a measure of SD nonequivalence and as a qualitative measure of the expected degree of nondynamic-correlation-related bias in FNDMC energy differences. This is explored on a set of small noncovalent complexes and covalent bond breaking of Si_{2} vs N_{2}. It turns out that NOON-based measures well reflect the magnitude and sign of the bias present in the data available, thus providing insights into the nature of bias cancellation in SD FNDMC energy differences.
Collapse
Affiliation(s)
- Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701 03 Ostrava, Czech Republic and ATRI, Faculty of Materials Science and Technology, Slovak University of Technology, Paulínska 16, 917 24 Trnava, Slovakia
| |
Collapse
|
35
|
Zhou X, Wang F. Barrier heights of hydrogen-transfer reactions with diffusion quantum monte carlo method. J Comput Chem 2017; 38:798-806. [DOI: 10.1002/jcc.24750] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2016] [Revised: 01/11/2017] [Accepted: 01/13/2017] [Indexed: 01/17/2023]
Affiliation(s)
- Xiaojun Zhou
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University; Chengdu 610065 People's Republic of China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University; Chengdu 610065 People's Republic of China
| |
Collapse
|
36
|
Hermann J, Alfè D, Tkatchenko A. Nanoscale π-π stacked molecules are bound by collective charge fluctuations. Nat Commun 2017; 8:14052. [PMID: 28169280 PMCID: PMC5309697 DOI: 10.1038/ncomms14052] [Citation(s) in RCA: 60] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2016] [Accepted: 11/15/2016] [Indexed: 11/24/2022] Open
Abstract
Non-covalent π−π interactions are central to chemical and biological processes, yet the full understanding of their origin that would unite the simplicity of empirical approaches with the accuracy of quantum calculations is still missing. Here we employ a quantum-mechanical Hamiltonian model for van der Waals interactions, to demonstrate that intermolecular electron correlation in large supramolecular complexes at equilibrium distances is appropriately described by collective charge fluctuations. We visualize these fluctuations and provide connections both to orbital-based approaches to electron correlation, as well as to the simple London pairwise picture. The reported binding energies of ten supramolecular complexes obtained from the quantum-mechanical fluctuation model joined with density functional calculations are within 5% of the reference energies calculated with the diffusion quantum Monte-Carlo method. Our analysis suggests that π−π stacking in supramolecular complexes can be characterized by strong contributions to the binding energy from delocalized, collective charge fluctuations—in contrast to complexes with other types of bonding. Attractive, non-covalent interactions between aromatic rings—termed π−π stacking—is common in chemistry but difficult to model. Here the authors report a quantum-mechanical model to show the importance of collective charge fluctuations for understanding pi-stacked supramolecular systems.
Collapse
Affiliation(s)
- Jan Hermann
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Dario Alfè
- Department of Earth Sciences, University College London, London WC1E 6BT, UK.,Department of Physics and Astronomy, University College London, London WC1E 6BT, UK.,London Centre for Nanotechnology and Thomas Young Centre@UCL, University College London, London WC1E 6BT, UK
| | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany.,Physics and Materials Science Research Unit, University of Luxembourg, 162a Avenue de la Faiencerie, Luxembourg L-1511, Luxembourg
| |
Collapse
|
37
|
Hoja J, Reilly AM, Tkatchenko A. First-principles modeling of molecular crystals: structures and stabilities, temperature and pressure. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2016. [DOI: 10.1002/wcms.1294] [Citation(s) in RCA: 91] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Johannes Hoja
- Fritz-Haber-Institut der Max-Planck-Gesellschaft; Berlin Germany
- Physics and Materials Science Research Unit; University of Luxembourg; Luxembourg City Luxembourg
| | | | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft; Berlin Germany
- Physics and Materials Science Research Unit; University of Luxembourg; Luxembourg City Luxembourg
| |
Collapse
|
38
|
|
39
|
Hostaš J, Sigwalt D, Šekutor M, Ajani H, Dubecký M, Řezáč J, Zavalij PY, Cao L, Wohlschlager C, Mlinarić-Majerski K, Isaacs L, Glaser R, Hobza P. A Nexus between Theory and Experiment: Non-Empirical Quantum Mechanical Computational Methodology Applied to Cucurbit[n
]uril⋅Guest Binding Interactions. Chemistry 2016; 22:17226-17238. [DOI: 10.1002/chem.201601833] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2016] [Indexed: 12/25/2022]
Affiliation(s)
- Jiří Hostaš
- Institute of Organic Chemistry and Biochemistry; Flemingovo náměstí 2 Prague 6 160 10 Czech Republic
- Department of Physical and Macromolecular Chemistry; Faculty of Science; Charles University in Prague; Albertov 6 12843 Prague Czech Republic
| | - David Sigwalt
- Department of Chemistry and Biochemistry; University of Maryland; College Park MD 20742 USA
| | - Marina Šekutor
- Department of Organic Chemistry and Biochemistry; Ruđer Bošković Institute; Bijenička cesta 54 10000 Zagreb Croatia
| | - Haresh Ajani
- Institute of Organic Chemistry and Biochemistry; Flemingovo náměstí 2 Prague 6 160 10 Czech Republic
- Regional Centre of Advanced Technologies and Materials; Department of Physical Chemistry; Palacký University; 77146 Olomouc Czech Republic
| | - Matúš Dubecký
- Regional Centre of Advanced Technologies and Materials; Department of Physical Chemistry; Palacký University; 77146 Olomouc Czech Republic
| | - Jan Řezáč
- Institute of Organic Chemistry and Biochemistry; Flemingovo náměstí 2 Prague 6 160 10 Czech Republic
| | - Peter Y. Zavalij
- Department of Chemistry and Biochemistry; University of Maryland; College Park MD 20742 USA
| | - Liping Cao
- Department of Chemistry and Biochemistry; University of Maryland; College Park MD 20742 USA
| | | | - Kata Mlinarić-Majerski
- Department of Organic Chemistry and Biochemistry; Ruđer Bošković Institute; Bijenička cesta 54 10000 Zagreb Croatia
| | - Lyle Isaacs
- Department of Chemistry and Biochemistry; University of Maryland; College Park MD 20742 USA
| | - Robert Glaser
- Institut für Organische Chemie; Johannes Kepler Universität Linz; 4040 Linz Austria), Department of Chemistry, Ben-Gurion University of the Negev, Beer-Sheva 84105 (Israel
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry; Flemingovo náměstí 2 Prague 6 160 10 Czech Republic
- Regional Centre of Advanced Technologies and Materials; Department of Physical Chemistry; Palacký University; 77146 Olomouc Czech Republic
| |
Collapse
|
40
|
Reeves KG, Yao Y, Kanai Y. Diffusion quantum Monte Carlo study of martensitic phase transition energetics: The case of phosphorene. J Chem Phys 2016; 145:124705. [DOI: 10.1063/1.4962759] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Kyle G. Reeves
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| |
Collapse
|
41
|
Azadi S, Cohen RE. Chemical accuracy from quantum Monte Carlo for the benzene dimer. J Chem Phys 2016; 143:104301. [PMID: 26374029 DOI: 10.1063/1.4930137] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report an accurate study of interactions between benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory using different van der Waals functionals. In our quantum Monte Carlo (QMC) calculations, we use accurate correlated trial wave functions including three-body Jastrow factors and backflow transformations. We consider two benzene molecules in the parallel displaced geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the coupled-cluster theory through perturbative triplets/complete basis set limit is -2.65(2) kcal/mol [Miliordos et al., J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.
Collapse
Affiliation(s)
- Sam Azadi
- Department of Earth Science and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - R E Cohen
- London Centre for Nanotechnology, University College London, London WC1E 6BT, United Kingdom and Extreme Materials Initiative, Geophysical Laboratory, Carnegie Institution of Washington, Washington, D.C. 20015, USA
| |
Collapse
|
42
|
Affiliation(s)
- Matúš Dubecký
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lubos Mitas
- Department
of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Petr Jurečka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
43
|
Frey JA, Holzer C, Klopper W, Leutwyler S. Experimental and Theoretical Determination of Dissociation Energies of Dispersion-Dominated Aromatic Molecular Complexes. Chem Rev 2016; 116:5614-41. [DOI: 10.1021/acs.chemrev.5b00652] [Citation(s) in RCA: 53] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jann A. Frey
- Departement
für Chemie und Biochemie, Universität Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Christof Holzer
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, D-76131 Karlsruhe, Germany
| | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, D-76131 Karlsruhe, Germany
| | - Samuel Leutwyler
- Departement
für Chemie und Biochemie, Universität Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| |
Collapse
|
44
|
Gillan MJ, Alfè D, Michaelides A. Perspective: How good is DFT for water? J Chem Phys 2016; 144:130901. [DOI: 10.1063/1.4944633] [Citation(s) in RCA: 478] [Impact Index Per Article: 59.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Michael J. Gillan
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Dario Alfè
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
- Department of Earth Sciences, University College London, London WC1E 6BT, United Kingdom
| | - Angelos Michaelides
- London Centre for Nanotechnology, Gordon St., London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, London WC1H 0AH, United Kingdom
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
45
|
Trouillas P, Sancho-García JC, De Freitas V, Gierschner J, Otyepka M, Dangles O. Stabilizing and Modulating Color by Copigmentation: Insights from Theory and Experiment. Chem Rev 2016; 116:4937-82. [PMID: 26959943 DOI: 10.1021/acs.chemrev.5b00507] [Citation(s) in RCA: 286] [Impact Index Per Article: 35.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Natural anthocyanin pigments/dyes and phenolic copigments/co-dyes form noncovalent complexes, which stabilize and modulate (in particular blue, violet, and red) colors in flowers, berries, and food products derived from them (including wines, jams, purees, and syrups). This noncovalent association and their electronic and optical implications constitute the copigmentation phenomenon. Over the past decade, experimental and theoretical studies have enabled a molecular understanding of copigmentation. This review revisits this phenomenon to provide a comprehensive description of the nature of binding (the dispersion and electrostatic components of π-π stacking, the hydrophobic effect, and possible hydrogen-bonding between pigment and copigment) and of spectral modifications occurring in copigmentation complexes, in which charge transfer plays an important role. Particular attention is paid to applications of copigmentation in food chemistry.
Collapse
Affiliation(s)
- Patrick Trouillas
- INSERM UMR 850, Univ. Limoges , Faculty of Pharmacy, 2 rue du Dr. Marcland, F-87025 Limoges, France.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc , tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Juan C Sancho-García
- Departamento de Química Física, Universidad de Alicante , Apartado de Correos 99, E-03080 Alicante, Spain
| | - Victor De Freitas
- REQUIMTE/LAQV - Research Unit, Faculty of Science, Porto University , Rua do Campo Alegre, 4169-007 Porto, Portugal
| | - Johannes Gierschner
- Madrid Institute for Advanced Studies - IMDEA Nanoscience , C/Faraday 9, Ciudad Universitaria de Cantoblanco, E-28049 Madrid, Spain
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc , tr. 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Olivier Dangles
- University of Avignon, INRA, UMR408 SQPOV , F-84000 Avignon, France
| |
Collapse
|
46
|
Řezáč J, Hobza P. Benchmark Calculations of Interaction Energies in Noncovalent Complexes and Their Applications. Chem Rev 2016; 116:5038-71. [DOI: 10.1021/acs.chemrev.5b00526] [Citation(s) in RCA: 281] [Impact Index Per Article: 35.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jan Řezáč
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
| | - Pavel Hobza
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Palacký University, 771 46 Olomouc, Czech Republic
| |
Collapse
|
47
|
Benali A, Shulenburger L, Krogel JT, Zhong X, Kent PRC, Heinonen O. Quantum Monte Carlo analysis of a charge ordered insulating antiferromagnet: the Ti4O7 Magnéli phase. Phys Chem Chem Phys 2016; 18:18323-35. [DOI: 10.1039/c6cp02067d] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Spin density surfaces in the low temperature phase of Ti4O7: the ferromagnetic state (left panel), and the lowest-energy antiferromagnetic state (right panel). The figures were generated using self-interaction corrected density functional theory.
Collapse
Affiliation(s)
- Anouar Benali
- Argonne Leadership Computing Facility
- Argonne National Laboratory
- Argonne
- USA
| | | | - Jaron T. Krogel
- Materials Science and Technology Division
- Oak Ridge National Laboratory
- Oak Ridge
- USA
| | - Xiaoliang Zhong
- Material Science Division
- Argonne National Laboratory
- Argonne
- USA
| | - Paul R. C. Kent
- Center for Nanophase Materials Sciences and Computer Science and Mathematics Division
- Oak Ridge National Laboratory
- Oak Ridge
- USA
| | - Olle Heinonen
- Material Science Division
- Argonne National Laboratory
- Argonne
- USA
- Northwestern-Argonne Institute for Science and Engineering
| |
Collapse
|
48
|
Pykal M, Jurečka P, Karlický F, Otyepka M. Modelling of graphene functionalization. Phys Chem Chem Phys 2016; 18:6351-72. [DOI: 10.1039/c5cp03599f] [Citation(s) in RCA: 168] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
This perspective describes the available theoretical methods and models for simulating graphene functionalization based on quantum and classical mechanics.
Collapse
Affiliation(s)
- Martin Pykal
- Regional Centre of Advanced Technologies and Materials
- Department of Physical Chemistry
- Faculty of Science
- Palacký University Olomouc
- 771 46 Olomouc
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials
- Department of Physical Chemistry
- Faculty of Science
- Palacký University Olomouc
- 771 46 Olomouc
| | - František Karlický
- Regional Centre of Advanced Technologies and Materials
- Department of Physical Chemistry
- Faculty of Science
- Palacký University Olomouc
- 771 46 Olomouc
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials
- Department of Physical Chemistry
- Faculty of Science
- Palacký University Olomouc
- 771 46 Olomouc
| |
Collapse
|
49
|
Gasperich K, Jordan KD. Diffusion Monte Carlo Study of the Parallel Displaced Form of the Benzene Dimer. ACS SYMPOSIUM SERIES 2016. [DOI: 10.1021/bk-2016-1234.ch007] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Affiliation(s)
- Kevin Gasperich
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Kenneth D. Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
50
|
Shulenburger L, Baczewski AD, Zhu Z, Guan J, Tománek D. The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus. NANO LETTERS 2015; 15:8170-5. [PMID: 26523860 DOI: 10.1021/acs.nanolett.5b03615] [Citation(s) in RCA: 120] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/14/2023]
Abstract
Sensitive dependence of the electronic structure on the number of layers in few-layer phosphorene raises a question about the true nature of the interlayer interaction in so-called "van der Waals (vdW) solids". We performed quantum Monte Carlo calculations and found that the interlayer interaction in bulk black phosphorus and related few-layer phosphorene is associated with a significant charge redistribution that is incompatible with purely dispersive forces and not captured by density functional theory calculations with different vdW corrected functionals. These findings confirm the necessity of more sophisticated treatment of nonlocal electron correlation in total energy calculations.
Collapse
Affiliation(s)
- L Shulenburger
- Sandia National Laboratories , Albuquerque, New Mexico 87185, United States
| | - A D Baczewski
- Sandia National Laboratories , Albuquerque, New Mexico 87185, United States
| | - Z Zhu
- Physics and Astronomy Department, Michigan State University , East Lansing, Michigan 48824, United States
| | - J Guan
- Physics and Astronomy Department, Michigan State University , East Lansing, Michigan 48824, United States
| | - D Tománek
- Physics and Astronomy Department, Michigan State University , East Lansing, Michigan 48824, United States
| |
Collapse
|