Orbital based electronic structural signatures of the guanine keto G-7H/G-9H tautomer pair as studied using dual space analysis

Electronic structural signatures of the guanine-7H and guanine-9H tautomers have been investigated on an orbital by orbital basis using dual space analysis. A combination of density functional theory (B3LYP/TZVP), the statistical average of model orbital potentials (SAOP/TZ2P) method and outer valence Green's function theory (OVGF/TZVP) has been used to generate optimal tautomer geometries and accurate ionization energy spectra for the guanine tautomer pair. The present work found that the non-planar form for both of the guanine keto pair possesses lower energies than their corresponding planar counterparts, and that the canonical form of the guanine-7H tautomer has slightly lower total energy than guanine-9H. This latter result is in agreement with previous experimental and theoretical findings. In the planar guanine pair the geometric parameters and anisotropic molecular properties are compared, focusing on changes caused by the mobile proton transfer. It is demonstrated that the mobile proton only causes limited disturbance to isotropic properties, such as geometry and the energetics, of the guanine keto tautomer pair. The exception to this general statement is for related local changes such as the N((7))-C((8)) and C((8))-N((9)) bond length resonance between the single and double bonds, reflecting the nitrogen atom being bonded with the mobile proton in the tautomers. The mobile proton distorts the electron distribution of the tautomers, which leads to significant changes in the molecular anisotropic properties. The dipole moment of guanine-7H is altered by about a factor of three, from 2.23 to 7.05 D (guanine-9H), and the molecular electrostatic potentials also reflect significant electron charge distortion. The outer valence orbital momentum distributions, which were obtained using the plane wave impulse approximation (PWIA), have demonstrated quantitatively that the outer valence orbitals of the tautomer pair can be divided into three groups. That is orbitals 1a'-7a' and 18a', which do not have visible alternations in the tautomeric process (which consist of either pi orbitals or are close to the inner valence shell); a second group comprising orbitals 19a'-22a', 25a', 26a', 28a', 29a' and 31a', which show small perturbations as a result of the mobile hydrogen locations; and group three, orbitals 23a', 24a', 27a', 30a' and 32a', which demonstrate significant changes due to the mobile proton transfer and are therefore considered as signature orbitals of the G-7H/G-9H keto tautomeric process.     

Ab initio study of electron transport in 4-(3-nitro-4-tetrafluorophenylthiolate-ethynyl, phenylethynyl) benzenethiolate

The electron transport of the 4-(3-nitro-4-tetrafluorophenylthiolate-ethynyl, phenylethynyl) benzenethiolate (S-FNPPB-o) molecule assembled in two Au (111) electrodes, was studied using two approaches: in the first approximate approach an electric field was applied to the pure molecule attached to two thiolate ends fixed, and in the second approach we used the nonequilibrium Green′s function formalism (NEGF) coupled to DFT to calculate the I-V curve and the voltage dependence of the transmission function in the extended system, molecule plus electrodes. By applying an electric field to the pure molecule plus thiolate ends fixed, and visualizing the changes in the spatial distribution of the frontier molecular orbitals, we can expect based on the continuity of the conduction pathway in electron transport, that if electron transport occurs through the frontier orbitals, only the LUMO orbital would create an open channel for electron transport due to its delocalized nature and large orbital density at the thiolate groups. The NEGF calculations indicate that at applied voltages lower than ±0.8 V, the current is related to transmission values through the tails of the broad LUMO orbital, and since this orbital is the one closer to the Fermi energy, and we observed very low current values in this region, higher current values at positive bias than at negative bias. As the voltage exceeds ±0.8 V the current increases from the contribution of more states from the broadened part of the transmission function from the LUMO orbital, and when the voltage approaches ±2 V, the LUMO + 1 orbital enters into the bias window and the current increases again.     

Relativistically parameterized extended Hückel calculations. 10. Lanthanide trihalides

The REX relativistically parameterized extended Hückel method is used to study the electronic structure of lanthanide trihalide molecules. All valence orbitals are described in terms of double-zeta Slater functions, with the atomic orbital parameters being determined by a least-squares fitting to published relativistic (Dirac- Fock) radial densities. Comparisons of orbital energies to experimental values are made and various trends are discussed. Ab initio all-electron calculations at the self-consistent field level and as a function of molecular geometry are reported for LaH₃, LaF₃, and LaCl₃. While LaH₃ and LaF₃ are calculated to be pyramidal, LaCl₃ is calculated to be planar.     

UV photoelectron studies of biological pyrimidines: the electronic structure of uracil

The photoelectron spectrum of uracil has been measured in the gas phase. A qualitative interpretation of the spectrum indicates that the first and third highest occupied orbitals in uracil are π orbitals, while the second and fourth highest orbitals are lone-pair orbitals associated with oxygen atoms. This ordering is consistent with that predicted by CNDO and INDO molecular orbital calculations.     

Electron Transport and Nonlinear Optical Properties of Substituted Aryldimesityl Boranes: A DFT Study

A comprehensive theoretical study was carried out on a series of aryldimesityl borane (DMB) derivatives using Density Functional theory. Optimized geometries and electronic parameters like electron affinity, reorganization energy, frontiers molecular contours, polarizability and hyperpolarizability have been calculated by employing B3PW91/6-311++G (d, p) level of theory. Our results show that the Hammett function and geometrical parameters correlates well with the reorganization energies and hyperpolarizability for the series of DMB derivatives studied in this work. The orbital energy study reveals that the electron releasing substituents increase the LUMO energies and electron withdrawing substituents decrease the LUMO energies, reflecting the electron transport character of aryldimesityl borane derivatives. From frontier molecular orbitals diagram it is evident that mesityl rings act as the donor, while the phenylene and Boron atom appear as acceptors in these systems. The calculated hyperpolarizability of secondary amine derivative of DMB is 40 times higher than DMB (1). The electronic excitation contributions to the hyperpolarizability studied by using TDDFT calculation shows that hyperpolarizability correlates well with dipole moment in ground and excited state and excitation energy in terms of the two-level model. Thus the results of these calculations can be helpful in designing the DMB derivatives for efficient electron transport and nonlinear optical material by appropriate substitution with electron releasing or withdrawing substituents on phenyl ring of DMB system.     

Polyhedral approximation approach to molecular orbital graphics

Interactive molecular orbital display and manipulation system is developed under the concept of polyhedral approximation of equivalued surfaces of orbital and density functions. Two new methods for generating polyhedral data, lscell-based’ and ‘propagation’ methods, are presented. The use of stored polyhedral data largely simplifies the shaded-image synthesis, direct manipulation, and animation display of molecular orbitals. Manual overlay of frontier orbitais is useful to estimate the geometries of chemically reacting molecules.     

In-silico analysis of substituent effect on the static first order hyperpolarizability of electron donating mono substituted Chalcone derivatives

Noncentrosymmetric π conjugated systems with suitable electron donor acceptor groups play a crucial role in material NLO activity. The influence of an electron donating mono substituent at the para position of the phenylene ring of chalcone was investigated as a resource for second harmonic generation. The geometrical optimization of 11 electron donating group substituted chalcones were performed using density functional theory at the B3LYP/6-311G(d,p) level and compared with experimental geometrical parameters of five reported chalcones. All the derivatives are transparent to visible radiation as shown by the electronic absorption spectra investigated by the TDDFT-CAM B3LYP/6-311G(d,p) method, and the maximum absorption wavelength was due to the π_PhB?→?π* transition. The first order hyperpolarizability β_tot, calculated using the CAM B3LYP/6-311G(d,p) method, increases with the electron donating ability of the substituent, and the largest β_tot was observed for dimethylamino substituent. The Hammett substituent constant (σ_p) shows good linear correlation with β, λ_max, and E_gap in the ground state. The Brown constant (σ_p⁺) was better correlated indicating the polarization of carbonyl group in the excited state. Frontier molecular orbitals also reveal the valence electron excitation. Correlation of σ_p with various parameters was analyzed to assess the property interrelationship with electronic reorganization in the molecule. The electronic structures of molecular fragments were described in terms of natural bond orbital analysis, which shows intramolecular interactions.     

An insight into the general relationship between the three dimensional structures of enzymes and their electronic wave functions: Implication for the prediction of functional sites of enzymes

In this study, we explored the general relationship between the three-dimensional (3D) structures of enzymes and their electronic wave functions. Furthermore, we developed a method for the prediction of their functionally important sites. For this purpose, we first performed linear-scaling molecular orbital calculations for 112 nonredundant, non-homologous enzymes with known structure and function. In consequence, we showed that the canonical molecular orbitals (MOs) of the enzymes could be classified into three groups according to the degree of electron delocalization: highly localized orbitals (Group A), highly delocalized orbitals whose electrons are distributed over almost the whole molecule (Group B), and moderately delocalized orbitals (Group C). The MOs belonging to Group A are located near the HOMO-LUMO band gap, and thereby include the frontier orbitals of a given enzyme. We inferred that the MOs of Group B play a role in stabilizing the 3D structure of the enzyme, while those of Group C contribute to constructing the covalent bond framework of the enzyme. Next, we investigated whether the frontier orbitals of enzymes could be used for identifying their potential functional sites. As a result, we found that the frontier orbitals of the 112 enzymes have a high propensity to be colocalized with the known functional sites, especially when the enzymes are hydrated. Such a propensity is shown to be remarkable when Glu or Asp is a functional site residue. On the basis of these results, we finally propose a protocol for the prediction of functional sites of enzymes.     

Spatial electron distribution and population analysis of amides, carboxylic acid, and peptides, and their relation to empirical potential functions.

The properties of the electron distribution in amides, peptides, and carboxylic acids, obtained from ab-initio molecular orbital calculations using both minimal and extended basis sets have been studied. These properties are discussed in terms of some of the common assumptions made in empirical conformational calculations of biomolecules. In particular, population analyses of 15 compounds in these families were carried out with both the minimal and extended basis sets, and compared with results of CNDO/2 calculations. It is suggested that population analysis is a useful tool for recognizing patterns of charge distributions, and investigating the transferability of parameters of different functional groups. However, its use for providing partial charges for conformational analysis is a questionable procedure. A more detailed analysis of the charge distribution was carried out by calculating the spatial electron distribution in the four compounds, N-methylacetamide, acetic acid, diketopiperazine, and N-acetyl-N′-methylalanine. Both total electron-density maps and differencedensity maps are presented. The properties of the overall shape of the molecule and the atoms in the molecule, are discussed in terms of the former along with three-dimensional shape plots of the total density. The distortion accompanying molecular formation, resulting in such features as the lone pair orbital and “bonding deensities” is discussed in terms of the difference maps. Semiquantitative estimates of the bonding and orbital densities resulting from the integration of the densities are also presented. Finally, one of the novel features of the study is the presentation of three-dimensional surfaces of constant difference densities from which the shapes of the orbitals and bonding densities emerge.     

DFT and the electromerism in complexes of iron with diatomic ligands

A reliable procedure is proposed for assigning the electronic structures for large biologically-relevant systems, where the size of the model confines one to the use of density functional theory (DFT) methods, and where the risk of over-interpreting DFT-derived molecular orbitals and spin densities still exists. The proposed approach focuses on the use of the only DFT-derived parameter that is unanimously recognized to be reliable: the geometry. We examine DFT-derived O-O bond lengths in formally ferrous-dioxygen models, and compare them to bond lengths in free, non metal-bound, dioxygen, superoxide and peroxide moieties. Likewise, we compare the N-O bond lengths within ferrous-nitrosyl {FeNO}7 models, with the same parameter in free NO+, NO*, and HNO species. This allows a calibrated, straightforward way of assigning the electronic structure in systems where electromerism makes detailed single-reference molecular orbital analysis unreliable.     

Theoretical investigation of auxiliary electronic acceptors in modifying D-D-π-A sensitizers for dye-sensitized solar cells

In light of the performance of the SD2 pigments in DSSC, in order to expand the absorption spectral scope, decrease the energy difference between the highest occupied and the lowest unoccupied molecular orbitals, with SD2 dye molecular electron donor and electron acceptor as the fundamental framework, the indole fragment and thiophene derivative in the prototype dye molecule were replaced by the two π-bridges (labeled PA, PB, respectively) and the four auxiliary electron acceptors (labeled A1, A2, A3, A4, respectively). For the sake of characterizing dye molecules as thoroughly as possible in DSSC, the frontier orbital energy levels, ultraviolet absorption spectra, natural bond orbital analysis, intramolecular charge transfer, charge and hole reorganization energies, parameters influencing the short-circuit current density and the open-circuit photovoltage for these eight individual dye molecules are carried out to try to fully characterize the properties of these dye molecules. According to these computational results of physical quantities and based on the performance of these dye molecules in the above aspects, in this paper, six free molecular models were picked out to combine with titanium dioxide cluster to calculate their geometrical structures, frontier orbital distributions, electron excitation energies, ultraviolet absorption spectra and the composition of the electronic transitions in chloroform solvent with polarizable continuum model. The results of these calculations show that the PA-A2 and PB-A4 dye molecule has better properties in electron transfer and spectral absorption range before and after the adsorption on the titanium dioxide.     

Quantum Chemical Studies on Proteins in the Reaction Center of Rhodobacter sphaeroides

The electronic structure of protein chains L and M in photosynthetic reaction center (PRC) of Rhodobacter sphaeroides (Van Niel) Imhoff, Truper et Pfennig) was studied by using the Overlapping Dimer Approximation method and the Extended Negative Factor Counter method at ab initio level. The result indicated that: (1) Amino acid residues, the molecular orbitals of which composed the main components of frontier orbitals of protein chain L (M ), are located at the random coil areas of chain L (α helix areas of chain M ). Since the random coil is flexible and more easy to change its conformation in the electron transfer process and to reduce the energy of the system, and the structure of the α helix is reletively stable, this difference might be one of the causes for the electron transfer in photosynthetic reaction center (PRC) only takes place along the L branch. (2) The His residues which axially coordinated to the “special pair” P and accessory chlorophyll molecules （ABChls） are essentially important for the ELUMO levels of P and ABChl. But, the corresponding molecular orbitals of these His residues do not appear in the composition of frontier orbitals of protein chains. It means that the interaction between pigment molecules and protein chains do not influence the contribution to the frontier orbitals of protein chains explicitly, but influences the corresponding ELUMO levels significantly.     

Electronic structure and stabilities of Ni-doped germanium nanoclusters: a density functional modeling study

The present study reports the geometry, electronic structure, growth behavior and stability of neutral and ionized nickel encapsulated germanium clusters containing 1–20 germanium atoms within the framework of a linear combination of atomic orbital density functional theory (DFT) under a spin polarized generalized gradient approximation. In the growth pattern, Ni-capped Ge_n and Ni-encapsulated Ge_n clusters appear mostly as theoretical ground state at a particular size. To explain the relative stability of the ground state clusters, variation of different parameters, such as average binding energy per atom (BE), embedding energy (EE) and fragmentation energy (FE) of the clusters, were studied together with the size of the cluster. To explain the chemical stability of the clusters, different parameters, e.g., energy gap between the highest occupied and lowest unoccupied molecular orbitals (HOMO–LUMO gap), ionization energy (IP), electron affinity (EA), chemical potential (μ), chemical hardness (η), and polarizability etc. were calculated and are discussed. Finally, natural bond orbital (NBO) analysis was applied to understand the electron counting rule applied in the most stable Ge₁₀Ni cluster. The importance of the calculated results in the design of Ge-based superatoms is discussed.

Structure and electronic properties of Alq3 derivatives with electron acceptor/donor groups at the C4 positions of the quinolate ligands: a theoretical study

The molecular structures of the ground (S(0)) and first singlet excited (S(1)) states of Alq3 derivatives in which pyrazolyl and 3-methylpyrazolyl groups are substituted at the C4 positions of the 8-hydroxyquinolate ligands as electron acceptors, and piperidinyl and N-methylpiperazinyl groups are substituted at the same positions as electron donors, have been optimized using the B3LYP/6-31G and CIS/6-31G methods, respectively. In order to analyze the electronic transitions in these derivatives, the frontier molecular orbital characteristics were analyzed systematically, and it was found that the highest occupied molecular orbital is localized on the A ligand while the lowest unoccupied molecular orbital is localized on the B ligand in their ground states, similar to what is seen for mer-Alq3. The absorption and emission spectra were evaluated at the TD-PBE0/6-31G level, and it was observed that electron acceptor substitution causes a red-shift in the emission spectra, which is also seen experimentally. The reorganization energies were calculated at the B3LYP/6-31G level and the results show that acceptor/donor substitution has a significant effect on the intrinsic charge mobilities of these derivatives as compared to mer-Alq3.     

A molecular orbital study on the oxidation of hydrogen donor molecules by peroxidase compound II

The electronic characteristics of some hydrogen donor substrate (phenol and aniline derivatives) for peroxidase reaction were calculated with the aid of the CNDO/2 and other methods. These results were compared with the experimental data concerning the rate of oxidation of these compounds by peroxidase and hydrogen peroxide. No simple relationship between the total or frontier electron densities on the oxygen or nitrogen atoms, or the lowest unoccupied orbital energies, and the rate of oxidation was found. It was, however, found that the logarithm of the rate of oxidation for the compounds studied correlates linearly with the highest occupied orbital energies. Based on these results, the mechanism of electron transfer from the substrate to compound II is discussed.     

Chlorophyll ring deformation modulates Qy electronic energy in chlorophyll-protein complexes and generates spectral forms

The possibility that the chlorophyll (chl) ring distortions observed in the crystal structures of chl-protein complexes are involved in the transition energy modulation, giving rise to the spectral forms, is investigated. The out-of-plane chl-macrocycle distortions are described using an orthonormal set of deformations, defined by the displacements along the six lowest-frequency, out-of-plane normal coordinates. The total chl-ring deformation is the linear combination of these six deformations. The two higher occupied and the two lower unoccupied chl molecular orbitals, which define the Q(y) electronic transition, have the same symmetry as four of the six out-of-plane lowest frequency modes. We assume that a deformation along the normal-coordinate having the same symmetry as a given molecular orbital will perturb that orbital and modify its energy. The changes in the chl Q(y) transition energies are evaluated in the Peridinin-Chl-Protein complex and in light harvesting complex II (LHCII), using crystallographic data. The macrocycle deformations induce a distribution of the chl Q(y) electronic energy transitions which, for LHCII, is broader for chla than for chlb. This provides the physical mechanism to explain the long-held view that the chla spectral forms in LHCII are both more numerous and cover a wider energy range than those of chlb.     

Electronic properties of neuroleptics: ionization energies of benzodiazepines

Vertical ionization energies (VIEs) of medazepam, nordazepam and their molecular subunits have been calculated using the electron propagator method in the P3/CEP-31G* approximation. Vertical electron affinities (VEAs) have been obtained with a ∆SCF procedure at the DFT-B3LYP/6-31+G* level of theory. Excellent correlations have been achieved between IE_calc and IE_exp, allowing reliable assignment of the ionization processes. Our proposed assignment differs in many instances from that previously reported in the literature. The electronic structure of the frontier Dyson orbitals shows that the IE and EA values of the benzodiazepines can be modulated by substitution at the benzene rings. Hardness values, evaluated as (IE − EA)/2, follow the trend of the experimental singlet transition energies. Medazepam is a less hard (i.e., less stable) compound than nordazepam.     

Equivalent potential of water for the electronic structure of glycine

First-principles, all-electron, ab initio calculations have been performed to construct an equivalent potential of water for the electronic structure of glycine (Gly) in solution. The calculation involved three steps. The first step was to search for the minimum-energy geometric structure of the Gly + nH₂O system. The second step was to calculate the electronic structure of Gly with the potential of water molecules via the self-consistent cluster-embedding method (SCCE), based on the result obtained in the first step. The last step was to calculate the electronic structure of Gly with the potential of dipoles after replacing the water molecules with dipoles. The results show that the occupied molecular orbitals of Gly are raised by about 0.0524 Ry on average due to the effect of water. The effect of water can be simulated well using the dipole potential. The equivalent potential obtained can be applied directly to electronic structure calculations of proteins in solution using the SCCE method.