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.
This content will become publicly available on August 5, 2025
First Principles Simulations of Optical Rotation of Chiral Molecular Crystals
In this work, we present simulations of the optical rotation (OR) for five molecular crystals at density functional theory level with periodic boundary conditions (DFT‐PBC). Calculations are compared with experimental measurements and show semi‐quantitative agreement with experimental data for three of the crystals: tartatic acid, benzil, and pentaerythritol. For the other two crystals, aspartic acid and glutamic acid, the calculated data are in qualitative agreement with, but two orders of magnitude smaller than, the experimental data. We provide some arguments that support the theoretical predictions and suggest that the experiments should be revisited. We also find that the position of H centers provided in experimental X‐ray data is not sufficiently reliable for simulating OR, and better results are obtained when H atoms are allowed to relax while keeping heavier elements fixed at the experimental positions. Comparison with molecular cluster calculations with a better functional and a larger basis set indicate that the role of intermolecular interactions (reproduced with the PBC technique) is as or more important than the choice of model chemistry. Despite the current limitations in the level of theory that can be employed, these simulations provide a promising avenue to investigate the effect of intermolecular interactions on this sensitive electronic property of molecules and materials.
more » « less- PAR ID:
- 10531149
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Chirality
- Volume:
- 36
- Issue:
- 8
- ISSN:
- 0899-0042
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
null (Ed.)9-Anthracene carboxylic acid derivatives comprise a family of thermally reversible photomechanical molecular crystals. The photomechanical response relies on a [4 + 4] photodimerization followed by dissociation that occurs on timescales of seconds to minutes. A combined theoretical and experimental investigation is undertaken to better understand how chemical modification of the anthracene core influences energetics of both the isolated molecule and the crystal lattice. We use both density functional theory and dispersion-corrected Moller–Plesset perturbation theory computational methods to establish orbital energies, photodimerization reaction energies, and lattice energies for a set of substituted 9-anthracene carboxylic acid molecules. The calculations reveal that steric interactions play a dominant role in the ability to form photodimers and indicate an energetic threshold of 80–90 kJ per mole for the dimerization reaction. Examination of intermolecular bonding in a subset of fluorinated 9ACs revealed the absence of H⋯F intermolecular bond formation and energy differences that can explain observed trends in the dissociation kinetics and mechanical reset times. Fluorescence recovery after photobleaching experiments shows that the photodimer dissociation kinetics depend on the amount of initial photodimer, preventing a straightforward correlation between halogen atom substitution and dissociation rates using the Bell–Evans–Polanyi principle. The results clarify how molecular structure affects intermolecular interactions and photoreactivity in this family of molecular crystals, but the origin of the complex photodimer dissociation dynamics remains an open question.more » « less
-
Organic chromophores initiate much of daytime aqueous phase chemistry in the environment. Thus, studying the absorption spectra of commonly used organic photosensitizers is paramount to fully understand their relevance in environmental processes. In this work, we combined UV-Vis spectroscopy, 1 H-NMR spectroscopy, quantum chemical calculations, and molecular dynamics simulations to investigate the absorption spectra of 4-benzoyl benzoic acid (4BBA), a widely used photosensitizer and a common proxy of environmentally relevant chromophores. Solutions of 4BBA at different pH values show that protonated and deprotonated species have an effect on its absorbance spectra. Theoretical calculations of these species in water clusters provide physical and chemical insights into the spectra. Quantum chemical calculations were conducted to analyze the UV-Vis absorbance spectra of 4BBA species using various cluster sizes, such as C 6 H 5 COC 6 H 4 COOH·(H 2 O) n , where n = 8 for relatively small clusters and n = 30 for larger clusters. While relatively small clusters have been successfully used for smaller chromophores, our results indicate that simulations of protonated species of 4BBA require relatively larger clusters of n = 30. A comparison between the experimental and theoretical results shows good agreement in the pH-dependent spectral shift between the hydrated cluster model and the experimental data. Overall, the theoretical and empirical results indicate that the experimental optical spectra of aqueous phase 4BBA can be represented by the acid–base equilibrium of the keto-forms, with a spectroscopically measured p K a of 3.41 ± 0.04. The results summarized here contribute to a molecular-level understanding of solvated organic molecules through calculations restricted to cluster models, and thereby, broader insight into environmentally relevant chromophores.more » « less
-
The absorption spectra of molecular organic chromophores in aqueous media are of considerable importance in environmental chemistry. In this work, the UV-vis spectra of benzoic acid (BA), the simplest aromatic carboxylic acid, in aqueous solutions at varying pH and in the presence of salts are measured experimentally. The solutions of different pH provide insights into the contributions from both the non-dissociated acid molecule and the deprotonated anionic species. The microscopic interpretation of these spectra is then provided by quantum chemical calculations for small cluster models of benzoic species (benzoic acid and benzoate anion) with water molecules. Calculations of the UV-vis absorbance spectra are then carried out for different clusters such as C 6 H 5 COOH·(H 2 O) n and C 6 H 5 COO − ·(H 2 O) n , where n = 0–8. The following main conclusions from these calculations and the comparison to experimental results can be made: (i) the small water cluster yields good quantitative agreement with observed solution experiments; (ii) the main peak position is found to be very similar at different levels of theory and is in excellent agreement with the experimental value, however, a weaker feature about 1 eV to lower energy (red shift) of the main peak is correctly reproduced only by using high level of theory, such as Algebraic Diagrammatic Construction (ADC); (iii) dissociation of the BA into ions is found to occur with a minimum of water molecules of n = 8; (iv) the deprotonation of BA has an influence on the computed spectrum and the energetics of the lowest energy electronic transitions; (v) the effect of the water on the spectra is much larger for the deprotonated species than for the non-dissociated acid. It was found that to reproduce experimental spectrum at pH 8.0, additional continuum representation for the extended solvent environment must be included in combination with explicit solvent molecules ( n ≥ 3); (vi) salts (NaCl and CaCl 2 ) have minimal effect on the absorption spectrum and; (vii) experimental results showed that B-band of neutral BA is not sensitive to the solvent effects whereas the effect of the water on the C-band is significant. The water effects blue-shift this band up to ∼0.2 eV. Overall, the results demonstrate the ability to further our understanding of the microscopic interpretation of the electronic structure and absorption spectra of BA in aqueous media through calculations restricted to small cluster models.more » « less
-
The 3d transition metal insulator Bi2CrAl3O9 forms with a quasi-one-dimensional structure characterized by linear chains of edge-sharing, Cr-and Al-centered, distorted octahedra. The UV/Vis spectrum of high-quality single crystals is marked by broad absorption edges corresponding to direct transitions across a 1.36-eV insulating gap. Measurements of dc magnetic susceptibility χ reveal a fluctuating moment of 2.60±0.01μB/Cr—reduced from the 3.87μB/Cr expected for Cr3+, while the Weiss temperature ΘW=−21±1 K implies that the prevailing local moment interactions are weakly antiferromagnetic in nature. Some 10% of the fluctuating moment is quenched, presumably due to the onset of an antiferromagnetic or spin glass phase at temperature T★=98±3 K, while measurements of magnetization versus field H at T≤10 K scale as H/T0.68(4), suggesting the presence of quantum fluctuations associated with a disordered phase. Density functional theory calculations carried out within the generalized gradient approximation are in excellent agreement with experimental results, asserting that short-range magnetic interactions remnant above T★ stabilize the insulating statemore » « less
