Confinement on the optical response in h-BNCs: Towards highly efficient SERS-active 2D substrates.

Several experimental and theoretical studies have shown that 2D hybrid structures formed by boron, nitrogen and carbon atoms (h-BNCs) possess a highly tunable linear and non-linear optical responses. Recent advances towards the controlled synthesis of these unique structures have motivated an important number of experimental and theoretical work. In this work, the confinement on the optical response induced by boron-nitride (BN) strings in h-BNC 2D structures is investigated using time-dependent density functional theory (TDDFT) and electron density response properties. The number of surrounding BN strings (NBN) necessary to "isolate" the optical modes of a carbon nanoisland (nanographene) from the remaining substrate has been characterized in two different nanoisland models: benzene and pyrene. It was found that for NBN ≥ 3 the excitation wavelengths of the optically active modes remain constant and the changes in the transition densities, the ground to excited state density differences and their associated electron deformation orbitals are negligible and strongly confined within the carbon nanoisland. Using a water molecule as model system, Raman enhancement factors of 10 [6] for the water vibrational modes are obtained when these electromagnetic "hot spots" are activated by an external electromagnetic field. The high tunability of the optical absorption bands of nanographenes through changes in size and morphology makes h-BNCs be perfect materials to construct platforms for surface enhancement Raman spectroscopy (SERS) for a wide range of laser sources.

On the Permeation of Polychlorinated Dibenzodioxins and Dibenzofurans through Lipid Membranes: Classical MD and Hybrid QM/MM-EDA Analysis.

The permeation of dioxin-like pollutants, namely, chlorinated dibenzodioxins and dibenzofurans, through lipid membranes has been simulated using classic molecular dynamics (CMD) combined with the umbrella sampling approach. The most toxic forms of chlorinated dibenzodioxin and dibenzofuran, 2,3,7,8-tetrachloro-p-dibenzodioxin (TCDD) and 2,3,7,8-tetrachlorodibenzofuran (TCDF), and a dioleyl-phosphatidylcholine (DOPC) lipid membrane of 50 Å wide have been chosen for our study. The free energy profile shows the penetration process is largely favoured thermodynamically (ΔG ≈ -12 kcal/mol), with a progressively decrease of the free energy until reaching the energy minima at distances of 8 Å and 9.5 Å from the centre of the membrane for, respectively, TCDD and TCDF. At the centre of the membrane, both molecules display subtle local maxima with free energy differences of 0.5 and 1 kcal/mol with respect to the energy minima for TCDD and TCDF, respectively. Furthermore, the intermolecular interactions between the molecules and the lipid membrane have been characterized at the minima and the local maxima using hybrid quantum mechanics/molecular mechanics energy decomposition analysis (QM/MM-EDA). Total interaction energies of -17.5 and -16.5 kcal/mol have been found at the energy minima for TCDD and TCDF, respectively. In both cases, the dispersion forces govern the molecule-membrane interactions, no significant changes have been found at the local maxima, in agreement with the classical free energy profile. The small differences found in the results obtained for TCDD and TCDF point out that the adsorption and diffusion processes through the cell membrane are not related to the different toxicity shown by these pollutants.

The effect of spin polarization on the electron transport of molecular wires with diradical character.

Some of the most promising materials for application in molecular electronics and spintronics are based on diradical chains. Herein, the proposed relation between increasing conductance with length and diradical character is revisited using ab initio methods that account for the static electron correlation effects. Electron transmission was previously obtained from restricted single determinant wavefuntions or tight-binding approximations, which are unable to account for static correlation. Broken Symmetry Unrestricted Kohn-Sham Density Functional Theory (BS-UKS-DFT) in combination with electron transport analysis based on electron deformation orbitals (EDOs) reflects an exponential decay of the electrical conductance with length. Also, other important effects such as quantum interference are correctly accounted for, leading to a decrease of the conductance as the diradical character increases. As a proof-of-concept, the electrical conductance obtained from BS-UKS-DFT and CASSCF(2,2) wavefunctions were compared in diradical graphene strips in the frame of the pseudo-π approach, obtaining very similar results.

Deciphering the Chemical Basis of Fluorescence of a Selenium-Labeled Uracil Probe when Bound at the Bacterial Ribosomal A-Site.

We unveil in this work the main factors that govern the turn-on/off fluorescence of a Se-modified uracil probe at the ribosomal RNA A-site. Whereas the constraint into an "in-plane" conformation of the two rings of the fluorophore is the main driver for the observed turn-on fluorescence emission in the presence of the antibiotic paromomycin, the electrostatics of the environment plays a minor role during the emission process. Our computational strategy clearly indicates that, in the absence of paromomycin, the probe prefers conformations that show a dark S1 electronic state with participation of nπ* electronic transition contributions between the selenium atom and the π-system of the uracil moiety.

Assessing the Reversed Exponential Decay of the Electrical Conductance in Molecular Wires: The Undeniable Effect of Static Electron Correlation.

An extraordinary new family of molecular junctions, inaccurately referred to as "anti-Ohmic" wires in the recent literature, has been proposed based on theoretical predictions. The unusual electron transport observed for these systems, characterized by a reversed exponential decay of their electrical conductance, might revolutionize the design of molecular electronic devices. This behavior, which has been associated with intrinsic diradical nature, is reexamined in this work. Since the diradical character arises from a near-degeneracy of the frontier orbitals, the employment of a multireference approach is mandatory. CASSCF calculations on a set of nanowires based on polycyclic aromatic hydrocarbons (PAHs) demonstrate that, in the frame of an appropriate multireference treatment, the ground state of these systems shows the expected exponential decay of the conductance. Interestingly, these calculations do evidence a reversed exponential decay of the conductance, although now in several excited states. Similar results have been obtained for other recently proposed candidates to "anti-Ohmic" wires. These findings open new horizons for possible applications in molecular electronics of these promising systems.

Anti-ohmic single molecule electron transport: is it feasible?

Hitherto, only molecular wires with a regular ohmic behavior in which the electric conductance decreases with the wire length have been synthesized. Implementation of molecular conductors with reversed conductance/length trend (anti-ohmic) might revolutionize the field of molecular electronics, allowing the development of electronic devices with extraordinary properties. It is for this reason that, recently, theoretical efforts have been focused on this topic and different structures have been proposed to show reversed conductance/length behavior on the basis of density functional theory non-equilibrium Green function approach (DFT-NEGF) and topological models. From the previous works, it can be stated that an anti-ohmic molecular wire must display a very small HOMO-LUMO gap and a reversed bond alternation pattern in the case of polyenes and related conjugated systems. In this work, the pursuit of a mechanism by which the anti-ohmic electron transport may arise was carried out by studying the paradigmatic anti-ohmic p-xylylene chain (pX2) at the DFT level in combination with topological models. It has been found that the electron transport in the anti-ohmic regime is favored by a long-range superexchange mechanism, which, contrary to what is expected, is reinforced by the increase in the length of the chain. Moreover, strong links between anti-ohmic character in molecular wires and one-dimensional topological insulator models have been established. Due to the small HOMO-LUMO gap predicted at DFT level, the anti-ohmic character has been put to the proof using a multireference scenario. Preliminary results point out to the presence of different ohmic and anti-ohmic electronic states. In the particular case of pX2 the anti-ohmic states do not correspond to the ground state. These findings require a reconsideration of previous studies on the reversed conductance/length behavior using single reference methodologies.

A computational study of the interaction of graphene structures with biomolecular units.

Due to the great interest that biochemical sensors constructed from graphene nanostructures have raised recently, in this work we analyse in detail the electronic factors responsible for the large affinity of biomolecular units for graphene surfaces using ab initio quantum chemical tools based on density functional theory. Both finite and periodic graphene structures have been employed in our study. Whereas the former allows the analysis of the different energy components contributing to the interaction energy separately, the periodic structure provides a more realistic calculation of the total adsorption energy in an extended graphene surface and serves to validate the results obtained using the finite model. In addition, qualitative relations between interaction energy and electron polarization upon adsorption have been established using the finite model. In this work, we have analysed thermodynamically stable adsorption complexes formed by glycine, melamine, pyronin cation, porphine, tetrabenzoporphine and phthalocyanine with a 2D structure of ninety six carbons and periodic structures formed by cells of fifty and seventy two carbons. Differences in the electrostatic, Pauli repulsion, induction and dispersion energies among aromatic and non-aromatic molecules, charged and non-charged molecules and H-π and stacking interactions have been thoroughly analysed in this work.

Tunable aromaticity in bicalicenes.

The unusual aromatic stability of cyclic bicalicene has been suggested to come from a tetraionic structure, where positive and negative charges are located on the cyclopropene and cyclopentadiene rings, respectively. Energetic, magnetic, geometric and electron delocalization analysis performed on a series of bicalicene derivatives, incorporating different electron donating and withdrawing groups, and electrically perturbed bicalicene structures provide additional proof of the role played by this tetraionic structure in the aromatic stability of bicalicene. In this work the aromatic stabilization is chemically and electrically tuned, enhancing or disrupting the electron delocalization and aromatic stability of the cyclopropene and cyclopentadiene rings by increasing or decreasing their corresponding charges. It is shown how the electron delocalization within these rings is similar to that of cyclopropene cation and cyclopentadiene anion for a perfect polarization of one electron.

Theoretical study of the adsorption of aromatic units on carbon allotropes including explicit (empirical) DFT dispersion corrections and implicitly dispersion-corrected functionals: the pyridine case.

The suitability of implicitly dispersion-corrected functionals, namely the M06-2X, for the determination of interaction energies and electron polarization densities in adsorption studies of aromatic molecules on carbon allotropes surfaces is analysed by comparing the results with those obtained using explicit dispersion through Grimme's empirical corrections. Several models of increasing size for the graphene sheet together with one-dimensional curved carbon structures, (5,5), (6,6) and (7,7) armchair single-walled nanotubes, and two-dimensional curved carbon structures, C60 fullerene, have been considered as substrates in this work, whereas pyridine has been chosen as an example for the adsorbed aromatic molecule. Comparison with recent experimental estimations of the adsorption energy and calculations using periodic boundary conditions on a supercell of 72 carbon atoms indicates that a finite model containing ninety six carbon atoms (C96) approaches quite well the adsorption on a graphene sheet. Analysis of the interaction energy components reveals that the M06-2X functional accounts for most of the dispersion energy implicitly, followed far by wB97X and B3LYP, whereas B97 and BLYP do not differ too much from HF. It has been found that M06-2X corrects only the energy component associated to dispersion and leaves the rest, electrostatic, Pauli and induction "unaltered" with respect to the other DFT functionals investigated. Moreover, only the M06-2X functional reflects the effect of dispersion on the electron polarization density, whereas for the remaining functionals the polarization density does not differ too much from the HF density. This makes the former functional more suitable a priori for the calculation of electron density related properties in these adsorption complexes.

Revisiting the Calculation of I/V Profiles in Molecular Junctions Using the Uncertainty Principle.

Ortiz and Seminario (J. Chem. Phys. 2007, 127, 111106/1-3) proposed some years ago a simple and direct approach to obtain I/V profiles from the combination of ab initio equilibrium electronic structure calculations and the uncertainty principle as an alternative or complementary tool to more sophisticated nonequilibrium Green's functions methods. In this work, we revisit the fundamentals of this approach and reformulate accordingly the expression of the electric current. By analogy to the spontaneous electron decay process in electron transitions, in our revision, the current is calculated upon the relaxing process from the "polarized" state induced by the external electric field to the electronic ground state. The electric current is obtained from the total charge transferred through the molecule and the corresponding electronic energy relaxation. The electric current expression proposed is more general compared with the previous expression employed by Ortiz and Seminario, where the charge variation must be tested among different slabs of atoms at the contact. This new approach has been tested on benzene-1,4-dithiolate attached to different gold clusters that represent the contact with the electrodes. Analysis of the total electron deformation density induced by the external electric voltage and properties associated with the electron deformation orbitals supports the conclusions obtained from the I/V profiles.

Analyzing the electric response of molecular conductors using "electron deformation" orbitals and occupied-virtual electron transfer.

The concept of "electron deformation orbitals" (EDOs) is used to investigate the electric response of conducting metals and oligophenyl chains. These orbitals and their eigenvalues are obtained by diagonalization of the deformation density matrix (difference between the density matrices of the perturbed and unperturbed systems) and can be constructed as linear combinations of the unperturbed molecular orbitals within "frozen geometry" conditions. This form of the EDOs allows calculating the part of the electron deformation density associated to an effective electron transfer from occupied to virtual orbitals (valence to conduction band electron transfer in the band model of conductivity). It is found that the "electron deformation" orbitals pair off, displaying the same eigenvalue but opposite sign. Each pair represents an amount of accumulation/depletion of electron charge at different molecular regions. In the oligophenyl systems investigated only one pair contributes effectively to the charge flow between molecular ends, resulting from the promotion of electrons from occupied orbitals to close in energy virtual orbitals of appropriate symmetry and overlapping. Analysis of this pair along explains the differences in conductance of olygophenyl chains based on phenyl units.

Aromaticity of closed-shell charged polybenzenoid hydrocarbons.

The aromatic stabilization of closed-shell charged polybenzenoid hydrocarbons (PBHs) has been scrutinized by means of energetic and magnetic aromaticity criteria and by direct measures of electron delocalization. Thus, topological resonance energies and their circuit contributions, ring current maps, and multicenter delocalization indices have been calculated for a series of 18 polybenzenoid cations containing from 3 to 10 benzene rings. All calculations indicate that the closed-shell cations have a similar degree of aromaticity compared to that of the corresponding closed-shell neutral PBHs. All cations investigated display a large degree of electronic delocalization in the ring, accompanied by significant aromatic stabilization and a strong diatropic peripheral electron current. Graph theoretical models describe perfectly the aromatic features of these hydrocarbon fragments, showing how they can be understood as a superposition of specific neutral PBHs. The large aromatic character of these systems suggests they may be relatively stable upon formation at combustion conditions, like those given in the interstellar medium. It has been postulated that closed-shell fragments of PBHs may play an important role in the photoluminescent phenomenon known as extended red emission.

Electronic properties of p-xylylene and p-phenylene chains subjected to finite bias voltages: a new highly conducting oligophenyl structure.

Recently, experimental and theoretical determination of electric currents induced by finite bias voltages in p-xylylene chains attached to gold contacts revealed higher conductance of these systems in comparison with p-phenylene homologous chains. To gain more insight into the conducting properties of these oligophenyl structures, ab initio studies were carried out on the electronic properties of two different p-xylylene-like chains (pX1 and pX2) and the p-phenylene (pP) chain attached to gold contacts, with molecular formulas AuCH2 (C6 H4 )n CH2 Au (n=1-5), Au2 C(C6 H4 )n CAu2 (n=1-5), and Au(C6 H4 )n Au (n=1-5), respectively. The molecules were subjected to finite bias voltages ranging from 0 to 5 V. Analysis of the intramolecular electron transfer and electron delocalization revealed a completely opposite response to electric perturbation of pX2 in comparison with pX1 and pP. Thus, in pX2 the applied voltage causes an increase in the electron delocalization within the rings together with a large electron transfer and energetic stabilization. On the contrary, the same voltages partially destroy the electron delocalization in pX1 and pP, produce a large local electron polarization in the benzene rings, and a smaller energetic stabilization. These differences can be rationalized in terms of the role played by polarized valence bond structures in the total wave function. Theoretical estimation of the I/V profiles indicates that pX2 chains are much better electronic conductors than pX1 and pP.