Torsional barrier, ionization potential, and electron affinity of biphenyl—A theoretical study

Bioaccumulative chemicals are either too hard or too soft: Conceptual density functional theory as a screening tool for emerging pollutants

Conceptual density functional theory (CDFT) descriptors were computed to predict the environmental fate of approximately 6,000 widely used industrial chemicals. CDFT descriptors aligned with a molecule's possible bioaccumulation mechanism, i.e., soft chemicals are lipophilic, whereas hard chemicals may bioaccumulate by other mechanisms such as protein binding. The results have provided us with a new "rule of thumb" to guide risk assessment of chemical hazards: suspected persistent organic pollutants are either too hard (η > 0.40 hartree) or too soft (α > 200 Å3). This offers a novel approach to environmental risk assessment using two fundamental properties of a molecule.

Planar in Brooker's mode and twisted in Reichardt's mode: defying the steric forces in biphenyl types of zwitterionic systems through metameric resonance stabilizations

To be planar or to be twisted at the bridge junctions in biphenyls or biaryl types of molecular systems depends on two conflicting forces: (1) steric repulsions (destabilizations) and (2) conjugation assisted electron delocalizations (resonance stabilizations). This work reports an unfamiliar kind of behaviour shown by metamers of a zwitterionic biphenyl type of system, where the Reichardt's metamer was found to be in an usual twisted conformation (delicate balance of conflicting forces), but the Brooker's metamer was found to be in a fully planar conformation. Interestingly, at the ωB97xD/aug-cc-pVDZ level, energetically (ΔE) the planar Brooker's metamer was found to be 16.7 kcal mol^-1 lower (22.9 kcal mol^-1 lower in the CASSCF method) in energy (more stable) than the isoelectronic twisted Reichardt's metamer, and also thermodynamic ΔG values were found to be close to ΔE values for various methods (for example, 15.6 kcal mol^-1 in the above case using the ωB97xD method). When the steric repulsions are in their full potentials at the ring junction site, attainment of a conformational planarity by any biaryl type of system has not been reported previously. Without reducing the steric constraints or even without inducing any attractive forces, determining what other factors were responsible for defying the steric forces is the main focus of this investigation. Using the results of quantum mechanical computations of NBO, rotational barriers, and other saddle points (metastable conformations in singlet and triplet surfaces) in the potential energy surfaces, the dominant contribution of the resonance stabilized quinonoid form to the ground state was delineated as the possible reason for this unusual behaviour.

Imaging Charge Localization in a Conjugated Oligophenylene

Polaron formation in conjugated polymers has a major impact on their optical and electronic properties. In polyphenylene, the molecular conformation is determined by a delicate interplay between electron delocalization and steric effects. Injection of excess charges is expected to increase the degree of conjugation, leading to structural distortions of the chain. Here we investigated at the single-molecule level the role of an excess charge in an individual oligophenylene deposited on sodium chloride films. By combining sub-molecular-resolved atomic force microscopy with redox-state-selective orbital imaging, we characterize both structural and electronical changes occurring upon hole injection. While the neutral molecule exhibits a delocalized frontier orbital, for the cationic radical the excess charge is observed to localize, inducing a partial planarization of the molecule. These results provide direct evidence for self-trapping of the excess charge in oligophenylenes, shedding light on the interplay of charge localization and structural distortion.

Absolute Configuration Assignment from Optical Rotation Data by Means of Biphenyl Chiroptical Probes

A non-empirical approach for the assignment of the absolute configuration of chiral 2-alkyl-substituted carboxylic acids and primary amines by [α]_D measurements has been developed. The method requires the conversion of the chiral acids or amines into the corresponding 4,4'-disubstituted biphenylamides or biphenylazepines, respectively. In these derivatives a central-to-axial chirality transfer induces a preferred torsion in the biphenyl moiety revealed by the sign of the biphenyl A band in the ECD spectrum. By 4,4'-substitution on the biphenyl moiety a redshift of the A band is obtained, leading to an increase of its relative contribution to optical rotation. This allows to reliably establish a direct correlation between the [α]_D sign, the biphenyl twist and, then, the substrate absolute configuration. This approach thus constitutes a really practical and reliable method to assign the absolute configuration of chiral carboxylic acids and primary amines by simple and straightforward [α]_D measurement, readily obtainable by a routine instrumentation like the polarimeter.

A relative energy gradient (REG) study of the planar and perpendicular torsional energy barriers in biphenyl

Biphenyl is a prototype molecule, the study of which is important for a proper understanding of stereo-electronic effects. In the gas phase it has an equilibrium central torsion angle of ~ 45° and shows both a planar (0°) and a perpendicular (90°) torsional energy barrier. The latter is analysed for the first time. We use the newly proposed REG method, which is an exhaustive procedure that automatically ranks atomic energy contributions according to their importance in explaining the energy profile of a total system. Here, the REG method operates on energy contributions computed by the interacting quantum atoms method. This method is minimal in architecture and provides a crisp picture of well-defined and well-separated electrostatic, steric and exchange (covalent) energies at atomistic level. It is shown that the bond critical point occurring between the ortho-hydrogens in the planar geometry has been wrongly interpreted as a sign of repulsive interaction. A convenient metaphor of analysing football matches is introduced to clarify the role of a REG analysis.

Torsional Potentials of Glyoxal, Oxalyl Halides, and Their Thiocarbonyl Derivatives: Challenges for Popular Density Functional Approximations

The reliability of popular density functionals was studied for the description of torsional profiles of 36 molecules: glyoxal, oxalyl halides, and their thiocarbonyl derivatives. HF and 18 functionals of varying complexity, from local density to range-separated hybrid approximations and double-hybrid, have been considered and benchmarked against CCSD(T)-level rotational profiles. For molecules containing heavy halogens, most functionals fail to reproduce barrier heights accurately and a number of functionals introduce spurious minima. Dispersion corrections show no improvement. Calibrated torsion-corrected atom-centered potentials rectify the shortcomings of PBE and also improve on σ-hole based intermolecular binding in dimers and crystals.

Stereoelectronic Effect-Induced Conductance Switching in Aromatic Chain Single-Molecule Junctions

Biphenyl, as the elementary unit of organic functional materials, has been widely used in electronic and optoelectronic devices. However, over decades little has been fundamentally understood regarding how the intramolecular conformation of biphenyl dynamically affects its transport properties at the single-molecule level. Here, we establish the stereoelectronic effect of biphenyl on its electrical conductance based on the platform of graphene-molecule single-molecule junctions, where a specifically designed hexaphenyl aromatic chain molecule is covalently sandwiched between nanogapped graphene point contacts to create stable single-molecule junctions. Both theoretical and temperature-dependent experimental results consistently demonstrate that phenyl twisting in the aromatic chain molecule produces different microstates with different degrees of conjugation, thus leading to stochastic switching between high- and low-conductance states. These investigations offer new molecular design insights into building functional single-molecule electrical devices.

Fabrication of carbon nanomembranes by helium ion beam lithography

The irradiation-induced cross-linking of aromatic self-assembled monolayers (SAMs) is a universal method for the fabrication of ultrathin carbon nanomembranes (CNMs). Here we demonstrate the cross-linking of aromatic SAMs due to exposure to helium ions. The distinction of cross-linked from non-cross-linked regions in the SAM was facilitated by transferring the irradiated SAM to a new substrate, which allowed for an ex situ observation of the cross-linking process by helium ion microscopy (HIM). In this way, three growth regimes of cross-linked areas were identified: formation of nuclei, one-dimensional (1D) and two-dimensional (2D) growth. The evaluation of the corresponding HIM images revealed the dose-dependent coverage, i.e., the relative monolayer area, whose density of cross-links surpassed a certain threshold value, as a function of the exposure dose. A complete cross-linking of aromatic SAMs by He(+) ion irradiation requires an exposure dose of about 850 µC/cm(2), which is roughly 60 times smaller than the corresponding electron irradiation dose. Most likely, this is due to the energy distribution of secondary electrons shifted to lower energies, which results in a more efficient dissociative electron attachment (DEA) process.

Mono- and dinuclear osmium N,N'-di- and tetraphenylbipyridyls and extended bipyridyls. Synthesis, structure and electrochemistry

The efficient synthesis of mono- and dinuclear Os(IV) bipyridyl complexes is reported. These compounds show a two-step oxidation process leading to notable structural changes, which are reflected in their emission properties. During the second oxidation process a tetracation with a hydride-dihydrogen structure (instead of a trihydride) is formed. This results in a significant bathochromic shift of the emission band, accompanied by a moderate increase in intensity.

Effect of ring torsion on intramolecular vibrational redistribution dynamics of 1,1'-binaphthyl and 2,2'-binaphthyl

The role of ring torsion in the enhancement of intramolecular vibrational energy redistribution (IVR) in aromatic molecules was investigated by conducting excitation and dispersed fluorescence spectroscopy of 1,1'-binaphthyl (1,1'-BN) and 2,2'-BN. The dispersed fluorescence spectra of 1,1'-BN in the origin region of S(1)-S(0) were well resolved, which presented 25-27 cm(-1) gaps of torsional mode in the ground state. The overall profile of the dispersed spectra of 1,1'-BN is similar to that of naphthalene. In contrast, the spectra of 2,2'-BN were not resolved due to the multitude of the active torsional modes. In both cases, dissipative IVR was observed to take place with a relatively small excess vibrational energy: 237.5 cm(-1) for 1,1'-BN and 658 cm(-1) for 2,2'-BN, which clearly shows that ring torsion efficiently enhances the IVR rate. Ab initio and density functional theory calculations with medium-sized basis sets showed that the torsional potential of 1,1'-BN has a very flat minimum over the range of torsional angles from ca. 60° to 120°, whereas that of 2,2'-BN showed two well-defined potential minima at ca. 40° and 140°, in resemblance to the case of biphenyl. In this work, we propose that aromatic molecules be classified into "strong" and "weak" torsional hindrance cases: molecules with strong hindrance case show shorter torsional progressions and more effective IVR dynamics than do those with weak hindrance.

Modeling of epitaxial graphene functionalization

A new model for graphene epitaxially grown on silicon carbide is proposed. Density functional theory modeling of epitaxial graphene functionalization by hydrogen, fluorine, methyl and phenyl groups has been performed, with hydrogen and fluorine showing a high probability of cluster formation in high adatom concentration. It has also been shown that the clusterization of fluorine adatoms provides midgap states in formation, due to significant flat distortion of graphene. The functionalization of epitaxial graphene using larger species (methyl and phenyl groups) renders cluster formation impossible, due to the steric effect, and results in uniform coverage with the energy gap opening.

Assessment of quantum-chemical methods for electronic properties and geometry of signaling biomolecules

A reasonable balance between accuracy and feasibility of quantum-chemical methods depends on the complexity of the molecular system and the scientific goals. Six series of indole-, naphthalene-, phenol-, benzoic-, phenoxy-, other auxin-derivatives, and a test set of similar organic molecules have been chosen for an assessment of 13 density functional and semi-empirical molecular orbital methods with respect to electronic and structural properties. The accuracy and precision of HOMO/LUMO calculations are determined by comparison with experimental ionization potentials and electron affinities. Further comparison was performed at atomic level by covariance analysis. The methods KMLYP, MSINDO, and PM3 are precise and accurate for the whole set of molecules. The method AM1 offers comparable accuracy with the exception of electron affinities of indole derivatives, where significant deviations from experiment were observed. Geometrical properties were best reproduced with the semi-empirical method MSINDO.

Structures and thermodynamics of biphenyl dihydrodiol stereoisomers and their metabolites in the enzymatic degradation of arene xenobiotics

A key step in the metabolic degradation of biphenyl xenobiotics is catechol formation upon dehydrogenation of cis- and trans-dihydrodiols in prokaryotic and eukaryotic pathways, respectively. Structure and thermodynamics of stereoisomers of cis-, trans-2,3-biphenyl-dihydrodiols (I) and their dehydrogenation products (hydroxyketones, II), as well as final catechol (2,3-biphenyldiol, III) are studied by means of ab initio MP2/6-311++G(2df,2p)//MP2/6-311G(d,p) calculations. Formation of stereoisomers I and II is exothermic and endergonic, whereas III is enthalpically and entropically driven. Dehydrogenations are endothermic (DeltaHR0 approximately 1.5-4 kcal mol(-1)) and exergonic (DeltaGR0 approximately -5 to -7.5 kcal mol(-1)) without noticeable differences between cis and trans pathways, although the same keto stereoisomer II-(2S) is found to be the more favored product from both cis- and trans-I. The final II --> III tautomerization is thermodynamically enhanced (DeltaHR0 approximately -27, DeltaGR0 approximately -28 kcal mol(-1)) but the process is shown to have a large activation energy if it had to occur via unimolecular path. Although this tautomerization is generally assumed to be a nonenzymatic process as it involves rearomatization of an oxygenated ring, proton transfer with an anionic intermediate might be a more probable process.

Decabromobiphenyl (PBB-209) activates the aryl hydrocarbon receptor while decachlorobiphenyl (PCB-209) is inactive: experimental evidence and computational rationalization of the different behavior of some halogenated biphenyls

In rat H4IIE cells permanently transfected with a luciferase gene under the control of AhR, incubation with PBB-209 led to a statistically significant increase of luminescence. In this system, PCB-209 only caused a small induction of luciferase activity. In a fish cell line, only PBB-209 was able to provoke an induction of ethoxyresorufin- O-deethylase activity. Ligand binding to the AhR was studied by means of a cell-free in vitro system in which the activation of AhR is very unlikely to occur without ligand binding. None of the biphenyls studied provoked any activation of AhR in this system. To rationalize the results and to get insight into the molecular mechanism of activation of AhR by PBB-209 as compared with PCB-209, a comprehensive computational study was carried out on these congeners as well as on PCB-126 and PCB-169, two potent AhR activators through ligand binding. The calculations include (i) conformational analysis and dipole moments of each conformer, (ii) aromaticity indices, (iii) molecular electrostatic potentials, (iv) quadrupole moments, (v) electronic and reactivity descriptors, and (vi) dissociation energies of C-Cl and C-Br bonds in model aromatic compounds. It was found that some molecular features of PBB-209, such as the electrostatic potential (EP) and EP-derived descriptors (Politzer's parameters), indicate that PBB-209 is more similar to PCB-126 and PCB-169 than to PCB-209, which share quite similar geometries based on the substitution pattern. The similarity between PBB-209, PCB-126, and PCB-169 seems to hint that these three compounds can share, at least partially, similar mechanisms of activation of AhR. It is unquestionable that PCB-126 and PCB-169 directly bind AhR and PBB-209 does not. We hypothesize that there are several simultaneous mechanisms for activation of AhR, and the most active compounds act for more than one mechanism.

Dependence of Optical Properties of Oligo-para-phenylenes on Torsional Modes and Chain Length

A systematic characterization of excited-state properties of para-phenylene oligomers constructed from two to eight aromatic rings is presented using density functional theory (DFT) and the coupled-cluster singles and doubles (CC2) method. Geometry optimizations have been performed for the ground state and for the electronically excited state. Vertical excitations and the fluorescence transitions have been calculated. Time-dependent DFT (TDDFT) method underestimates excitation and fluorescence energies systematically in comparison with experimental results. The computed TDDFT lifetime for the polymer limit (0.43 ns) is in agreement with the experimental value of 0.55 ns. The TDDFT torsional potential curves were investigated for biphenyl, terphenyl, and quarterphenyl oligomers in their electronic ground and excited states. Our calculations show an increase in the separation of the lowest excited state (S1) to the next higher one with increasing molecular size. No indication is found for state crossings of the S1 state with higher ones from planar structures up to torsional angles of 60 degrees to 70 degrees. Thus, an adiabatic description of the dynamics of the S1 state might significantly simplify any dynamics simulations of torsional broadenings.

Assessing a new nonempirical density functional: Difficulties in treating π-conjugation effects

The reliability of the Tao-Perdew-Staroverov-Scuseria (TPSS) exchange-correlation functional for the description of conjugation effects in model pi-conjugated systems has been thoroughly assessed through the calculation of torsion energy profiles. The functional reproduces qualitatively the shape of torsional potentials but, interestingly, the mixing of TPSS and exact exchange governs the quantitative results: thus, well-defined hybrid extensions of the functional are consistently employed to improve the results. The hybrid approaches led to more accurate descriptions of conjugation effects but, however, the finest performance along the whole range of dihedral angles was obtained by a customized mixing of pure or hybrid TPSS functionals and wave function methods in a multicoefficient fashion. Despite the successful construction of this nonempirical functional, higher rungs of the ladder of methods in which TPSS is based are hoped to reduce the errors with respect to reference data for pi-conjugated systems.

Inner Shell Excitation Spectroscopy of Biphenyl and Substituted Biphenyls: Probing Ring−Ring Delocalization

Quantitative optical oscillator strength spectra for C 1s excitation and ionization of gas-phase biphenyl, decafluorobiphenyl, and 2,2'-bis(bromomethyl)-1,1'-biphenyl have been derived from electron energy loss spectroscopy recorded under electric dipole dominated conditions. The C 1s X-ray absorption spectrum of hexaphenylbenzene has been recorded in the solid state. The C 1s spectral features are interpreted with the aid of ab initio calculations for core excitation of benzene, biphenyl, hexafluorobenzene, and decafluorobiphenyl. A weak feature at 287.7 eV in biphenyl is identified as a C 1s --> pi(deloc) transition, characteristic of ring-ring delocalization. Its intensity and position are shown to be related to the average torsion angle and thus the extent of pi-pi-interaction between adjacent aromatic rings. The effects of perfluoro substitution on core excitation spectra are also characterized and discussed.

The Fate of Dicationic States in Molecular Clusters of Benzene and Related Compounds

Calculations employing density functional theory indicate that, rather than undergoing fragmentation, dicationic clusters of benzene, hexafluorobenzene, and naphthalene produced by sequential one-electron or sudden double-ionization experiments on the neutrals can relax via the formation of inter-ring covalent C-C bonds, along with a series of proton transfers that enable a substantial reduction of inter- and intramolecular Coulomb repulsions. The theoretically predicted chemically bound structures correspond to deep local energy minima on the potential energy surface pertaining to the lowest electronic state of the dications and can therefore be regarded as metastable (kinetically long-lived) species. This discovery invalidates on theoretical grounds the liquid-droplet model of multiply charged clusters and sheds very unexpected light on possible consequences in chemistry of the intermolecular Coulombic decay (ICD) mechanism [Cederbaum, L. S.; et al. Phys. Rev. Lett. 1997, 79, 4778; Jahnke, T.; et al. Phys. Rev. Lett. 2004, 93, 163401] for deep inner-valence ionized states. Propagation of charge rearrangement reactions and proton transfers to several monomers may eventually lead to the formation of rather extended dicationic assemblies.

Probing the Shape and Stereochemistry of Molecular Orbitals in Locally Flexible Aromatic Chains: A Penning Ionization Electron Spectroscopy and Green's Function Study of the Electronic Structure of Biphenyl

We report on the results of an exhaustive study of the interplay between the valence electronic structure, the topology and reactivity of orbitals, and the molecular structure of biphenyl by means of Penning ionization electron spectroscopy in the gas phase upon collision with metastable He*(2(3)S) atoms. The measurements are compared with one-particle Green's function calculations of one-electron and shake-up valence ionization spectra employing the third-order algebraic diagrammatic construction scheme [ADC(3)]. Penning ionization intensities are also analyzed by means of the exterior electron-density model and comparison with photoelectron spectra: in contrast with the lines originating from sigma orbitals, ionization lines belonging to the pi-band system have large Penning ionization cross sections due to their greater extent outside the molecular van der Waals surface. The involved chemi-ionization processes are further experimentally investigated using collision-energy-resolved Penning ionization electron spectroscopy. The cross sections of pi-ionization bands exhibit a markedly negative collision-energy dependence and indicate that the interaction potential that prevails between the molecule and the He*(2(3)S) atom is strongly attractive in the pi-orbital region. On the other hand, the partial ionization cross sections pertaining to sigma-ionization channels are characterized by more limited collision-energy dependencies, as a consequence of rather repulsive interactions within the sigma-orbital region. A comparison of ADC(3) simulations with the Penning ionization electron spectra and UV photoelectron spectra measured by Kubota et al. [Chem. Phys. Lett. 1980, 74, 409] on thin films of biphenyl deposited at 170 and 109 K on copper demonstrates that biphenyl molecules lying at the surface of polycrystalline layers adopt predominantly a planar configuration, whereas within an amorphous sample most molecules have twisted structures similar to those prevailing in the gas phase.

Spectral broadening and diffusion by torsional motion in biphenyl

We have studied biphenyl by time-dependent density-functional theory. In particular, we have analyzed the dependence of singlet excitation energies and transition dipoles on the torsional angle between the phenyl groups. The torsional spectrum has been computed quantum mechanically as well as semiclassically in order to understand how this influences the broadening of absorption and luminescence spectra. Our results are in best agreement with supersonic jet spectroscopy data, but also fit astonishingly well to spectra of biphenyl in condensed phase. Furthermore, we compare the torsional and vibrational relaxation and discuss qualitatively the general consequences for poly-para-phenylenes and related conjugated polymers as poly-thiophenes, considering, in particular, how side chains and solvents may affect the optical spectra.