Hydrogen Bonded Complexes of Acetylacetone and Methanol: HF and DFT level Study

Summary. Several possible hydrogen-bonded complexes between the tautomeric forms of acetylacetone and methanol were studied by ab initio methods using 6-311G** and D95** basis functions at the HF and DFT (B3LYP) levels of theory. The calculations were carried out for isolated molecules and solvent assisted complexes by means of the isodensity polarized model (IPCM). The theoretical frequencies were compared with the experimental IR spectrum of an equimolar mixture of acetylacetone and methanol. It was proved that the most stable H-bonded complex acetylacetone–methanol is formed between O-bonded methanol- and the enol molecule.     

autoCAS: A Program for Fully Automated Multiconfigurational Calculations

We present our implementation autoCAS for fully automated multiconfigurational calculations, which we also make available free of charge on our webpages. The graphical user interface of autoCAS connects a general electronic structure program with a density-matrix renormalization group program to carry out our recently introduced automated active space selection protocol for multiconfigurational calculations (Stein and Reiher, J. Chem. Theory Comput., 2016, 12, 1760). Next to this active space selection, autoCAS carries out several steps of multiconfigurational calculations so that only a minimal input is required to start them, comparable to that of a standard Kohn–Sham density-functional theory calculation, so that black-box multiconfigurational calculations become feasible. Furthermore, we introduce a new extension to the selection algorithm that facilitates automated selections for molecules with large valence orbital spaces consisting of several hundred orbitals. © 2019 Wiley Periodicals, Inc.     

Environmental Effects in Computational Spectroscopy: Accuracy and Interpretation

Spectroscopic techniques are valuable tools for understanding the structure and dynamics of complex systems, such as biomolecules or nanomaterials. Most of the current research is devoted to the development of new experimental techniques for improving the intrinsic resolution of different spectra. However, the subtle interplay of several different effects acting at different length and time scales still makes the interpretation and analysis of such spectra a very difficult task. In this respect, computational spectroscopy is becoming a needful and versatile tool for the assignment and interpretation of experimental spectra. It is in fact possible nowadays to model with relatively high accuracy the physical–chemical properties of complex molecules in different environments, and to link spectroscopic evidence directly to the structural and dynamical properties of optically or magnetically active solvated probes. In this Review, significant steps toward the simulation of entire spectra in condensed phases are presented together with some basic aspects of computational spectroscopy, which highlight how intramolecular and intermolecular degrees of freedom influence several spectroscopic parameters.     

Raman and photoacoustic infrared spectra of fluorene derivatives: Experiment and calculations

Raman and photoacoustic (PA) infrared spectra of fluorene and four derivatives (2,3-benzofluorene, 2-methylfluorene, 2-ethylfluorene and 1,8-dimethylfluorene) were recorded and analyzed in this investigation. Mid- and far-infrared PA spectra were examined from about 2000 to 100 cm⁻¹. The Raman spectra spanned the same wavenumber range. Observed bands in both PA and Raman spectra were compared with DFT (harmonic) and variational (anharmonic) calculations, and with published spectra. The DFT calculations provided single-molecule frequencies, whereas the variational method yielded results for both monomeric and dimeric species. Many previously unknown bands, including numerous features due to combination and overtone transitions, were identified and assigned in this work.     

Evaluation of the Electrophilicites and Band Gaps of Conductive Polymers: A New Model

It is shown that the obtained hardness and electrophilicity values of the short length oligomers for a heterocyclic conductive polymer from density functional theory method are quadraticly correlated to those obtained from semiempirical method. Therefore these quantities for all oligomers (up to 15 repeating units) are predicted using the obtained quadratic relations. Both of these predicted quantities for different oligomers are fitted to a new exponential model (Y = Y_∞ e^d/n) to estimate the electronic properties for the considered conjugated polymers. The calculated band gaps from this model show better agreement with the reported experimental data than those predicted by the previous models. Because of a wide range of variations, the electrophilicity could be a better index for investigating the doping effect in polymers than the hardness.

     

The Si-B chromophore: A joint experimental and theoretical investigation

Silylboranes with aromatic substituents linked to boron and silicon exhibit an unexpected absorption band in the UV-Vis spectral region. When polar groups were introduced, a marked solvatochromic effect was observed in their fluorescence emission spectra, revealing a strong excited state dipole moment. Semi-empirical MNDO/d and AM1 calculations showed that, upon UV excitation, the polarity of the Si-B bond increased and the aromatic π-electrons migrated toward the Si-B bond, consistent with experimental observations.     

Front Cover: Twisting versus Delocalization in CAAC- and NHC-Stabilized Boron-Based Biradicals: The Roles of Sterics and Electronics (Chem. Eur. J. 16/2021)

Invited for the cover of this issue is Bernd Engels, Holger Braunschweig, Volker Engel and their coworkers at University of Würzburg. The image depicts bridged boron compounds which possess fascinating relationships between their composition and their geometrical and electronic structures, the latter ranging from closed-shell to biradical triplet or singlet ground state. Read the full text of the article at 10.1002/chem.202004619 .     

Interaction of triols with formaldehyde and acetone: Experimental and theoretical study

Experimental and theoretical aspects of the condensation of glycerol and its homologs (1,2,3- and 1,2,4-butanetriols) with formaldehyde and acetone are studied under conditions of acid catalysis. Calculation of the thermodynamic parameters of the resulting products by the composite method CBS-QB3 shows that the six-membered heterocycles, the products of the interaction of triols with formaldehyde, are thermodynamically more stable than the five-membered acetals, while the reaction of the same triols with acetone is preferable for the formation of the five-membered acetals. This is due to the fact that the regioselectivity of the studied reactions is determined by the structural features and reactivity of the carbocations formed in a condensed medium during the course of the reaction. According to the theoretical data obtained experimentally, during the condensation of glycerol and 1,2,4-butanetriol with formaldehyde in the most stable form of the six-membered cyclic carbocation, intramolecular hydrogen bonding and anomeric stabilization due to the axially oriented hydroxyl group take place. As a result, cation 1b–1 is 1.2–1.6 kJ/mol more stable than its five-membered isomers ( 1a–1 and 1b–2 ). It leads to the predominant formation of 1,3-dioxane ( 3b ). However, upon condensation of butanetriol-1,2,3 with formaldehyde, the intermediate cation 4a–1 turns out to be significantly more stable than the other isomers due to the strong intramolecular hydrogen bond in the six-membered ring with the participation of the hydroxyl group of the substituent and the hydroxyl group of the cationic center, leading to the predominant formation of the dioxolane 6a .     

Theoretical Characterisation of Phosphinyl Radicals and Their Magnetic Properties: g Matrix

The g matrices (g tensors) of various phosphinyl radicals (R₂P^.) were calculated using the DFT and multireference configuration interaction (MRCI) methods. The g matrices were distinctly dependent on the molecular structure of the radical. To thoroughly examine this dependence, the contributions from individual atoms and excited states were calculated. The former revealed the gain from the phosphorus atom to be preeminent unless P?O or P?S bonds are present in the radical molecule. The contributions owing to excited states arising from electronic transitions between doubly occupied molecular orbitals and the SOMO were clearly positive, as in the case of semiquinone and niroxide radicals. The transitions from the phosphorus lone pair were of paramount importance. Surprisingly, unlike for semiquinones and nitroxides, a significant negative contribution was observed from excitations from the SOMO to unoccupied molecular orbitals. For radicals with P?O bonds, this contribution to the g₂ component was dominant.     

A hybrid ab initio/free electron computational model for conjugated dye molecules: Simple cyanines and oxonols

Justifications developed for the application the free electron model to the π‐orbitals of conjugated molecules suggest that the optical properties of these molecules would be well described by a one‐dimensional free electron model with a potential chosen to reproduce the energy level spacing of the ground state occupied π‐orbitals. Such a hybrid ab initio/free electron modeling approach, where the free electron potential parameters are optimized on a molecule‐by‐molecule basis, is developed, and applied to a series of simple cyanine and oxonol dyes. The ensuing predictions for λ_max, oscillator strengths, and redox properties compare well to available experimental information. Two important strengths of this approach are that no explicit calculations of the excited electronic state are required, and that the ab initio determination of the occupied π‐orbital level spacing considers all the electrons (π and σ) of the entire molecule in a specified geometry, environment, etc. This second characteristic gives the ability to efficiently model modifications of the optical properties of conjugated molecules resulting from chemical and/or physical modifications occuring within and remote to the conjugated region of the molecule. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 943–953, 2000     

Excited States of Xanthene Analogues: Photofragmentation and Calculations by CC2 and Time‐Dependent Density Functional Theory

Action spectroscopy has emerged as an analytical tool to probe excited states in the gas phase. Although comparison of gas‐phase absorption properties with quantum‐chemical calculations is, in principle, straightforward, popular methods often fail to describe many molecules of interest—such as xanthene analogues. We, therefore, face their nano‐ and picosecond laser‐induced photofragmentation with excited‐state computations by using the CC2 method and time‐dependent density functional theory (TDDFT). Whereas the extracted absorption maxima agree with CC2 predictions, the TDDFT excitation energies are blueshifted. Lowering the amount of Hartree–Fock exchange in the DFT functional can reduce this shift but at the cost of changing the nature of the excited state. Additional bandwidth observed in the photofragmentation spectra is rationalized in terms of multiphoton processes. Observed fragmentation from higher‐lying excited states conforms to intense excited‐to‐excited state transitions calculated with CC2. The CC2 method is thus suitable for the comparison with photofragmentation in xanthene analogues.     

Force fields and vibrational spectroscopy: A discussion on status and future prospects

A start is given for a discussion about the use of force fields for describing and interpreting molecular geometries, relative energies and vibrational spectra. The development of modern, generic force fields, e.g. Class II force fields, has diminished the necessity for specific force fields, and they should be applied in far more cases than they are being used now.     

Electronic Structures of LixV2O5 (x=0.5 and 1): A Theoretical Study

The electronic properties of α‐Li_xV₂O₅ (x=0.5 and 1) are investigated using first principle calculations based on density functional theory with local density approximation. Different intercalation sites for Li in the V₂O₅ lattices are considered, showing different influences on the electronic structures of Li_xV₂O₅. The lowest total energy is found when Li is only intercalated along the c axis between two bridging oxygen ions of sequential V₂O₅ layers. The intercalation of Li into V₂O₅ does not change the electron transition property of V₂O₅, which is an indirect band gap semiconductor, but leads to a reduction of vanadium ions and an increase of the Fermi level of Li_xV₂O₅ arising from the electron transfer from the Li 2 s orbital to the initially empty conduction band of the V₂O₅ host.     

Local structure and thermodynamics of a core-softened potential fluid: theory and simulation.

Phase behavior and structural properties of homogeneous and inhomogeneous core-softened (CS) fluid consisting of particles interacting via the potential, which combines the hard-core repulsion and double attractive well interaction, are investigated. The vapour-liquid coexistence curves and critical points for various interaction ranges of the potential are determined by discrete molecular dynamics simulations to provide guidance for the choice of the bulk density and potential parameters for the study of homogeneous and inhomogeneous structures. Spatial correlations in the homogeneous CS system are studied by the Ornstein-Zernike integral equation in combination with the modified hypernetted chain (MHNC) approximation. The local structure of CS fluid subjected to diverse external fields maintaining the equilibrium with the bulk CS fluid are studied on the basis of a recently proposed third order+second order perturbation density functional approximation (DFA). The accuracy of DFA predictions is tested against the results of a grand canonical ensemble Monte Carlo simulation. Reasonable agreement between the results of both methods proves that the DFA theory applied in this work is a convenient theoretical tool for the investigation of the CS fluid, which is practically applicable for modeling numerous real systems.     

Molecular Modeling of Bioorganometallic Compounds: Thermodynamic Properties of Molybdocene–Glutathione Complexes and Mechanism of Peptide Hydrolysis

The computational study of bioinorganic complexes between transition metals and flexible ligands is still challenging, given that, besides requiring extensive conformational searches, the treatment of metal–ligand bonds demands the application of quantum chemical methods. Herein, the adducts formed between molybdocene, which exhibits antitumor activity and reacts with thiol groups to give stable water‐soluble complexes, and the tripeptide glutathione, which is a major source of biological thiols, are studied. Conformational searches are performed using the semiempirical PM6 method followed by geometry optimizations and single‐point calculations using density functional theory methods. In addition, molecular dynamics simulations of the molybdocene–glutathione complex involved in the regioselective hydrolysis of the Cys–Gly linkage are performed in explicit solvent. The reactive process is also studied theoretically on cluster models of both the molybdocene‐bound and the free peptide.     

Ab initio and DFT studies on van der Waals trimers: the OCS.(CO2)2 complexes

Ab initio calculations [MP2, MP4SDTQ, and QCISD(T)] using different basis sets [6-31G(d,p), cc-pVXZ (X = D, T, Q), and aug-cc-pVDZ] and density functional theory [B3LYP/6-31G(d,p)] calculations were carried out to study the OCS.(CO2)2 van der Waals trimer. The DFT has proved inappropriate to the study of this type of systems where the dispersion forces are expected to play a relevant role. Three minima isomers (two noncyclic and one cyclic) were located and characterized. The most stable isomer exhibits a noncyclic barrel-like structure whose bond lengths, angles, rotational constants, and dipole moment agree quite well with the corresponding experimental values of the only structure observed in recent microwave spectroscopic studies. The energetic proximity of the three isomers, with stabilization energies of 1442, 1371, and 1307 cm-1, respectively, at the CBS-MP2/cc-pVXZ (X = D, T, Q) level, strongly suggests that the two unobserved structures should also be detected as in the case of the (CO2)3 trimer where both noncyclic and cyclic isomers have been reported to exist. The many-body symmetry-adapted perturbation theory is employed to analyze the nature of the interactions leading to the formation of the different structures. The three-body contributions are small and stabilizing for the two most stable structures and almost negligible for the cyclic isomer.     

C-H stretching vibrational shift of benzene dimer: consistency of experiment and calculation.

Three low-energy structures of the benzene dimer are investigated by several theoretical procedures (RI-MP2, CCSD(T), RI-DFT-D, DFT/BH&H) covering London dispersion energy. The RI-DFT-D and CCSD(T) calculations are used to verify the DFT/BH&H dimer characteristics, as only at this level can anharmonic calculations be performed. It is ascertained that the T-shaped (C(2v)) structure, in which the C-H stretching frequency of the proton donor shows a significant blue shift, is not stable at any level of theory. It is either a transition structure or a minimum which is easily transformed into a parallel-displaced structure or a T-shaped (C(s)) structure, even at low temperature. Consequently, no blue shift can be detected. On the other hand, the calculated anharmonic IR spectra of the two most stable structures of benzene dimer, namely, the T-shaped (C(s)) and the parallel-displaced ones, give rise to a small red (and no blue) shift of the C-H stretching vibration. This finding is fully consistent with the experimental results.     

Magnetic Properties and Electronic Structure of Neptunyl(VI) Complexes: Wavefunctions,Orbitals, and Crystal‐Field Models

The electronic structure and magnetic properties of neptunyl(VI), NpO₂²⁺, and two neptunyl complexes, [NpO₂(NO₃)₃]^? and [NpO₂Cl₄]^2?, were studied with a combination of theoretical methods: ab initio relativistic wavefunction methods and density functional theory (DFT), as well as crystal‐field (CF) models with parameters extracted from the ab initio calculations. Natural orbitals for electron density and spin magnetization from wavefunctions including spin–orbit coupling were employed to analyze the connection between the electronic structure and magnetic properties, and to link the results from CF models to the ab initio data. Free complex ions and systems embedded in a crystal environment were studied. Of prime interest were the electron paramagnetic resonance g‐factors and their relation to the complex geometry, ligand coordination, and nature of the nonbonding 5f orbitals. The g‐factors were calculated for the ground and excited states. For [NpO₂Cl₄]^2?, a strong influence of the environment of the complex on its magnetic behavior was demonstrated. Kohn–Sham DFT with standard functionals can produce reasonable g‐factors as long as the calculation converges to a solution resembling the electronic state of interest. However, this is not always straightforward.