Density functional theory applied to calculating optical and spectroscopic properties of metal complexes: NMR and optical activity

Synthesis, characterization and NLO properties of 1,4-phenylenediamine-based Schiff bases: a combined theoretical and experimental approach

In the current study, three novel 1,4-phenylenediamine-based chromophores (3a-3c) were synthesized and characterized and then their nonlinear optical (NLO) characteristics were explored theoretically. The characterization was done by spectroscopic analysis, i.e. FT-IR, UV-Visible, and NMR spectroscopy, and elemental analysis. Notably, these chromophores exhibited UV-Visible absorption within the range of 378.635-384.757 nm in acetonitrile solvent. Additionally, the FMO findings for 3a-3c revealed the narrowest band gap (4.129 eV) for 3c. The GRPs for these chromophores were derived from HOMO-LUMO energy values, which showed correspondence with FMO results by depicting a minimum hardness (2.065 eV) for 3c. Among these compounds, 3c displayed the highest nonlinear behavior with maximum μtot, βtot and γtot values of 4.79 D, 8.00 × 10-30 and 8.13 × 10-34 a.u., respectively. Our findings disclosed that the synthesized 1,4-phenylenediamine chromophores may be considered promising candidates for nonlinear optical materials, showing potential applications in the realm of optoelectronic devices.

Systematic Evaluation of Modern Density Functional Methods for the Computation of NMR Shifts of 3d Transition-Metal Nuclei

A wide range of density functionals from all rungs of Jacob's ladder have been evaluated systematically for a set of experimental 3d transition-metal NMR shifts of 70 complexes encompassing 12 × ⁴⁹Ti, 10 × ⁵¹V, 10 × ⁵³Cr, 11 × ⁵⁵Mn, 9 × ⁵⁷Fe, 9 × ⁵⁹Co, and 9 × ⁶¹Ni shift values, as well as a diverse range of electronic structure characteristics. The overall 39 functionals evaluated include one LDA, eight GGAs, seven meta-GGAs (including their current-density-functional─CDFT─versions), nine global hybrids, four range-separated hybrids, eight local hybrids, and two double hybrids, and we also include Hartree-Fock and MP2 calculations. While recent evaluations of the same functionals for a very large coupled-cluster-based benchmark of main-group shieldings and shifts achieved in some cases aggregate percentage mean absolute errors clearly below 2%, the best results for the present 3d-nuclei set are in the range between 4 and 5%. Strikingly, the overall best-performing functionals are the recently implemented CDFT versions of two meta-GGAs, namely cM06-L (4.0%) and cVSXC (4.3%), followed by cLH14t-calPBE (4.9%), B3LYP (5.0%), and cLH07t-SVWN (5.1%), i.e., the previously best-performing global hybrid and two local hybrids. A number of further functionals achieve aggregate deviations in the range 5-6%. Range-separated hybrids offer no particular advantage over global hybrids. Due to the overall poor performance of Hartree-Fock theory for all systems except the titanium complexes, MP2 and double-hybrid functionals are unsuitable for these 3d-nucleus shifts and provide large errors. Global hybrid functionals with larger EXX admixtures, such as BHLYP or M06-2X, also perform poorly, and some other highly parametrized global hybrids also are unsuitable. For many functionals depending on local kinetic energy τ, their CDFT variants perform much better than their "non-CDFT" versions. This holds notably also for the above-mentioned M06-L and VSXC, while the effect is small for τ-dependent local hybrids and can even be somewhat detrimental to the agreement with experiment for a few other cases. The separation between well-performing and more poorly performing functionals is mainly determined by their results for the most critical nuclei ⁵⁵Mn, ⁵⁷Fe, and ⁵⁹Co. Here either moderate exact-exchange admixtures or CDFT versions of meta-GGAs are beneficial for the accuracy. The overall deviations of the better-performing global or local hybrids are then typically dominated by the ⁵³Cr shifts, where triplet instabilities appear to disfavor exact-exchange admixture. Further detailed analyses help to pinpoint specific nuclei and specific types of complexes that are challenges for a given functional.

Polyoxometalates in dye-sensitized solar cells

Dye-sensitized solar cells (DSSCs) are the third generation of photovoltaic cells developed by Grätzel and O'Regan. They have the characteristics of low cost, simple manufacturing process, tunable optical properties, and higher photoelectric conversion efficiency (PCE). With an ever increasing energy crisis, there is an urgent need to develop highly efficient, environmentally benign, and energy-saving cell materials. Polyoxometalates (POMs), a kind of molecular inorganic quasi-semiconductor, are promising candidates for use in different parts of DSSCs due to their excellent photosensitivity, redox, and catalytic properties, as well as their relative stability. Following a brief introduction to the development of DSSCs and the potential virtues of POMs in DSSCs, we attempt to make some generalizations about the energy level regulation of POMs that is the underlying theoretical basis for their application in DSSCs, and then we summarize the research progress of POMs in DSSCs in recent years. This is organized in terms of the properties of POMs, namely, electron acceptor, photosensitivity, redox and catalysis, based on the accumulation of our research into POMs over many years. Meanwhile, in view of the fact that the properties of POMs depend primarily on their electronic structural diversity, we keep this point in mind throughout the article with a view to revealing their structure-property relationships. Finally we provide a short summary and remarks on the future outlook. This review may be of interest to synthetic chemists devoted to designing POMs with specific structures, and researchers engaged in the extension of POMs to photoelectric materials.

Relativity-Induced Bonding Pattern Change in Coinage Metal Dimers M2 (M = Cu, Ag, Au, Rg)

The periodic table provides a fundamental protocol for qualitatively classifying and predicting chemical properties based on periodicity. While the periodic law of chemical elements had already been rationalized within the framework of the nonrelativistic description of chemistry with quantum mechanics, this law was later known to be affected significantly by relativity. We here report a systematic theoretical study on the chemical bonding pattern change in the coinage metal dimers (Cu₂, Ag₂, Au₂, Rg₂) due to the relativistic effect on the superheavy elements. Unlike the lighter congeners basically demonstrating ns- ns bonding character and a 0_g⁺ ground state, Rg₂ shows unique 6d-6d bonding induced by strong relativity. Because of relativistic spin-orbit (SO) coupling effect in Rg₂, two nearly degenerate SO states, 0_g⁺ and 2_u, exist as candidate of the ground state. This relativity-induced change of bonding mechanism gives rise to various unique alteration of chemical properties compared with the lighter dimers, including higher intrinsic bond energy, force constant, and nuclear shielding. Our work thus provides a rather simple but clear-cut example, where the chemical bonding picture is significantly changed by relativistic effect, demonstrating the modified periodic law in heavy-element chemistry.

Absolute configuration determination through the unique intramolecular excitonic coupling in the circular dichroisms of o,p′-DDT and o,p′-DDD. A combined experimental and theoretical study

Circular dichroisms (CDs) of theo,p′-isomers of 1,1,1-trichloro- and 1,1-dichloro-2,2-bis(chlorophenyl)ethanes (DDTandDDD) were investigated experimentally and theoretically.

Elucidation of molecular interactions between a DBU based protic ionic liquid and organic solvents: thermophysical and computational studies

Energy profile of 1,8-diazabicyclo[5.4.0]undec-7-en-8-ium trifluoroacetate [DBUTFA].

Geometries of Third-Row Transition-Metal Complexes from Density-Functional Theory

A set of 41 metal-ligand bond distances in 25 third-row transition-metal complexes, for which precise structural data are known in the gas phase, is used to assess optimized and zero-point averaged geometries obtained from DFT computations with various exchange-correlation functionals and basis sets. For a given functional (except LSDA) Stuttgart-type quasi-relativistic effective core potentials and an all-electron scalar relativistic approach (ZORA) tend to produce very similar geometries. In contrast to the lighter congeners, LSDA affords reasonably accurate geometries of 5d-metal complexes, as it is among the functionals with the lowest mean and standard deviations from experiment. For this set the ranking of some other popular density functionals, ordered according to decreasing standard deviation, is BLYP > VSXC > BP86 ≈ BPW91 ≈ TPSS ≈ B3LYP ≈ PBE > TPSSh > B3PW91 ≈ B3P86 ≈ PBE hybrid. In this case hybrid functionals are superior to their nonhybrid variants. In addition, we have reinvestigated the previous test sets for 3d- (Bühl M.; Kabrede, H. J. Chem. Theory Comput. 2006, 2, 1282-1290) and 4d- (Waller, M. P.; Bühl, M. J. Comput. Chem. 2007, 28, 1531-1537) transition-metal complexes using all-electron scalar relativistic DFT calculations in addition to the published nonrelativistic and ECP results. For this combined test set comprising first-, second-, and third-row metal complexes, B3P86 and PBE hybrid are indicated to perform best. A remarkably consistent standard deviation of around 2 pm in metal-ligand bond distances is achieved over the entire set of d-block elements.

The fundamental flaw of the HSAB principle is revealed by a complete speciation analysis of the [PtCl(6-n)Br(n)](2-) (n = 0-6) system

Bjerrum's model of step-wise ligand exchange is extended to compute a complete speciation diagram for the [PtCl6-nBrn](2-) (n = 0-6) system including all 17 equilibrium constants concerning the Pt(IV) chlorido-bromido exchange reaction network (HERN). In contrast to what the hard soft acid base (HSAB) principle "predicts", the thermodynamic driving force for the replacement of chloride by bromide in an aqueous matrix, for each individual ligand exchange reaction present in the Pt(IV) HERN, is due to the difference in halide hydration energy and not bonding interactions present in the acid-base complex. A generalized thermodynamic test calculation was developed to illustrate that the HSAB classified class (b) metal cations Ag(+), Au(+), Au(3+), Rh(3+), Cd(2+), Pt(2+), Pt(4+), Fe(3+), Cd(2+), Sn(2+) and Zn(2+) all form thermodynamically stable halido complexes in the order F(-) ≫ Cl(-) > Br(-) > I(-) irrespective of the sample matrix. The bonding interactions in the acid-base complex, e.g. ionic-covalent σ-bonding, Π-bonding and electron correlation effects, play no actual role in the classification of these metal cations using the HSAB principle. Instead, it turns out that the hydration/solvation energy of halides is the reason why metal cations are categorized into two classes using the HSAB principle which highlights the fundamental flaw of the HSAB principle.

A novel iron complex containing an N,O-type bidentate oxazoline ligand: Synthesis, X-ray studies, DFT calculations and catalytic activity

A five-coordinated Fe(III) complex with the distorted trigonal bipyramidal configuration was synthesized by reactions of FeCl3⋅6H2O and 2-(2'-hydroxyphenyl)oxazoline (Hphox) as a bidentate ON donor oxazoline ligand. Complex [Fe(phox)2Cl] was fully characterized, including by single-crystal X-ray structure analysis. DFT calculations were accompanied with experimental results in order to obtain a deeper insight into the electronic structure and vibrational normal modes of complex. Oxidation of sulfides to sulfoxides in one-step was conducted by this complex as catalyst using urea hydrogen peroxide (UHP) in mixture of CH2Cl2/CH3OH (1:1) under air at room temperature. The results show that using this system in oxidation of sulfides, sulfoxides are obtained as the main products, together with variable amounts of sulfones (≤13%), depending on the nature of the substrate.

Spectroscopic and structural investigations of iron(III) isothiocyanates. A comparative theoretical and experimental study

A combined experimental and theoretical study on the molecular structure and vibrational spectra of [Fe(NCS)](2+) complex in the aqueous solution at the pH∼2 ± 0.1 have been performed. Experimental Raman spectra of the iron(III) isothiocyanate with higher coordination number in the acidic aqueous solution have been analyzed. Molecular modeling of the iron(III) monoisothiocyanate complex was accomplished by the density functional theory (DFT) method using B3LYP and PBE1PBE functionals. Theoretical vibrational spectra of the iron(III) monoisothiocyanate were interpreted by means of the potential energy distributions (PEDs). The influence of different solvation models and position of SO4(2)(-) ligand vs. NCS(-) ligand upon its geometry and vibrational frequencies have been evaluated. The effect of H2O/D2O isotopic substitution on the experimental and calculated Raman spectra of iron(III) isothiocyanates has been examined. Procedures of Raman spectra subtraction have been applied for the extractions of weak and/or obscured Raman signals. As a result, the presence of bound SO4(2)(-) ion and water molecules in the first coordination sphere in the acidic aqueous iron(III) isothiocyanate solution was confirmed. The vibrational assignments for the investigated iron(III) isothiocyanates were proposed here for the first time.

An efficient implementation of two-component relativistic exact-decoupling methods for large molecules

We present an efficient algorithm for one- and two-component relativistic exact-decoupling calculations. Spin-orbit coupling is thus taken into account for the evaluation of relativistically transformed (one-electron) Hamiltonian. As the relativistic decoupling transformation has to be evaluated with primitive functions, the construction of the relativistic one-electron Hamiltonian becomes the bottleneck of the whole calculation for large molecules. For the established exact-decoupling protocols, a minimal matrix operation count is established and discussed in detail. Furthermore, we apply our recently developed local DLU scheme [D. Peng and M. Reiher, J. Chem. Phys. 136, 244108 (2012)] to accelerate this step. With our new implementation two-component relativistic density functional calculations can be performed invoking the resolution-of-identity density-fitting approximation and (Abelian as well as non-Abelian) point group symmetry to accelerate both the exact-decoupling and the two-electron part. The capability of our implementation is illustrated at the example of silver clusters with up to 309 atoms, for which the cohesive energy is calculated and extrapolated to the bulk.

[Pt@Pb12]2– — A challenging system for relativistic density functional theory calculations of 195Pt and 207Pb NMR parameters

We report computations of NMR chemical shifts and indirect spin-spin coupling constants (J couplings) for the [Pt@Pb12]2– “superatom”. The system is strongly influenced by relativistic effects. The Pt–Pb coupling constant is predicted to be negative, with its magnitude being in reasonable agreement with experiment. Pt and Pb chemical shifts also agree reasonably well with experiment. The Pb shielding tensor is strongly anisotropic, with a large deshielding principal component dominated by magnetic coupling between frontier orbitals of the cluster that resemble atomic g orbitals. The NMR parameters are sensitive to approximations made in the computations and require the inclusion of spin-orbit coupling in the theoretical model to achieve reliable results. Computing the NMR parameters of the compact [Pt@Pb12]2– system with its many electrons proves to be a challenging test case for relativistic density functional methods.

A comparison of experimental and DFT calculations of ¹⁹⁵Pt NMR shielding trends for [PtX(n)Y(6-n)](2-) (X, Y = Cl, Br, F and I) anions

A comparison between experimental and calculated gas-phase as well as the conductor-like screening model DFT (195)Pt chemical shifts of a series of octahedral [PtX(6-n)Y(n)](2-) complexes for X = Cl, Br, F, I was carried out to assess the accuracy of computed NMR shielding and to gain insight into the dominant σ(dia), σ(para) and σ(SO) shielding contributions. The discrepancies between the experimental and the DFT-calculated (195)Pt chemical shifts vary for these complexes as a function of the coordinated halide ions, the largest being obtained for the fluorido-chlorido and fluorido-bromido complexes, while negligible discrepancies are found for the [PtCl(6-n)Br(n)](2-) series; the discrepancies are somewhat larger where a significant deviation from the ideal octahedral symmetry such as for the geometric cis/trans or fac/mer isomers of [PtF(6-n)Cl(n)](2-) and [PtF(6-n)Br(n)](2-) may be expected. The discrepancies generally increase in the order [PtCl(6-n)Br(n)](2-) < [PtBr(6-n)I(n)](2-) < [PtCl(6-n)I(n)](2-) < [PtF(6-n)Br(n)](2-) ≈ [PtF(6-n)Cl(n)](2-), and show a striking correlation with the increase in electronegativity difference Δχ between the two halide ligands (X(-) and Y(-)) bound to Pt(IV) for these anions: 0.09 < 0.52 < 0.63 < 1.36 ≈ 1.27, respectively. The computed (195)Pt sensitivity to Pt-X bond displacement, ∂(δ(195)Pt)/∂(ΔPt-X), of these complexes is very large and depends on the halide ion, decreasing from 24 800, 18 300, 15 700 to 12 000 ppm/Å for [PtF(6)](2-), [PtCl(6)](2-), [PtBr(6)](2-) and [PtI(6)](2-), respectively.

Modeling of heavy-atom-ligand NMR spin-spin coupling in solution: molecular dynamics study and natural bond orbital analysis of Hg-C coupling constants

Ab initio molecular dynamics (MD) and relativistic density functional NMR methods were applied to calculate the one-bond Hg-C NMR indirect nuclear spin-spin coupling constants (J) of [Hg(CN)(2) ] and [CH(3) HgCl] in solution. The MD averages were obtained as J((199) Hg-(13) C)=3200 and 1575 Hz, respectively. The experimental Hg-C spin-spin coupling constants of [Hg(CN)(2) ] in methanol and [CH(3) HgCl] in DMSO are 3143 and 1674 Hz, respectively. To deal with solvent effects in the calculations, finite "droplet" models of the two systems were set up. Solvent effects in both systems lead to a strong increase of the Hg-C coupling constant. From a relativistic natural localized molecular orbital (NLMO) analysis, it was found that the degree of delocalization of the Hg 5d(σ) nonbonding orbital and of the HgC bonding orbital between the two coupled atoms, the nature of the trans Hg-C/Cl bonding orbital, and the s character of these orbitals, exhibit trends upon solvation of the complexes that, when combined, lead to the strong increase of J(Hg-C).

NMR line shapes from AB spin systems in solids — The role of antisymmetric spin–spin coupling

NMR parameters such as indirect nuclear spin–spin coupling (J), nuclear magnetic shielding (σ), direct dipolar coupling (D), and electric field gradient (V) are properly described by second-rank tensors. Each may be decomposed into isotropic, symmetric, and antisymmetric components; the number of these three components which may be nonzero is a distinguishing attribute of each interaction tensor. The rank-1 antisymmetric portion of J (Janti) holds the distinction of remaining the only nonzero part of these fundamental NMR interaction tensors which has never been observed experimentally. Accordingly, effects from Janti are usually ignored, but it is important to consider when this is valid. An experimental strategy for observing Janti in powdered samples of tightly coupled homonuclear spin pairs, based on ideas originally presented by Andrew and Farnell ( Mol. Phys. 1968, 15, 157 ), is described. The theory of Andrew and Farnell is extended to powder samples, and methods for analyzing NMR spectra from powdered samples are presented. It is found that, in certain rare cases, Janti has the potential to affect the NMR line shapes from AB spin systems, but that even in these systems, the most intense features of the spectra are not affected and may be analyzed independently of Janti. Furthermore, Janti will only have an observable effect on the NMR spectra when its magnitude is comparable with that of Jiso and with the difference in chemical shifts (in Hz) between the two sites. Finally, the first experimental attempts to measure Janti are reported, and experimental proof that no elements of Janti(119Sn,119Sn) in hexa(p-tolyl)ditin are larger than 2900 Hz is given. The benefits of modern double-quantum filtering NMR pulse sequences in isolating effects from Janti are also illustrated.

Optical signatures of molecular dissymmetry: combining theory with experiments to address stereochemical puzzles

Modern chemistry emerged from the quest to describe the three-dimensional structure of molecules: van't Hoff's tetravalent carbon placed symmetry and dissymmetry at the heart of chemistry. In this Account, we explore how modern theory, synthesis, and spectroscopy can be used in concert to elucidate the symmetry and dissymmetry of molecules and their assemblies. Chiroptical spectroscopy, including optical rotatory dispersion (ORD), electronic circular dichroism (ECD), vibrational circular dichroism (VCD), and Raman optical activity (ROA), measures the response of dissymmetric structures to electromagnetic radiation. This response can in turn reveal the arrangement of atoms in space, but deciphering the molecular information encoded in chiroptical spectra requires an effective theoretical approach. Although important correlations between ECD and molecular stereochemistry have existed for some time, a battery of accurate new theoretical methods that link a much wider range of chiroptical spectroscopies to structure have emerged over the past decade. The promise of this field is considerable: theory and spectroscopy can assist in assigning the relative and absolute configurations of complex products, revealing the structure of noncovalent aggregates, defining metrics for molecular diversity based on polarization response, and designing chirally imprinted nanomaterials. The physical organic chemistry of chirality is fascinating in its own right: defining atomic and group contributions to optical rotation (OR) is now possible. Although the common expectation is that chiroptical response is determined solely by a chiral solute's electronic structure in a given environment, chiral imprinting effects on the surrounding medium and molecular assembly can, in fact, dominate the chiroptical signatures. The theoretical interpretation of chiroptical markers is challenging because the optical properties are subtle, resulting from the strong electric dipole and the weaker electric quadrupole and magnetic dipole perturbations by the electromagnetic field. Moreover, OR arises from a combination of nearly canceling contributions to the electronic response. Indeed, the challenge posed by the chiroptical properties delayed the advent of even qualitatively accurate descriptions for some chiroptical signatures until the past decade when, for example, prediction of the observed sign of experimental OR became accessible to theory. The computation of chiroptical signatures, in close coordination with synthesis and spectroscopy, provides a powerful framework to diagnose and interpret the dissymmetry of chemical structures and molecular assemblies. Chiroptical theory now produces new schemes to elucidate structure, to describe the specific molecular sources of chiroptical signatures, and to assist in our understanding of how dissymmetry is templated and propagated in the condensed phase.