Accurate computations of the structures and binding energies of the imidazole ... benzene and pyrrole ... benzene complexes

Using explicitly-correlated coupled-cluster theory with single and double excitations, the intermolecular distances and interaction energies of the T-shaped imidazole⋯⋯benzene and pyrrole⋯⋯benzene complexes have been computed in a large augmented correlation-consistent quadruple-zeta basis set, adding also corrections for connected triple excitations and remaining basis-set-superposition errors. The results of these computations are used to assess other methods such as Møller–Plesset perturbation theory (MP2), spin-component-scaled MP2 theory, dispersion-weighted MP2 theory, interference-corrected explicitly-correlated MP2 theory, dispersion-corrected double-hybrid density-functional theory (DFT), DFT-based symmetry-adapted perturbation theory, the random-phase approximation, explicitly-correlated ring-coupled-cluster-doubles theory, and double-hybrid DFT with a correlation energy computed in the random-phase approximation.

Highly-effective gating of single-molecule junctions: an electrochemical approach

We report an electrochemical gating approach with [similar]100% efficiency to tune the conductance of single-molecule 4,4′-bipyridine junctions using scanning-tunnelling-microscopy break junction technique. Density functional theory calculation suggests that electrochemical gating aligns molecular frontier orbitals relative to the electrode Fermi-level, switching the molecule from an off resonance state to “partial” resonance.

Toward Cove-Edged Low Band Gap Graphene Nanoribbons

Graphene nanoribbons (GNRs), defined as nanometer-wide strips of graphene, have attracted increasing attention as promising candidates for next-generation semiconductors. Here, we demonstrate a bottom-up strategy toward novel low band gap GNRs (E-g = 1.70 eV) with a well-defined cove-type periphery both in solution and on a solid substrate surface with chrysene as the key monomer. Corresponding cyclized chrysene-based oligornerS consisting of the dimer and tetramer are obtained via an Ullmann Coupling followed by oxidative intramolecular cyclodehydrogenation in solution, and much higher GNR homologues via on-surface synthesis. These oligomers adopt nonplanar structures due to the isteric repulsion between the two C-H bonds at the inner cove position. Characterizations by single crystal X-ray analysis, UV-vis absorption spectroscopy, NMR spectroscopy, and scanning tunneling microscopy (STM) are described. The interpretation is assisted by density functional theory (DFT) calculations.

Competition and cooperation between antiferrodistortive and ferroelectric instabilities in the model perovskite SrTiO₃

We use density functional theory to explore the interplay between octahedral rotations and ferroelectricity in the model compound SrTiO3. We find that over the experimentally relevant range octahedral rotations suppress ferroelectricity as is generally assumed in the literature. Somewhat surprisingly, we observe that at larger angles the previously weakened ferroelectric instability strengthens significantly. By analyzing geometry changes, energetics, force constants and charges, we explain the mechanisms behind this transition from competition to cooperation with increasing octahedral rotation angle.

Adsorption of biomedical coating molecules, amino acids, and short peptides on magnetite (110)

Superparamagnetic iron oxide nanoparticles for biomedical applications are usually coated with organic molecules to form a steric barrier against agglomeration. The stability of these coatings is well established in the synthesis medium but is more difficult to assess in physiological environment. To obtain a first theoretical estimate of their stability in such an environment, we perform density functional theory calculations of the adsorption of water, polyvinyl alcohol (PVA) and polyethylene glycol (PEG) coating molecules, as well as the monomer and dimer of glycine as a prototype short peptide, on the (110) surface of magnetite (Fe3O4) in vacuo. Our results show that PVA binds significantly stronger to the surface than both PEG and glycine, while the difference between the latter two is quite small. Depending on the coverage, the wateradsorption strength is intermediate between PVA and glycine. Due to its strongly interacting OH side groups, PVA is likely to remain bound to the surface in the presence of short peptides. This stability will have to be further assessed by molecular dynamics in the solvated state for which the present work forms the basis.

Distributed Atomic Polarizabilities of Amino Acids and their Hydrogen-Bonded Aggregates

With the purpose of rational design of optical materials, distributed atomic polarizabilities of amino acid molecules and their hydrogen-bonded aggregates are calculated in order to identify the most efficient functional groups, able to buildup larger electric susceptibilities in crystals. Moreover, we carefully analyze how the atomic polarizabilities depend on the one-electron basis set or the many-electron Hamiltonian, including both wave function and density functional theory methods. This is useful for selecting the level of theory that best combines high accuracy and low computational costs, very important in particular when using the cluster method to estimate susceptibilities of molecular-based materials.

Reaction of Copper(II) Chloroacetate with Pyrazole. Synthesis of a One-Dimensional Coordination Polymer and Unexpected Dehydrochlorination Reaction

The reaction of copper(II) chloroacetate (1d) with pyrazole (Hpz) mainly yielded the mononuclear compound [Cu(μ-ClCH2COO)2(Hpz)2] (2m), which self-assembled generating a one-dimensional coordination polymer. Moreover, the concomitant isolation of the tetranuclear [{Cu2(μ-pz)(μ-OCH2COO)(Hpz)(MeOH)}2(μ-ClCH2COO)2] (3t) and hexanuclear [{Cu3(μ3-OH)(μ-pz)3(Hpz)2}2(μ-ClCH2COO)2](Cl)2 (4h) species evidenced the occurrence of a peculiar, previously unreported, dehydrochlorination reaction and the formation of the trinuclear triangular moiety [Cu3(μ3-OH)(μ-pz)3]. Theoretical calculations based on density functional theory including solvation effects indicate a possible pathway for the formation of 3t. Interestingly, besides the energy minimum corresponding to 3t, a further relative energy minimum is found for a species which can be considered a possible reaction intermediate.

On-Surface Synthesis of BN-Substituted Heteroaromatic Networks

We report on the bottom-up fabrication of BN-substituted heteroaromatic networks achieved by surface-assisted polymerization and subsequent cyclodehydrogenation of specifically designed BN-substituted precursor monomers based on a borazine core structural element. To get insight into the cyclodehydrogenation pathway and the influence of molecular flexibility on network quality, two closely related precursor monomers with different degrees of internal cyclodehydrogenation have been employed. Scanning tunneling microscopy shows that, for both monomers, surface-assisted cyclodehydrogenation allows for complete monomer cyclization and the formation of covalently interlinked BN-substituted polyaromatic hydrocarbon networks on the Ag(111) surface. In agreement with experimental observations, density functional theory calculations reveal a significantly lower energy barrier for the cyclodehydrogenation of the conformationally more rigid precursor monomer, which is also reflected in a higher degree of long-range order of the obtained heteroaromatic network. Our proof-of-concept study will allow for the fabrication of atomically precise substitution patterns within BNC heterostructures.

Experimental and Calculated Spectra of π-Stacked Mild Charge-Transfer Complexes: Jet-Cooled Perylene·(Tetrachloroethene) n , n = 1,2

The S0 ↔ S1 spectra of the mild charge-transfer (CT) complexes perylene·tetrachloroethene (P·4ClE) and perylene·(tetrachloroethene)2 (P·(4ClE)2) are investigated by two-color resonant two-photon ionization (2C-R2PI) and dispersed fluorescence spectroscopy in supersonic jets. The S0 → S1 vibrationless transitions of P·4ClE and P·(4ClE)2 are shifted by δν = −451 and −858 cm–1 relative to perylene, translating to excited-state dissociation energy increases of 5.4 and 10.3 kJ/mol, respectively. The red shift is ∼30% larger than that of perylene·trans-1,2-dichloroethene; therefore, the increase in chlorination increases the excited-state stabilization and CT character of the interaction, but the electronic excitation remains largely confined to the perylene moiety. The 2C-R2PI and fluorescence spectra of P·4ClE exhibit strong progressions in the perylene intramolecular twist (1au) vibration (42 cm–1 in S0 and 55 cm–1 in S1), signaling that perylene deforms along its twist coordinate upon electronic excitation. The intermolecular stretching (Tz) and internal rotation (Rc) vibrations are weak; therefore, the P·4ClE intermolecular potential energy surface (IPES) changes little during the S0 ↔ S1 transition. The minimum-energy structures and inter- and intramolecular vibrational frequencies of P·4ClE and P·(4ClE)2 are calculated with the dispersion-corrected density functional theory (DFT) methods B97-D3, ωB97X-D, M06, and M06-2X and the spin-consistent-scaled (SCS) variant of the approximate second-order coupled-cluster method, SCS-CC2. All methods predict the global minima to be π-stacked centered coplanar structures with the long axis of tetrachloroethene rotated by τ ≈ 60° relative to the perylene long axis. The calculated binding energies are in the range of −D0 = 28–35 kJ/mol. A second minimum is predicted with τ ≈ 25°, with ∼1 kJ/mol smaller binding energy. Although both monomers are achiral, both the P·4ClE and P·(4ClE)2 complexes are chiral. The best agreement for adiabatic excitation energies and vibrational frequencies is observed for the ωB97X-D and M06-2X DFT methods.

Oxygen Reduction Reaction on Pt Overlayers Deposited onto a Gold Film: Ligand, Strain, and Ensemble Effect

We study the oxygen reduction reaction (ORR), the catalytic process occurring at the cathode in fuel cells, on Pt layers prepared by electrodeposition onto an Au substrate. Using a nominal Pt layer by layer deposition method previously proposed, imperfect layers of Pt on Au are obtained. The ORR on deposited Pt layers decreases with increasing Pt thickness. In the submonolayer region, however, the ORR activity is superior to that of bulk Pt. Using density functional theory (DFT) calculations, we correlate the observed activity trend to strain, ligand, and ensemble effects. At submonolayer coverage certain atom configurations weaken the binding energies of reaction intermediates due to a ligand and ensemble effect, thus effectively increasing the ORR activity. At higher Pt coverage the activity is governed by a strain effect, which lowers the activity by decreasing the oxidation potential of water. This study is a nice example of how the influence of strain, ligand, and ensemble effects on the ORR can be deconvoluted.

Correlation between Accurate Electron Density and Linear Optical Properties in Amino Acid Derivatives: L-Histidinium Hydrogen Oxalate

The accurate electron density and linear optical properties of L-histidinium hydrogen oxalate are discussed. Two high-resolution single crystal X-ray diffraction experiments were performed and compared with density functional calculations in the solid state as well as in the gas phase. The crystal packing and the hydrogen bond network are accurately investigated using topological analysis based on quantum theory of atoms in molecules, Hirshfeld surface analysis, and electrostatic potential mapping. The refractive indices are computed from couple perturbed Kohn-Sham calculations and measured experimentally. Moreover, distributed atomic polarizabilities are used to analyze the origin of the linear susceptibility in the crystal, in order to separate molecular and intermolecular causes. The optical properties are also correlated with the electron density distribution. This compound also offers the possibility to test the electron density building block approach for material science and different refinement schemes for accurate positions and displacement parameters of hydrogen atoms, in the absence of neutron diffraction data.

Strain-induced structural instability in FeRh

We perform density functional calculations to investigate the structure of the intermetallic alloy FeRh under epitaxial strain. Bulk FeRh exhibits a metamagnetic transition from a low-temperature antiferromagnetic (AFM) phase to a ferromagnetic phase at 350 K, and its strain dependence is of interest for tuning the transition temperature to the room-temperature operating conditions of typical memory devices. We find an unusually strong dependence of the structural energetics on the choice of exchange-correlation functional, with the usual local density approximation yielding the wrong ground-state structure, and generalized gradient (GGA) extensions being in better agreement with the bulk experimental structure. Using the GGA we show the existence of a metastable face-centered-cubic-like AFM structure that is reached from the ground-state body-centered-cubic-like AFM structure by following the epitaxial Bain path. We show that the behavior is well described using nonlinear elasticity theory, which captures the softening and eventual sign change of the orthorhombic shear modulus under compressive strain, consistent with this structural instability. Finally, we predict the existence of an additional unit-cell-doubling lattice instability, which should be observable at low temperature.