Molecular structures and properties of the complete series of bromophenols: Density functional theory calculations

The complete series of 19 bromophenols have been studied by density functional theory (DFT) calculations at the B3LYP/6-311G++(d,p) level. The molecular structures and properties of bromophenols are strongly influenced by intramolecular hydrogen bonding of ortho-bromine, steric and inductive effects of substituted bromine, and other intramolecular electrostatic interactions. Systematic trends in several structural parameters and molecular properties of bromophenols have been found with the increasing number of bromine substitutions, including increase in O-H bond length, decrease in C-O bond length, red shift in O-H stretching frequency, and blue shift in O-H torsional frequency. Correlations among several key molecular parameters as well as those with available aqueous pKa values are examined. Comparisons with chlorophenols have indicated that the inductive effect of substituted bromine appears larger and bromophenols are slightly stronger acids than chlorophenols.     

Structural and energetic consequences of the formation of intramolecular hydrogen bonds

Formation of intramolecular hydrogen bonds leads to structural modifications in the whole molecule, which are discussed on the basis of B3LYP/6-31G(d,p) calculations. The energy and the structure of various hydrogen-bonded and open conformers are considered for two groups of ortho-substituted phenols–N-dimethylaminomethylphenols (Mannich bases) and N-methylbenzylideneamines (Schiff bases). The energy of intramolecular hydrogen bond formation in Mannich bases was corrected for non-bonded interactions within the molecules, based on a thermodynamic cycle. Structural data were used to estimate the fraction of the ortho-quinoid (keto) form in particular tautomers. It is shown that proton transfer in Schiff bases leads to an increase of this fraction to about 40%, while opening of the hydrogen bond in the proton transferred form increases the keto fraction to 70%.     

Revealing Intra- and Intermolecular Interactions Determining Physico-Chemical Features of Selected Quinolone Carboxylic Acid Derivatives

The intra- and intermolecular interactions of selected quinolone carboxylic acid derivatives were studied in monomers, dimers and crystals. The investigated compounds are well-recognized as medicines or as bases for further studies in drug design. We employed density functional theory (DFT) in its classical formulation to develop gas-phase and solvent reaction field (PCM) models describing geometric, energetic and electronic structure parameters for monomers and dimers. The electronic structure was investigated based on the atoms in molecules (AIM) and natural bond orbital (NBO) theories. Special attention was devoted to the intramolecular hydrogen bonds (HB) present in the investigated compounds. The characterization of energy components was performed using symmetry-adapted perturbation theory (SAPT). Finally, the time-evolution methods of Car–Parrinello molecular dynamics (CPMD) and path integral molecular dynamics (PIMD) were employed to describe the hydrogen bond dynamics as well as the spectroscopic signatures. The vibrational features of the O-H stretching were studied using Fourier transformation of the autocorrelation function of atomic velocity. The inclusion of quantum nuclear effects provided an accurate depiction of the bridged proton delocalization. The CPMD and PIMD simulations were carried out in the gas and crystalline phases. It was found that the polar environment enhances the strength of the intramolecular hydrogen bonds. The SAPT analysis revealed that the dispersive forces are decisive factors in the intermolecular interactions. In the electronic ground state, the proton-transfer phenomena are not favourable. The CPMD results showed generally that the bridged proton is localized at the donor side, with possible proton-sharing events in the solid-phase simulation of stronger hydrogen bridges. However, the PIMD enabled the quantitative estimation of the quantum effects inclusion—the proton position was moved towards the bridge midpoint, but no qualitative changes were detected. It was found that the interatomic distance between the donor and acceptor atoms was shortened and that the bridged proton was strongly delocalized.     

Spectral investigation of the intramolecular charge-transfer in some aminotriazole Schiff bases

3-Amino-1,2,4-triazole Schiff bases were reported to contain intramolecular charge-transfer. The enhancing and depressing effects were remarkable as the substituent was changed from electron-donating to electron-withdrawing groups. The path of the resonating delocalization was reversed in the case of the p-NO2 group. To validate these results we effectively used Weinhold et al's natural bond orbital analysis to assess the UV and FT-IR spectrophotometric monitoring of the change reflected in this phenomenon when the substituent in the benzene ring is altered. The NBO analysis was simulated by ab inito computations at the HF/6-31G(d) level of theory, in order to properly detect any possible presence of a hydrogen bond association. The changes occurring in electron occupancies of double-centered bonds, antibonding orbitals and in lone-pair orbitals appraised the results, as did the s and p character listings of the two-centered bonds and the simultaneous changes occurring in the geometric parameters of the molecules in question. Contrary to its normal preference, in these molecules the nitrogen used sp2 hybrid orbitals for its interaction, housing its electron lone-pair in the third p hybrid orbital. Furthermore, NBO analysis reflected the presence of a very soft intramolecular hydrogen association (C-H?π), labelled by UV and FT-IR assignments, between the benzene and triazole rings in all Schiff bases but p-N(Me)2. The n-π* stabilization energy decreased in the order: p-OH>p-OCH3>p-Cl>p-CH3>H>p-NO2>o-OH. The relation between the band position and Hammett substitution constant is interpreted in relation to the molecular structure.     

NMR studies of solvent-assisted proton transfer in a biologically relevant Schiff base: toward a distinction of geometric and equilibrium H-bond isotope effects

The tautomeric equilibrium in a Schiff base, N-(3,5-dibromosalicylidene)-methylamine 1, a model for the hydrogen bonded structure of the cofactor pyridoxal-5'-phosphate PLP which is located in the active site of the enzyme, was measured by means of 1H and 15N NMR and deuterium isotope effects on 15N chemical shifts at variable temperature and in different organic solvents. The position of the equilibrium was estimated using the one-bond 1J(OHN) and vicinal 3J(H(alpha)CNH) scalar coupling constants. Additionally, DFT calculations of a series of Schiff bases, N-(R1-salicylidene)-alkyl(R2)amines, were performed to obtain the hydrogen bond geometries. The latter made it possible to investigate a broad range of equilibrium positions. The increase of the polarity of the aprotic solvent shifts the proton in the intramolecular OHN hydrogen bond closer to the nitrogen. The addition of methanol and of hexafluoro-2-propanol to 1 in aprotic solvents models the PLP-water interaction in the enzymatic active site. The alcohols, which vary in acidity and change the polarity around the hydrogen bond, also stabilize the equilibrium, so that the proton is shifted to the nitrogen.     

Effect of cooperative hydrogen bonding in azo-hydrazone tautomerism of azo dyes

Azo-hydrazone tautomerism in azo dyes has been modeled by using density functional theory (DFT) at the B3LYP/6-31+G(d,p) level of theory. The most stable tautomer was determined both for model compounds and for azo dyes Acid Orange 7 and Solvent Yellow 14. The effects of the sulfonate group substitution and the replacement of the phenyl group with naphthyl on the tautomer stability and on the behavior in solvent have been discussed. Intramolecular hydrogen bond energies have been estimated for the azo and hydrazone tautomers to derive a relationship between the tautomer stability and the hydrogen bond strength. The transition structures for proton transfer displayed resonance assisted strong hydrogen bonding properties within the framework of the electrostatic-covalent hydrogen bond model (ECHBM). Evolution of the intramolecular hydrogen bond with changing structural and environmental factors during the tautomeric conversion process has been studied extensively by means of the atoms-in-molecules (AIM) analysis of the electron density. The bulk solvent effect was examined using the self-consistent reaction field model. Special solute-solvent interactions were further investigated by means of quantum mechanical calculations after defining the first-solvation shell by molecular dynamics simulations. The effect of cooperative hydrogen bonding with solvent molecules on the tautomer stability has been discussed.     

Note: Accounting for pressure effects on the calculated equilibrium structure of glassy GeSe(2)

First-principles molecular dynamics (FPMD) data on the structural properties of glassy GeSe(2) under ambient conditions are produced by carefully considering and minimizing the effect of a residual pressure on the periodic system. When compared to previous FPMD results, this strategy leads to an improved agreement between theory and neutron diffraction experiments.     

Intramolecular Hydrogen Bonds in Low‐Molecular‐Weight Polyethylene Glycol

We used static DFT calculations to analyze, in detail, the intramolecular hydrogen bonds formed in low‐molecular‐weight polyethylene glycol (PEG) with two to five repeat subunits. Both red‐shifted O?H???O and blue‐shifting C?H???O hydrogen bonds, which control the structural flexibility of PEG, were detected. To estimate the strength of these hydrogen bonds, the quantum theory of atoms in molecules was used. Car–Parrinello molecular dynamics simulations were used to mimic the structural rearrangements and hydrogen‐bond breaking/formation in the PEG molecule at 300 K. The time evolution of the H???O bond length and valence angles of the formed hydrogen bonds were fully analyzed. The characteristic hydrogen‐bonding patterns of low‐molecular‐weight PEG were described with an estimation of their lifetime. The theoretical results obtained, in particular the presence of weak C?H???O hydrogen bonds, could serve as an explanation of the PEG structural stability in the experimental investigation.     

Deuterium isotope effects on 13C NMR chemical shifts of some α-(2-hydroxyaryl)-N-phenylnitrones (Schiff base N-oxides)

The deuterium isotope effect on the ¹³C NMR chemical shifts of some α-2-hydroxyaryl-N-phenylnitrones (Schiff base N-oxides) was studied. The existence of an intramolecular hydrogen bond with the proton localized on the phenolic oxygen atom was evidenced. Exceptionally large isotope effects ΔC-2(D) and ΔC-α(D) suggest that the substitution of the proton of the OH group by deuterium leads to a weakening of the hydrogen bond and some conformational changes in the molecule. This conclusion was drawn on the basis of a comparison of the deuterium isotope effects of Schiff base N-oxides and parent Schiff bases. Copyright © 2001 John Wiley & Sons, Ltd.     

Car-Parrinello and path integral molecular dynamics study of the hydrogen bond in the chloroacetic acid dimer system

We have studied the double proton transfer (DPT) reaction in the cyclic dimer of chloroacetic acid using both classical and path integral Car-Parrinello molecular dynamics. We also attempt to quantify the errors in the potential energy surface that arise from the use of a pure density functional. In the classical dynamics a clear reaction mechanism can be identified, where asynchronized DPT arises due to coupling between the O-H stretching oscillator and several low energy intermolecular vibrational modes. This mechanism is considerably altered when quantum tunneling is permitted in the simulation. The introduction of path integrals leads to considerable changes in the thermally averaged molecular geometry, leading to shorter and more centered hydrogen bond linkages.     

气相水杨酸激发态分子内质子转移特性分析

采用密度泛函理论(DFT)和含时密度泛函理论(TD-DFT)研究了气相水杨酸(SA)分子的激发态氢键动力学过程。通过对水杨酸分子基态和激发态结构的优化,以及对其稳态吸收和发射光谱特性、前线分子轨道、红外振动光谱和势能曲线的计算分析,阐明水杨酸分子内质子转移可在激发态下自发地发生,导致其激发态可存在烯醇式和酮式两种异构体结构,并揭示了这种质子转移源于分子内电荷转移的激发态氢键的加强机制。     

Deuterium isotope effect on 13C chemical shifts of tetrabutylammonium salts of Schiff bases amino acids

Deuterium isotope effects on 13C chemical shift of tetrabutylammonium salts of Schiff bases, derivatives of amino acids (glycine, L-alanine, L-phenylalanine, L-valine, L-leucine, L-isoleucine and L-methionine) and various ortho-hydroxyaldehydes in CDCl3 have been measured. The results have shown that the tetrabutylammonium salts of the Schiff bases amino acids, being derivatives of 2-hydroxynaphthaldehyde and 3,5-dibromosalicylaldehyde, exist in the NH-form, while in the derivatives of salicylaldehyde and 5-bromosalicylaldehyde a proton transfer takes place. The interactions between COO- and NH groups stabilize the proton-transferred form through a bifurcated intramolecular hydrogen bond.     

Car-Parrinello molecular dynamics study of anharmonic systems: a Mannich base in solution

A Car-Parrinello molecular dynamics study was performed for 4,5-dimethyl-2-(N,N-dimethylaminomethyl)phenol, a Mannich base, to investigate the vibrational properties in solution of its intramolecular hydrogen bond. The dynamic behavior of this hydrogen-bonded system was investigated using an explicit solvent model. Addition of a nonpolar solvent permitted inclusion of delicate environmental effects on the strongly anharmonic system which was studied from first principles. Molecular dynamics and a posteriori quantization of the O-H motion were applied to reproduce the vibrational features of the O-H stretching mode. Consistent application of Car-Parrinello dynamics based on the density functional theory with subsequent solution of the vibrational Schr?dinger equation for the O-H stretching motion offers an effective method for strongly anharmonic systems, and this is supported by the comparison of the results with experimental spectra. As a further element of the intramolecular hydrogen bond study, the effects of deuteration were taken into account and a successful application of the O-H stretching mode quantization technique to the liquid phase is demonstrated. This provides a valuable computational methodology for investigations incorporating nuclear quantum effects in the liquid phase and enzyme active centers and can be used to investigate numerous systems that are not readily susceptible to experimental analysis.     

Theoretical study on the intramolecular hydrogen bond in chloro-substituted N,N-dimethylaminomethylphenols. I. Structural effects

Ab initio and density functional calculations are applied to study the influence of an increasing number of chlorine substituents on the properties of the intramolecular hydrogen bond in substituted Mannich bases. It is shown, that not only the acidity of the proton donor, which depends on the number of chlorine atoms at the aromatic ring, but also steric interactions modify the geometry of the hydrogen bond. Specific interactions of O–ClH–O hydrogen-bonding in some derivatives are estimated by calculations on related chlorophenols.     

Simulation of intramolecular hydrogen bond dynamics in manzamine A as a sensitive test for charge distribution quality

Subtle balance of inter- and intramolecular hydrogen bond strength in aqueous solutions often governs the structure and dynamics of molecular species used as potential drugs and in supramolecular applications. In silico molecular dynamics study of water solution of manzamine A has been performed with different atomic charges in order to investigate the influence of charge distribution choice on predicting qualitative and quantitative features of the simulated systems. Various well known charge schemes (MK-ESP, RESP, Mulliken, AMI-BCC, Gasteiger-Hückel, Gasteiger-Marsili, MMFF94, and Dynamic Electronegativity Relaxation - DENR) led to qualitatively different pictures of dynamic behavior of the intramolecular hydrogen bond. The reported calculation framework represents a relatively rare case where differences in charge distributions lead to noticeable differences in simulated properties, thus providing a useful test case for force field and charge distribution development, provided high quality experiments are conducted to use as references.     

Improving the convergence of closed and open path integral molecular dynamics via higher order Trotter factorization schemes

Higher order factorization schemes are developed for path integral molecular dynamics in order to improve the convergence of estimators for physical observables as a function of the Trotter number. The methods are based on the Takahashi-Imada and Susuki decompositions of the Boltzmann operator. The methods introduced improve the averages of the estimators by using the classical forces needed to carry out the dynamics to construct a posteriori weighting factors for standard path integral molecular dynamics. The new approaches are straightforward to implement in existing path integral codes and carry no significant overhead. The Suzuki higher order factorization was also used to improve the end-to-end distance estimator in open path integral molecular dynamics. The new schemes are tested in various model systems, including an ab initio path integral molecular dynamics calculation on the hydrogen molecule and a quantum water model. The proposed algorithms have potential utility for reducing the cost of path integral molecular dynamics calculations of bulk systems.     

Brønsted basicities of diamines in the gas phase, acetonitrile, and tetrahydrofuran

A comprehensive basicity study of alpha,omega-alkanediamines and related bases has been carried out. Basicities in acetonitrile (AN, pK(a) values), tetrahydrofuran (THF, pK(alpha) values), and gas phase (GP, GB values), were measured for 16, 14, and 9 diamine bases and for several related monoamines. In addition the gas-phase basicities and equilibrium geometries were computed for 19 diamino bases and several related monoamines at the DFT B3LYP 6-311+G** level. The effects of the different factors (intrinsic basicity of the amino groups, formation of intramolecular hydrogen bonds, and molecular strain) determining the diamine basicities were estimated by using the method of isodesmic reactions. The results are discussed in terms of molecular structure and solvation effects. The GP basicity is determined by the molecular size and polarizability, the extent of alkylation, and the energy effect of intramolecular hydrogen bond formation in the protonated base. The basicity trends in the solvents differ very much from those in GP: 1) The solvents severely compress the basicity range of the bases studied (3.5 times for the 1,3-propanediamine family in AN, and 7 times in THF), and 2) while stepwise alkylation of the basicity center leads to a steady basicity increase in the gas phase, the picture is complex in the solvents. Significant differences are also evident between THF and AN. The high hydrogen bond acceptor strength of THF leads to this solvent favoring the bases with "naked" protonation centers. In particular, the basicity order of N-methylated 1,3-propanediamines is practically inverse to that in the gas phase. The picture in AN is intermediate between that of GP and THF.