A theoretical study of the photochemistry of 1,3-cyclopentadiene and its cyano derivatives bound to a water dimer: Assessing reactivity of ionized clusters and possible photoproducts.

We theoretically investigate the photoionization scenarios of molecular complexes involving cyclopentadiene and cyanocyclopentadiene bound to water dimers. Using electronic structure calculations within density-functional theory (DFT) and time dependent DFT (TD-DFT), we explore the potential photochemical pathways following ionization, and determine the charge transfer excitations related to the possible subsequent reactions. Our findings suggest that the investigated photochemical pathways of the hydrated complexes take place in two well-defined ultraviolet regions: (i) 8.2-9.5 eV for the cyclic compounds and (ii) 11.2-11.4 eV for the bound water dimer. We quantify how H-bonding effects can influence the photoionization channels. Before forming possible photoproducts, we also examine the regiospecificity of OH addition to 1,3-cyclopentadiene and its cyano derivatives We analyze our results in light of photoionization studies of jet-cooled molecular complexes and possible implications in astrochemical environments.

Influence of Solvents and Halogenation on ESIPT of Benzimidazole Derivatives for Designing Turn-on Fluorescence Probes.

This work reports a theoretical investigation of the solvent polarity as well as the halogenation of benzimidazole derivatives during excited state intramolecular proton transfer (ESIPT). It details how the environment and halogen substitution may contribute to the efficiency of ESIPT upon keto-enol tautomerism and exploits this effect to design fluorescence sensing. For this purpose, we first examine the conformational equilibrium of benzimidazole derivatives containing different halogen atoms, which results in intramolecular proton transfer, using density-functional theory (DFT) combined with the polarizable continuum model (PCM). Then we evaluate the fluorescence of the benzimidazole derivatives in different dielectric constants within time-dependent DFT (TD-DFT) approaches. Our results quantitatively allow the determination of large Stokes shifts in nonpolar solvents around 100 nm. These theoretical results are in agreement with experimental solvatochromism studies of benzimidazoles. The effect of halogenation, with fluorine, chlorine, and bromine, is less important than solvent polarization when ESIPT takes place. Thus, halogenation can be properly chosen depending on the interest of the synthesis of benzimidazole-based turn-on fluorescence in appropriate solvents.

Predicting pressure-dependent rate constants for the furan + OH reactions and their impact under tropospheric conditions.

In this work, we have evaluated the influence of temperature and pressure on the mechanism of furan oxidation by the OH radical. The stationary points on the potential energy surface were described at the M06-2X/aug-cc-pVTZ level of theory. In the kinetic treatment at the high-pressure limit (HPL), we have combined the multistructural canonical variational theory with multidimensional small-curvature tunneling corrections and long-range transition state theory. The system-specific quantum Rice-Ramsperger-Kassel theory was employed to estimate the pressure-dependent rate. In the HPL, the OH addition on the α carbon is the dominant pathway in the mechanism, producing a product via the ring-opening process, also confirmed by the product branching ratio calculations. The overall rate constant, obtained by a kinetic Monte Carlo simulation, reads the form koverall=5.22×10-13T/3001.10⁡exp1247(K/T) and indicates that the furan oxidation by OH radicals is a pressure-independent reaction under tropospheric conditions.

Regioselectivity in the Nitration of Eugenol Is Independent of Inorganic Reagents: An Experimental and Theoretical Investigation with Synthetic and Mechanistic Implications.

In this study, we reinvestigated the straightforward nitration of eugenol using traditional reagents and bismuth nitrate. NMR analysis of the obtained products revealed that the regioselectivity of eugenol nitration was independent of the inorganic nitrating reagent used, consistently resulting in the formation of 6-nitroeugenol. This contradicts previous literature reports because the elusive synthesis of 5-nitroeugenol using Bi(NO3)3·5H2O was not achievable through straightforward methods; instead, this isomer could only be prepared via the well-established three-step synthesis. Theoretical investigations using DFT calculations, considering both the dielectric constant of the medium and explicit water molecules, substantiated this regioselectivity. It was found that hydration water played a critical role in the formation of a Zundel cation, shifting the thermodynamic equilibrium toward the exclusive production of 6-nitroeugenol. These results imply that all biological studies involving eugenol derivatives synthesized via direct nitration with Bi(NO3)3·5H2O should be reviewed, as they dealt with 6-substituted eugenol derivatives rather than the previously assumed 5-substituted eugenol.

Photoabsorption of Microhydrated Naphthalene and Its Cyano-Substituted Derivatives: Probing Prereactive Models for Photodissociation in Molecular Clouds.

We investigate the photoionization pathways of naphthalene, 1-cyanonaphthalene, and 2-cyanonaphthalene upon complexation with the water dimer, aiming to understand the photodissociation process under conditions of the interstellar medium (ISM). We analyze the intermolecular bonding pattern, equilibrium rotational properties, energy complexation, far-IR spectra, and ionic trends of the possible photoproducts using dispersion-corrected density functional theory (DFT-D) and time-dependent DFT (TD-DFT). For the different configurations, we evaluate the possible charge-transfer (CT) excitations near the photoionization limit. Our results indicate that, in high-radiation regions of the ISM (>8.0 eV), CT excitations occur from localized occupied molecular orbitals (MOs) in the aromatic molecules to mixed unoccupied MOs in the complexes, favoring cationic aromatic species in these conditions. We notice that the photoabsorption spectra depend on the type of intermolecular interaction (H-bonds or O-H···π bonds) in the complexes, as well as the presence and position (1 or 2) of the cyano-functional group in naphthalene. For hydrated naphthalene, the O-H···π complexes assume a more relevant role for photodissociation. In the case of the cyano-substituted derivatives, the H-bonded structures are more favorable to be considered as prereactive models. However, the cyano group at position 2 indicates that CT excitations toward the water dimer are more likely to occur.

Optical properties of organosilicon compounds containing sigma-electron delocalization by quasiparticle self-consistent GW calculations.

We investigate theoretically the electronic and optical absorption properties of two sub-classes of oligosilanes: (i) Si(CH3)4, Si4(CH3)8, and Si8(CH3)8 that contain Si dot, ring and cage, respectively, and exhibit typical SiC and SiSi bonds; and (ii) persilastaffanes Si7H6(CH3)6 and Si12H6(CH3)12, which contain extended delocalized σ-electrons in SiSi bonds over three-dimensional Si frameworks. Our modeling is performed within the GW approach up to the partially self-consistent GW0 approximation, which is more adequate for reliably predicting the optical band gaps of materials. We examine how the optical properties of these organosilicon compounds depend on their size, geometric features, and Si/C composition. Our results indicate that the present methodology offers a viable way of describing the optical excitations of tailored functional Si-C-based clusters and molecular optical tags with potential use as efficient light absorbers/emitters in molecular optical devices.

A Theoretical Assessment of Spin and Charge States in Binuclear Cobalt-Ruthenium Complexes: Implications for a Creutz-Taube Model Ion Separated by a C60-Derivative Bridging Ligand.

We investigate the spin-state energetics and the role of ionic charges in the electronic configuration of binuclear complexes of the form [(NH3)5Co(py)-X-(py)Ru(NH3)5]q+. In these compounds with q = 4-6, py = pyridine, and X = C≡C and C60, the Co-Ru distance varies from ∼1.4 to ∼2.1 nm. We carry out a systematic electronic structure calculation using different exchange-correlation (xc) approaches within spin-density functional theory, which are largely employed to investigate the properties of a variety of coordination complexes. To evaluate the effects of spin states and type of spacer in the bridging ligand on the valence tautomerism between Co2+/3+ and Ru2+/3+, we examine in more detail the case of Creutz-Taube-type ions [(NH3)5Co(py)-X-(py)Ru(NH3)5]5+. Our analysis shows that the stabilization of low- and high-spin states critically depends on the total charge of the complex, type of X-bridged ligand, and employed xc approach to calculate the electron spin density. Importantly, the C60-bridged group may result in a blockage of the valence tautomerism of the Creutz-Taube complex, inducing bistable charge configurations. Overall, our results also show that an adiabatic description in terms of the frontier molecular spin-orbitals for analyzing the distinct spin-charge states of these complexes may dramatically depend on the density-functional description.

Modeling the Field Emission Enhancement Factor for Capped Carbon Nanotubes Using the Induced Electron Density.

In many field electron emission experiments on single-walled carbon nanotubes (SWCNTs), the SWCNT stands on one of two well-separated parallel plane plates, with a macroscopic field FM applied between them. For any given location "L" on the SWCNT surface, a field enhancement factor (FEF) is defined as FL/FM, where FL is a local field defined at "L". The best emission measurements from small-radii capped SWCNTs exhibit characteristic FEFs that are constant (i.e., independent of FM). This paper discusses how to retrieve this result in quantum-mechanical (as opposed to classical electrostatic) calculations. Density functional theory (DFT) is used to analyze the properties of two short, floating SWCNTs, capped at both ends, namely, a (6,6) and a (10,0) structure. Both have effectively the same height (∼5.46 nm) and radius (∼0.42 nm). It is found that apex values of local induced FEF are similar for the two SWCNTs, are independent of FM, and are similar to FEF values found from classical conductor models. It is suggested that these induced-FEF values are related to the SWCNT longitudinal system polarizabilities, which are presumed similar. The DFT calculations also generate "real", as opposed to "induced", potential-energy (PE) barriers for the two SWCNTs, for FM values from 3 V/µm to 2 V/nm. PE profiles along the SWCNT axis and along a parallel "observation line" through one of the topmost atoms are similar. At low macroscopic fields, the details of barrier shape differ for the two SWCNT types. Even for FM = 0, there are distinct PE structures present at the emitter apex (different for the two SWCNTs); this suggests the presence of structure-specific chemically induced charge transfers and related patch-field distributions.

Close Relationships between NMR J-Coupling Alternation (JCA) and Molecular Properties of Carbon Chains.

We propose a J-coupling alternation (JCA) value that is demonstrated to be a suitable parameter to evaluate the nuclear magnetic resonance (NMR) indirect spin-spin coupling constants (SSCCs) as a function of molecular properties of chains by increasing their length. As an application, we report a theoretical study of the SSCCs for the interactions between neighbor nuclei in increasingly patterned carbon chains within density functional theory. First, we examine the J-coupling constants between 1H and 13C nuclei ( n JHC) considering the separation distance, as well as between two adjacent 13C nuclei (1 JCC) considering their relative positions in polyynes and cumulenes. Further, we define and determine JCA in terms of the differences of 1 JCC, which is investigated as a function of several molecular properties, e.g., cohesive energy, characteristic vibrational frequency, average polarizability, and energy gap of the systems. We also determine JCA for other types of carbon chains, such as diphenyl-capped polyynes, polyacetylene and polythiophene. The behavior of JCA as a function of the energy gap may be related to highly π-conjugated low-band-gap carbon chains. Overall, JCA correlates very well with the electronic properties of these chains, especially with their energy gap, exhibiting positive values for pristine polyyne and polythiophene and negative values for pristine cumulene and plyacetylene. These findings indicate an alternative way to establish an appropriate SSCC descriptor that characterizes the electronic nature of the system, such as the proposed JCA value averaging the whole system, instead of using only the individual J-coupling values to give insights into the properties of large and extended systems.

Self-Healing in Carbon Nitride Evidenced As Material Inflation and Superlubric Behavior.

All known materials wear under extended mechanical contacting. Superlubricity may present solutions, but is an expressed mystery in C-based materials. We report negative wear of carbon nitride films; a wear-less condition with mechanically induced material inflation at the nanoscale and friction coefficient approaching ultralow values (0.06). Superlubricity in carbon nitride is expressed as C-N bond breaking for reduced coupling between graphitic-like sheets and eventual N2 desorption. The transforming surface layer acts as a solid lubricant, whereas the film bulk retains its high elasticity. The present findings offer new means for materials design at the atomic level, and for property optimization in wear-critical applications like magnetic reading devices or nanomachines.

Probing Lewis Acid-Base Interactions with Born-Oppenheimer Molecular Dynamics: The Electronic Absorption Spectrum of p-Nitroaniline in Supercritical CO2.

The structure and dynamics of p-nitroaniline (PNA) in supercritical CO2 (scCO2) at T = 315 K and ρ = 0.81 g cm(-3) are investigated by carrying out Born-Oppenheimer molecular dynamics, and the electronic absorption spectrum in scCO2 is determined by time dependent density functional theory. The structure of the PNA-scCO2 solution illustrates the role played by Lewis acid-base (LA-LB) interactions. In comparison with isolated PNA, the ν(N-O) symmetric and asymmetric stretching modes of PNA in scCO2 are red-shifted by -17 and -29 cm(-1), respectively. The maximum of the charge transfer (CT) absorption band of PNA in scSCO2 is at 3.9 eV, and the predicted red-shift of the π → π* electronic transition relative to the isolated gas-phase PNA molecule reproduces the experimental value of -0.35 eV. An analysis of the relationship between geometry distortions and excitation energies of PNA in scCO2 shows that the π → π* CT transition is very sensitive to changes of the N-O bond distance, strongly indicating a correlation between vibrational and electronic solvatochromism driven by LA-LB interactions. Despite the importance of LA-LB interactions to explain the solvation of PNA in scCO2, the red-shift of the CT band is mainly determined by electrostatic interactions.

A first principles approach to the electronic properties of liquid and supercritical CO2.

The electronic absorption spectra of liquid and supercritical CO2 (scCO2) are investigated by coupling a many-body energy decomposition scheme to configurations generated by Born-Oppenheimer molecular dynamics. A Frenkel exciton Hamiltonian formalism was adopted and the excitation energies were calculated with time dependent density functional theory. A red-shift of ∼ 0.2 eV relative to the gas-phase monomer is observed for the first electronic absorption maximum in liquid and scCO2. The origin of this shift, which is not very dependent on deviations from the linearity of the CO2 molecule, is mainly related to polarization effects. However, the geometry changes of the CO2 monomer induced by thermal effects and intermolecular interactions in condensed phase lead to the appearance of an average monomeric electric dipole moment〈µã€‰= 0.26 ± 0.04 D that is practically the same at liquid and supercritical conditions. The predicted average quadrupole moment for both liquid and scCO2 is〈Θ〉= - 5.5 D Å, which is increased by ∼ -0.9 D Å relative to its gas-phase value. The importance of investigating the electronic properties for a better understanding of the role played by CO2 in supercritical solvation is stressed.

Including thermal disorder of hydrogen bonding to describe the vibrational circular dichroism spectrum of zwitterionic L-alanine in water.

The vibrational circular dichroism (VCD) spectrum of l-alanine amino acid in aqueous solution in ambient conditions has been studied. The emphasis has been placed on the inclusion of the thermal disorder of the solute-solvent hydrogen bonds that characterize the aqueous solution condition. A combined and sequential use of molecular mechanics and quantum mechanics was adopted. To calculate the average VCD spectrum, the DFT B3LYP/6-311++G(d,p) level of calculation was employed, over one-hundred configurations composed of the solute plus all water molecules making hydrogen bonds with the solute. Simplified considerations including only four explicit solvent molecules and the polarizable continuum model were also made for comparison. Considering the large number of vibration frequencies with only limited experimental results a direct comparison is presented, when possible, and in addition a statistical analysis of the calculated values was performed. The results are found to be in line with the experiment, leading to the conclusion that including thermal disorder may improve the agreement of the vibrational frequencies with experimental results, but the thermal effects may be of greater value in the calculations of the rotational strengths.

Unusual electronic properties and transmission in hexagonal SiB monolayers.

After the success of graphene, several two-dimensional (2D) layers have been proposed and investigated both theoretically and experimentally in order to evaluate their structural stability and possible applications of these unusual materials in electronics. Except for graphene, only silicon and germanium were predicted to form semi-metallic honeycomb monolayers, while most of the binary graphene-like compounds are all semiconductors. These predictions have been corroborated for several 2D structures experimentally synthesized. Considering the possibility of finding other candidates in this realm, exhibiting exceptional electron mobility, we have explored low-dimensional silicon-boron compounds containing planar sp(2)-bonding silicon atoms, through first-principles density-functional theory calculations. We have demonstrated that the so-called h-SiB sheet, which is a structural analogue of 2D honeycomb binary compounds, exhibits good structural stability, compared to the structure of silicene, for example, and predicted that this structure is also able to roll up into thermally stable single-walled silicon-boron nanotubes. The h-SiB sheet exhibits a delocalized charge density like in graphene, but the partially filled π band and two highest occupied σ bands are above the Fermi level, leading to the metallic behaviour of this SiB sheet. In this sense, we perform first-principles electron transport calculations, based on the nonequilibrium Green's function formalism, which has demonstrated that h-SiB exhibits higher transmission around the Fermi energy than the transmission in graphene. Our results indicate the unusual conductivity of this new material and open up new possibilities for the realization of metallic graphene-like systems for electronic transport in low dimensions.

Origin of the Red Shift for the Lowest Singlet π → π* Charge-Transfer Absorption of p-Nitroaniline in Supercritical CO2.

The origin of the unusual solvatochromic shift of p-nitroaniline (PNA) in supercritical carbon dioxide (SCCO2) is theoretically investigated on the basis of experimental data. Ab initio quantum chemistry calculations have been employed to unveil the interaction of CO2 with this archetypical molecule. It is demonstrated that the nitro group of PNA works as an electron-donating site binding to the electron-deficient carbon atom of CO2, most probably via a Lewis acid-base interaction. Moreover, a cooperative C-H···O hydrogen bond seems to act as an additional stabilizing source during the solvation process of PNA in SCCO2. To support the influence of solute-solvent specific interactions on the lowest singlet π → π* charge-transfer excitation, we perform a sequential Monte Carlo time-dependent density functional theory simulation to evaluate the excited states of PNA in SCCO2 (T = 315 K, ρ = 0.81 g/cm(3)). A critical assessment of this simulation, compared to calculations carried out within the polarized continuum model, gives strong evidence that our proposed complexes are important in describing the solvatochromic shift of PNA in SCCO2. The calculated red shift from the gas phase accounts for 66% to 80% (depending on the degree of complexation) of the experimental data. Finally, these results also alleviate possible failures commonly attributed to long-range corrected functionals in reproducing the solvatochromism of PNA.

Exploring hydrogenation and fluorination in curved 2D carbon systems: a density functional theory study on corannulene.

Corannulene has been a useful prototype for studying C-based nanostructures as well as surface chemistry and reactivity of sp(2)-hybridized carbon-based materials. We have investigated fluorination and hydrogenation of corannulene carrying out density functional theory calculations. In general, the fluorination is energetically more favorable than hydrogenation of corannulene. The substitution of the peripheral H atoms in the corannulene molecule by F atoms leads to a larger cohesive energy gain than when F (or H) atoms are bonded to the hub carbon and bridge carbon sites of this molecule. As expected for doped C-based nanostructures, the hydrogenation or fluorination significantly changes the HOMO-LUMO gap of the system. We have obtained HOMO-LUMO gap variations of 0.13-3.46 eV for F-doped and 0.38-1.52 eV for H-doped systems. These variations strongly depend on the concentration and position of the incorporated F/H atoms, instead of the structural stability of the doped systems. Considering these calculations, we avoid practical difficulties associated with the addition/substitution reactions of larger curved two-dimensional (2D) carbon nanostructures, and we obtain a comprehensive and systematic understanding of a variety of F/H 2D doped systems.

Optical response of liquid acetonitrile at ambient conditions: the dynamical dielectric behavior from ab initio calculations.

We probe the linear optical properties of the neat liquid acetonitrile (CH(3)CN) at ambient conditions using ab initio density functional theory. Uncorrelated structures extracted from Monte Carlo simulation are employed to efficiently calculate average electronic properties. It becomes evident that condensation leads to a conduction band with a large degree of dispersion, which is consistent with the description of dipolar liquids. This allows an interpretation of the dielectric spectrum based on the electronic structure of liquid CH(3)CN, and clearly shows the influence of intermolecular interactions in the absorption features. We find that the lowest-lying excitation of the condensed phase occurs at 7.8 eV, which is reasonable as compared to the 8-9.5 eV absorption region measured in the gas phase.

Note on the free energy of transfer of fullerene C60 simulated by using classical potentials.

Diverse atomistic parameters of C60 have been developed and utilized to simulate fullerene solutions in biological environments. However, no thermodynamic assessment and validation of these parameters have been so far realized. Here, we employ extensive molecular dynamics simulations with the thermodynamic integration method in the isothermal-isobaric ensemble to investigate the transfer of a single fullerene C60 between different solvent environments using different potential models. A detailed analysis is performed on the structure and standard Gibbs free energy of transfer of C60 from benzene to ethanol. All of the interactions concerned in the transfer process are included via atomistic models. We notice that having only structural and dynamical properties is not decisive to validate reliable atomic parameters capable of describing a more realistic thermodynamic process. Thus, we employ the calculated free energy of transfer to validate more accurate atomic parameters for the solvation thermodynamics of fullerenes by direct comparison with the solubility experimental data.

Lewis acid-base interactions in weakly bound formaldehyde complexes with CO2, HCN, and FCN: considerations on the cooperative H-bonding effects.

Ab initio quantum chemistry calculations reveal that HCN and mainly FCN can form Lewis acid-base complexes with formaldehyde associated with cooperative H bonds, as first noticed by Wallen et al. (Blatchford, M. A.; Raveendran, P.; Wallen, S. L. J. Am. Chem. Soc. 2002, 124, 14818-14819) for CO2-philic materials under supercritical conditions. The present results, obtained with MP2(Full)/aug-cc-pVDZ calculations, show that the degeneracy of the nu(2) mode in free HCN or FCN is removed upon complexation in the same fashion as that of CO2. The splitting of these bands along with the electron structure analysis provides substantial evidence of the interaction of electron lone pairs of the carbonyl oxygen with the electron-deficient carbon atom of the cyanides. Also, this work investigates the role of H bonds acting as additional stabilizing interactions in the complexes by performing the energetic and geometric characterization.

Effects of hydroxyl group distribution on the reactivity, stability and optical properties of fullerenols.

We investigate the impact of hydroxyl groups on the properties of C(60)(OH)(n) systems, with n = 1, 2, 3, 4, 8, 10, 16, 18, 24, 32 and 36 by means of first-principles density functional theory calculations. A detailed analysis from the local density of states has shown that adsorbed OH groups can induce dangling bonds in specific carbon atoms around the adsorption site. This increases the tendency to form polyhydroxylated fullerenes (fullerenols). The structural stability is analyzed in terms of the calculated formation enthalpy of each species. Also, a careful examination of the electron density of states for different fullerenols shows the possibility of synthesizing single molecules with tunable optical properties.