Sebastian A Ramos, Leonard J Mueller, Gregory J. O. Beran
Ab initio chemical shift prediction plays a central role in nuclear magnetic resonance (NMR) crystallography, and the accuracy with which chemical shifts can be predicted relative to experiment impacts the confidence with which structures can be assigned. For organic crystals, periodic density functional theory calculations with the gauge-including projector augmented wave (GIPAW) approximation and the PBE functional are widely used at present. Many previous studies have examined how using more advanced density functionals can increase the accuracy of predicted chemical shifts relative to experiment, but nearly all of those studies employed crystal structures that were optimized with generalized-gradient approximation (GGA) functionals. Here, we investigate how the accuracy of the predicted chemical shifts in organic crystals is affected by replacing GGA-level PBE-D3(BJ) crystal geometries with more accurate hybrid functional PBE0-D3(BJ) ones. Based on benchmark data sets containing 132 13C and 35 15N chemical shifts plus case studies on testosterone, acetaminophen, and phenobarbital, we find that switching from GGA-level geometries and chemical shifts to hybrid functional ones reduces 13C and 15N chemical shift errors by ~40-60% versus experiment. However, most of the improvement stems from the use of the hybrid functional for the chemical shift calculations, rather than from the refined geometries. In addition, even with the improved geometries, we find that double-hybrid functionals still do not systematically increase chemical shift agreement with experiment beyond what hybrid functionals provide. In the end, these results suggest that the combination of GGA-level crystal structures and hybrid-functional chemical shifts represents a particularly cost-effective combination for NMR crystallography in organic systems.
{"title":"The interplay of density functional selection and crystal structure for accurate NMR chemical shift predictions","authors":"Sebastian A Ramos, Leonard J Mueller, Gregory J. O. Beran","doi":"10.1039/d4fd00072b","DOIUrl":"https://doi.org/10.1039/d4fd00072b","url":null,"abstract":"Ab initio chemical shift prediction plays a central role in nuclear magnetic resonance (NMR) crystallography, and the accuracy with which chemical shifts can be predicted relative to experiment impacts the confidence with which structures can be assigned. For organic crystals, periodic density functional theory calculations with the gauge-including projector augmented wave (GIPAW) approximation and the PBE functional are widely used at present. Many previous studies have examined how using more advanced density functionals can increase the accuracy of predicted chemical shifts relative to experiment, but nearly all of those studies employed crystal structures that were optimized with generalized-gradient approximation (GGA) functionals. Here, we investigate how the accuracy of the predicted chemical shifts in organic crystals is affected by replacing GGA-level PBE-D3(BJ) crystal geometries with more accurate hybrid functional PBE0-D3(BJ) ones. Based on benchmark data sets containing 132 <small><sup>13</sup></small>C and 35 <small><sup>15</sup></small>N chemical shifts plus case studies on testosterone, acetaminophen, and phenobarbital, we find that switching from GGA-level geometries and chemical shifts to hybrid functional ones reduces <small><sup>13</sup></small>C and <small><sup>15</sup></small>N chemical shift errors by ~40-60% versus experiment. However, most of the improvement stems from the use of the hybrid functional for the chemical shift calculations, rather than from the refined geometries. In addition, even with the improved geometries, we find that double-hybrid functionals still do not systematically increase chemical shift agreement with experiment beyond what hybrid functionals provide. In the end, these results suggest that the combination of GGA-level crystal structures and hybrid-functional chemical shifts represents a particularly cost-effective combination for NMR crystallography in organic systems.","PeriodicalId":76,"journal":{"name":"Faraday Discussions","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141061921","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Javier Valenzuela Reina, Federico Civaia, Angela F Harper, Christoph Scheurer, Simone Köcher
We present a comprehensive study on the best practices for integrating first principles simulations in experimental quadrupolar solid-state nuclear magnetic resonance (SS-NMR), exploiting the synergies between theory and experiment for achieving the optimal interpretation of both. Most high performance materials (HPM), such as battery electrodes, exhibit complex SS-NMR spectra due to dynamic effects or amorphous phases. NMR crystallography for such challenging materials requires reliable, accurate, efficient computational methods for calculating NMR observables from first principles for the transfer between material structure models and the interpretation of their experimental SS-NMR spectra. NMR-active nuclei within HPMs are routinely probed by their chemical shielding anisotropy (CSA). However, several nuclear isotopes of interest, e.g. 7Li and 27Al, have a nuclear quadrupole and experience additional interactions with the surrounding electric field gradient (EFG). The quadrupolar interaction is a valuable source of information about atomistic structure and in particular local symmetry, complementing the CSA. As such, there is a range of different methods and codes to choose from for calculating EFGs, from all-electron to plane wave methods. We benchmark the accuracy of different simulation strategies for computing the EFG tensor of quadrupolar nuclei with plane wave density functional theory (DFT) and study the impact of the material structure, and the details of the simulation strategy. Especially for small nuclei with few electrons, such as 7Li, we can show that the choice of physical approximations and simulation parameters has a large effect on the transferability of the simulations results. To the best of our knowledge, we present the first comprehensive reference scale and literature survey for 7Li quadrupolar couplings. The results allow us to establish practical guidelines for developing the best simulation strategy for correlating DFT to experimental data extracting the maximum benefit and information from both, thereby advancing further research into HPMs.
{"title":"The EFG Rosetta Stone: Translating between DFT calculations and solid state NMR experiments","authors":"Javier Valenzuela Reina, Federico Civaia, Angela F Harper, Christoph Scheurer, Simone Köcher","doi":"10.1039/d4fd00075g","DOIUrl":"https://doi.org/10.1039/d4fd00075g","url":null,"abstract":"We present a comprehensive study on the best practices for integrating first principles simulations in experimental quadrupolar solid-state nuclear magnetic resonance (SS-NMR), exploiting the synergies between theory and experiment for achieving the optimal interpretation of both. Most high performance materials (HPM), such as battery electrodes, exhibit complex SS-NMR spectra due to dynamic effects or amorphous phases. NMR crystallography for such challenging materials requires reliable, accurate, efficient computational methods for calculating NMR observables from first principles for the transfer between material structure models and the interpretation of their experimental SS-NMR spectra. NMR-active nuclei within HPMs are routinely probed by their chemical shielding anisotropy (CSA). However, several nuclear isotopes of interest, e.g. <small><sup>7</sup></small>Li and <small><sup>27</sup></small>Al, have a nuclear quadrupole and experience additional interactions with the surrounding electric field gradient (EFG). The quadrupolar interaction is a valuable source of information about atomistic structure and in particular local symmetry, complementing the CSA. As such, there is a range of different methods and codes to choose from for calculating EFGs, from all-electron to plane wave methods. We benchmark the accuracy of different simulation strategies for computing the EFG tensor of quadrupolar nuclei with plane wave density functional theory (DFT) and study the impact of the material structure, and the details of the simulation strategy. Especially for small nuclei with few electrons, such as <small><sup>7</sup></small>Li, we can show that the choice of physical approximations and simulation parameters has a large effect on the transferability of the simulations results. To the best of our knowledge, we present the first comprehensive reference scale and literature survey for <small><sup>7</sup></small>Li quadrupolar couplings. The results allow us to establish practical guidelines for developing the best simulation strategy for correlating DFT to experimental data extracting the maximum benefit and information from both, thereby advancing further research into HPMs.","PeriodicalId":76,"journal":{"name":"Faraday Discussions","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-05-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140935460","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Angela F Harper, Tabea Huss, Simone Koecher, Christoph Scheurer
We present for the first time a multiscale machine learning approach to jointly simulate atomic structure and dynamics with the corresponding solid state Nuclear Magnetic Resonance (ssNMR) observables. We study the use-case of spin-alignment echo (SAE) NMR for exploring Li-ion diffusion within the solid state electrolyte material Li3PS4 (LPS) by calculating quadrupolar frequencies of 7Li. SAE NMR probes long-range dynamics down to microsecond-timescale hopping processes. Therefore only a few machine learning force field schemes are able to capture the time- and length scales required for accurate comparison with experimental results. By using a new class of machine learning interatomic potentials, known as ultra-fast potentials (UFPs), we are able to efficiently access timescales beyond the microsecond regime. In tandem, we have developed a machine learning model for predicting the full 7Li electric field gradient (EFG) tensors in LPS. By combining the long timescale trajectories from the UFP with our model for 7Li EFG tensors, we are able to extract the autocorrelation function (ACF) for 7Li quadrupolar frequencies during Li diffusion. We extract the decay constants from the ACF for both crystalline β-LPS and amorphous LPS, and find that the predicted Li hopping rates are on the same order of magnitude as those predicted from the Li dynamics. This demonstrates the potential for machine learning to finally make predictions on experimentally relevant timescales and temperatures, and opens a new avenue of NMR crystallography: using machine learning dynamical NMR simulations for accessing polycrystalline and glass ceramic materials.
{"title":"Tracking Li atoms in real-time with ultra-fast NMR simulations","authors":"Angela F Harper, Tabea Huss, Simone Koecher, Christoph Scheurer","doi":"10.1039/d4fd00074a","DOIUrl":"https://doi.org/10.1039/d4fd00074a","url":null,"abstract":"We present for the first time a multiscale machine learning approach to jointly simulate atomic structure and dynamics with the corresponding solid state Nuclear Magnetic Resonance (ssNMR) observables. We study the use-case of spin-alignment echo (SAE) NMR for exploring Li-ion diffusion within the solid state electrolyte material Li<small><sub>3</sub></small>PS<small><sub>4</sub></small> (LPS) by calculating quadrupolar frequencies of <small><sup>7</sup></small>Li. SAE NMR probes long-range dynamics down to microsecond-timescale hopping processes. Therefore only a few machine learning force field schemes are able to capture the time- and length scales required for accurate comparison with experimental results. By using a new class of machine learning interatomic potentials, known as ultra-fast potentials (UFPs), we are able to efficiently access timescales beyond the microsecond regime. In tandem, we have developed a machine learning model for predicting the full <small><sup>7</sup></small>Li electric field gradient (EFG) tensors in LPS. By combining the long timescale trajectories from the UFP with our model for <small><sup>7</sup></small>Li EFG tensors, we are able to extract the autocorrelation function (ACF) for <small><sup>7</sup></small>Li quadrupolar frequencies during Li diffusion. We extract the decay constants from the ACF for both crystalline β-LPS and amorphous LPS, and find that the predicted Li hopping rates are on the same order of magnitude as those predicted from the Li dynamics. This demonstrates the potential for machine learning to finally make predictions on experimentally relevant timescales and temperatures, and opens a new avenue of NMR crystallography: using machine learning dynamical NMR simulations for accessing polycrystalline and glass ceramic materials.","PeriodicalId":76,"journal":{"name":"Faraday Discussions","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-05-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140833918","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Kalman Szenes, Maximilian Moerchen, Paul Fischill, Markus Reiher
Multi-configurational electronic structure theory delivers the most versatile approximations to many-electron wavefunctions, flexible enough to deal with all sorts of transformations, ranging from electronic excitations, to open-shell molecules and chemical reactions. Multi-configurational models are therefore essential to establish universally applicable, predictive ab initio methods for chemistry. Here, we present a discussion of explicit correlation approaches which address the nagging problem of dealing with static and dynamic electron correlation in multi-configurational approaches. We review the latest developments and then point to their key obstacles. Our discussion is supported by new data obtained with tensor network methods. We argue in favor of simple electrons-only correlator expressions that may allow one to define transcorrelated models, in which the correlator does not bear a dependence on molecular structure.
{"title":"Striking the right balance of encoding electron correlation in the Hamiltonian and wavefunction ansatz","authors":"Kalman Szenes, Maximilian Moerchen, Paul Fischill, Markus Reiher","doi":"10.1039/d4fd00060a","DOIUrl":"https://doi.org/10.1039/d4fd00060a","url":null,"abstract":"Multi-configurational electronic structure theory delivers the most versatile approximations to many-electron wavefunctions, flexible enough to deal with all sorts of transformations, ranging from electronic excitations, to open-shell molecules and chemical reactions. Multi-configurational models are therefore essential to establish universally applicable, predictive ab initio methods for chemistry. Here, we present a discussion of explicit correlation approaches which address the nagging problem of dealing with static and dynamic electron correlation in multi-configurational approaches. We review the latest developments and then point to their key obstacles. Our discussion is supported by new data obtained with tensor network methods. We argue in favor of simple electrons-only correlator expressions that may allow one to define transcorrelated models, in which the correlator does not bear a dependence on molecular structure.","PeriodicalId":76,"journal":{"name":"Faraday Discussions","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-05-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140833916","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Zijun Zhao, Maria-Andreea Filip, Alexander James William Thom
The multi-reference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multi-reference problem. These include a logarithmically scaling metric-tree based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants at essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.
{"title":"Rapidly convergent quantum Monte Carlo using a Chebyshev projector","authors":"Zijun Zhao, Maria-Andreea Filip, Alexander James William Thom","doi":"10.1039/d4fd00035h","DOIUrl":"https://doi.org/10.1039/d4fd00035h","url":null,"abstract":"The multi-reference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multi-reference problem. These include a logarithmically scaling metric-tree based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants at essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.","PeriodicalId":76,"journal":{"name":"Faraday Discussions","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-05-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140834088","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Kemal Atalar, Yannic Rath, Rachel Crespo-Otero, George Booth
We build on the concept of eigenvector continuation to develop an efficient multi-state method for the rigorous and smooth interpolation of a small training set of many-body wavefunctions through chemical space at mean-field cost. The inferred states are represented as variationally optimal linear combinations of the training states transferred between the many-body basis of different nuclear geometries. We show that analytic multi-state forces and nonadiabatic couplings from the model enable application to nonadiabatic molecular dynamics, developing an active learning scheme to ensure a compact and systematically improvable training set. This culminates in application to the nonadiabatic molecular dynamics of a photoexcited 28-atom hydrogen chain, with surprising complexity in the resulting nuclear motion. With just 22 DMRG calculations of training states from the low-energy correlated electronic structure at different geometries, we infer the multi-state energies, forces and nonadiabatic coupling vectors at 12,000 geometries with provable convergence to high accuracy along an ensemble of molecular trajectories, which would not be feasible with a brute force approach. This opens up a route to bridge the timescales between accurate single-point correlated electronic structure methods and timescales of relevance for photo-induced molecular dynamics.
{"title":"Fast and accurate nonadiabatic molecular dynamics enabled through variational interpolation of correlated electron wavefunctions","authors":"Kemal Atalar, Yannic Rath, Rachel Crespo-Otero, George Booth","doi":"10.1039/d4fd00062e","DOIUrl":"https://doi.org/10.1039/d4fd00062e","url":null,"abstract":"We build on the concept of eigenvector continuation to develop an efficient multi-state method for the rigorous and smooth interpolation of a small training set of many-body wavefunctions through chemical space at mean-field cost. The inferred states are represented as variationally optimal linear combinations of the training states transferred between the many-body basis of different nuclear geometries. We show that analytic multi-state forces and nonadiabatic couplings from the model enable application to nonadiabatic molecular dynamics, developing an active learning scheme to ensure a compact and systematically improvable training set. This culminates in application to the nonadiabatic molecular dynamics of a photoexcited 28-atom hydrogen chain, with surprising complexity in the resulting nuclear motion. With just 22 DMRG calculations of training states from the low-energy correlated electronic structure at different geometries, we infer the multi-state energies, forces and nonadiabatic coupling vectors at 12,000 geometries with provable convergence to high accuracy along an ensemble of molecular trajectories, which would not be feasible with a brute force approach. This opens up a route to bridge the timescales between accurate single-point correlated electronic structure methods and timescales of relevance for photo-induced molecular dynamics.","PeriodicalId":76,"journal":{"name":"Faraday Discussions","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-05-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140842277","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Lujo Matasovic, Hugo Bronstein, Richard Friend, Felix Plasser
Diradicals are of high current interest as emerging materials for next generation optoelectronic applications. To tune their excited-state properties it would be greatly beneficial to have a detailed understanding of the wave functions of the different states involved but this endeavour is hampered by formal and practical barriers. To tackle these challenges, we present a formal analysis as well as concrete results on diradical excited states. We start with a detailed investigation of the avail- able states of a two-orbital two-electron model viewed from, both, the valence-bond and molecular orbital perspectives. We highlight the presence of diradical and zwitterionic states and illustrate their connections to the states found in closed-shell molecules. Subsequently, we introduce practi- cal protocols for analysing states from realistic multireference computations applying these to the paraquinodimethane (pQDM) molecule. The analysis reveals four different categories of states – diradical, zwitterionic, HOMO-SOMO as well as biexciton – while also providing insight into their energetics and optical properties. Twisting the CH2 groups allows us to interconvert between the closed and open-shell forms of pQDM illustrating the connection between the states in both forms. More generally, we hope that this work will lay the foundations for a more powerful rational design approach to diradicals for photophysical applications.
{"title":"Classification and quantitative characterisation of the excited states of π-conjugated diradicals","authors":"Lujo Matasovic, Hugo Bronstein, Richard Friend, Felix Plasser","doi":"10.1039/d4fd00055b","DOIUrl":"https://doi.org/10.1039/d4fd00055b","url":null,"abstract":"Diradicals are of high current interest as emerging materials for next generation optoelectronic applications. To tune their excited-state properties it would be greatly beneficial to have a detailed understanding of the wave functions of the different states involved but this endeavour is hampered by formal and practical barriers. To tackle these challenges, we present a formal analysis as well as concrete results on diradical excited states. We start with a detailed investigation of the avail- able states of a two-orbital two-electron model viewed from, both, the valence-bond and molecular orbital perspectives. We highlight the presence of diradical and zwitterionic states and illustrate their connections to the states found in closed-shell molecules. Subsequently, we introduce practi- cal protocols for analysing states from realistic multireference computations applying these to the paraquinodimethane (pQDM) molecule. The analysis reveals four different categories of states – diradical, zwitterionic, HOMO-SOMO as well as biexciton – while also providing insight into their energetics and optical properties. Twisting the CH2 groups allows us to interconvert between the closed and open-shell forms of pQDM illustrating the connection between the states in both forms. More generally, we hope that this work will lay the foundations for a more powerful rational design approach to diradicals for photophysical applications.","PeriodicalId":76,"journal":{"name":"Faraday Discussions","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-04-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140797788","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We present the first application to real molecular systems of the recently proposed linear-response theory for the density-based basis-set correction method [J. Chem. Phys. {bf 158}, 234107 (2023)]. We apply this approach to accelerate the basis-set convergence of excitation energies in the equation-of-motion coupled-cluster singles doubles (EOM-CCSD) method. We use an approximate linear-response framework which neglects the second-order derivative of the basis-set correction density functional and consists in simply adding to the usual Hamiltonian the one-electron potential generated by the first-order derivative of the functional. This additional basis-set correction potential is evaluated at the Hartree-Fock density, leading to a very computationally cheap basis-set correction. We tested this approach over a set of about 30 excitation energies computed for five small molecular systems and found that the excitation energies from the ground state to Rydberg states are the main source of basis-set error. These excitation energies systematically increase when the size of the basis set is increased, suggesting a biased description in favour of the excited state. Despite the simplicity of the present approach, the results obtained with the basis-set corrected EOM-CCSD method are encouraging as they yield to a mean absolute deviation of 0.02 eV for the aug-cc-pVTZ basis set, while it is of 0.04 eV using the standard EOM-CCSD method. This might open the path to an alternative to explicitly correlated approaches to accelerate the basis-set convergence of excitation energies.
{"title":"Accelerated basis-set convergence of coupled-cluster excitation energies using the density-based basis-set correction method","authors":"Diata Traoré, Julien Toulouse, Emmanuel Giner","doi":"10.1039/d4fd00033a","DOIUrl":"https://doi.org/10.1039/d4fd00033a","url":null,"abstract":"We present the first application to real molecular systems of the recently proposed linear-response theory for the density-based basis-set correction method [J. Chem. Phys. {bf 158}, 234107 (2023)]. We apply this approach to accelerate the basis-set convergence of excitation energies in the equation-of-motion coupled-cluster singles doubles (EOM-CCSD) method. We use an approximate linear-response framework which neglects the second-order derivative of the basis-set correction density functional and consists in simply adding to the usual Hamiltonian the one-electron potential generated by the first-order derivative of the functional. This additional basis-set correction potential is evaluated at the Hartree-Fock density, leading to a very computationally cheap basis-set correction. We tested this approach over a set of about 30 excitation energies computed for five small molecular systems and found that the excitation energies from the ground state to Rydberg states are the main source of basis-set error. These excitation energies systematically increase when the size of the basis set is increased, suggesting a biased description in favour of the excited state. Despite the simplicity of the present approach, the results obtained with the basis-set corrected EOM-CCSD method are encouraging as they yield to a mean absolute deviation of 0.02 eV for the aug-cc-pVTZ basis set, while it is of 0.04 eV using the standard EOM-CCSD method. This might open the path to an alternative to explicitly correlated approaches to accelerate the basis-set convergence of excitation energies.","PeriodicalId":76,"journal":{"name":"Faraday Discussions","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-04-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140797755","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Nicolas Scaglione, Jocasta Avila, Agilio A.H. Padua, Margarida Costa Gomes
We used a library of thermally stable tetraalkylphosphonium carboxylate ILs that were easily pre- pared from available carboxylic acids. Depending on the 𝑝𝐾𝑎 of the precursor acids, the resulting ionic liquids either dissolve or reversibly chemically absorb CO2, with some exhibiting notable gas capacities, reaching up to 0.2 mole fraction of CO2 at 1 bar and 343 K. While equilibrium con- stants and ionic liquid capacities generally correlate with the 𝑝𝐾𝑎 of the acids, certain exceptions underscore the influence of liquid structure and physical properties, elucidated through molecular dynamics simulations and density functional theory (DFT) calculations. Unlike trends observed in other CO2-absorbing ILs, phosphonium carboxylates do not experience increased viscosity upon CO2 absorption; instead, enhanced diffusivities are observed, facilitating efficient gas-liquid trans- fer and reducing energy requirements in carbon capture processes.
我们使用了一个热稳定性四烷基膦羧酸盐离子液体库,它很容易从现有的羧酸中预制出来。根据前体酸𝑝𝑎的不同,得到的离子液体可以溶解或可逆化学吸收二氧化碳,其中一些离子液体表现出显著的气体容量,在 1 巴和 343 K 条件下,二氧化碳的摩尔分数高达 0.2。虽然平衡态和离子液体容量通常与酸的𝑝𝑎相关,但某些例外情况强调了液体结构和物理性质的影响,分子动力学模拟和密度泛函理论(DFT)计算阐明了这一点。与在其他吸收二氧化碳的液体中观察到的趋势不同,膦羧酸盐在吸收二氧化碳时粘度并没有增加;相反,观察到的扩散性增强,从而促进了气液的高效传输,降低了碳捕获过程中的能量需求。
{"title":"Tailored Carbon Dioxide Capacity in Carboxylate-Based Ionic Liquids†","authors":"Nicolas Scaglione, Jocasta Avila, Agilio A.H. Padua, Margarida Costa Gomes","doi":"10.1039/d4fd00052h","DOIUrl":"https://doi.org/10.1039/d4fd00052h","url":null,"abstract":"We used a library of thermally stable tetraalkylphosphonium carboxylate ILs that were easily pre- pared from available carboxylic acids. Depending on the 𝑝𝐾𝑎 of the precursor acids, the resulting ionic liquids either dissolve or reversibly chemically absorb CO2, with some exhibiting notable gas capacities, reaching up to 0.2 mole fraction of CO2 at 1 bar and 343 K. While equilibrium con- stants and ionic liquid capacities generally correlate with the 𝑝𝐾𝑎 of the acids, certain exceptions underscore the influence of liquid structure and physical properties, elucidated through molecular dynamics simulations and density functional theory (DFT) calculations. Unlike trends observed in other CO2-absorbing ILs, phosphonium carboxylates do not experience increased viscosity upon CO2 absorption; instead, enhanced diffusivities are observed, facilitating efficient gas-liquid trans- fer and reducing energy requirements in carbon capture processes.","PeriodicalId":76,"journal":{"name":"Faraday Discussions","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-04-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140797730","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In the domain of exchange-coupled PNTM clusters, local emergent symmetries exist which can be exploited to greatly increase the sparsity of the CI eigensolutions of such systems. Sparsity of the CI secular problem is revealed by exploring the site permutation space within spin-adapted many-body bases, and highly compressed wave functions may arise by finding optimal site orderings. However, the factorial cost of searching through the permutation space remains a bottleneck for clusters with a large number of metal centers. In this work, we explore ways to reduce the factorial scaling, by combining permutation and point group symmetry arguments, and using commutation relations between cumulative partial spin and the Hamiltonian operators, [(S(n))2, H]. Certain site orderings lead to commuting operators, from which more sparse wave functions arise. Two graphical strategies will be discussed, one to rapidly evaluate the commutators of interest, and one in the form of a tree search algorithm to predict how many and which distinct site permutations are to be analyzed, eliminating redundancies in the permutation space. Particularly interesting is the case of the singlet spin states for which an additional reversal symmetry can be utilized to further reduce the distinct site permutations.
在交换耦合 PNTM 簇领域,存在着局部新兴对称性,利用这些对称性可以大大提高此类系统的 CI 特征解的稀疏性。通过探索自旋适配多体基内的位点排列空间,可以揭示 CI 世俗问题的稀疏性;通过寻找最佳位点排序,可以产生高度压缩的波函数。然而,对于具有大量金属中心的集群来说,搜索排列空间的因子成本仍然是一个瓶颈。在这项工作中,我们结合了置换和点群对称性论证,并利用累积部分自旋和哈密顿算子[(S(n))2, H]之间的换向关系,探索了减少因子缩放的方法。某些位点排序会导致换向算子,从而产生更稀疏的波函数。我们将讨论两种图形策略,一种是快速评估感兴趣的换向算子,另一种是以树形搜索算法的形式预测有多少种不同的位点排列需要分析,从而消除排列空间中的冗余。特别有趣的是单子自旋态,可以利用额外的反转对称性来进一步减少不同位点的排列。
{"title":"Permutation Symmetry in Spin Adapted Many-Body Wave Functions","authors":"Maru Song, Ali Alavi, Giovanni Li Manni","doi":"10.1039/d4fd00061g","DOIUrl":"https://doi.org/10.1039/d4fd00061g","url":null,"abstract":"In the domain of exchange-coupled PNTM clusters, local emergent symmetries exist which can be exploited to greatly increase the sparsity of the CI eigensolutions of such systems. Sparsity of the CI secular problem is revealed by exploring the site permutation space within spin-adapted many-body bases, and highly compressed wave functions may arise by finding optimal site orderings. However, the factorial cost of searching through the permutation space remains a bottleneck for clusters with a large number of metal centers. In this work, we explore ways to reduce the factorial scaling, by combining permutation and point group symmetry arguments, and using commutation relations between cumulative partial spin and the Hamiltonian operators, [(S<small><sup>(n)</sup></small>)<small><sup>2</sup></small>, H]. Certain site orderings lead to commuting operators, from which more sparse wave functions arise. Two graphical strategies will be discussed, one to rapidly evaluate the commutators of interest, and one in the form of a tree search algorithm to predict how many and which distinct site permutations are to be analyzed, eliminating redundancies in the permutation space. Particularly interesting is the case of the singlet spin states for which an additional reversal symmetry can be utilized to further reduce the distinct site permutations.","PeriodicalId":76,"journal":{"name":"Faraday Discussions","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-04-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140797771","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}