Multiscale Modeling of Broadband Perfect Absorbers Based on Gold Metallic Molecules

The modeling of functional photonic devices that rely on molecular materials continues to be a major contemporary challenge. It is a challenge because, in the Maxwell equations, which govern the light-matter interaction, material properties are primarily introduced on phenomenological grounds and not from first principles. To overcome such limitations, we outline a multiscale modeling approach that bridges multiple length scales. We can predict with our approach the optical response of a photonic device that exploits in its design molecular materials whose properties were determined using time-dependent density functional theory. The specifically considered device is a broadband perfect absorber that uses in part a thin film comprising gold molecules made from 144 atoms. Our methodology discloses various chemical and physical effects that define such a device’s response. Our methodology is versatile, and a larger number of applications will profit from this development

Robust relativistic many-body Green’s function based approaches for assessing core ionized and excited states

A two-component contour deformation (CD) based GW method that employs frequency sampling to drastically reduce the computational effort when assessing quasiparticle states far away from the Fermi level is outlined. Compared to the canonical CD-GW method, computational scaling is reduced by an order of magnitude without sacrificing accuracy. This allows for an efficient calculation of core ionization energies. The improved computational efficiency is used to provide benchmarks for core ionized states, comparing the performance of 15 density functional approximations as Kohn–Sham starting points for GW calculations on a set of 65 core ionization energies of 32 small molecules. Contrary to valence states, GW calculations on core states prefer functionals with only a moderate amount of Hartree–Fock exchange. Moreover, modern ab initio local hybrid functionals are also shown to provide excellent generalized Kohn–Sham references for core GW calculations. Furthermore, the core–valence separated Bethe–Salpeter equation (CVS-BSE) is outlined. CVS-BSE is a convenient tool to probe core excited states. The latter is tested on a set of 40 core excitations of eight small inorganic molecules. Results from the CVS-BSE method for excitation energies and the corresponding absorption cross sections are found to be in excellent agreement with those of reference damped response BSE calculations

The GW/BSE Method in Magnetic Fields

The GW approximation and the Bethe–Salpeter equation have been implemented into the Turbomole program package for computations of molecular systems in a strong, finite magnetic field. Complex-valued London orbitals are used as basis functions to ensure gauge-invariant computational results. The implementation has been benchmarked against triplet excitation energies of 36 small to medium-sized molecules against reference values obtained at the approximate coupled-cluster level (CC2 approximation). Finally, a spectacular change of colour from orange to green of the tetracene molecule is induced by applying magnetic fields between 0 and 9,000 T perpendicular to the molecular plane

Current density functional framework for spin–orbit coupling

Relativistic two-component density functional calculations are carried out in a non-collinear formalism to describe spin–orbit interactions, where the exchange–correlation functional is constructed as a generalization of the non-relativistic density functional approximation. Contrary to non-relativistic density functional theory (DFT), spin–orbit coupling, however, leads to a non-vanishing paramagnetic current density. Density functionals depending on the kinetic energy density, such as meta-generalized gradient approximations, should therefore be constructed in the framework of current DFT (CDFT). The latter has previously exclusively been used in the regime of strong magnetic fields. Herein, we present a consistent CDFT approach for relativistic DFT, including spin–orbit coupling. Furthermore, we assess the importance of the current density terms for ground-state energies, excitation energies, nuclear magnetic resonance shielding, and spin–spin coupling constants, as well as hyperfine coupling constants, $\Delta$g-shifts, and the nuclear quadrupole interaction tensor in electron paramagnetic resonance (EPR) spectroscopy. The most notable changes are found for EPR properties. The impact of the current-dependent terms rises with the number of unpaired electrons, and consequently, the EPR properties are more sensitive toward CDFT. Considerable changes are observed for the strongly constrained and appropriately normed functionals, as well as the B97M family and TASK. The current density terms are less important when exact exchange is incorporated. At the same time, the current-dependent kernel ensures the stability of response calculations in all cases. We, therefore, strongly recommend to use the framework of CDFT for self-consistent spin–orbit calculations

Synthesis of New Donor‐Substituted Biphenyls: Pre‐ligands for Highly Luminescent (C^C^D) Gold(III) Pincer Complexes

We herein report on new synthetic strategies for the preparation of pyridine and imidazole substituted 2,2’‐dihalo biphenyls. These structures are pre‐ligands suitable for the preparation of respective stannoles. The latter can successfully be transmetalated to K[AuCl$_{4}$] forming non‐palindromic [(C^C^D)Au$^{III}$] pincer complexes featuring a lateral pyridine (D=N) or N‐heterocyclic carbene (NHC, D=C’) donor. The latter is the first report on a pincer complex with two formally anionic sp$^{2}$ and one carbenic carbon donor. The [(C^C^D)Au$^{III}$] complexes show intense phosphorescence in solution at room temperature. We discuss the developed multistep strategy and touch upon synthetic challenges. The prepared complexes have been fully characterized including X‐ray diffraction analysis. The gold(III) complexes’ photophysical properties have been investigated by absorption and emission spectroscopy as well as quantum chemical calculations on the quasi‐relativistic two‐component TD‐DFT and GW/Bethe–Salpeter level including spin–orbit coupling. Thus, we shed light on the electronic influence of the non‐palindromic pincer ligand and reveal non‐radiative relaxation pathways of the different ligands employed

A Multi-Scale Approach to Simulate the Nonlinear Optical Response of Molecular Nanomaterials

Nonlinear optics is essential for many recent photonic technologies. Here, we introduce a novel multi-scale approach to simulate the nonlinear optical response of molecular nanomaterials combining ab initio quantum-chemical and classical Maxwell-scattering computations. In this approach, the first hyperpolarizability tensor is computed with time-dependent density-functional theory and translated into a multi-scattering formalism that considers the optical interaction between neighboring molecules. A novel object is introduce to perform this transition from quantum-chemistry to classical scattering theory: the Hyper-Transition(T)-matrix. With this object at hand, the nonlinear optical response from single molecules and also from entire photonic devices can be computed, incorporating the full tensorial and dispersive nature of the optical response of the molecules. To demonstrate the applicability of our novel approach, the generation of a second-harmonic signal from a thin film of a Urea molecular crystal is computed and compared to more traditional simulations. Furthermore, an optical cavity is designed, which enhances the second-harmonic response of the molecular film by more than two orders of magnitude. Our approach is highly versatile and accurate and can be the working horse for the future exploration of nonlinear photonic molecular materials in structured photonic environments

GW quasiparticle energies of atoms in strong magnetic fields

Quasiparticle energies of the atoms H–Ne have been computed in the GW approximation in the presence of strong magnetic fields with field strengths varying from 0 to 0.25 atomic units (0.25 B 0 =0.25 ℏe −1 a −2 0 ≈58 763 0.25 B0=0.25 ℏe−1a0−2≈58 763 T). The GW quasiparticle energies are compared with equation-of-motion ionization-potential (EOM-IP) coupled-cluster singles-and-doubles (CCSD) calculations of the first ionization energies. The best results are obtained with the evGW@PBE0 method, which agrees with the EOM-IP-CCSD model to within about 0.20 eV. Ionization potentials have been calculated for all atoms in the series, representing the first systematic study of ionization potentials for the first-row atoms at field strengths characteristic of magnetic white dwarf stars. Under these conditions, the ionization potentials increase in a near-linear fashion with the field strength, reflecting the linear field dependence of the Landau energy of the ionized electron. The calculated ionization potentials agree well with the best available literature data for He, Li, and Be

Supplementary data for article: Peschel, L. M.; Holzer, C.; Mihajlovic-Lalic, L.; Belaj, F.; Mösch-Zanetti, N. C. Coordinative Flexibility of a Thiophenolate Oxazoline Ligand in Nickel(II), Palladium(II), and Platinum(II) Complexes. European Journal of Inorganic Chemistry 2015, 2015 (9), 1569–1578. https://doi.org/10.1002/ejic.201403108

Supporting information for: [https://doi.org/10.1002/ejic.201403108]Related to published version: [http://cherry.chem.bg.ac.rs/handle/123456789/1677

A T‐Matrix Based Approach to Homogenize Artificial Materials

The accurate and efficient computation of the electromagnetic response of objects made from artificial materials is crucial for designing photonic functionalities and interpreting experiments. Advanced fabrication techniques can nowadays produce new materials as 3D lattices of scattering unit cells. Computing the response of objects of arbitrary shape made from such materials is typically computationally prohibitive unless an effective homogeneous medium approximates the discrete material. In here, a homogenization method based on the effective transition (T-)matrix, $T_{eff}$ is introduced. Such a matrix captures the exact response of the discrete material, is determined by the T-matrix of the isolated unit cell and the material lattice vectors, and is free of spatial dispersion. The truncation of $T_{eff}$ to dipolar order determines the common bi-anisotropic constitutive relations. When combined with quantum-chemical and Maxwell solvers, the method allows one to compute the response of arbitrarily-shaped volumetric patchworks of structured molecular materials and metamaterials