RESUMEN
In this work, the isotope effect in optical rotation (OR) is examined by exploring structure-property relationships for H â D substitutions in chiral molecules. While electronic effects serve as the dominant source of optical activity, there is a non-negligible contribution from nuclear vibrations, which changes with isotopic substitution. We employ a test set of 50 small organic molecules: three-membered rings with varying heteroatoms (PCl, PH, S, NCl, NH, O, and NBr) and functional groups (Me, F), and simulations were run at the B3LYP/aug-cc-pVDZ level of theory. The objectives of this work are to determine locations of isotopic substitution that result in significant changes in the vibrational correction to the OR and to evaluate which vibrational modes and electronic response are the major contributors to the isotope effect. Molecules with more polarizable heteroatoms in the ring (e.g., S and P) have the largest change in the vibrational correction compared to the unsubstituted parent molecules. In many cases, isotopic substitution made to the hydrogens on the opposite side of the ring from the functional group provides the largest change in the OR. H/D wagging modes and C vibrations (for D-C centers) are the largest contributors to the isotope effect. This is explained with a molecular orbital decomposition analysis of the OR. The relevant vibrational modes affect the orbital transitions that are already significant at the equilibrium geometry. However, this effect is only large when polarizable heteroatoms are involved because the electron density surrounding them is diffuse enough to feel the subtle effect of change in mass due to isotopic substitution on the relevant vibrational modes.
RESUMEN
This paper is dedicated to the quantum chemical package Jaguar, which is commercial software developed and distributed by Schrödinger, Inc. We discuss Jaguar's scientific features that are relevant to chemical research as well as describe those aspects of the program that are pertinent to the user interface, the organization of the computer code, and its maintenance and testing. Among the scientific topics that feature prominently in this paper are the quantum chemical methods grounded in the pseudospectral approach. A number of multistep workflows dependent on Jaguar are covered: prediction of protonation equilibria in aqueous solutions (particularly calculations of tautomeric stability and pKa), reactivity predictions based on automated transition state search, assembly of Boltzmann-averaged spectra such as vibrational and electronic circular dichroism, as well as nuclear magnetic resonance. Discussed also are quantum chemical calculations that are oriented toward materials science applications, in particular, prediction of properties of optoelectronic materials and organic semiconductors, and molecular catalyst design. The topic of treatment of conformations inevitably comes up in real world research projects and is considered as part of all the workflows mentioned above. In addition, we examine the role of machine learning methods in quantum chemical calculations performed by Jaguar, from auxiliary functions that return the approximate calculation runtime in a user interface, to prediction of actual molecular properties. The current work is second in a series of reviews of Jaguar, the first having been published more than ten years ago. Thus, this paper serves as a rare milestone on the path that is being traversed by Jaguar's development in more than thirty years of its existence.
RESUMEN
This work reports a computational investigation of the effect of ancillary ligands on the activity of an Rh catalyst for hydrogen evolution based on the [Cp*Rh] motif (Cp* = η5-pentamethylcyclopentadienyl). Specifically, we investigate why a bipyridyl (bpy) ligand leads to H2 generation but diphenylphosphino-based (dpp) ligands do not. We compare the full ligands to simplified models and systematically vary structural features to ascertain their effect on the reaction energy of each catalytic step. The calculations based on density functional theory show that the main effect on reactivity is the choice of linker atom, followed by its coordination. In particular, P stabilizes the intermediate Rh-hydride species by donating electron density to the Rh, thus inhibiting the reaction toward H2 generation. Conversely, N, a more electron-withdrawing center, favors H2 generation at the price of destabilizing the hydride intermediate, which cannot be isolated experimentally and makes determining the mechanism of this reaction more difficult. We also find that the steric effects of bulky substituents on the main ligand scaffold can lead to large effects on the reactivity, which may be challenging to fine-tune. On the other hand, structural features like the bite angle of the bidentate ligand have a much smaller impact on reactivity. Therefore, we propose that the choice of linker atom is key for the catalytic activity of this species, which can be further fine-tuned by a proper choice of electron-directing groups on the ligand scaffold.
RESUMEN
In this work, we explore the issue of origin dependence in optical rotation (OR) calculations in the length dipole gauge (LG) using standard approximate methods belonging to density functional theory (DFT) and coupled cluster (CC) theory. We use the origin-invariant LG approach, LG(OI), that we recently proposed as reference for the calculations, and we study whether a proper choice of coordinate origin and molecular orientation can be made such that diagonal elements of the LG-OR tensor match those of the LG(OI) tensor. Using a numerical search algorithm, we show that multiple spatial orientations can be found where the LG and LG(OI) results match. However, a simple analytical procedure provides a spatial orientation where the origin of the coordinate system is close to the center of mass of the molecule. At the same time, we also show that putting the origin at the center of mass is not an ideal choice for every molecule (relative errors in the OR up to 70% can be obtained in out test set). Finally, we show that the choice of coordinate origin based on the analytical procedure is transferable across different methods and it is superior to putting the origin in the center of mass or center of nuclear charge. This is important because the LG(OI) approach is trivial to implement for DFT, but not necessarily for nonvariational methods in the CC family. Therefore, one can determine an optimal coordinate origin at DFT level and use it for standard LG-CC response calculations.
RESUMEN
This paper reports the derivation and implementation of the electric dipole-magnetic dipole and electric dipole-electric quadrupole polarizability tensors at the density functional theory level with periodic boundary conditions (DFT-PBC). These tensors are combined to evaluate the Buckingham/Dunn tensor that describes the optical rotation (OR) in oriented chiral systems. We describe several aspects of the derivation of the equations and present test calculations that verify the correctness of the tensor formulation and their implementation. The results show that the full OR tensor is completely origin invariant as for molecules and that PBC calculations match molecular cluster calculations on 1D chains. A preliminary investigation on the choice of density functional, basis set, and gauge indicates a similar dependence as for molecules: the functional is the primary factor that determines the OR magnitude, followed by the basis set and to a much smaller extent the choice of gauge. However, diffuse functions may be problematic for PBC calculations even if they are necessary for the molecular case. A comparison with experimental data of OR for the tartaric acid crystal shows reasonable agreement given the level of theory employed. The development presented in this paper offers the opportunity to simulate the OR of chiral crystalline materials with general-purpose DFT-PBC methods, which, in turn, may help to understand the role of intermolecular interactions on this sensitive electronic property.
RESUMEN
In this work, the basis set dependence of optical rotation (OR) calculations is examined for various choices of gauge/level of theory. The OR is calculated for a set of 50 molecules using B3LYP and CAM-B3LYP and 17 molecules using coupled cluster with single and double excitations (CCSD). The calculations employ the correlation-consistent basis sets, aug-cc-pVζZ with ζ = D, T, Q. An inverse-power extrapolation formula is then utilized to obtain OR values at the complete basis set (CBS) limit. We investigate the basis set convergence for these methods and three choices of gauge: length gauge (with gauge-including atomic orbitals, LG(GIAOs), for DFT), the origin-invariant length gauge [LG(OI)], and the modified velocity gauge (MVG). The results show that all methods converge smoothly to the CBS limit and that the LG(OI) approach has a slightly faster convergence rate than the other choices of gauge. While the DFT methods reach gauge invariance at the CBS limit, CCSD does not. The significant difference between the MVG and LG(OI) results at the CBS limit, 26%, indicates that CCSD is not quite at convergence in the description of electron correlation for this property. On the other hand, gauge invariance at the CBS limit for DFT does not lead to the same OR values for the two density functionals, which is also due to electron correlation incompleteness. A limited comparison to gas-phase experimental OR values for the DFT methods shows that CAM-B3LYP seems more accurate than B3LYP. Overall, this study shows that the LG(OI) approach with the aug-cc-pVTZ basis set for DFT, and with the CBS(DT) extrapolation for CCSD, provides a good cost/accuracy balance.
RESUMEN
We present an origin-invariant approach to compute the full optical rotation tensor (Buckingham/Dunn tensor) in the length dipole gauge without recourse to London atomic orbitals, called LG(OI). The LG(OI) approach is simpler and less computationally demanding than the more common length gauge (LG)-London and modified velocity gauge (MVG) approaches, and it can be used with any approximate wave function or density functional method. We report an implementation at the coupled cluster with single and double excitations level (CCSD), for which we present the first simulations of the origin-invariant Buckingham/Dunn tensor in the LG. We compare LG(OI) and MVG results on a series of 22 organic molecules, showing good linear correlation between the approaches, although for small tensor elements, they provide values of opposite sign. We also attempt to decouple the effects of electron correlation and basis set incompleteness on the choice of gauge for specific rotation calculations on simple test systems. The simulations show a smooth convergence of the LG(OI) and MVG results with the basis set size toward the complete basis set limit. However, these preliminary results indicate that CCSD may not be close to a complete description of the electron correlation effects on this property even for small molecules and that basis set incompleteness may be a less important cause of discrepancy between choices of gauge than electron correlation incompleteness.
RESUMEN
The measurement of optical rotation (OR) is a foundational technique for the detection and characterization of chiral molecules, but it is poorly understood how the observed property relates to the structure of the molecule. Over the years, several schemes have been developed to decompose the OR into more chemically intuitive contributions. In this paper, we introduce two alternative formulations of our previously developed SÌ molecular orbital space decomposition. These new expressions use the modified velocity gauge-magnetic (MVG-M) and -electric (MVG-E) definitions of OR, rather than the length gauge-magnetic (LG-M) definition used in the original paper. Comparing these formulations across a small set of previously studied chiral molecules, we find that these different definitions produce consistent physical interpretations of the OR. These results demonstrate that the SÌ methodology for investigations of structure-property relationships in chiral molecules is insensitive to the choice of gauge.
RESUMEN
This work presents the first simulations of the full optical rotation (OR) tensor at coupled cluster with single and double excitations (CCSD) level in the modified velocity gauge (MVG) formalism. The CCSD-MVG OR tensor is origin independent, and each tensor element can in principle be related directly to experimental measurements on oriented systems. We compare the CCSD results with those from two density functionals, B3LYP and CAM-B3LYP, on a test set of 22 chiral molecules. The results show that the functionals consistently overestimate the CCSD results for the individual tensor components and for the trace (which is related to the isotropic OR), by 10%-20% with CAM-B3LYP and 20%-30% with B3LYP. The data show that the contribution of the electric dipole-magnetic dipole polarizability tensor to the OR tensor is on average twice as large as that of the electric dipole-electric quadrupole polarizability tensor. The difficult case of (1S,4S)-(-)-norbornenone also reveals that the evaluation of the former polarizability tensor is more sensitive than the latter. We attribute the better agreement of CAM-B3LYP with CCSD to the ability of this functional to better reproduce electron delocalization compared with B3LYP, consistent with previous reports on isotropic OR. The CCSD-MVG approach allows the computation of reference data of the full OR tensor, which may be used to test more computationally efficient approximate methods that can be employed to study realistic models of optically active materials.