Cooperativity effects in linear formaldehyde oligomers using density functional theory calculations

This work reports hydrogen bonding interaction in linear formaldehyde oligomers using density functional theory method. Many-body analysis technique has been used to study the various interactions in these oligomers and to obtain % contributions from individual many-body energy terms to the binding energies of these oligomers. Co-operativity effects are studied using different indicators viz. hydrogen bond strength, inter- and intramolecular distances, dissociation energy, dipole co-operativity, energy per hydrogen bond, excess energy and non-additive energy. All these indicators show strong positive hydrogen bond co-operativity in linear formaldehyde oligomers. The dipole moment changes from 2.51 D in monomer to 20.92 D in formaldehyde heptamer.     

Conformational study of erythritol and threitol in the gas state by density functional theory calculations

Density functional theory calculations using the B3LYP functional and the 6-311++G(d,p) basis set were carried out on the isolated molecules of erythritol and L-threitol. For the meso isomer, a relatively large number of conformers have to be considered to describe the gas state structure. The lowest energy conformer is characterized by the establishment of a strong intramolecular H-bond between the two terminal hydroxyl groups, giving rise to a seven-membered ring and two additional weaker H-bonds between vicinal OH groups. In the case of L-threitol, two conformers are predominant in the gas state, and both are stabilized by the formation of a cyclic system of four intramolecular hydrogen bonds involving all OH groups. The conformational stability in both diastereomers is discussed in terms of the electronic energy and of the Gibbs energy. The weighted mean enthalpy of both diastereomers in the gas state at 298.15 K was obtained from the thermodynamic data and Boltzmann populations of the low-energy conformers.     

A comparative study of the dissolubility of pure and silicon substituted hydroxyapatite from density functional theory calculations

Journal of Molecular Modeling - We have investigated the mechanism of the nucleation of acetaminophen on poly(methyl-methacrylate) and poly(vinyl-acetate) utilizing a combination of quantum...     

Electronic structure study using density functional theory in organic dendrimers

The electronic and structural properties of pyrrolic ring derivatives were studied using density functional theory (DFT) in terms of their application as organic semiconductor materials in photovoltaic devices. The B3LYP hybrid functional in combination with Pople type 6-31G(d) basis set with a polarization function was used in order to determine the optimized geometries and the electronic properties of the ground state, while transition energies and excited state properties were obtained from time-dependent (TD)-DFT with B3LYP/6-31G(d) calculation. The investigation of pyrrolic derivatives formed by the arrangement of several monomeric units revealed that three-dimensional (3D) conjugated architectures in which the combination of a triphenylamine (TPA) core with π-conjugated rings attached to the core, present the best geometric and electronic characteristics for use as an organic semiconductor material. The highest occupied molecular orbital (HOMO) − lowest unoccupied molecular orbital (LUMO) energy gap was decreased in 3D-structures that extend the absorption spectrum toward longer wavelengths, revealing a feasible intramolecular charge transfer process in these systems. All calculations in this work were performed using the Gaussian 03 W software package.     

Time-dependent density functional theory calculations of the solvatochromism of some azo sulfonamide fluorochromes

The absorption and emission spectra of three azo sulfonamide compounds in different solvents were investigated theoretically by using response functions combined with density functional theory (DFT), while the solvent effect on the structure and the electronic transitions was determined using the integral equation formalism for the polarizable continuum model (IEF-PCM). The results show that the applied different exchange-correlation functionals can reproduce the experimental values well. DFT calculations of the title compounds showed that the H-bond formed between the solute and solvent molecules is one of the major causes of the reversible solvatochromism observed in measured spectra. This is due to a better stabilization of the neutral form than the zwitterionic form in the polar protic solvents, which is characteristic of the hypsochromic shift. On the other hand, the molecules considered exhibit a monotonic behavior regarding the polarity of the low-lying excited state (Δμ_g–CT) as a function of the solvent polarity. This dependence occurs in the case of the positive solvatochromism and confirms the thesis regarding the H-bond solute–solvent interactions. Theoretically determined values of the two-photon cross section revealed that the (σ _OF ⁽²⁾ ) shows similar trends with changes in λ _abs, in contrast to 〈δ ^OF〉 values. In conclusion, the results demonstrate that the investigated molecules can be used successfully as fluorochromes in bioimaging.     

Natural polyprenylated benzophenone: keto-enol tautomerism from density functional calculations and the AIM theory

The quantum theory of atoms in molecules (QTAIM) and density functional theory (DFT) calculations were employed to investigate the structure and tautomeric equilibrium of epiclusianone, a polyisoprenylated benzophenone with interesting biological activities. Two different exchange-correlation functionals were employed, namely ωB97x-D and M06-2x, including implicit solvent models (benzene and DMSO). Our results for the thermodynamic properties show that the isomer in which the H atom is bonded to the oxygen away from the benzene ring is the most stable tautomer form of the epiclusianone, thus confirming previous charge density analysis from X-ray diffraction data (Martins et al. J Braz Chem Soc 18(8):1515–1523, 22).     

Theoretical study on the adsorption of phenol on activated carbon using density functional theory

Density functional theory (DFT) calculations performed at the PBE/DZP level using the DFT-D2 method were utilized to investigate the adsorption of phenol on pristine activated carbon (AC) and on activated carbon functionalized with OH, CHO, or COOH groups. Over the pristine AC, the phenol molecule undergoes weak physical adsorption due to van der Waals interactions between the aromatic part of the phenol and the basal planes of the AC. Among the three functional groups used to functionalize the AC, the carboxylic group was found to interact most strongly with the hydroxyl group of phenol. These results suggest that functionalized AC-COOH has great potential for use in environmental applications as an adsorbent of phenol molecules in aqueous phases.     

Modification of the phylloquinone in the A1 binding site in photosystem I studied using time-resolved FTIR difference spectroscopy and density functional theory

A phylloquinone molecule (2-methyl-3-phytyl-1,4-naphthoquinone) occupies the A1 binding site in photosystem I. Previously, we have obtained A1(-)/A1 FTIR difference spectra using labeled and unlabeled photosystem I particles and proposed assignments for many of the bands in the spectra [Sivakumar, V., Wang, R., and Hastings, G. (2005) Biochemistry 44, 1880-1893]. In particular, we suggested that a negative/positive band at 1654/1495 cm(-1) in A1(-)/A1 FTIR DS is due to a C=O/C-:O mode of the neutral/anionic phylloquinone, respectively. To test this hypothesis, we have obtained A1(-)/A1 FTIR DS for menG mutant PS I particles. In menG mutant PS I, phylloquinone in the A1 binding site is replaced with an analogue in which the methyl group at position 2 of the quinone ring is replaced with a hydrogen atom (2-phytyl-1,4-naphthoquinone). In A1(-)/A1 FTIR DS obtained using menG mutant PS I particles, we find that the 1654/1495 cm(-1) bands are upshifted by approximately 6 cm(-1). To test if such upshifts are likely for C=O/C-:O modes of neutral/anionic phylloquinone, we have used density functional theory to calculate the "anion minus neutral" infrared difference spectra for both phylloquinone and its methyl-less analogue. We have also undertaken calculations in which the C4=O carbonyl group of phylloquinone and its methyl-less analogue are hydrogen bonded (to a water or leucine molecule). We find that, irrespective of the hydrogen bonding state of the C4=O group, the C=O/C-:O modes of neutral/reduced phylloquinone are indeed expected to be upshifted by at least 6 cm(-1) upon replacement of the methyl group at position 2 with hydrogen. The calculations also suggest that certain C=C/C-:C modes of neutral/reduced phylloquinone do not shift upon replacement of the methyl group. On the basis of these calculated results, we suggest which bands in the A1(-)/A1 FTIR DS may be associated with C=C/C-:C modes of neutral/reduced phylloquinone, respectively.     

Computational characterization of sodium selenite using density functional theory

In this theoretical study we used density functional theory to calculate the molecular and crystalline structures of sodium selenite. Our structural results were compared with experimental data. From the molecular structure we determined the ionization potential, electronic affinity, and global reactivity parameters like electronegativity, hardness, softness and global electrophilic index. A significant difference in the IP and EA values was observed, and this difference was dependent on the calculation method used (employing either vertical or adiabatic energies). Thus, values obtained for the electrophilic index (2.186 eV from vertical energies and 2.188 eV from adiabatic energies) were not significantly different. Selectivity was calculated using the Fukui functions. Since the Mulliken charge study predicted a negative value, it is recommended that AIM should be used in selectivity characterization. It was evident from the selectivity index that sodium atoms are the most sensitive sites to nucleophilic attack. The results obtained in this work provide data that will aid the characterization of compounds used in crop biofortification.     

A density functional theory study of the decomposition mechanism of nitroglycerin

The detailed decomposition mechanism of nitroglycerin (NG) in the gas phase was studied by examining reaction pathways using density functional theory (DFT) and canonical variational transition state theory combined with a small-curvature tunneling correction (CVT/SCT). The mechanism of NG autocatalytic decomposition was investigated at the B3LYP/6-31G(d,p) level of theory. Five possible decomposition pathways involving NG were identified and the rate constants for the pathways at temperatures ranging from 200 to 1000 K were calculated using CVT/SCT. There was found to be a lower energy barrier to the β-H abstraction reaction than to the α-H abstraction reaction during the initial step in the autocatalytic decomposition of NG. The decomposition pathways for CHOCOCHONO₂ (a product obtained following the abstraction of three H atoms from NG by NO₂) include O–NO₂ cleavage or isomer production, meaning that the autocatalytic decomposition of NG has two reaction pathways, both of which are exothermic. The rate constants for these two reaction pathways are greater than the rate constants for the three pathways corresponding to unimolecular NG decomposition. The overall process of NG decomposition can be divided into two stages based on the NO₂ concentration, which affects the decomposition products and reactions. In the first stage, the reaction pathway corresponding to O–NO₂ cleavage is the main pathway, but the rates of the two autocatalytic decomposition pathways increase with increasing NO₂ concentration. However, when a threshold NO₂ concentration is reached, the NG decomposition process enters its second stage, with the two pathways for NG autocatalytic decomposition becoming the main and secondary reaction pathways.     

A molecular orbital selection approach for fast calculations of specific rotation with density functional theory

In this work, we describe a simple approach to select the most important molecular orbitals (MOs) to compute the optical rotation tensor through linear response (LR) Kohn-Sham density functional theory (KS-DFT). Taking advantage of the iterative nature of the algorithms commonly used to solve the LR equations, we select the MOs with contributions to the guess perturbed density that are larger than a certain threshold and solve the LR equations with the selected MOs only. We propose two criteria for the selection, and two definitions of the selection threshold. We then test the approach with two functionals (B3LYP and CAM-B3LYP) and two basis sets (aug-cc-pVDZ and aug-cc-pVTZ) on a set of 51 organic molecules with specific rotation spanning five orders of magnitude, 10⁰–10⁴ deg (dm⁻¹ (g/mL)⁻¹). We show that this approach indeed can provide very accurate values of specific rotation with estimated speedup that ranges from 2 to 8× with the most conservative selection criterion, and up to 20 to 30× with the intermediate criterion.     

Exchange interactions in a Fe5 complex: A theoretical study using density functional theory

Theoretical methods based on density functional theory have been employed to analyze the exchange interactions in an Fe₅ complex. The calculated exchange coupling constants are consistent with an S = 15/2 ground state and agree well with those reported previously for other Fe^III polynuclear complexes. The strongest antiferromagnetic interactions are those through a double oxo-benzotriazole bridging ligand, where the benzotriazole ligand seems to play a minor role and the exchange coupling constants present values typical of single oxo bridging systems.     

Syntheses, characterization, density functional theory calculations, and activity of tridentate SNS zinc pincer complexes

A series of tridentate SNS ligand precursors were metallated with ZnCl₂ to give new tridentate SNS pincer zinc complexes. The zinc complexes serve as models for the zinc active site in liver alcohol dehydrogenase (LADH) and were characterized with single crystal X-ray diffraction, ¹H, ¹³C, and HSQC NMR spectroscopies and electrospray mass spectrometry. The bond lengths and bond angles of the zinc complexes correlate well to those in horse LADH. The zinc complexes feature SNS donor atoms and pseudotetrahedral geometry about the zinc center, as is seen for liver alcohol dehydrogenase. The SNS ligand precursors were characterized with ¹H, ¹³C, and HSQC NMR spectroscopies and cyclic voltammetry, and were found to be redox active. Gaussian calculations were performed and agree quite well with the experimentally observed oxidation potential for the pincer ligand. The zinc complexes were screened for the reduction of electron poor aldehydes in the presence of a hydrogen donor, 1-benzyl-1,4-dihydronicotinamide (BNAH). The zinc complexes enhance the reduction of electron poor aldehydes. Density functional theory calculations were performed to better understand why the geometry about the zinc center is pseudo-tetrahedral rather than pseudo-square planar, which is seen for most pincer complexes. For the SNS tridentate pincer complexes, the data indicate that the pseudo-tetrahedral geometry was 43.8 kcal/mol more stable than the pseudo-square planar geometry. Density functional theory calculations were also performed on zinc complexes with monodentate ligands and the data indicate that the pseudo-tetrahedral geometry was 30.6 kcal/mol more stable than pseudo-square planar geometry. Overall, the relative stabilities of the pseudo-tetrahedral and pseudo-square planar systems are the same for this coordination environment whether the ligand set is a single tridentate SNS system or is broken into three separate units. The preference of a d¹⁰ Zn center to attain a tetrahedral local environment trumps any stabilization gained by removal of constraints within the ligand set.     

Insight into the copper coordination environment in the prion protein through density functional theory calculations of EPR parameters

Density functional theory (DFT) calculations of Cu(II) electron paramagnetic resonance (EPR) parameters for the octarepeat unit of the prion protein were conducted. Model complexes were constructed and optimized using the crystal structure of the octarepeat unit of the prion protein. Copper g and A tensors and nitrogen hyperfine and quadrupole coupling constants were calculated using DFT. Solvent effects were incorporated using the conductor-like screening model as well as through the inclusion of explicit water molecules. Calculations using the model with an additional axial water molecule added to the coordination sphere of the Cu(II) metal center give the best qualitative agreement for the copper g and A tensors. The S-band experimental EPR spectra were interpreted in light of the DFT calculations of the directly coordinated nitrogen hyperfine coupling constants which indicate that the three directly coordinated nitrogen atoms in the octarepeat unit are not equivalent. These results demonstrate that DFT calculations of EPR parameters can provide important insight with respect to the structural interpretation of experimental EPR data. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Revisiting the Holy Grail: using plant functional traits to understand ecological processes

One of ecology's grand challenges is developing general rules to explain and predict highly complex systems. Understanding and predicting ecological processes from species' traits has been considered a ‘Holy Grail’ in ecology. Plant functional traits are increasingly being used to develop mechanistic models that can predict how ecological communities will respond to abiotic and biotic perturbations and how species will affect ecosystem function and services in a rapidly changing world; however, significant challenges remain. In this review, we highlight recent work and outstanding questions in three areas: (i) selecting relevant traits; (ii) describing intraspecific trait variation and incorporating this variation into models; and (iii) scaling trait data to community‐ and ecosystem‐level processes. Over the past decade, there have been significant advances in the characterization of plant strategies based on traits and trait relationships, and the integration of traits into multivariate indices and models of community and ecosystem function. However, the utility of trait‐based approaches in ecology will benefit from efforts that demonstrate how these traits and indices influence organismal, community, and ecosystem processes across vegetation types, which may be achieved through meta‐analysis and enhancement of trait databases. Additionally, intraspecific trait variation and species interactions need to be incorporated into predictive models using tools such as Bayesian hierarchical modelling. Finally, existing models linking traits to community and ecosystem processes need to be empirically tested for their applicability to be realized.     

Hydrogen bonds in galactopyranoside and glucopyranoside: a density functional theory study

Density functional theory calculations on two glycosides, namely, n-octyl-β-D-glucopyranoside (C₈O-β-Glc) and n-octyl-β-D-galactopyranoside (C₈O-β-Gal) were performed for geometry optimization at the B3LYP/6-31G level. Both molecules are stereoisomers (epimers) differing only in the orientation of the hydroxyl group at the C4 position. Thus it is interesting to investigate electronically the effect of the direction (axial/equatorial) of the hydroxyl group at the C4 position. The structure parameters of X-H???Y intramolecular hydrogen bonds were analyzed, while the nature of these bonds and the intramolecular interactions were considered using the atoms in molecules (AIM) approach. Natural bond orbital analysis (NBO) was used to determine bond orders, charge and lone pair electrons on each atom and effective non-bonding interactions. We have also reported electronic energy and dipole moment in gas and solution phases. Further, the electronic properties such as the highest occupied molecular orbital, lowest unoccupied molecular orbital, ionization energy, electron affinity, electronic chemical potential, chemical hardness, softness and electrophilicity index, are also presented here for both C₈O-β-Glc and C₈O-β-Gal. These results show that, while C₈O-β-Glc possess– only one hydrogen bond, C₈O-β-Gal has two intramolecular hydrogen bonds, which further confirms the anomalous stability of the latter in self-assembly phenomena.