Triplet excited electronic state switching induced by hydrogen bonding: A transient absorption spectroscopy and time-dependent DFT study

The solvent plays a decisive role in the photochemistry and photophysics of aromatic ketones. Xanthone (XT) is one such aromatic ketone and its triplet-triplet (T-T) absorption spectra show intriguing solvatochromic behavior. Also, the reactivity of XT towards H-atom abstraction shows an unprecedented decrease in protic solvents relative to aprotic solvents. Therefore, a comprehensive solvatochromic analysis of the triplet-triplet absorption spectra of XT was carried out in conjunction with time dependent density functional theory using the ad hoc explicit solvent model approach. A detailed solvatochromic analysis of the T-T absorption bands of XT suggests that the hydrogen bonding interactions are different in the corresponding triplet excited states. Furthermore, the contributions of non-specific and hydrogen bonding interactions towards differential solvation of the triplet states in protic solvents were found to be of equal magnitude. The frontier molecular orbital and electron density difference analysis of the T-1 and T-2 states of XT indicates that the charge redistribution in these states leads to intermolecular hydrogen bond strengthening and weakening, respectively, relative to the S-0 state. This is further supported by the vertical excitation energy calculations of the XT-methanol supra-molecular complex. The intermolecular hydrogen bonding potential energy curves obtained for this complex in the S-0, T-1, and T-2 states support the model. In summary, we propose that the different hydrogen bonding mechanisms exhibited by the two lowest triplet excited states of XT result in a decreasing role of the n pi* triplet state, and are thus responsible for its reduced reactivity towards H-atom abstraction in protic solvents. (C) 2016 AIP Publishing LLC.

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density �� A or �� B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

Variational calculations of rovibrational states: a precise high-energy potential surface for HCN [and discussion]

We report the results of variational calculations of the rovibrational energy levels of HCN for J = 0, 1 and 2, where we reproduce all the ca. 100 observed vibrational states for all observed isotopic species, with energies up to 18000 cm$^{-1}$, to about $\pm $1 cm$^{-1}$, and the corresponding rotational constants to about $\pm $0.001 cm$^{-1}$. We use a hamiltonian expressed in internal coordinates r$_{1}$, r$_{2}$ and $\theta $, using the exact expression for the kinetic energy operator T obtained by direct transformation from the cartesian representation. The potential energy V is expressed as a polynomial expansion in the Morse coordinates y$_{i}$ for the bond stretches and the interbond angle $\theta $. The basis functions are built as products of appropriately scaled Morse functions in the bond-stretches and Legendre or associated Legendre polynomials of cos $\theta $ in the angle bend, and we evaluate matrix elements by Gauss quadrature. The hamiltonian matripx is factorized using the full rovibrational symmetry, and the basis is contracted to an optimized form; the dimensions of the final hamiltonian matrix vary from 240 $\times $ 240 to 1000 $\times $ 1000.We believe that our calculation is converged to better than 1 cm$^{-1}$ at 18 000 cm$^{-1}$. Our potential surface is expressed in terms of 31 parameters, about half of which have been refined by least squares to optimize the fit to the experimental data. The advantages and disadvantages and the future potential of calculations of this type are discussed.

C-H center dot center dot center dot O and N-H center dot center dot center dot O hydrogen bonds in liquid amides investigated by Monte Carlo simulation

Monte Carlo simulations of liquid formamide, N-methylformamide (MF), and N,N-dimethylformamide (DMF) have been performed in the isothermal and isobaric ensemble at 298 K and 1 atm, aiming to investigate the C-H ... O and N-H ... O hydrogen bonds. The interaction energy was calculated using the classical 6-12 Lennard-Jones pairwise potential plus a Coulomb term on a rigid six-site molecular model with the potential parameters being optimized in this work. Theoretical values obtained for heat of vaporization and liquid densities are in good agreement with the experimental data. The radial distribution function [RDF, g(r)] obtained compare well with R-X diffraction data available. The RDF and molecular mechanics (MM2) minimization show that the C-H ... O interaction has a significant role in the structure of the three liquids. These results are supported by ab initio calculations. This Interaction is particularly important in the structure of MF. The intensity of the N-H ... O hydrogen bond is greater in the MF than formamide. This could explain some anomalous properties verified in MF. (C) 1997 John Wiley & Sons, Inc.

C - H ��� 0 and N - H ��� 0 hydrogen bonds in liquid amides investigated by monte carlo simulation

Monte Carlo simulations of liquid formamide, N-methylformamide (MF), and N, N-dimethytformamide (DMF) have been performed in the isothermal and isobaric ensemble at 298 K and 1 atm, aiming to investigate the C-H ��� O and N-H ��� O hydrogen bonds. The interaction energy was calculated using the classical 6-12 Lennard-Jones pairwise potential plus a Coulomb term on a rigid six-site molecular model with the potential parameters being optimized in this work. Theoretical values obtained for heat of vaporization and liquid densities are in good agreement with the experimental data. The radial distribution function [RDF, g(r)] obtained compare well with R-X diffraction data available. The RDF and molecular mechanics (MM2) minimization show that the C-H ��� O interaction has a significant role in the structure of the three liquids. These results are supported by ab initio calculations. This interaction is particularly important in the structure of MF. The intensity of the N - H ��� O hydrogen bond is greater in the MF than formamide. This could explain some anomalous properties verified in MF. �� 1997 John Wiley & Sons, Inc.

Global Search for Minimum Energy (H2O)n Clusters, n = 3���5

The Gaussian-3 (G3) model chemistry method has been used to calculate the relative ��G�� values for all possible conformers of neutral clusters of water, (H2O)n, where n = 3���5. A complete 12-fold conformational search around each hydrogen bond produced 144, 1728, and 20���736 initial starting structures of the water trimer, tetramer, and pentamer. These structures were optimized with PM3, followed by HF/6-31G* optimization, and then with the G3 model chemistry. Only two trimers are present on the G3 potential energy hypersurface. We identified 5 tetramers and 10 pentamers on the potential energy and free-energy hypersurfaces at 298 K. None of these 17 structures were linear; all linear starting models folded into cyclic or three-dimensional structures. The cyclic pentamer is the most stable isomer at 298 K. On the basis of this and previous studies, we expect the cyclic tetramers and pentamers to be the most significant cyclic water clusters in the atmosphere.

Theoretical study of two states reactivity of methane activation on iron atom and iron dimer

Density functional theory (DFT) calculations have been carried out to explore the catalytic activation of C���H bonds in methane by the iron atom, Fe, and the iron dimer, Fe2. For methane activation on an Fe atom, the calculations suggest that the activation of the first C���H bond is mediated via the triplet excited-state potential energy surface (PES), with initial excitation of Fe to the triplet state being necessary for the reaction to be energetically feasible. Compared with the breaking of the first C���H bond, the cleavage of the second C���H bond is predicted to involve a significantly higher barrier, which could explain experimental observations of the HFeCH3 complex rather than CH2FeH2 in the activation of methane by an Fe atom. For methane activation on an iron dimer, the cleavage of the first C���H bond is quite facile with a barrier only 11.2, 15.8 and 8.4 kcal/mol on the septet state energy surface at the B3LYP/6-311+G(2df,2dp), BPW91/6-311+G(2df,2dp) and M06/B3LYP level, respectively. Cleavage of the second C���H bond from HFe2CH3 involves a barrier calculated respectively as 18.0, 10.7 and 12.4 kcal/mol at the three levels. The results suggest that the elimination of hydrogen from the dihydrogen complex is a rate-determining step. Overall, our results indicate that the iron dimer Fe2 has a stronger catalytic effect on the activation of methane than the iron atom.

Reaction of Aromatic Peroxyl Radicals with Alkynes: A Mass Spectrometric and Computational Study Using the Distonic Radical Ion Approach

The reaction of the aromatic distonic peroxyl radical cations N-methyl pyridinium-4-peroxyl (PyrOO center dot+) and 4-(N,N,N-trimethyl ammonium)-phenyl peroxyl (AnOO center dot+), with symmetrical dialkyl alkynes 10?ac was studied in the gas phase by mass spectrometry. PyrOO center dot+ and AnOO center dot+ were produced through reaction of the respective distonic aryl radical cations Pyr center dot+ and An center dot+ with oxygen, O2. For the reaction of Pyr center dot+ with O2 an absolute rate coefficient of k1=7.1X10-12 cm3 molecule-1 s-1 and a collision efficiency of 1.2?% was determined at 298 K. The strongly electrophilic PyrOO center dot+ reacts with 3-hexyne and 4-octyne with absolute rate coefficients of khexyne=1.5X10-10 cm3 molecule-1 s-1 and koctyne=2.8X10-10 cm3 molecule-1 s-1, respectively, at 298 K. The reaction of both PyrOO center dot+ and AnOO center dot+ proceeds by radical addition to the alkyne, whereas propargylic hydrogen abstraction was observed as a very minor pathway only in the reactions involving PyrOO center dot+. A major reaction pathway of the vinyl radicals 11 formed upon PyrOO center dot+ addition to the alkynes involves gamma-fragmentation of the peroxy O?O bond and formation of PyrO center dot+. The PyrO center dot+ is rapidly trapped by intermolecular hydrogen abstraction, presumably from a propargylic methylene group in the alkyne. The reaction of the less electrophilic AnOO center dot+ with alkynes is considerably slower and resulted in formation of AnO center dot+ as the only charged product. These findings suggest that electrophilic aromatic peroxyl radicals act as oxygen atom donors, which can be used to generate alpha-oxo carbenes 13 (or isomeric species) from alkynes in a single step. Besides gamma-fragmentation, a number of competing unimolecular dissociative reactions also occur in vinyl radicals 11. The potential energy diagrams of these reactions were explored with density functional theory and ab initio methods, which enabled identification of the chemical structures of the most important products.

Total energy shaping of a class of underactuated Port-Hamiltonian Systems using a new set of closed-loop potential shape variables

This paper proposes a method for designing set-point regulation controllers for a class of underactuated mechanical systems in Port-Hamiltonian System (PHS) form. A new set of potential shape variables in closed loop is proposed, which can replace the set of open loop shape variables-the configuration variables that appear in the kinetic energy. With this choice, the closed-loop potential energy contains free functions of the new variables. By expressing the regulation objective in terms of these new potential shape variables, the desired equilibrium can be assigned and there is freedom to reshape the potential energy to achieve performance whilst maintaining the PHS form in closed loop. This complements contemporary results in the literature, which preserve the open-loop shape variables. As a case study, we consider a robotic manipulator mounted on a flexible base and compensate for the motion of the base while positioning the end effector with respect to the ground reference. We compare the proposed control strategy with special cases that correspond to other energy shaping strategies previously proposed in the literature.

A comparative study of the energetics, structures, and mechanisms of the HCN H HNC and LiCN <-> LiNC isomerizations

The potential energy surfaces of the HCN<->HNC and LiCN<->LiNC isomerization processes were determined by ab initio theory using fully optimized triple-zeta double polarization types of basis sets. Both the MP2 corrections and the QCISD level of calculations were performed to correct for the electron correlation. Results show that electron correlation has a considerable influence on the energetics and structures. Analysis of the intramolecular bond rearrangement processes reveals that, in both cases, H (or Li+) migrates in an almost elliptic path in the plane of the molecule. In HCN<->HNC, the migrating hydrogen interacts with the in-plane pi,pi* orbitals of CN, leading to a decrease in the C-N bond order. In LiCN<->LiNC, Li+ does not interact with the corresponding pi,pi* orbitals of CN.

X-Pro peptides. A theoretical study of the hydrogen bonded conformations of (alpha-aminoisobutyryl-L -prolyl)n sequences

Intramolecularly hydrogen bonded conformations of (Aib-Pro)n sequences have been analysed theoretically. Both 4-1 (C10 and 3-1 (C7 hydrogen bonded regular structures are shown to be stereochemically feasible. Conformational energies for the helical structures have been estimated using classical potential energy methods. Both C10 and C7 conformations have very similar energies. Pyrrolidine ring puckering has a pronounced effect on the energies, and only Cv-endo puckered Pro residues can be accommodated. The theoretical calculations using spectroscopic data suggest that the recently proposed novel 310 helical conformation for benzyloxycarbonyl(Aib-Pro)4-methyl ester is in solution, is indeed energetically and stereochemically favourable.

Theoretical studies on penicillins, penicillin sulphones and their 1-oxa-1-dethia and 1-carba-1-dethia analogues: binding specificities of transeptidases and penicillinases

Empirical potential energy calculations have been carried out to determine the preferred conformations of penicillins and penicillin sulphones and their 1-oxa-1-dethia and 1-carba-1-dethia analogues. With the exception of 1-oxa-1-dethia penicillins, all the other compounds favour C2 and the C3 puckered conformations of their five-membered rings. Replacement of C2 methyl groups by hydrogen atoms as in bisnorpenicillin V or oxidation of sulphur in position 1 as in sulphones, makes the C3 puckered form much less favourable. Addition of an amino-acyl group at the C6 atom, however, makes the C3 puckered form more favoured in penicillin G or V and in 1-carba-1-dethia penicillins. Through the replacement of the sulphur atom at position 1 by an oxygen atom or by a -CH2 group increases the non-planarity of the lactam peptide bond, it significantly affects the relative disposition of the C3 carboxyl group with respect to the ��-lactam ring. These conformational differences have been correlated with the biological activities of these compounds. The present study suggests that the conformation of the bicyclic ring system may be more important for initial binding with the crosslinking enzyme(s) involved in the biosynthesis of bacterial cell-wall peptidoglycan and that the mode of binding is influenced by the nature of the side-group at the C6 atom. These studies predict, in agreement with experimental results, that the 1-oxa-1-dethia penicillin nulceus is an inhibitor of penicillianses. The study also suggests that the stereospecificities of the crosslinking enzyme(s) and penicillinases are very similar with regard to the nature of the side-group at the 6 atom and the confirmation of the bicyclic ring system. However, the confirmational requirement for the bicyclic ring system appears to be more specific in the former enzyme than in the latter.

Conformation of the LL and LD hairpin bends with internal hydrogen bonds in proteins and peptides

The conformation of three linked peptide units having an internal 4 ��� 1 type of hydrogen bond has been studied in detail, and the low energy conformations are listed. These conformations all lead to the reversal of the chain direction, and may therefore be called as ���hairpin bends��� or ���U-bends���. Since this bend can occur at the end of two chains hydrogen-bonded in the antiparallel ��-conformation, it is also known as the �����-bend���. Two types of conformation are possible when the residues at the second and third C�� atoms are both of type L (the LL bend), while only one type is possible for the LD and the DL bend. The LL bend can also accommodate the sequences LG, GL, GG (G = glycine), while the LD bend can accommodate the sequences LG, GD and GG. The conformations for the sequences DD and DL are exact inverses (or mirror images) of those for the sequences LL and LD, respectively, and have dihedral angles (phi2, ��2), (phi3, ��3) of the same magnitudes, but of opposite signs as those for the former types, which are listed, along with the characteristics (length, angle and energy) of the hydrogen bonds. A comparison of the theoretical predictions with experimental data (from X-ray diffraction and NMR studies) on proteins and peptides, show reasonably good agreement. However, a systematic trend is observable in the experimental data, slightly deviating from theory, which indicates that some deformations occur in the shapes of the peptide units forming the bend, differing from that of the standard planar peptide unit.

X-Pro peptides. A theoretical study of the hydrogen bonded conformations of (a-aminoisobutyryl-L -prolyl)n sequences

Intramolecularly hydrogen bonded conformations of (Aib-Pro)n sequences have been analysed theoretically. Both 4��������1 (C10 and 3��������1 (C7 hydrogen bonded regular structures are shown to be stereochemically feasible. Conformational energies for the helical structures have been estimated using classical potential energy methods. Both C10 and C7 conformations have very similar energies. Pyrrolidine ring puckering has a pronounced effect on the energies, and only C����-endo puckered Pro residues can be accommodated. The theoretical calculations using spectroscopic data suggest that the recently proposed novel 310 helical conformation for benzyloxycarbonyl(Aib-Pro)4-methyl ester is in solution, is indeed energetically and stereochemically favourable.

Studies On Hydrogen Bonds: Part IV. Proposed Working Criteria for Assessing Qualitative Strength of Hydrogen Bonds

The data obtained in the earlier parts of this series for the donor and acceptor end parameters of N-H. O and O-H. O hydrogen bonds have been utilised to obtain a qualitative working criterion to classify the hydrogen bonds into three categories: "very good" (VG), "moderately good" (MG) and weak (W). The general distribution curves for all the four parameters are found to be nearly of the Gaussian type. Assuming that the VG hydrogen bonds lie between 0 and �� la, MG hydrogen bonds between �� 1 and �� 2, W hydrogen bonds beyond �� 2 (where is the standard deviation), suitable cut-off limits for classifying the hydrogen bonds in the three categories have been derived. These limits are used to get VG and MG ranges for the four parameters 1 and �� (at the donor end) and �� and �� (at the acceptor end). The qualitative strength of a hydrogen bond is decided by the cumulative application of the criteria to all the four parameters. The criterion has been further applied to some practical examples in conformational studies such as ��-helix and can be used for obtaining suitable location of hydrogen atoms to form good hydrogen bonds. An empirical approach to the energy of hydrogen bonds in the three categories has also been presented.