1
|
Gustafson A, Sburlati S, Kahr B. Computed Gyration Tensors of Knotted Chiral and Achiral Topological Stereoisomers of C 60 Cyclocarbons. Chemphyschem 2024; 25:e202400277. [PMID: 38606486 DOI: 10.1002/cphc.202400277] [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: 03/12/2024] [Revised: 04/10/2024] [Accepted: 04/10/2024] [Indexed: 04/13/2024]
Abstract
The electronic origins of the computed optical rotations of the simplest chiral and achiral chemical knots with comparatively simple compositions and large, anticipated magnetoelectric polarizabilities are provided. Linear response theory (LRT) is used to calculate the gyration at 1064 nm of two knotted polyyne chains, topological stereoisomers of cyclo[60]carbon. One isomer is analogous to the trefoil knot with approximate D3 symmetry and the other to the figure eight knot with approximate S4 symmetry. The response in each case can be attributed largely to the magnetic dipole term that arises in a near degenerate E-like excited state. An oriented achiral figure eight knot is as optically active in some directions as the chiral knot in any direction, and its absolute eigenvalues are larger.
Collapse
Affiliation(s)
- Afton Gustafson
- Department of Chemistry and Molecular Design Institute, New York University, 100 Washington Square East, Room 1001, New York City, New York, 10003, USA
| | - Sophia Sburlati
- Department of Chemistry and Molecular Design Institute, New York University, 100 Washington Square East, Room 1001, New York City, New York, 10003, USA
| | - Bart Kahr
- Department of Chemistry and Molecular Design Institute, New York University, 100 Washington Square East, Room 1001, New York City, New York, 10003, USA
| |
Collapse
|
2
|
Zara A, Albu RM, Stoica I, Barzic AI, Dimitriu DG, Dorohoi DO. New Insights on Solvent-Induced Changes in Refractivity and Specific Rotation of Poly(propylene oxide) Systems Extracted from Channeled Spectra. Int J Mol Sci 2024; 25:4682. [PMID: 38731902 PMCID: PMC11083607 DOI: 10.3390/ijms25094682] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2024] [Revised: 04/23/2024] [Accepted: 04/24/2024] [Indexed: 05/13/2024] Open
Abstract
Investigation of chiroptical polymers in the solution phase is paramount for designing supramolecular architectures for photonic or biomedical devices. This work is devoted to the case study of poly(propylene oxide) (PPO) optical activity in several solvents: benzonitrile, carbon disulfide, chloroform, ethyl acetate, and p-dioxane. To attain information on the interactions in these systems, rheological testing was undertaken, showing distinct variations of the rheological parameters as a function of the solvent type. These aspects are also reflected in the refractive index dispersive behavior, from which linear and non-linear optical properties are extracted. To determine the circular birefringence and specific rotation of the PPO solutions, the alternative method of the channeled spectra was employed. The spectral data were correlated with the molecular modeling of the PPO structural unit in the selected solvents. Density functional theory (DFT) computational data indicated that the torsional potential energy-related to the O1-C2-C3-O4 dihedral angle from the polymer repeating unit-was hindered in solvation environments characterized by high polarity and the ability to interact via hydrogen bonding. This was in agreement with the optical characterization of the samples, which indicated a lower circular birefringence and specific rotation for the solutions of PPO in ethyl acetate and p-dioxane. Also, the shape of optical rotatory dispersion curves was slightly modified for PPO in these solvents compared with the other ones.
Collapse
Affiliation(s)
- Alexandru Zara
- Faculty of Physics, “Alexandru Ioan Cuza” University, 700506 Iasi, Romania; (A.Z.); (D.O.D.)
| | - Raluca Marinica Albu
- “Petru Poni” Institute of Macromolecular Chemistry, 700487 Iasi, Romania; (R.M.A.); (I.S.)
| | - Iuliana Stoica
- “Petru Poni” Institute of Macromolecular Chemistry, 700487 Iasi, Romania; (R.M.A.); (I.S.)
| | - Andreea Irina Barzic
- “Petru Poni” Institute of Macromolecular Chemistry, 700487 Iasi, Romania; (R.M.A.); (I.S.)
| | - Dan Gheorghe Dimitriu
- Faculty of Physics, “Alexandru Ioan Cuza” University, 700506 Iasi, Romania; (A.Z.); (D.O.D.)
| | - Dana Ortansa Dorohoi
- Faculty of Physics, “Alexandru Ioan Cuza” University, 700506 Iasi, Romania; (A.Z.); (D.O.D.)
| |
Collapse
|
3
|
Parsons T, Balduf T, Caricato M. On the choice of coordinate origin in length gauge optical rotation calculations. Chirality 2023; 35:708-717. [PMID: 37137811 DOI: 10.1002/chir.23575] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2023] [Accepted: 04/11/2023] [Indexed: 05/05/2023]
Abstract
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.
Collapse
Affiliation(s)
- Taylor Parsons
- Department of Chemistry, University of Kansas, Lawrence, Kansas, USA
| | - Ty Balduf
- Department of Chemistry, University of Kansas, Lawrence, Kansas, USA
| | - Marco Caricato
- Department of Chemistry, University of Kansas, Lawrence, Kansas, USA
| |
Collapse
|
4
|
Morgante P, Ludowieg HD, Autschbach J. Comparative Study of Vibrational Raman Optical Activity with Different Time-Dependent Density Functional Approximations: The VROA36 Database. J Phys Chem A 2022; 126:2909-2927. [PMID: 35512708 DOI: 10.1021/acs.jpca.2c00951] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
A new database, VROA36, is introduced to investigate the performance of computational approaches for vibrational Raman optical activity (VROA) calculations. The database is composed of 36 molecules with known experimental VROA spectra. It includes 93 conformers. Normal modes calculated with B3LYP-D3(BJ)/def2-TZVP are used to compute the VROA spectra with four functionals, B3LYP-D3(BJ), ωB97X-D, M11, and optimally tuned LC-PBE, as well as several basis sets. SimROA indices and frequency scaling factors are used to compare calculated spectra with each other and with experimental data. The four functionals perform equally well independently of the basis set and usually achieve good agreement with the experimental data. For molecules in near- or at-resonance conditions, the inclusion of a complex (damped) linear response approach is important to obtain physically meaningful VROA intensities. The use of any of the tested functional approximations with the def2-SVPD Gaussian-type basis set, or a basis of similar flexibility, can be recommended for efficient and reliable theoretical VROA studies.
Collapse
Affiliation(s)
- Pierpaolo Morgante
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Herbert D Ludowieg
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| |
Collapse
|
5
|
Parsons T, Balduf T, Cheeseman JR, Caricato M. Basis Set Dependence of Optical Rotation Calculations with Different Choices of Gauge. J Phys Chem A 2022; 126:1861-1870. [PMID: 35271772 DOI: 10.1021/acs.jpca.2c00201] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Taylor Parsons
- Department of Chemistry, University of Kansas, 1567 Irving Hill Road, Lawrence, Kansas 66045, United States
| | - Ty Balduf
- Department of Chemistry, University of Kansas, 1567 Irving Hill Road, Lawrence, Kansas 66045, United States
| | - James R Cheeseman
- Gaussian, Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Marco Caricato
- Department of Chemistry, University of Kansas, 1567 Irving Hill Road, Lawrence, Kansas 66045, United States
| |
Collapse
|
6
|
Pearce KC, Fuentes RG, Calderon S, Marolikar R, Kingston DGI, Crawford TD. Structure Elucidation and Confirmation of Phloroglucinols from the Roots of Garcinia dauphinensis by Comparison of Experimental and Calculated ECD Spectra and Specific Rotations. JOURNAL OF NATURAL PRODUCTS 2021; 84:1163-1174. [PMID: 33823109 DOI: 10.1021/acs.jnatprod.0c01208] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Eight phloroglucinols from Garcinia dauphinensis were recently reported to have good to moderate antiplasmodial and anticancer activities, consistent with other phloroglucinol derivatives isolated from natural sources. Chiroptical properties were previously calculated and compared to experimental data for compound 2 as a means to deduce its absolute configuration. Tentative assignments for the remaining compounds were also reported based on these data. In order to arrive at stereochemical assignments for phloroglucinols 1 and 3-8, ECD spectra and specific rotations were computed for all stereoisomers of each compound. Molecular orbital analyses were also carried out for the most energetically favorable conformers of each compound. Absolute configurations are reported for all eight phloroglucinols for the first time.
Collapse
Affiliation(s)
- Kirk C Pearce
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Rolly G Fuentes
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Virginia Tech Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
- Division of Natural Sciences and Mathematics, University of the Philippines Visayas Tacloban College, 6500 Tacloban City, Philippines
| | - Susana Calderon
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Rageshwari Marolikar
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - David G I Kingston
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Virginia Tech Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - T Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Molecular Sciences Software Institute, 1880 Pratt Drive, Suite 1100, Blacksburg, Virginia 24060, United States
| |
Collapse
|
7
|
D'Cunha R, Crawford TD. Modeling Complex Solvent Effects on the Optical Rotation of Chiral Molecules: A Combined Molecular Dynamics and Density Functional Theory Study. J Phys Chem A 2021; 125:3095-3108. [PMID: 33829790 DOI: 10.1021/acs.jpca.1c00803] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The challenge of assigning the absolute stereochemical configuration to a chiral compound can be overcome via accurate ab initio predictions of optical rotation, a sensitive molecular property that is further complicated by solvent effects. The solvent's "chiral imprint"-the transfer of the chirality from the solute to the surrounding achiral solvent-is explored here using conformational averaging and time-dependent density-functional theory. These complex solvent effects are taken into account via simple averaging over a molecular dynamics trajectory together with the explicit quantum mechanical consideration of the solvent molecules within the solute's cybotactic region and implicit modeling of the bulk solvent. We consider several axes along which the system's optical rotation varies, including the sampling of the dynamical trajectory, the quality of the one-electron basis set, and the use of continuum solvent models to account for bulk effects.
Collapse
Affiliation(s)
- Ruhee D'Cunha
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - T Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States.,Molecular Sciences Software Institute, 1880 Pratt Drive, Suite 1100, Blacksburg, Virginia 24060, United States
| |
Collapse
|
8
|
Tai TY, Liu YH, Lai CC, Peng SM, Chiu SH. Absolute Configurations of Topologically Chiral [2]Catenanes and the Acid/Base-Flippable Directions of Their Optical Rotations. Org Lett 2019; 21:5708-5712. [PMID: 31283244 DOI: 10.1021/acs.orglett.9b02062] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
The absolute configurations of the two enantiomers of a topologically chiral [2]catenane were determined unambiguously based on HPLC resolution and X-ray crystal analysis. Although structurally dissimilar to simple amino acids, the optical rotations of these separated [2]catenanes share the Clough-Lutz-Jirgensons behavior of amino acids: the optical rotation flips direction in the presence of acid and base, the first example of such behavior for a mechanically interlocked topologically chiral catenane.
Collapse
Affiliation(s)
- Tzu-Yi Tai
- Department of Chemistry , National Taiwan University , No. 1, Sec. 4, Roosevelt Road , Taipei 10617 , Taiwan , R.O.C
| | - Yi-Hung Liu
- Department of Chemistry , National Taiwan University , No. 1, Sec. 4, Roosevelt Road , Taipei 10617 , Taiwan , R.O.C
| | - Chien-Chen Lai
- Institute of Molecular Biology , National Chung Hsing University , Taichung City 402 , Taiwan , R.O.C
| | - Shie-Ming Peng
- Department of Chemistry , National Taiwan University , No. 1, Sec. 4, Roosevelt Road , Taipei 10617 , Taiwan , R.O.C
| | - Sheng-Hsien Chiu
- Department of Chemistry , National Taiwan University , No. 1, Sec. 4, Roosevelt Road , Taipei 10617 , Taiwan , R.O.C
| |
Collapse
|
9
|
Caricato M. Linear response coupled cluster theory with the polarizable continuum model within the singles approximation for the solvent response. J Chem Phys 2018; 148:134113. [DOI: 10.1063/1.5021781] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Affiliation(s)
- Marco Caricato
- Department of Chemistry, University of Kansas, 1251 Wescoe Hall Dr., Lawrence, Kansas 66045, USA
| |
Collapse
|
10
|
Aharon T, Lemler P, Vaccaro PH, Caricato M. Comparison of measured and predicted specific optical rotation in gas and solution phases: A test for the polarizable continuum model of solvation. Chirality 2018; 30:383-395. [PMID: 29419897 DOI: 10.1002/chir.22822] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2017] [Revised: 12/17/2017] [Accepted: 12/18/2017] [Indexed: 11/12/2022]
Abstract
A comparative theoretical and experimental study of dispersive optical activity is presented for a set of small, rigid organic molecules in gas and solution phases. Target species were chosen to facilitate wavelength-resolved measurements of specific rotation in rarefied vapors and in organic solvents having different polarities, while avoiding complications due to conformational flexibility. Calculations were performed with two density functionals (B3LYP and CAM-B3LYP) and with the coupled-cluster singles and doubles (CCSD) ansatz, and solvent effects were included through use of the polarizable continuum model (PCM). Across the various theoretical methods surveyed, CCSD with the modified velocity gauge provided the best overall performance for both isolated and solvated conditions. Zero-point vibrational corrections to equilibrium calculations of chiroptical response tended to improve agreement with gas-phase experiments, but the quality of performance realized for solutions varied markedly. Direct comparison of measured and predicted specific-rotation suggests that PCM, in general, is not able to reproduce attendant solvent shifts (neither between gas and solution phases nor among solvents) and fares better in estimating actual medium-dependent values of this property (although the error is rather system dependent). Thus, more elaborate solvation models seem necessary for a proper theoretical description of solvation in dispersive optical activity.
Collapse
Affiliation(s)
- Tal Aharon
- Department of Chemistry, University of Kansas, Lawrence, Kansas
| | - Paul Lemler
- Department of Chemistry, Yale University, New Haven, Connecticut
| | | | - Marco Caricato
- Department of Chemistry, University of Kansas, Lawrence, Kansas
| |
Collapse
|