Can model Hamiltonians describe the electron–electron interaction in π-conjugated systems?: PAH and graphene

Model Hamiltonians have been, and still are, a valuable tool for investigating the electronic structure of systems for which mean field theories work poorly. This review will concentrate on the application of Pariser–Parr–Pople (PPP) and Hubbard Hamiltonians to investigate some relevant properties of polycyclic aromatic hydrocarbons (PAH) and graphene. When presenting these two Hamiltonians we will resort to second quantisation which, although not the way chosen in its original proposal of the former, is much clearer. We will not attempt to be comprehensive, but rather our objective will be to try to provide the reader with information on what kinds of problems they will encounter and what tools they will need to solve them. One of the key issues concerning model Hamiltonians that will be treated in detail is the choice of model parameters. Although model Hamiltonians reduce the complexity of the original Hamiltonian, they cannot be solved in most cases exactly. So, we shall first consider the Hartree–Fock approximation, still the only tool for handling large systems, besides density functional theory (DFT) approaches. We proceed by discussing to what extent one may exactly solve model Hamiltonians and the Lanczos approach. We shall describe the configuration interaction (CI) method, a common technology in quantum chemistry but one rarely used to solve model Hamiltonians. In particular, we propose a variant of the Lanczos method, inspired by CI, that has the novelty of using as the seed of the Lanczos process a mean field (Hartree–Fock) determinant (the method will be named LCI). Two questions of interest related to model Hamiltonians will be discussed: (i) when including long-range interactions, how crucial is including in the Hamiltonian the electronic charge that compensates ion charges? (ii) Is it possible to reduce a Hamiltonian incorporating Coulomb interactions (PPP) to an 'effective' Hamiltonian including only on-site interactions (Hubbard)? The performance of CI will be checked on small molecules. The electronic structure of azulene and fused azulene will be used to illustrate several aspects of the method. As regards graphene, several questions will be considered: (i) paramagnetic versus antiferromagnetic solutions, (ii) forbidden gap versus dot size, (iii) graphene nano-ribbons, and (iv) optical properties.

Textural Characterization of Micro- and Mesoporous Carbons Using Combined Gas Adsorption and n-Nonane Preadsorption

Porous carbon and carbide materials with different structures were characterized using adsorption of nitrogen at 77.4 K before and after preadsorption of n-nonane. The selective blocking of the microporosity with n-nonane shows that ordered mesoporous silicon carbide material (OM-SiC) is almost exclusively mesoporous whereas the ordered mesoporous carbon CMK-3 contains a significant amount of micropores (25%). The insertion of micropores into OM-SiC using selective extraction of silicon by hot chlorine gas leads to the formation of ordered mesoporous carbide-derived carbon (OM-CDC) with a hierarchical pore structure and significantly higher micropore volume as compared to CMK-3, whereas a CDC material from a nonporous precursor is exclusively microporous. Volumes of narrow micropores, calculated by adsorption of carbon dioxide at 273 K, are in linear correlation with the volumes blocked by n-nonane. Argon adsorption measurements at 87.3 K allow for precise and reliable calculation of the pore size distribution of the materials using density functional theory (DFT) methods.

Polarized–unpolarized ground state of small polycyclic aromatic hydrocarbons

Do polyacenes, circumacenes, periacenes, nanographenes, and graphene nanoribbons show a spin polarized ground state? In this work, we present monodeterminantal (Hartree–Fock (HF) and density functional theory (DFT) types), and multideterminantal calculations (Møller–Plesset and Coupled Cluster), for several families of unsaturated organic molecules (n-Acenes, n-Periacenes and n-Circumacenes). All HF calculations and many DFT show a spin-polarized (antiferromagnetic) ground state, in agreement with previous calculations. Nevertheless, the multideterminantal calculations, carried out with perturbative and variational wavefunctions, show that the more stable state is obtained starting from the unpolarized HF wavefunction. The trend of the stabilization of wavefunctions (polarized or unpolarized) with respect to exchange and correlation potentials, and to the number of benzene rings, has been analyzed. A study of the spin (〈Ŝ2〉) and the spin density on the carbon atoms has also been carried out.

Spectroelectrochemical Study of the Photoinduced Catalytic Formation of 4,4′-Dimercaptoazobenzene from 4-Aminobenzenethiol Adsorbed on Nanostructured Copper

Surface-enhanced raman scattering (SERS) spectra of self-assembled monolayers of 4-aminobenzenethiol (4-ABT) on copper (Cu) and silver (Ag) surfaces decorated with Cu and Ag nanostructures, respectively, have been obtained with lasers at 532, 632.8, 785, and 1064 nm. Density functional theory (DFT) has been used to obtain calculated vibrational frequencies of the 4-ABT and 4,4′-dimercaptoazobenzene (4,4′-DMAB) molecules adsorbed on model Cu surfaces. The features of the SERS spectra depend on the electrode potential and the type and power density of the laser. SERS spectra showed the formation of the 4,4′-DMAB on the nanostructured Cu surface independently of the laser employed. For the sake of comparison SERS spectra of a self-assembled monolayer of the 4-ABT on Ag surfaces decorated with Ag nanostructures have been also obtained with the same four lasers. When using the 532 and 632.8 nm lasers, the 4,4′-DMAB is formed on Cu surface at electrode potentials as low as −1.0 V (AgCl/Ag) showing a different behavior with respect to Ag (and others metals such as Au and Pt). On the other hand, the surface-enhanced infrared reflection absorption (SEIRA) spectra showed that in the absence of the laser excitation the 4,4′-DMAB is not produced from the adsorbed 4-ABT on nanostructured Cu in the whole range of potentials studied. These results point out the prevalence of the role of electron–hole pairs through surface plasmon activity to explain the obtained SERS spectra.

Large spin splitting in the conduction band of transition metal dichalcogenide monolayers

We study the conduction band spin splitting that arises in transition metal dichalcogenide (TMD) semiconductor monolayers such as MoS2, MoSe2, WS2, and WSe2 due to the combination of spin-orbit coupling and lack of inversion symmetry. Two types of calculation are done. First, density functional theory (DFT) calculations based on plane waves that yield large splittings, between 3 and 30 meV. Second, we derive a tight-binding model that permits to address the atomic origin of the splitting. The basis set of the model is provided by the maximally localized Wannier orbitals, obtained from the DFT calculation, and formed by 11 atomiclike orbitals corresponding to d and p orbitals of the transition metal (W, Mo) and chalcogenide (S, Se) atoms respectively. In the resulting Hamiltonian, we can independently change the atomic spin-orbit coupling constant of the two atomic species at the unit cell, which permits to analyze their contribution to the spin splitting at the high symmetry points. We find that—in contrast to the valence band—both atoms give comparable contributions to the conduction band splittings. Given that these materials are most often n-doped, our findings are important for developments in TMD spintronics.

Electronic properties of transition metal atoms on Cu2N/Cu(100)

We study the nature of spin excitations of individual transition metal atoms (Ti, V, Cr, Mn, Fe, Co, and Ni) deposited on a Cu2N/Cu(100) surface using both spin-polarized density functional theory (DFT) and exact diagonalization of an Anderson model derived from DFT. We use DFT to compare the structural, electronic, and magnetic properties of different transition metal adatoms on the surface. We find that the average occupation of the transition metal d shell, main contributor to the magnetic moment, is not quantized, in contrast with the quantized spin in the model Hamiltonians that successfully describe spin excitations in this system. In order to reconcile these two pictures, we build a zero bandwidth multi-orbital Anderson Hamiltonian for the d shell of the transition metal hybridized with the p orbitals of the adjacent nitrogen atoms, by means of maximally localized Wannier function representation of the DFT Hamiltonian. The exact solutions of this model have quantized total spin, without quantized charge at the d shell. We propose that the quantized spin of the models actually belongs to many-body states with two different charge configurations in the d shell, hybridized with the p orbital of the adjacent nitrogen atoms. This scenario implies that the measured spin excitations are not fully localized at the transition metal.

Derivation of the spin Hamiltonians for Fe in MgO

A method to calculate the effective spin Hamiltonian for a transition metal impurity in a non-magnetic insulating host is presented and applied to the paradigmatic case of Fe in MgO. In the first step we calculate the electronic structure employing standard density functional theory (DFT), based on generalized gradient approximation (GGA), using plane waves as a basis set. The corresponding basis of atomic-like maximally localized Wannier functions is derived and used to represent the DFT Hamiltonian, resulting in a tight-binding model for the atomic orbitals of the magnetic impurity. The third step is to solve, by exact numerical diagonalization, the N electron problem in the open shell of the magnetic atom, including both effects of spin–orbit and Coulomb repulsion. Finally, the low energy sector of this multi-electron Hamiltonian is mapped into effective spin models that, in addition to the spin matrices S, can also include the orbital angular momentum L when appropriate. We successfully apply the method to Fe in MgO, considering both the undistorted and Jahn–Teller (JT) distorted cases. Implications for the influence of Fe impurities on the performance of magnetic tunnel junctions based on MgO are discussed.

Enantioselective Synthesis of Polysubstituted Spiro-nitroprolinates Mediated by a (R,R)-Me-DuPhos·AgF-Catalyzed 1,3-Dipolar Cycloaddition

The synthesis of constrained spirocycles is achieved effectively by means of 1,3-dipolar cyclodditions employing α-imino γ-lactones as azomethine ylide precursors and nitroalkenes as dipolarophiles. The complex formed by (R,R)-Me-DuPhos 18 and AgF is the most efficient bifunctional catalyst. Final spiro-nitroprolinates cycloadducts are obtained in good to moderate yields and both high diastereo- and enantioselectivities. Density functional theory (DFT) calculations supported the expected absolute configuration as well as other stereochemical parameters.

Systematic Improvement of Density Functionals through Parameter-Free Hybridization Schemes

It is suggested here that the ultimate accuracy of DFT methods arises from the type of hybridization scheme followed. This idea can be cast into a mathematical formulation utilizing an integrand connecting the noninteracting and the interacting particle system. We consider two previously developed models for it, dubbed as HYB0 and QIDH, and assess a large number of exchange-correlation functionals against the AE6, G2/148, and S22 reference data sets. An interesting consequence of these hybridization schemes is that the error bars, including the standard deviation, are found to markedly decrease with respect to the density-based (nonhybrid) case. This improvement is substantially better than variations due to the underlying density functional used. We thus finally hypothesize about the universal character of the HYB0 and QIDH models.

Isoprene-mediated lithiation of imidazole derivatives: mechanism considerations

The isoprene-mediated lithiation, with lithium metal, of different imidazole derivatives is an interesting methodology for their functionalization. Studies of different possible intermediates involved in the reaction employing density functional theory calculations, at the B3LYP/6-311++G(d,p) level are considered. A plausible mechanism is described, in which isoprene is reduced, to the corresponding radical anion, in the presence of Li(s), acting then as a base deprotonating N-methylimidazole (NMI) and producing the 1,1-dimethylallyl radical. This radical is further reduced by the excess of lithium proceeding once more as a base. This final step produces stable final products that compensate the previous equilibriums, making favourable the whole process.

Unraveling the performance of dispersion-corrected functionals for the accurate description of weakly bound natural polyphenols

Long-range non-covalent interactions play a key role in the chemistry of natural polyphenols. We have previously proposed a description of supramolecular polyphenol complexes by the B3P86 density functional coupled with some corrections for dispersion. We couple here the B3P86 functional with the D3 correction for dispersion, assessing systematically the accuracy of the new B3P86-D3 model using for that the well-known S66, HB23, NCCE31, and S12L datasets for non-covalent interactions. Furthermore, the association energies of these complexes were carefully compared to those obtained by other dispersion-corrected functionals, such as B(3)LYP-D3, BP86-D3 or B3P86-NL. Finally, this set of models were also applied to a database composed of seven non-covalent polyphenol complexes of the most interest.

Theoretical Rationalization of the Singlet–Triplet Gap in OLEDs Materials: Impact of Charge-Transfer Character

New materials for OLED applications with low singlet–triplet energy splitting have been recently synthesized in order to allow for the conversion of triplet into singlet excitons (emitting light) via a Thermally Activated Delayed Fluorescence (TADF) process, which involves excited-states with a non-negligible amount of Charge-Transfer (CT). The accurate modeling of these states with Time-Dependent Density Functional Theory (TD-DFT), the most used method so far because of the favorable trade-off between accuracy and computational cost, is however particularly challenging. We carefully address this issue here by considering materials with small (high) singlet–triplet gap acting as emitter (host) in OLEDs and by comparing the accuracy of TD-DFT and the corresponding Tamm-Dancoff Approximation (TDA), which is found to greatly reduce error bars with respect to experiments thanks to better estimates for the lowest singlet–triplet transition. Finally, we quantitatively correlate the singlet–triplet splitting values with the extent of CT, using for it a simple metric extracted from calculations with double-hybrid functionals, that might be applied in further molecular engineering studies.

Stretching dependence of the vibration modes of a single-molecule Pt-H2-Pt bridge

A conducting bridge of a single hydrogen molecule between Pt electrodes is formed in a break junction experiment. It has a conductance near the quantum unit, G0=2e2∕h, carried by a single channel. Using point-contact spectroscopy three vibration modes are observed and their variation upon isotope substitution is obtained. The stretching dependence for each of the modes allows uniquely classifying them as longitudinal or transversal modes. The interpretation of the experiment in terms of a Pt-H2-Pt bridge is verified by density-functional theory calculations for the stability, vibrational modes, and conductance of the structure.

Transport in magnetically ordered Pt nanocontacts

Pt nanocontacts, like those formed in mechanically controlled break junctions, are shown to develop spontaneous local magnetic order. Our density functional calculations predict that a robust local magnetic order exists in the atoms presenting low coordination, i.e., those forming the atom-sized neck. We thus find that the electronic transport can be spin polarized, although the net value of the conductance still agrees with available experimental information. Experimental implications of the formation of this new type of nanomagnet are discussed.

Mechanical annealing of metallic electrodes at the atomic scale

The process of creating an atomically defined and robust metallic tip is described and quantified using measurements of contact conductance between gold electrodes and numerical simulations. Our experiments show how the same conductance behavior can be obtained for hundreds of cycles of formation and rupture of the nanocontact by limiting the indentation depth between the two electrodes up to a conductance value of approximately 5G0 in the case of gold. This phenomenon is rationalized using molecular dynamics simulations together with density functional theory transport calculations which show how, after repeated indentations (mechanical annealing), the two metallic electrodes are shaped into tips of reproducible structure. These results provide a crucial insight into fundamental aspects relevant to nanotribology or scanning probe microscopies.