Different Levels Of Theory Research Articles

Laser Spectroscopic Characterization of Supersonic Jet-Cooled 2,6-Diazaindole (26DAI).

The article presents a comprehensive laser spectroscopic characterization of a nitrogen-rich indole derivative, namely, 2,6-diazaindole (26DAI), in the gas phase. A supersonic jet-cooled molecular beam of 26DAI was characterized using two-color resonant two-photon ionization (2C-R2PI) and laser-induced fluorescence spectroscopy (LIF) to investigate the electronic excitation. The S1 ← S0 origin transition was obtained at 33915 cm-1, which was red-shifted from that of one (indole) and two (7-azaindole) nitrogen-containing indole derivatives by 1317 and 713 cm-1, respectively. The molecular orbital and energy analysis for the S1 ← S0 transition shows the significant stabilization of LUMO on subsequent N-insertion, resulting in the lowering of the S1 ← S0 (ππ*) transition energy. The single vibronic level fluorescence spectrum from the vibrationless S1 state of the molecule was recorded. The spectrum displayed an extensive Franck-Condon activity until 2500 cm-1 for the vibrational modes of the S0 state of the 26DAI molecule. The experimental ground state vibrational frequencies were compared to the calculated ones obtained at three different levels of theories. More accurate results were found at DFT B3LYP-D4 than those at the wave function-based MP2 and CCSD levels of theories. Further, the N-H stretching frequency of 26DAI in the S0 state was measured at 3524 cm-1 using fluorescence-dip infrared (FDIR) spectroscopy. The stability of 26DAI against ionization radiation was probed by measuring the two-color photoionization energy (IE2P) of 26DAI at 71866 cm-1. The IE2P value is significantly higher than those of N-poor counterparts (indole and 7-azaindole). The NBO charges and spin density (SD) values of the 26DAI molecule have shown that electronegative N(6) makes the cationic ground state less stable due to the position of the positive centers on the N atom. The results provided insights into the stability of N-rich biomolecules against photodamage. The current investigation can shed light on nature's way of stabilizing biomolecules with a possible N-insertion mechanism.

Local potential energy density – A DFT analysis and the local binding energy in complexes with multiple interactions

A recently developed method, so-called local potential energy density, LPED, provides the binding energy density of intra/intermolecular interactions. The LPED cannot directly give the binding energy of intra/intermolecular interactions. However, it can indirectly give the binding energy through the linear equation between LPED and supramolecular binding energy, SME. In addition, the LPED can be used to obtain the SME of local or individual interactions indirectly for the case of complexes with multiple interactions, which cannot be obtained for any other method to our knowledge. The calculation of the LPED was evaluated with three different levels of theory using density functional methods. The linearity of LPED and SME between the reference level of theory (ωB97X-D/aug-cc-pVTZ) and the other levels of theory are similar among the studied levels of theory. In addition, LPED was used indirectly to obtain the local binding energy of intermolecular interactions of complexes with multiple interactions, such as the EDTA-Ca+2 and the fosfomycin-Ca+2.

Ab initio methods applied to carbon-containing defects in hexagonal boron nitride

The functionalities activated by defect centers in solids are constantly growing, opening new avenues for sustainable future technologies. These may extend to quantum optoelectronics if suitable defect centers are created and their properties understood. Recent progress in developing quantum emitters in hexagonal boron nitride (hBN) associated with carbon impurities enabled the realization of such concepts in atomically thin films, where the defect centers exhibit an unprecedented level of sensitivity toward the environment. The complexity of defects, together with new control knobs provided by van der Waals technology, poses a challenge for theory to accurately predict the properties of defect centers and to match them with experimental results. Here, we review the ab initio methods applied to carbon-containing defect centers in hBN, exploring the predictive capabilities of different levels of theory for their structural and optoelectronic properties.

In silico modelling of radiative efficiencies of anthropogenic greenhouse gases

Radiative efficiency (RE) is a climate metric adopted in international reports on climate change to quantify the greenhouse capacity of gases, and hence to guide decision-making processes and drive transitions in the production and utilization of chemicals in different application fields. Key quantities for the determination of the RE of a gas are the atmospheric irradiance profile and the infrared (IR) absorption cross section spectrum. The latter is usually measured experimentally, even though acquiring high-quality IR spectra can pose severe challenges, sometimes limiting the accuracy or the accessible spectral range. While computational quantum chemistry methods have emerged as valuable tools to simulate IR absorption properties, their application to REs estimation is still limited to the use of the double-harmonic approximation, which presents fundamental limitations. In this work, a cost-effective quantum chemical (QC) workflow including non-empirical anharmonic contributions to spectral properties and an automatic identification of conformer distribution is presented for the accurate evaluation of REs using a range of atmospheric irradiance profiles. Different levels of theory are considered, according to the current state-of-the-art, and the accuracy of the QC RE tool is demonstrated with reference to a number of representative halocarbons widely used in refrigeration, manufacturing, and pharmaceutical fields. The results show that REs can be computed with an average accuracy of 5% using double-hybrid functionals, which overshoot the widely used B3LYP method. Finally, the QC methodology is applied to determine the REs of selected halocarbons for which data is limited, or to address some contradictory results appeared in the literature for some species. The outcomes of this work demonstrate that QC anharmonic IR cross section spectra can be used to estimate REs with an accuracy on par with that of experimental measurements, hence applicable to challenging cases for providing data for policymakers as well for screening purposes when seeking new replacement compounds.

Density functional theory study of substituent effects on gas-phase heterolytic Fe-N bond energies of m-G-C6H4NHFe(CO)2(η5-C5H5) and m-G-C6H4N(COMe)Fe(CO)2(η5-C5H5)

The nature and strength of metal–ligand bonds in organotransition-metal complexes are crucial to the understanding of organometallic reactions and catalysis. Quantum chemical calculations at different levels of theory have been used to investigate heterolytic Fe-N bond energies of meta-substituted anilinyldicarbonyl(η5-cyclopentadienyl)iron [m-G-C6H4NH(η5-C5H5)Fe(CO)2, abbreviated as m-G-C6H4NHFp (1), where G=NO2, CN, COMe, CO2Me, CF3, Br, Cl, F, H, Me, MeO and NMe2] and meta-substituted α-acetylanilinyldicarbonyl(η5-cyclopentadienyl)iron [m-G-C6H4N(COMe)(η5-C5H5)Fe(CO)2, abbreviated as m-G-C6H4N(COMe)Fp (2)] complexes. The results show that BP86 and TPSSTPSS can provide the best price/performance ratio and more accurate predictions in the study of ΔHhet(Fe-N)′s. ΔΔHhet(Fe-N)′s (1 and 2) conform to the captodative principle. There are excellent linear free energy relations among ΔΔHhet(Fe-N)′s and the experimental and computational substituent effects on acidities of m-G-C6H4NH2, the differences of acidic dissociation constants (ΔpKas) of NH bonds for m-G-C6H4NH2 for series 1 or the substituent σm constants. The former correlations imply that the govering factors for these bond scissions are similar; the latter ones suggest that polar effects of meta-substituents show the dominant role to the magnitudes of ΔHhet(Fe-N)′s. And these correlations are in accordance with Hammett linear free energy relationships. The Fe-N bonds in series 1 are stronger than those in series 2. Therefore m-G-C6H4N(COMe)- are more stable than m-G-C6H4NH-. The absolute magnitudes of AEs, AEα-COMes and AEα-COMe,para-Gs are far larger than MEs, MEα-COMes and MEα-COMe,para-Gs. The influences of MEs, MEα-COMes and MEα-COMe,para-Gs are subtle. A better understanding of organometallic bond energies can suggest whether proposed new catalytic systems may work well.

Comparing thiol and selenol reactivity towards peroxynitrite by computer simulation

Peroxynitrite is a very reactive species implicated in a variety of pathophysiological cellular processes. Particularly, peroxynitrite-mediated oxidation of cellular thiol-containing compounds such as cysteine residues is a key process which has been extensively studied. Cysteine plays roles in many redox biochemistry pathways. In contrast, selenocysteine, the 21st amino acid, is only present in 25 human proteins. Investigating the molecular basis of selenocysteine's reactivity may provide insights into its unique role in these selenocysteine-containing proteins. The two-electron oxidation of thiols or selenols by peroxynitrite is a process that is carried out by the thiolate/selenate forms and peroxynitrous acid.In this work, we shed light on the molecular basis of the differential reactivity of both species towards peroxynitrite by means of state-of-the-art computer simulations. We performed electronic structure calculations of the reaction in the methanethiolate and methaneselenolate model systems with peroxynitrous acid at different levels of theory using an implicit solvent scheme. In addition, we employed a multi-scale quantum mechanics/molecular mechanics approach for obtaining free energy profiles of these chemical reactions in aqueous solution, which enabled the comparison between the simulations and the available experimental data. Our results suggest that the larger reactivity observed in the selenocysteine case at physiological pH is mainly due to the lower pKa, which affords a larger fraction of the reactive anionic species in these conditions, and in a second place to a slightly enhanced intrinsic reactivity of the selenate form due to its larger nucleophilicity.

The Simulation of Solvent Polarizabilities and Dipolarities with Polarizable Continuum Model.

The ability of polarizable continuum models (PCM) to simulate nonspecific solvent effects (dipolarity and polarizability) was evaluated by calculating the transition energies of 1,1,10,10-tetrabutyldecanonaene (ttbp9) and 2-N,N-dimethylamino-7-nitrofluorene (DMANF), basis of Catalán's polarizability (SP) and dipolarity (SdP) solvent scales, respectively. Time-dependent density-functional theory (TD-DFT) calculations were performed at different levels of theory, employing four basis sets in 10 different solvents, covering the full range of the normalized SP and SdP scales. Transition energies were calculated using linear response (LR) and corrected linear response (cLR2) schemes. Although these methods yielded variable mean absolute errors, the LR-PCM calculations reproduced medium polarizability and dipolarity trends. While calculated ttbp9 transition energies correlated with SP and Laurence's dispersion-induced (DI) scales, the DMANF transition energies correlated poorly with SdP or Laurence's ES dipolarity scales. This result agrees with the fact that DMANF solvatochromism is "contaminated" by solvent polarizability and HB acidity. The incorporation of SP or DI contributions led to much better (r2 > 0.95) correlations with the DMANF-calculated transitions. The results offer a clearer picture of the limitations of continuum models in simulating the behavior of solvatochromic dyes in solution by pointing out their poor performance when specific solvent effects, such as hydrogen-bond interactions, play a significant role in their solvatochromism.

Iodine Clusters in the Atmosphere I: Computational Benchmark and Dimer Formation of Oxyacids and Oxides.

The contribution of iodine-containing compounds to atmospheric new particle formation is still not fully understood, but iodic acid and iodous acid are thought to be significant contributors. While several quantum chemical studies have been carried out on clusters containing iodine, there is no comprehensive benchmark study quantifying the accuracy of the applied methods. Here, we present the first study in a series that investigate the role of iodine species in atmospheric cluster formation. In this work, we have studied the iodic acid, iodous acid, iodine tetroxide, and iodine pentoxide monomers and their dimers formed with common atmospheric precursors. We have tested the accuracy of commonly applied methods for calculating the geometry of the monomers, thermal corrections of monomers and dimers, the contribution of spin-orbit coupling to monomers and dimers, and finally, the accuracy of the electronic energy correction calculated at different levels of theory. We find that optimizing the structures either at the ωB97X-D3BJ/aug-cc-pVTZ-PP or the M06-2X/aug-cc-pVTZ-PP level achieves the best thermal contribution to the binding free energy. The electronic energy correction can then be calculated at the ZORA-DLPNO-CCSD(T0) level with the SARC-ZORA-TZVPP basis for iodine and ma-ZORA-def2-TZVPP for non-iodine atoms. We applied this methodology to calculate the binding free energies of iodine-containing dimer clusters, where we confirm the qualitative trends observed in previous studies. However, we identify that previous studies overestimate the stability of the clusters by several kcal/mol due to the neglect of relativistic effects. This means that their contributions to the currently studied nucleation pathways of new particle formation are likely overestimated.

A cadmium-thiocyanate coordination polymer with bridging 2-aminopyridine: Synthesis, characterization DFT studies, Hirshfeld Surface analysis and docking studies

A one-dimensional [Cd(SCN)2]n coordination polymer doubly bridged by 2-aminopyridine has been synthesized via slow evaporation method at low temperature and it was characterized by single-crystal X-ray diffraction, Fourier Transform Infrared (FT-IR) spectroscopy, electronic spectroscopy, thermogravimetric analysis (TGA), and elemental analysis. The polymeric complex crystallizes in the triclinic crystal system with space group P-1. An analysis of the crystal structure as well as Hirshfeld surfaces (HS) associated with 2D fingerprint plots, indicates the existence of C–H· · · π and π-stacking interactions alongside NH—N and NH—S H-bonding interactions contributing to the stabilization of the layered assembly of the compound. Theoretical Density Functional Theory (DFT) Calculations (using different levels of theory) were done on the asymmetric unit of the polymer and a polymer fragment that reflects the coordination environment of Cd2+ in the polymer, to gain more information about the coordination environment of Cd2+ in the polymer. This also permitted a quantification of the interaction energies associated with hydrogen bonding and non-covalent interactions, as well as other quantum chemical parameters. Good consistency was found between the calculated results and the experimental structure. All the levels of theory used predict that the polymer fragment was less stable compared to the asymmetric unit except at the B3LYP/SSD level of theory. Analysis of the HOMO and LUMO was used to explain the charge transfer within the polymer. Thermal analysis of the polymer indicates it decomposes in two steps resulting in a residue indexed to the hexagonal phase of CdS JCPDS [06–0314]. The activity of the metal complex against breast, lung, and liver cancer was measured. The interactions occurring in the docking calculation were examined in detail with PLIP (Protein-Ligand Interaction Profiler) analysis. Lastly, Swiss- ADME (absorption, distribution, metabolism and excretion) analysis was performed to examine the pharmacological properties of the metal complex.

Investigation of infra-red spectra of three chiral liquid crystalline compounds with different fluorosubstitution of benzene ring

The IR spectra of three liquid crystalline compounds, differing only by the number of fluorine atoms in the molecular core, are investigated. The DFT+D3/BLYP-def2SVP calculations for isolated molecules are applied for the band assignments of the experimental IR spectra collected in the crystal phase. Additional DFT calculations, including the B3LYP exchange-correlation functional and def2TZVPP basis set, are also performed for the compound without the fluorine atoms in the molecular core. The scaling factors between the experimental and calculated peak positions, obtained at different levels of theory, are compared. For all compounds, the absorption bands assigned to the C=O stretching vibrations are analyzed as a function of temperature in the isotropic liquid, smectic phases, and crystal phase. The formation of hydrogen bonding involving the C=O groups during crystallization is inferred. The differences in the crystallization temperature of three compounds and the presence of the hexatic smectic phase for only one compound are discussed.

Rovibrational analysis of AlCO3, OAlO2, and HOAlO2 for possible atmospheric detection.

The lack of observational data for the AlO molecule in the mesosphere/lower thermosphere may be due to ablated aluminum reacting quickly to form other species. Previously proposed reaction pathways show that aluminum could be ablated in the atmosphere from meteoritic activity, but there currently exist very limited spectroscopic data on the intermediates in these reactions, limiting the possible detection of said molecules. As such, rovibrational spectroscopic data are computed herein using quartic force field methodology at four different levels of theory for the neutral intermediates AlCO3, OAlO2, and HOAlO2. Each molecule exhibits multiple vibrational modes with large vibrational transition intensities. For instance, the C-O stretch (ν1) in AlCO3 has a harmonic intensity of 536 km mol-1, the Al-O stretch (ν2) in OAlO2 has an intensity of 678 km mol-1, and the out-of-plane torsion (ν9) in HOAlO2 has an intensity of 158 km mol-1. All three molecules have exceptionally large dipole moments of 6.27, 4.21, and 5.04D, respectively. These properties indicate that all three molecules are good candidates for potential atmospheric observation utilizing vibrational and/or rotational spectroscopic techniques.

Effective prediction of SnO2 conduction band edge potential: The key role of surface oxygen vacancies.

Several theoretical studies at different levels of theory have attempted to calculate the absolute position of the SnO2 conduction band, whose knowledge is key for its effective application in optoelectronic devices such us, for example, perovskite solar cells. However, the predicted band edges fall outside the experimentally measured range. In this work, we introduce a computational scheme designed to calculate the conduction band minimum values of SnO2, yielding results aligned with experiments. Our analysis points out the fundamental role of encompassing surface oxygen vacancies to properly describe the electronic profile of this material. We explore the impact of both bridge and in-plane oxygen vacancy defects on the structural and electronic properties of SnO2, explaining from an atomistic perspective the experimental observables. The results underscore the importance of simulating both types of defects to accurately predict SnO2 features and provide new fundamental insights that can guide future studies concerning design and optimization of SnO2-based materials and functional interfaces.

Reactivities of acrylamide warheads toward cysteine targets: a QM/ML approach to covalent inhibitor design.

Covalent inhibition offers many advantages over non-covalent inhibition, but covalent warhead reactivity must be carefully balanced to maintain potency while avoiding unwanted side effects. While warhead reactivities are commonly measured with assays, a computational model to predict warhead reactivities could be useful for several aspects of the covalent inhibitor design process. Studies have shown correlations between covalent warhead reactivities and quantum mechanic (QM) properties that describe important aspects of the covalent reaction mechanism. However, the models from these studies are often linear regression equations and can have limitations associated with their usage. Applications of machine learning (ML) models to predict covalent warhead reactivities with QM descriptors are not extensively seen in the literature. This study uses QM descriptors, calculated at different levels of theory, to train ML models to predict reactivities of covalent acrylamide warheads. The QM/ML models are compared with linear regression models built upon the same QM descriptors and with ML models trained on structure-based features like Morgan fingerprints and RDKit descriptors. Experiments show that the QM/ML models outperform the linear regression models and the structure-based ML models, and literature test sets demonstrate the power of the QM/ML models to predict reactivities of unseen acrylamide warhead scaffolds. Ultimately, these QM/ML models are effective, computationally feasible tools that can expedite the design of new covalent inhibitors.

Conformational Stability of 3-aminopropionitrile: DFT and Ab initio Calculations.

Many conformers of 3-aminopropionitrile are known. Due to the biomedical importance of 3-aminopropionitrile a full investigation of structural, vibrational, and other associated properties of all possible conformers was performed. The geometrical structures, relative stability, and vibrational frequencies of the gauche and trans 3-aminopropionitrile conformers have been studied using ab initio (CCSD/6-311+G(d,p)) and DFT (B3LYP and M06 functionals at 6-311+G(d,p) and aug-cc-pVDZ basis set) calculations. The conformational and vibrational studies of 3-aminopropionitrile molecule were presented here are in very good interpretation of the calculated data compared with very poor interpretation in previous studies. The results showed that the gauche 2 conformer is more stable by 0.19 kcal/mol than gauche 1, outlined as enthalpy change ΔH between the conformers, at CCSD/6-311+G(d,p). Additionally, the population analysis shows that the gauche conformers are more prevalent than the trans conformers in the gas phase, present at 72.8%, with gauche 2 being the dominating gauche conformer at 40.1%. These results are in good agreement with earlier experimental and theoretical conclusions. All minima conformers' thermodynamic characteristics have also been studied. The relevant bond lengths, bond angles, and dihedral angles were calculated at a different level of theory for all possible conformers. The geometrical outcomes of the conformers agree very well with the previous experimental results. Electrostatic potential surface (ESP) has been used to interpret the structure-activity relationship. The atomic charges are examined, together with the energy difference between HOMO (Highest Occupied Molecular Orbital) and LUMO (Lowest Unoccupied Molecular Orbital). Additionally, the HOMO-LUMO energy gap and other relevant molecular properties are computed. The most stable conformers' stabilization energy has been determined by the Natural Bond Orbital (NBO) analysis.

Evidence of Sharp Transitions between Octahedral and Capped Trigonal Prism States of the Solvation Shell of the Aqueous Fe3+ Ion.

The structure of the solvation shell of the aqueous Fe3+ ion has been a subject of controversy due to discrepancies between experiments and different levels of theory. We address this issue by performing simulations for a wide range of ion concentrations, using various potential energy functions, supplemented by density functional theory calculations of selected configurations. The solvation shell undergoes abrupt transitions between two states: a hexacoordinated octahedral (OH) state and a capped trigonal prism (CTP) state with 7-fold coordination. The lifetime of these states is dependent on concentration. In dilute FeCl3 solutions, the lifetimes of both are similar (≈1 ns). However, the lifetime of the OH state increases with ion concentration, while that of the CTP state decreases slightly. When a uniform negative background charge is used instead of explicit counterions, the lifetime of the OH state is greatly overestimated. These findings underscore the need for further experimental measurements and higher-level simulations.

Machine-Learned Kohn-Sham Hamiltonian Mapping for Nonadiabatic Molecular Dynamics.

In this work, we report a simple, efficient, and scalable machine-learning (ML) approach for mapping non-self-consistent Kohn-Sham Hamiltonians constructed with one kind of density functional to the nearly self-consistent Hamiltonians constructed with another kind of density functional. This approach is designed as a fast surrogate Hamiltonian calculator for use in long nonadiabatic dynamics simulations of large atomistic systems. In this approach, the input and output features are Hamiltonian matrices computed from different levels of theory. We demonstrate that the developed ML-based Hamiltonian mapping method (1) speeds up the calculations by several orders of magnitude, (2) is conceptually simpler than alternative ML approaches, (3) is applicable to different systems and sizes and can be used for mapping Hamiltonians constructed with arbitrary density functionals, (4) requires a modest training data, learns fast, and generates molecular orbitals and their energies with the accuracy nearly matching that of conventional calculations, and (5) when applied to nonadiabatic dynamics simulation of excitation energy relaxation in large systems yields the corresponding time scales within the margin of error of the conventional calculations. Using this approach, we explore the excitation energy relaxation in C60 fullerene and Si75H64 quantum dot structures and derive qualitative and quantitative insights into dynamics in these systems.

Crystal structure, Hirshfeld surface analysis and DFT study of N-(2-nitro-phen-yl)male-imide.

The title compound [systematic name: 1-(2-nitro-phen-yl)pyrrole-2,5-dione], C10H6N2O4, crystallizes in the monoclinic system (space group P21/n) with two mol-ecules in the asymmetric unit, which are linked by C-H⋯O hydrogen bonds. Hirshfeld surface analysis showed that the most significant contributions to the crystal packing are from H⋯O/O⋯H, H⋯C/C⋯H and H⋯H inter-actions, which contribute 54.7%, 15.2% and 15.6%, respectively. A DFT study was conducted using three different levels of theory [(B3LYP/6-311+G(d,p), wB97XD/Def2TZVPP and LC-wpbe/6-311(2 d,2p)] in order to determine the stability, structural and electronic properties of the title mol-ecule with a view to its potential applications and photochemical and copolymer properties.

Vibrational bands of formaldoxime isotopologue 13CD2NOH in the 280–4000 cm−1 region and rovibrational analysis of its [formula omitted] band by Fourier transform infrared (FTIR) spectroscopy

A total of 11 fundamental and 3 overtone bands of the formaldoxime isotopologue 13CD2NOH were identified using its Fourier transform infrared (FTIR) spectra which were recorded with a low resolution (0.50 cm−1) in the 500–4000 cm−1 region, and high resolution (0.00096 cm−1) in the 280–500 cm−1 region. Their relative infrared (IR) band intensities were also measured. Furthermore, a rovibrational analysis of the IR transitions of the ν12 band of 13CD2NOH was carried out using its high-resolution FTIR spectrum which was recorded at the Australian Synchrotron. A total of 1077 IR transitions of the C-type ν12 band were assigned and fitted using the Watson's A-reduced Hamiltonian in the Ir representation to derive its band center and the v12 = 1 state rovibrational constants up to all 5 quartic centrifugal distortion terms for the first time, with a root-mean-square (rms) deviation of 0.00044 cm−1. The band center of the ν12 band of 13CD2NOH were found to be 391.054446(36) cm−1. The ground state rovibrational constants up to all 5 quartic terms were determined for the first time by fitting 407 ground state combination differences (GSCDs) derived from the assigned IR transitions of the ν12 band of 13CD2NOH of this work. The rms deviation of the GSCD fit was 0.00040 cm−1. Additionally, all 3 rotational constants and 5 quartic centrifugal distortion terms of the ground state and 3 rotational constants of the v12 = 1 state of 13CD2NOH were computed from theoretical anharmonic calculations at two different levels of theory, B3LYP and MP2 with the cc-pVTZ basis set, for comparison with the experimental results. Close agreement was found for the calculated and experimental rovibrational constants of 13CD2NOH for both ground and v12 = 1 states. The vibrational anharmonic frequencies of the 12 fundamental bands of 13CD2NOH in the 280–4000 cm−1 region, and their IR band intensities were also calculated using B3LYP and MP2 with the cc-pVTZ basis set, and they were compared with the respective experimental data. Finally, the ground state rotational constants and the band center of the ν12 band of the cis conformer of 13CD2NOH were calculated and compared with those of the trans conformer of this work.

Exploring tunneling ESEEM beyond methyl groups in nitroxides at low temperatures.

Tunneling of methyl rotors coupled to an electron spin causes magnetic field independent electron spin echo envelope modulation (ESEEM) at low temperatures. For nitroxides containing alkyl substituents, we observe this effect as a contribution at the beginning of the Hahn echo decay signal occurring on a faster time scale than the matrix-induced decoherence. The tunneling ESEEM contribution includes information on the local environment of the methyl rotors, which manifests as a distribution of rotation barriers P(V3) when measuring the paramagnetic species in a glassy matrix. Here, we investigate the differences in tunneling behaviour of geminal methyl and ethyl group rotors in nitroxides while exploring different levels of theory in our previously introduced methyl quantum rotor (MQR) model. Moreover, we extend the MQR model to analyze the tunneling ESEEM originating from two different rotor types coupled to the same electron spin. We find that ethyl groups in nitroxides give rise to stronger tunneling ESEEM contributions than methyl groups because the difference between hyperfine couplings of their methyl protons better matches the tunneling frequency. The methyl rotors of both ethyl and propyl groups exhibit distributions at lower rotation barriers compared to geminal methyl groups. This is in good agreement with density functional theory (DFT) calculations of their rotation barriers and showcases that conformational flexibility impacts the hindrance of rotation. Using Monte-Carlo based fitting in combination with an identifiability analysis of the MQR model parameter space, we extract rotation barrier distributions for the individual rotor types in mixed-rotor nitroxides as well as identify which rotors dominate the observed tunneling contribution in the Hahn echo decay signal.

Performance of point charge embedding schemes for excited states in molecular organic crystals.

Modeling excited state processes in molecular crystals is relevant for several applications. A popular approach for studying excited state molecular crystals is to use cluster models embedded in point charges. In this paper, we compare the performance of several embedding models in predicting excited states and S1-S0 optical gaps for a set of crystals from the X23 molecular crystal database. The performance of atomic charges based on ground or excited states was examined for cluster models, Ewald embedding, and self-consistent approaches. We investigated the impact of various factors, such as the level of theory, basis sets, embedding models, and the level of localization of the excitation. We consider different levels of theory, including time-dependent density functional theory and Tamm-Dancoff approximation (TDA) (DFT functionals: ωB97X-D and PBE0), CC2, complete active space self-consistent field, and CASPT2. We also explore the impact of selection of the QM region, charge leakage, and level of theory for the description of different kinds of excited states. We implemented three schemes based on distance thresholds to overcome overpolarization and charge leakage in molecular crystals. Our findings are compared against experimental data, G0W0-BSE, periodic TDA, and optimally tuned screened range-separated functionals.