1
|
Mandalia R, Fauser S, Trushin E, Görling A. Assessment of RPA and σ-functional methods for the calculation of dipole moments and static polarizabilities and hyperpolarizabilities. J Chem Phys 2025; 162:184106. [PMID: 40341926 DOI: 10.1063/5.0267912] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2025] [Accepted: 04/16/2025] [Indexed: 05/11/2025] Open
Abstract
In the present paper, we assess the performance of methods based on the random phase approximation (RPA) and on σ-functionals for predicting static optical properties, i.e., dipole moment, polarizability, and first and second hyperpolarizability, of small- and medium-sized molecules, including chain-like systems. First, we provide accurate reference data by coupled-cluster singles, doubles, with perturbative triples calculations with sufficiently large basis sets. The RPA and σ-functional calculations are carried out post-self-consistently using input orbitals and eigenvalues from the hybrid density-functional calculation. The optimal fraction of exact non-local exchange in these calculations is found to be quite high, around 0.5-0.6 in RPA and around 0.8-1.0 in σ-functional methods. σ-functional methods, however, proved to be less sensitive than RPA methods with respect to the amount of exact non-local exchange used in the generation of their input data. σ-functional methods are shown to outperform in accuracy RPA methods and various other considered density-functional theory methods for static optical properties and, thus, are well-suited for the calculation of linear and non-linear optical properties.
Collapse
Affiliation(s)
- Raviraj Mandalia
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Steffen Fauser
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Egor Trushin
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
- Erlangen National High Performance Computing Center (NHR@FAU), Martensstr. 1, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
- Erlangen National High Performance Computing Center (NHR@FAU), Martensstr. 1, D-91058 Erlangen, Germany
| |
Collapse
|
2
|
Raposo-Hernández G, Pappalardo RR, Réal F, Vallet V, Sánchez Marcos E. Toward a realistic theoretical electronic spectra of metal aqua ions in solution: The case of Ce(H2O)n3+ using statistical methods and quantum chemistry calculations. J Chem Phys 2024; 161:144109. [PMID: 39387406 DOI: 10.1063/5.0228155] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2024] [Accepted: 09/24/2024] [Indexed: 10/15/2024] Open
Abstract
Accurately predicting spectra for heavy elements, often open-shell systems, is a significant challenge typically addressed using a single cluster approach with a fixed coordination number. Developing a realistic model that accounts for temperature effects, variable coordination numbers, and interprets experimental data is even more demanding due to the strong solute-solvent interactions present in solutions of heavy metal cations. This study addresses these challenges by combining multiple methodologies to accurately predict realistic spectra for highly charged metal cations in aqueous media, with a focus on the electronic absorption spectrum of Ce3+ in water. Utilizing highly correlated relativistic quantum mechanical (QM) wavefunctions and structures from molecular dynamics (MD) simulations, we show that the convolution of individual vertical transitions yields excellent agreement with experimental results without the introduction of empirical broadening. Good results are obtained for both the normalized spectrum and that of absolute intensity. The study incorporates a statistical machine learning algorithm, Gaussian Mixture Models-Nuclear Ensemble Approach (GMM-NEA), to convolute individual spectra. The microscopic distribution provided by MD simulations allows us to examine the contributions of the octa- and ennea-hydrate of Ce3+ in water to the final spectrum. In addition, the temperature dependence of the spectrum is theoretically captured by observing the changing population of these hydrate forms with temperature. We also explore an alternative method for obtaining statistically representative structures in a less demanding manner than MD simulations, derived from QM Wigner distributions. The combination of Wigner-sampling and GMM-NEA broadening shows promise for wide application in spectroscopic analysis and predictions, offering a computationally efficient alternative to traditional methods.
Collapse
Affiliation(s)
| | - Rafael R Pappalardo
- Department of Physical Chemistry, University of Seville, 41012 Seville, Spain
| | - Florent Réal
- Université de Lille, CNRS, UMR 8523-PhLAM, Physique des Lasers, Atomes et Molecules, F-59000 Lille, France
| | - Valérie Vallet
- Université de Lille, CNRS, UMR 8523-PhLAM, Physique des Lasers, Atomes et Molecules, F-59000 Lille, France
| | | |
Collapse
|
3
|
Nagy PR. State-of-the-art local correlation methods enable affordable gold standard quantum chemistry for up to hundreds of atoms. Chem Sci 2024:d4sc04755a. [PMID: 39246365 PMCID: PMC11376132 DOI: 10.1039/d4sc04755a] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2024] [Accepted: 07/30/2024] [Indexed: 09/10/2024] Open
Abstract
In this feature, we review the current capabilities of local electron correlation methods up to the coupled cluster model with single, double, and perturbative triple excitations [CCSD(T)], which is a gold standard in quantum chemistry. The main computational aspects of the local method types are assessed from the perspective of applications, but the focus is kept on how to achieve chemical accuracy (i.e., <1 kcal mol-1 uncertainty), as well as on the broad scope of chemical problems made accessible. The performance of state-of-the-art methods is also compared, including the most employed DLPNO and, in particular, our local natural orbital (LNO) CCSD(T) approach. The high accuracy and efficiency of the LNO method makes chemically accurate CCSD(T) computations accessible for molecules of hundreds of atoms with resources affordable to a broad computational community (days on a single CPU and 10-100 GB of memory). Recent developments in LNO-CCSD(T) enable systematic convergence and robust error estimates even for systems of complicated electronic structure or larger size (up to 1000 atoms). The predictive power of current local CCSD(T) methods, usually at about 1-2 order of magnitude higher cost than hybrid density functional theory (DFT), has become outstanding on the palette of computational chemistry applicable for molecules of practical interest. We also review more than 50 LNO-based and other advanced local-CCSD(T) applications for realistic, large systems across molecular interactions as well as main group, transition metal, bio-, and surface chemistry. The examples show that properly executed local-CCSD(T) can contribute to binding, reaction equilibrium, rate constants, etc. which are able to match measurements within the error estimates. These applications demonstrate that modern, open-access, and broadly affordable local methods, such as LNO-CCSD(T), already enable predictive computations and atomistic insight for complicated, real-life molecular processes in realistic environments.
Collapse
Affiliation(s)
- Péter R Nagy
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics Műegyetem rkp. 3. H-1111 Budapest Hungary
- HUN-REN-BME Quantum Chemistry Research Group Műegyetem rkp. 3. H-1111 Budapest Hungary
- MTA-BME Lendület Quantum Chemistry Research Group Műegyetem rkp. 3. H-1111 Budapest Hungary
| |
Collapse
|
4
|
Kaka KS, Castet F, Champagne B. On the third-order nonlinear optical responses of cis and trans stilbenes - a quantum chemistry investigation. Phys Chem Chem Phys 2024; 26:14808-14824. [PMID: 38717796 DOI: 10.1039/d4cp00522h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/23/2024]
Abstract
The second hyperpolarizabilities (γ) of the stilbene molecular switch in its trans and cis forms have been calculated using quantum chemistry methods to address their third-order nonlinear optical contrasts, to assess the reliability of lower-cost DFT methods, and to make comparisons with experiments. First, the reference CCSD(T) method shows that trans-stilbene presents a γ‖ value twice larger than its cis isomer (its γTHS value is 2.7 times larger). Among more cost-effective methods, reliable results are obtained at MP2 as well as with DFT, provided the CAM-B3LYP or ωB97X-D XCFs are employed. Supplementary DFT calculations have investigated the relationships between the accuracy of the exchange-correlation functionals, the fulfillment of Koopmans' theorem, and the delocalization error, and they demonstrated that satisfying Koopmans' theorem is not the condition for the best accuracy but that functionals with small delocalization errors are generally efficient. Using the selected CAM-B3LYP, large γ enhancements by about 70% (trans-stilbene) and 50% (cis-stilbene) have been evidenced when accounting for solvent effects using an implicit solvation model (IEFPCM), even for apolar solvents. Then, the frequency dispersion of the γ responses has been described using Bishop polynomial expansions, allowing comparisons with a broad set of experimental data. To a certain extent, no systematic agreement between the calculations and the measured values was found. On the one hand, the agreement is satisfactory for the γ(-ω;ω,-ω,ω) quantities, provided that the dominant vibrational contribution is taken into account. On the other hand, the agreement is poor for the γ(-2ω;ω,ω,0) and γ(-3ω;ω,ω,ω) quantities, while some inconsistencies between experimental values are also highlighted.
Collapse
Affiliation(s)
- Komlanvi Sèvi Kaka
- Theoretical Chemistry Laboratory, Unit of Theoretical and Structural Physical Chemistry, NISM (Namur Institute of Structured Matter), University of Namur (UNamur), B-5000 Namur, Belgium.
| | - Frédéric Castet
- University of Bordeaux, CNRS, Bordeaux INP, ISM, UMR 5255, F-33405 Cedex Talence, France
| | - Benoît Champagne
- Theoretical Chemistry Laboratory, Unit of Theoretical and Structural Physical Chemistry, NISM (Namur Institute of Structured Matter), University of Namur (UNamur), B-5000 Namur, Belgium.
| |
Collapse
|
5
|
Rodríguez-Mayorga M, Besalú-Sala P, Pérez-Jiménez ÁJ, Sancho-García JC. Application to nonlinear optical properties of the RSX-QIDH double-hybrid range-separated functional. J Comput Chem 2024; 45:995-1001. [PMID: 38206899 DOI: 10.1002/jcc.27302] [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: 11/10/2023] [Revised: 12/19/2023] [Accepted: 12/21/2023] [Indexed: 01/13/2024]
Abstract
The effective calculation of static nonlinear optical properties requires a considerably high accuracy at a reasonable computational cost, to tackle challenging organic and inorganic systems acting as precursors and/or active layers of materials in (nano-)devices. That trade-off implies to obtain very accurate electronic energies in the presence of externally applied electric fields to consequently obtain static polarizabilities (α i j ) and hyper-polarizabilities (β i j k andγ i j k l ). Density functional theory is known to provide an excellent compromise between accuracy and computational cost, which is however largely impeded for these properties without introducing range-separation techniques. We thus explore here the ability of a modern (double-hybrid and range-separated) Range-Separated eXchange Quadratic Integrand Double-Hybrid exchange-correlation functional to compete in accuracy with more costly and/or tuned methods, thanks to its robust and parameter-free nature.
Collapse
Affiliation(s)
- M Rodríguez-Mayorga
- Department of Physical Chemistry, University of Alicante, Alicante, Spain
- Université Grenoble Alpes, CNRS, Inst. NÉEL, Grenoble, France
| | - P Besalú-Sala
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute for Molecular and Life Sciences (AIMMS), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Á J Pérez-Jiménez
- Department of Physical Chemistry, University of Alicante, Alicante, Spain
| | - J C Sancho-García
- Department of Physical Chemistry, University of Alicante, Alicante, Spain
| |
Collapse
|
6
|
Zhao D, Zhao Y, He X, Li Y, Ayers PW, Liu S. Accurate and Efficient Prediction of Post-Hartree-Fock Polarizabilities of Condensed-Phase Systems. J Chem Theory Comput 2023; 19:6461-6470. [PMID: 37676647 DOI: 10.1021/acs.jctc.3c00646] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/08/2023]
Abstract
To accurately and efficiently predict the molecular response properties (such as polarizability) at post-Hartree-Fock levels for condensed-phase systems under periodic boundary conditions (PBC) is still an unaccomplished and ongoing task. We demonstrate that static isotropic polarizabilities can be cost-effectively predicted at post-Hartree-Fock levels by combining the linear-scaling generalized energy-based fragmentation (GEBF) and information-theoretic approach (ITA) quantities. In PBC-GEBF, the total molecular polarizability of an extended system is obtained as a linear combination of the corresponding quantities of a series of small embedded subsystems of several monomers. Here, we show that in the PBC-GEBF-ITA framework, one can obtain the molecular polarizabilities and establish linear relations to ITA quantities. Once these relations are established for smaller subsystems, one can predict the polarizabilities of larger subsystems directly from the molecular wavefunction (or electron density) via ITA quantities. Alternatively, one can determine the total molecular polarizability via a linear combination equation in PBC-GEBF. We have corroborated that this newly proposed PBC-GEBF-ITA protocol is much more efficient than the original PBC-GEBF approach but is not much less accurate and that this conclusion holds for both many-body perturbation theory and the coupled cluster calculations. Good efficiency and transferability of the PBC-GEBF-ITA protocol are demonstrated for periodic systems with several hundred atoms in a unit cell.
Collapse
Affiliation(s)
- Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming 650500, P. R. China
| | - Yilin Zhao
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON L8S 4M1, Canada
| | - Xin He
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao 266237, P. R. China
| | - Yunzhi Li
- School of Chemistry and Chemical Engineering, Linyi University, Linyi 276000, P. R. China
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON L8S 4M1, Canada
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| |
Collapse
|
7
|
Besalú-Sala P, Bruneval F, Pérez-Jiménez ÁJ, Sancho-García JC, Rodríguez-Mayorga M. RPA, an Accurate and Fast Method for the Computation of Static Nonlinear Optical Properties. J Chem Theory Comput 2023; 19:6062-6069. [PMID: 37696751 PMCID: PMC10861135 DOI: 10.1021/acs.jctc.3c00674] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Indexed: 09/13/2023]
Abstract
The accurate computation of static nonlinear optical properties (SNLOPs) in large polymers requires accounting for electronic correlation effects with a reasonable computational cost. The Random Phase Approximation (RPA) used in the adiabatic connection fluctuation theorem is known to be a reliable and cost-effective method to render electronic correlation effects when combined with density-fitting techniques and integration over imaginary frequencies. We explore the ability of the RPA energy expression to predict SNLOPs by evaluating RPA electronic energies in the presence of finite electric fields to obtain (using the finite difference method) static polarizabilities and hyperpolarizabilities. We show that the RPA based on hybrid functional self-consistent field calculations yields accurate SNLOPs as the best-tuned double-hybrid functionals developed today, with the additional advantage that the RPA avoids any system-specific adjustment.
Collapse
Affiliation(s)
- Pau Besalú-Sala
- Department
of Chemistry and Pharmaceutical Sciences, Amsterdam Institute for
Molecular and Life Sciences (AIMMS), Vrije
Universiteit Amsterdam, De Boelelaan 1083, HV Amsterdam 1081, The Netherlands
- Institut
de Química Computacional i Catàlisi and Departament
de Química, Universitat de Girona, Girona 17003, Spain
| | - Fabien Bruneval
- Université
Paris-Saclay, CEA, Service de recherche en Corrosion et Comportement
des Matériaux, SRMP, Gif-sur-Yvette 91191, France
| | | | | | | |
Collapse
|