1
|
Batista PR, Ducati LC, Autschbach J. Dynamic and relativistic effects on Pt-Pt indirect spin-spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations. J Chem Phys 2024; 160:114307. [PMID: 38497474 DOI: 10.1063/5.0196853] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2024] [Accepted: 02/22/2024] [Indexed: 03/19/2024] Open
Abstract
Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent-solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin-orbit Kohn-Sham (2c-ZKS) and four-component Dirac-Kohn-Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1-Pt2(NH3)4(Am)2-L2]n+ (Am = α-pyrrolidonate and pivalamidate; L = H2O, Cl-, and Br-) are studied. The computed Pt-Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N-C-O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X-) and aquahalo (L1 = H2O; L2 = X-) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.
Collapse
Affiliation(s)
- Patrick R Batista
- Department of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, Av. Prof. Lineu Prestes, 748, 05508-000 São Paulo, SP, Brazil
| | - Lucas C Ducati
- Department of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, Av. Prof. Lineu Prestes, 748, 05508-000 São Paulo, SP, Brazil
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, USA
| |
Collapse
|
2
|
Schenberg LA, Ducati LC, Autschbach J. Inquiring 199Hg NMR Parameters by Combining Ab Initio Molecular Dynamics and Relativistic NMR Calculations. Inorg Chem 2024; 63:2082-2089. [PMID: 38207278 DOI: 10.1021/acs.inorgchem.3c03878] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
Ab initio molecular dynamics (AIMD) sampling followed by relativistic density functional theory (DFT) 199Hg NMR calculations were performed for Hg organometallic complexes in water, dimethyl sulfoxide, and chloroform. The spin-orbit coupling, a relativistic effect, is a key factor for predicting δ(Hg) and 1J(Hg-C) accurately, in conjunction with a dynamic treatment of the systems. Good agreement between the theoretical and experimental results is reached by adopting implicit (based on a continuum model) and explicit (solvent molecules treated quantum mechanically) solvation models. Broader trends appearing in the experimental data available in the literature are reproduced by the calculations, and therefore, quantum chemistry is able to assist in the assignment and interpretation of 199Hg NMR data. Less pronounced trends, such as changes in the 199Hg chemical shift in different systems with the same atom types bound to Hg, are too weak to be predicted reliably by the current state-of-the-art theoretical methods based on AIMD sampling and relativistic DFT with hybrid functionals for NMR calculations.
Collapse
Affiliation(s)
- Leonardo Araujo Schenberg
- Department of Fundamental Chemistry, Institute of Chemistry University of Sao Paulo, Sao Paulo, Sao Paulo 05508-000, Brazil
| | - Lucas Colucci Ducati
- Department of Fundamental Chemistry, Institute of Chemistry University of Sao Paulo, Sao Paulo, Sao Paulo 05508-000, Brazil
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo State University of New York, Buffalo, New York 14260-3000, United States
| |
Collapse
|
3
|
Biswas R, Chen Y, Vela J, Rossini AJ. Relativistic DFT Calculations of Cadmium and Selenium Solid-State NMR Spectra of CdSe Nanocrystal Surfaces. ACS OMEGA 2023; 8:44362-44371. [PMID: 38027327 PMCID: PMC10666156 DOI: 10.1021/acsomega.3c07680] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/03/2023] [Revised: 10/11/2023] [Accepted: 10/16/2023] [Indexed: 12/01/2023]
Abstract
Solid-state NMR spectra have been used to probe the structure of CdSe nanocrystals and propose detailed models of their surface structures. Density functional theory (DFT)-optimized cluster models that represent probable molecular structures of carboxylate-coordinated surface sites have been proposed. However, to the best of our knowledge, 113Cd and 77Se chemical shifts have not been calculated for these surface models. We performed relativistic DFT calculations of cadmium and selenium magnetic shielding tensors on model compounds with previously measured solid-state NMR spectra with (i) the four-component Dirac-Kohn-Sham (DKS) Hamiltonian and (ii) the scalar and (iii) spin-orbit levels within the ZORA Hamiltonian. Molecular clusters with Cd and Se sites in varying bonding environments were used to model CdSe (100) and CdSe(111) surfaces capped with carboxylic acid ligands. Our calculations identify the observed 113Cd isotropic chemical shifts δ(iso) of -465, -318, and -146 ppm arising from CdSeO3, CdSe2O2, and CdSe3O surface groups, respectively, with very good agreement with experimental measurements. The 113Cd chemical shifts linearly decrease with the number of O-neighbors. The calculated spans (δ11 - δ33) encompass the experimental values for CdSe3O and CdSe2O2 clusters but are slightly larger than the measured value for CdSeO3 clusters. Relativistic DFT calculations predicted a one-bond 113Cd-77Se scalar coupling of 258 Hz, which is in good agreement with the experimental values of 250 Hz. With a dense coverage of carboxylic acid ligands, the CdSe (100) surface shows a distribution of Cd-Se bond lengths and J-couplings. Relativistic DFT simulations thus aid in interpretation of NMR spectra of CdSe nanocrystals and related nanomaterials.
Collapse
Affiliation(s)
- Rana Biswas
- U.S.
Department of Energy Ames National Laboratory, Ames, Iowa 50011. United States
- Department
of Physics and Astronomy; Electrical & Computer Engineering; Microelectronics
Research Center, Iowa State University, Ames, Iowa 50011. United States
| | - Yunhua Chen
- U.S.
Department of Energy Ames National Laboratory, Ames, Iowa 50011. United States
- Department
of Chemistry, Iowa State University, Ames, Iowa 50011. United States
| | - Javier Vela
- U.S.
Department of Energy Ames National Laboratory, Ames, Iowa 50011. United States
- Department
of Chemistry, Iowa State University, Ames, Iowa 50011. United States
| | - Aaron J. Rossini
- U.S.
Department of Energy Ames National Laboratory, Ames, Iowa 50011. United States
- Department
of Chemistry, Iowa State University, Ames, Iowa 50011. United States
| |
Collapse
|
4
|
Fernández-Alarcón A, Autschbach J. Relativistic Density Functional NMR Tensors Analyzed with Spin-free Localized Molecular Orbitals. Chemphyschem 2023; 24:e202200667. [PMID: 36169984 DOI: 10.1002/cphc.202200667] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Revised: 09/21/2022] [Indexed: 01/07/2023]
Abstract
The implementation of fast relativistic methods based on density functional theory, in conjunction with localized molecular orbital (LMO) based analysis, allows straightforward interpretations of NMR parameters in terms of contributions from core shells, lone pairs, and bonds, for compounds containing elements from across the periodic table. We present a conceptual review of a frequently used LMO analysis of NMR parameters calculated in the presence of spin-orbit interactions and other relativistic effects. An accompanying example focuses on the 15 N shielding in a heavy metal complex.
Collapse
Affiliation(s)
- Alberto Fernández-Alarcón
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, NY 14260-3000, USA
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, NY 14260-3000, USA
| |
Collapse
|
5
|
Batista PR, Penna TC, Ducati LC, Correra TC. p-Aminobenzoic acid protonation dynamics in an evaporating droplet by ab initio molecular dynamics. Phys Chem Chem Phys 2021; 23:19659-19672. [PMID: 34524295 DOI: 10.1039/d1cp01495a] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Protonation equilibria are known to vary from the bulk to microdroplet conditions, which could induce many chemical and physical phenomena. Protonated p-aminobenzoic acid (PABA + H+) can be considered a model for probing the protonation dynamics in an evaporating droplet, as its protonation equilibrium is highly dependent on the formation conditions from solution via atmospheric pressure ionization sources. Experiments using diverse experimental techniques have shown that protic solvents allow formation of the O-protomer (PABA protonated in the carboxylic acid group) stable in the gas phase, while aprotic solvents yield the N-protomer (protonated in the amino group) that is the most stable protomer in solution. In this work, we explore the protonation equilibrium of PABA solvated by different numbers of water molecules (n = 0 to 32) using ab initio molecular dynamics. For n = 8-32, the protonation is either at the NH2 group or in the solvent network. The solvent network interacts with the carboxylic acid group, but there is no complete proton transfer to form the O-protomer. For smaller clusters, however, solvent-mediated proton transfers to the carboxylic acid were observed, both via the Grotthuss mechanism and the vehicle or shuttle mechanism (for n = 1 and 2). Thermodynamic considerations allowed a description of the origins of the kinetic trapping effect, which explains the observation of the solution structure in the gas phase. This effect likely occurs in the final evaporation steps, which are outside the droplet size range covered by previous classical molecular dynamics simulations of charged droplets. These results may be considered relevant in determining the nature of the species observed in the ubiquitous ESI based mass spectrometry analysis, and in general for droplet chemistry, explaining how protonation equilibria are drastically changed from bulk to microdroplet conditions.
Collapse
Affiliation(s)
- Patrick R Batista
- Department of Fundamental Chemistry, Institute of Chemistry - University of São Paulo, Av. Prof. Lineu Prestes, 748, Cidade Universitária, São Paulo, SP, Brazil.
| | - Tatiana C Penna
- Department of Fundamental Chemistry, Institute of Chemistry - University of São Paulo, Av. Prof. Lineu Prestes, 748, Cidade Universitária, São Paulo, SP, Brazil.
| | - Lucas C Ducati
- Department of Fundamental Chemistry, Institute of Chemistry - University of São Paulo, Av. Prof. Lineu Prestes, 748, Cidade Universitária, São Paulo, SP, Brazil.
| | - Thiago C Correra
- Department of Fundamental Chemistry, Institute of Chemistry - University of São Paulo, Av. Prof. Lineu Prestes, 748, Cidade Universitária, São Paulo, SP, Brazil.
| |
Collapse
|
6
|
Batista PR, Ducati LC, Autschbach J. Solvent effect on the 195Pt NMR properties in pyridonate-bridged Pt III dinuclear complex derivatives investigated by ab initio molecular dynamics and localized orbital analysis. Phys Chem Chem Phys 2021; 23:12864-12880. [PMID: 34075921 DOI: 10.1039/d0cp05849a] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
An ab initio molecular dynamics investigation of the solvent effect (water) on the structural parameters, 195Pt NMR spin-spin coupling constants (SSCCs) and chemical shifts of a series of pyridonate-bridged PtIII dinuclear complexes is performed using Kohn-Sham (KS) Car-Parrinello molecular dynamics (CPMD) and relativistic hybrid KS NMR calculations. The indirect solvent effect (via structural changes) has a dramatic effect on the 1JPtPt SSCCs. The complexes exhibit a strong trans influence in solution, where the Pt-Pt bond lengthens with increasing axial ligand σ-donor strength. In the diaqua complex, where the solvent effect is more pronounced, the SSCCs averaged for CPMD configurations with explicit plus implicit solvation agree much better with the experimental data, while the calculations for static geometry and CPMD unsolvated configurations show large deviations with respect to experiment. The combination of CPMD with hybrid KS NMR calculations provides a much more realistic computational model that reproduces the large magnitudes of 1JPtPt and 195Pt chemical shifts. An analysis of 1JPtPt in terms of localized and canonical orbitals shows that the SSCCs are driven by changes in the s-character of the natural atomic orbitals of Pt atoms, which affect the 'Fermi contact' mechanism.
Collapse
Affiliation(s)
- Patrick R Batista
- Department of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, Av. Prof. Lineu Prestes, 748, 05508-000, São Paulo, SP, Brazil.
| | | | | |
Collapse
|
7
|
Cuesta Calvo PV, Rodrigues Batista P, Rodrigues de Oliveira Silva R, Converti A, Al Arni S, Solisio C, Ducati LC, Alves Palma MS. Flow Synthesis of 2‐[Methyl(pyridin‐2‐yl)amino]ethanol: An Experimental and Computational Study. Chem Eng Technol 2020. [DOI: 10.1002/ceat.202000423] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Paulo Victor Cuesta Calvo
- University of São Paulo School of Pharmaceutical Sciences Av. Prof. Lineu Prestes, 580 – Butantã 05508-000 São Paulo SP Brazil
| | - Patrick Rodrigues Batista
- University of São Paulo Institute of Chemistry Av. Prof. Lineu Prestes, 748 – Butantã 05508-000 São Paulo SP Brazil
| | | | - Attilio Converti
- University of Genoa, Pole of Chemical Engineering Department of Civil, Chemical and Environmental Engineering via Opera Pia 15 16145 Genoa Italy
| | - Saleh Al Arni
- University of Ha'il Department of Chemical Engineering College of Engineering P.O. Box 2440 81441 Ha'il Saudi Arabia
| | - Carlo Solisio
- University of Genoa, Pole of Chemical Engineering Department of Civil, Chemical and Environmental Engineering via Opera Pia 15 16145 Genoa Italy
| | - Lucas C. Ducati
- University of São Paulo Institute of Chemistry Av. Prof. Lineu Prestes, 748 – Butantã 05508-000 São Paulo SP Brazil
| | - Mauri Sergio Alves Palma
- University of São Paulo School of Pharmaceutical Sciences Av. Prof. Lineu Prestes, 580 – Butantã 05508-000 São Paulo SP Brazil
| |
Collapse
|
8
|
Philips A, Autschbach J. Proton NMR relaxation from molecular dynamics: intramolecular and intermolecular contributions in water and acetonitrile. Phys Chem Chem Phys 2019; 21:26621-26629. [DOI: 10.1039/c9cp04976b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Accurate 1H NMR relaxation rates for protons in pure water and acetonitrile are computed via ab initio and force field molecular dynamics. Dipole–dipole and spin-rotation mechanisms are considered.
Collapse
Affiliation(s)
- Adam Philips
- Department of Chemistry
- University at Buffalo
- State University of New York
- Buffalo
- USA
| | - Jochen Autschbach
- Department of Chemistry
- University at Buffalo
- State University of New York
- Buffalo
- USA
| |
Collapse
|
9
|
Saielli G. Computational Spectroscopy of Ionic Liquids for Bulk Structure Elucidation. ADVANCED THEORY AND SIMULATIONS 2018. [DOI: 10.1002/adts.201800084] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Giacomo Saielli
- CNR Institute on Membrane Technology; Unit of Padova; Via Marzolo 1-35131 Padova Italy
- Department of Chemical Sciences; University of Padova; Via Marzolo 1-35131 Padova Italy
| |
Collapse
|
10
|
Lino JBR, Rocha EP, Ramalho TC. Value of NMR Parameters and DFT Calculations for Quantum Information Processing Utilizing Phosphorus Heterocycles. J Phys Chem A 2017; 121:4486-4495. [PMID: 28441482 DOI: 10.1021/acs.jpca.6b12728] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Quantum computing is the field of science that uses quantum-mechanical phenomena, such as superposition and entanglement, to perform operations on data. The fundamental information unit used in quantum computing is the quantum bit or qubit. It is well-known that quantum computers could theoretically be able to solve problems much more quickly than any classical computers. Currently, the first and still the most successful implementations of quantum information processing (QIP) have been based on nuclear spins in liquids. However, molecules that enable many qubits NMR QIP implementations should meet some conditions: have large chemical shifts and be appropriately dispersed for qubit addressability, appreciable spin-spin coupling between any pair of spins, and a long relaxation time. In this line, benzyldene-2,3-dihydro-1H-[1,3]diphosphole (BDF) derivatives have been theoretically tested for maximizing large chemical shifts, spin-spin coupling, and minimizing the hyperfine coupling constant. Thus, the structures were optimized at the B3LYP/6-311G(d,p) level and showed a significant similarity with the experimental geometrical parameters. The NMR spectroscopic parameters (δ and J) were calculated with six different DFT functionals. The τ-HCTH/6-31G(2d) level is in better agreement with the experimental data of 31P and 13C chemical shifts, while PCM-B3LYP/cc-pVDZ level shows a decrease on deviation between calculated and experimental values for P-P and P-C SSCC. The surface response technique was employed to rationalize how the hyperfine constant varies with the chemical shifts and coupling constants values. From our findings, BDF-NO2 was the best candidate for NMR quantum computations (NMR-QC) among the studied series.
Collapse
Affiliation(s)
- Jéssica B R Lino
- Department of Chemistry, Federal University of Lavras , P.O. Box 3037, 37200-000, Lavras, MG, Brazil
| | - Eduardo P Rocha
- Federal Institute of Science, Education and Technology of Southwest MG , Campus Rio, Pomba, MG, Brazil
| | - Teodorico C Ramalho
- Department of Chemistry, Federal University of Lavras , P.O. Box 3037, 37200-000, Lavras, MG, Brazil.,Center for Basic and Applied Research, Faculty of Informatics and Management, University of Hradec Kralove , Hradec Kralove, Czech Republic
| |
Collapse
|
11
|
Marchenko A, Truflandier LA, Autschbach J. Uranyl Carbonate Complexes in Aqueous Solution and Their Ligand NMR Chemical Shifts and 17O Quadrupolar Relaxation Studied by ab Initio Molecular Dynamics. Inorg Chem 2017; 56:7384-7396. [PMID: 28598146 DOI: 10.1021/acs.inorgchem.7b00396] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Dynamic structural effects, NMR ligand chemical shifts, and 17O NMR quadrupolar relaxation rates are investigated in the series of complexes UO22+, UO2(CO3)34-, and (UO2)3(CO3)66-. Car-Parrinello molecular dynamics (CPMD) is used to simulate the dynamics of the complexes in water. NMR properties are computed on clusters extracted from the CPMD trajectories. In the UO22+ complex, coordination at the uranium center by water molecules causes a decrease of around 300 ppm for the uranyl 17O chemical shift. The final value of this chemical shift is within 40 ppm of the experimental range. The UO2(CO3)34- and (UO2)3(CO3)66- complexes show a solvent dependence of the terminal carbonate 17O and 13C chemical shifts that is less pronounced than that for the uranyl oxygen atom. Corrections to the chemical shift from hybrid functionals and spin-orbit coupling improve the accuracy of chemical shifts if the sensitivity of the uranyl chemical shift to the uranyl bond length (estimated at 140 ppm per 0.1 Å from trajectory data) is taken into consideration. The experimentally reported trend in the two unique 13C chemical shifts is correctly reproduced for (UO2)3(CO3)66-. NMR relaxation rate data support large 17O peak widths, but remain below those noted in the experimental literature. Comparison of relaxation data for solvent-including versus solvent-free models suggest that carbonate ligand motion overshadows explicit solvent effects.
Collapse
Affiliation(s)
- Alex Marchenko
- Department of Chemistry, University at Buffalo, State University of New York , Buffalo, New York 14260-3000, United States
| | - Lionel A Truflandier
- Institut des Sciences Moleculaires, Universite Bordeaux , CNRS UMR 5255, 33405 Talence cedex, France
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York , Buffalo, New York 14260-3000, United States
| |
Collapse
|