Theoretical Study of the Gaseous Hydrolysis of NO2 in the Presence of Amines.

The effects on the hydrolysis of NO2 in the presence of methylamine and dimethylamine molecules were investigated by theoretical calculations of a series of the molecular clusters 2NO2-mH2O-CH3NH2 (m = 1-3) and 2NO2-mH2O-(CH3)2NH (m = 1, 2). With methylamine included in the clusters, the energy barrier is reduced by 3.2 kcal/mol from that with ammonia, and the corresponding products may form without an energy barrier. The results show that amines have larger effects than ammonia in promoting the hydrolysis of NO2 on thermodynamics. The additional water molecules can stabilize the transition states and the product complexes, and we infer that adding more water molecules in the reactions mainly act as solvent and promoting to form the methylamine nitrate (CH3NH3+NO3-). In addition, the interactions of CH3NH2 and (CH3)2NH on the hydration of HNO3 are also more effective than NH3, and the NH4NO3, CH3NH3NO3, and (CH3)2NH2NO3 complexes tend to form the larger aerosols with the increasing of water molecules. The equilibrium geometries, harmonic vibrational frequencies, and intensities of both HONO-CH3NH2 and HONO-NH3 complexes were investigated. Calculations predict that the binding energies of both HONO-CH3NH2 complexes are larger than HONO-NH3 complexes, and the OH stretching vibrational frequencies and intensities are most affected. The natural bond orbital analysis was performed to describe the donor-acceptor interactions on a series of complexes in the reactions 2NO2 + H2O + CH3NH2 and 2NO2 + H2O + (CH3)2NH, as well as the complexes of HONO-NH3 and HONO-CH3NH2. The results show that the interactions with amines are relatively larger, and the higher stabilization energies between CH3NH2 and HONO are found.

Hydrolysis of Sulfur Dioxide in Small Clusters of Sulfuric Acid: Mechanistic and Kinetic Study.

The deposition and hydrolysis reaction of SO2 + H2O in small clusters of sulfuric acid and water are studied by theoretical calculations of the molecular clusters SO2-(H2SO4)n-(H2O)m (m = 1,2; n = 1,2). Sulfuric acid exhibits a dramatic catalytic effect on the hydrolysis reaction of SO2 as it lowers the energy barrier by over 20 kcal/mol. The reaction with monohydrated sulfuric acid (SO2 + H2O + H2SO4 - H2O) has the lowest energy barrier of 3.83 kcal/mol, in which the cluster H2SO4-(H2O)2 forms initially at the entrance channel. The energy barriers for the three hydrolysis reactions are in the order SO2 + (H2SO4)-H2O > SO2 + (H2SO4)2-H2O > SO2 + H2SO4-H2O. Furthermore, sulfurous acid is more strongly bonded to the hydrated sulfuric acid (or dimer) clusters than the corresponding reactant (monohydrated SO2). Consequently, sulfuric acid promotes the hydrolysis of SO2 both kinetically and thermodynamically. Kinetics simulations have been performed to study the importance of these reactions in the reduction of atmospheric SO2. The results will give a new insight on how the pre-existing aerosols catalyze the hydrolysis of SO2, leading to the formation and growth of new particles.

Mechanism of the gaseous hydrolysis reaction of SO2: Effects of NH3 versus H2O.

Effects of ammonia and water molecules on the hydrolysis of sulfur dioxide are investigated by theoretical calculations of two series of the molecular clusters SO2-(H2O)n (n = 1-5) and SO2-(H2O)n-NH3 (n = 1-3). The reaction in pure water clusters is thermodynamically unfavorable. The additional water in the clusters reduces the energy barrier for the reaction, and the effect of each water decreases with the increasing number of water molecules in the clusters. There is a considerable energy barrier for reaction in SO2-(H2O)5, 5.69 kcal/mol. With ammonia included in the cluster, SO2-(H2O)n-NH3, the energy barrier is dramatically reduced, to 1.89 kcal/mol with n = 3, and the corresponding product of hydrated ammonium bisulfate NH4HSO3-(H2O)2 is also stabilized thermodynamically. The present study shows that ammonia has larger kinetic and thermodynamic effects than water in promoting the hydrolysis reaction of SO2 in small clusters favorable in the atmosphere.

Theoretical study of the hydrogen abstraction of substituted phenols by nitrogen dioxide as a source of HONO.

The mild yet promiscuous reactions of nitrogen dioxide (NO2) and phenolic derivatives to produce nitrous acid (HONO) have been explored with density functional theory calculations. The reaction is found to occur via four distinct pathways with both proton coupled electron transfer (PCET) and hydrogen atom transfer (HAT) mechanisms available. While the parent reaction with phenol may not be significant in the gas phase, electron donating groups in the ortho and para positions facilitate the reduction of nitrogen dioxide by electronically stabilizing the product phenoxy radical. Hydrogen bonding groups in the ortho position may additionally stabilize the nascent resonantly stabilized radical product, thus enhancing the reaction. Catechol (ortho-hydroxy phenol) has a predicted overall free energy change ΔG(0) = -0.8 kcal mol(-1) and electronic activation energy Ea = 7.0 kcal mol(-1). Free amines at the ortho and para positions have ΔG(0) = -3.8 and -1.5 kcal mol(-1); Ea = 2.3 and 2.1 kcal mol(-1), respectively. The results indicate that the hydrogen abstraction reactions of these substituted phenols by NO2 are fast and spontaneous. Hammett constants produce a linear correlation with bond dissociation energy (BDE) demonstrating that the BDE is the main parameter controlling the dark abstraction reaction. The implications for atmospheric chemistry and ground-level nitrous acid production are discussed.

A combined experimental and theoretical study of dinitrosyl iron complexes containing chelating bis(diphenyl)phosphinoX (X = benzene, propane and ethylene): X-ray crystal structures and properties influenced by the presence or absence of π-bonds in chelating ligands.

Recent discoveries involving the roles of nitric oxide in humans have stimulated intense interest in transition metal nitrosyl complexes. A series of dinitrosyl iron complexes with the formula [(DPPX)Fe(NO)2], {DPPX = 1,2-bis(diphenylphosphino)benzene (1), 1,3-bis(diphenylphosphino)propane (2), and cis-1,2-bis(diphenylphosphino)ethylene (3)} has been prepared and characterized through a combination of FT-IR, NMR, UV-vis, X-ray crystallography, and electrochemical techniques. Infrared spectroscopy showed NO shifts to the region of 1723 and 1674 cm-1 for complexes 1 and 3, and 1708 and 1660 cm-1 for 2, indicating that ligand 2 acts as a stronger σ-donor. The X-ray crystallographic data showed that 1 and 3 possess the rare repulso conformation while 2 has the attracto conformation. CV studies on compounds 1, 2 and 3 display two quasi-reversible oxidations with the E°1/2 values at 0.101 and 0.186 V, 0.121 and 0.184 V, and 0.019 and 0.342 V, respectively. The larger ΔE value for compound 2 compared with that of 1 and 3 is attributed to the lack of π-bonds between the two phosphorus atoms. Theoretical calculations using density functional theory were carried out on the synthesized compounds and model compounds and the results are consistent with the experimental data. The calculated HOMO-LUMO gaps for compounds 1, 2 and 3 are 3.736, 4.060, and 3.669 eV, respectively, which supports the stronger back-donation for compound 2 than that of compounds 1 and 3.

Theoretical study on stable small clusters of oxalic acid with ammonia and water.

Thermodynamically stable small clusters of oxalic acid (CO2H)2, ammonia (NH3), and water (H2O) are studied through quantum chemical calculations. The (CO2H)2-NH3 core system with up to three waters of hydration was examined by B3LYP density functional theory and MP2 molecular orbital theory with the aug-cc-pVDZ basis set. The (CO2H)2-NH3 core complexes are observed to hydrogen bond strongly and should be found in appreciably significant concentrations in the atmosphere. Subsequent hydration of the (CO2H)2-NH3 core, however, is found to be somewhat prohibitive under ambient conditions. Relative populations of the examined clusters are predicted and the binding patterns detailed. Atmospheric implications related to new particle formations are discussed.

Ab initio calculations of the Ar-ethane intermolecular potential energy surface using bond function basis sets.

The intermolecular potential energy surface (PES) of argon with ethane has been studied by ab initio calculations at the levels of second-order Møller-Plesset perturbation (MP2) theory and coupled-cluster theory with single, double, and noniterative triple configurations (CCSD(T)) using a series of augmented correlation-consistent basis sets. Two sets of bond functions, bf1 (3s3p2d) and bf2 (6s6p4d2f), have been added to the basis sets to show a dramatic and systematic improvement in the convergence of the entire PES. The PES of Ar-ethane is characterized by a global minimum at a near T-shaped configuration with a well depth of 0.611 kcal mol(-1), a second minimum at a collinear configuration with a well depth of 0.456 kcal mol(-1), and a saddle point connecting the two minima. It is shown that an augmented correlation-consistent basis set with a set of bond functions, either bf1 or bf2, can effectively produce results equivalent to the next larger augmented correlation-consistent basis set, that is, aug-cc-pVDZ-bf1 ≈ aug-cc-pVTZ, aug-cc-pVTZ-bf1 ≈ aug-cc-pVQZ. Very importantly, the use of bond functions improves the PES globally, resulting accurate potential anisotropy. Finally, MP2 method is inadequate for accurate calculations, because it gives a potentially overestimated well depth and, more seriously, a poor potential anisotropy.

Theoretical study on the structure and stabilities of molecular clusters of oxalic acid with water.

The importance of aerosols to humankind is well-known, playing an integral role in determining Earth's climate and influencing human health. Despite this fact, much remains unknown about the initial events of nucleation. In this work, the molecular properties of common organic atmospheric pollutant oxalic acid and its gas phase interactions with water have been thoroughly examined. Local minima single-point energies for the monomer conformations were calculated at the B3LYP and MP2 level of theory with both 6-311++G(d,p) and aug-cc-pVDZ basis sets and are compared with previous works. Optimized geometries, relative energies, and free energy changes for the stable clusters of oxalic acid conformers with up to six waters were then obtained from B3LYP calculations with 6-31+G(d) and 6-311++G(d,p) basis sets. Initially, cooperative binding is predicted to be the most important factor in nucleation, but as the clusters grow, dipole cancellations are found to play a pivotal role. The clusters of oxalic acid hydrated purely with water tend to produce extremely stable and neutral core systems. Free energies of formation and atmospheric implications are discussed.

Titanium-mediated rearrangement of cyclopropenylmethyl acetates to (E)-halodienes.

TiCl(4) and TiBr(4) rapidly transform cyclopropenylmethyl acetates to (E)-halodienes via ring-opening to allyl-vinyl cations. DFT calculations suggest that the regioselectivity of the halogenation of this cationic intermediate by [TiX(4)OAc](-) is under thermodynamic control, while the stereoselectivity is governed by kinetics.

Interpreting vibrational sum-frequency spectra of sulfur dioxide at the air/water interface: a comprehensive molecular dynamics study.

We investigated the solvation and spectroscopic properties of SO(2) at the air/water interface using molecular simulation. Molecular interactions from both Kohn-Sham (KS) density functional theory (DFT) and classical polarizable models were used to understand the properties of SO(2):(H(2)O)(x) complexes in the vicinity of the air/water interface. The KS-DFT was included to allow comparisons with vibrational sum-frequency spectroscopy through the identification of surface SO(2):(H(2)O)(x) complexes. Using our simulation results, we were able to develop a much more detailed picture of the surface structure of SO(2) consistent with spectroscopic data obtained by Richmond and co-workers (J. Am. Chem. Soc. 2005, 127, 16806). We also found many similarities and differences between the two interaction potentials, including a noticeable weakness of the classical potential model in reproducing the asymmetric hydrogen bonding of water with SO(2) due to its inability to account for SO(2) resonance structures.

Theoretical study of the quantitative structure-activity relationships for the toxicity of dibenzo-p-dioxins.

Density functional theory calculations of polychlorinated dibenzo-p-dioxins (PCDDs) were carried out to obtain the electronic descriptors, polarizabilities, and traceless quadrupole moments of 76 PCDD congeners. No simple relationships were found for the electronic descriptors with the relevant aryl hydrocarbon receptor (AhR) binding affinities of PCDDs, which suggests that they alone may not be sufficient to explain the variation in toxicity. However, quantitative structure-activity relationships (QSARs) were developed with the polarizabilities and traceless quadrupole moments, explaining about 74% and 59% of variation in AhR binding affinities of PCDDs, respectively. To explain the nature of toxic interaction, a mathematical model based on the ligand-receptor binding and solute-solvent interaction was presented, and then a multiple regression analysis of all the above parameters was performed to evaluate the contributions of the parameters to the bonding affinities. Other data for PCDFs found in the literature were also included in the regression analysis to minimize data over-fitting. The results suggest that dispersion and electrostatic interactions are equally important for the interaction of PCDD/Fs with the AhR.

Correlations and predictions of carboxylic acid pKa values using intermolecular structure and properties of hydrogen-bonded complexes.

Density functional theory calculations have been preformed on a series of hydrogen-bonded complexes of substituted aliphatic and aromatic carboxylic acids with ammonia. Molecular properties, particularly those related to hydrogen bonding, have been carefully examined for their interdependence as well as dependence on the acidity of the acid. The bond length and stretching frequency of the hydroxyl group and the hydrogen-bond length and energy of the complex are shown to be highly correlated with each other and are linearly correlated with available literature pKa values of the carboxylic acids. The linear correlations resulting from the fit to the available pKa values can be used to predict the pKa values of similar carboxylic acids. The pKa values so predicted using the different molecular properties are highly consistent and in good agreement with the literature values. This study suggests that calculated molecular properties of hydrogen-bonded complexes allow effective and systematic prediction of pKa values for a large range of organic acids using the established linear correlations. This approach is unique in its capability to determine the acidity of a particular functional group or the local acidity within a large molecular system such as a protein.

Theoretical study on the chemical properties of polybrominated diphenyl ethers.

Density functional theory calculations at the B3LYP/6-31+G(d) and B3LYP/aug-cc-pVDZ levels were performed to obtain the equilibrium structures, thermodynamic properties, and electron affinities (EA) of 14 polybrominated diphenyl ether (PBDE) congeners in the gas phase. All congeners except for those of symmetric BDE are found to have two or more conformational isomers. The optimized geometries of the most stable conformational isomers are in agreement with recently published X-ray crystallographic data. The thermodynamic properties of the congeners with a given number of bromine substitutions are strongly dependent on the substitution pattern, whereas the EA values also depend on the number of bromine substitutions. The vertical electron affinities (EA(Ver)) calculated for the selected BDE congeners at the B3LYP/aug-cc-pVDZ level are all positive except for di-BDEs, and are correlated with the initial reductive debromination rate constants obtained recently [Keum, Y.-S., Li, Q.X., 2005. Reductive debromination of polybrominated diphenyl ethers by zerovalent iron. Environ. Sci. Technol. 39, 2280]. All adiabatic electron affinities (EA(Ada)) are positive, and suggest that the BDE congeners act as electron acceptors when reacting with receptors in living cells. The calculated EA(Ada) values differ considerably from those of EA(Ver) because of the large geometrical relaxation from the neutral to the anionic BDE congeners, highlighted by the lengthening of a C-Br bond. The elongated C-Br bond, which occurs at the alpha position, is directly involved in the debromination of n-bromodiphenyl to (n-1)-bromodiphenyl ethers in the reductive debromination experiments.

Efficient bond function basis set for pi-pi interaction energies.

Systematic study has been carried out to investigate the accuracy of mid-bond functions in describing pi-pi interactions in benzene dimer. Potential energy curves are calculated for the sandwich, T-shaped, and parallel-displaced configurations of benzene dimer by adding bond functions in MP2 (second-order Møller-Plesset perturbation theory) calculations with a wide range of split-valence and augmented, correlation-consistent basis sizes. At MP2 level, the largest basis set used with a bond function (denoted aug-cc-pVDZf-6s6p4d2f) differs by only approximately 0.1 kcal/mol relative to the result obtained from the standard aug-cc-pVQZ basis calculation (without the bond function). The calculated potential energy curves from the bond function basis aug-cc-pVDZf-6s6p4d2f and the larger standard basis aug-cc-pVTZ are in excellent agreement with each other for all three configurations. The number of bond function basis aug-cc-pVDZf-6s6p4d2f is 526 compared to 828 of aug-cc-pVTZ and 1512 of aug-cc-pVQZ. Current study shows that bond functions can be effectively employed to give accurate description of pi-pi interactions with the addition of only a minimal number of bond functions.

Structures, reductive dechlorination, and electron affinities of selected polychlorinated dibenzo-p-dioxins: density functional theory study.

Density functional theory calculations were performed to obtain the structures, vertical electron affinities, and adiabatic affinities of 15 polychlorinated dibenzo-p-dioxins (PCDDs), including several extremely toxic congeners. A three-parameter hybrid density functional, B3LYP, was utilized with two different basis sets, 6-311G(d,p) and 6-311+G(2d,2p). The optimized structures of all PCDDs under consideration were planar, while all corresponding anions attained nonplanar geometries. One of the C-Cl bonds on each PCDD anion was considerably elongated, and the dechlorination of PCDDs occurred as the departing chlorine bent off the aromatic ring plane for effective pi-sigma orbital mixing. The characteristic electron energy-dependent regioselective chloride ion loss channels for 1,2,3,7,8-pentaCDD were elucidated by transition-state theory calculations. The relative low-energy barrier for the dechlorination of 1,2,3,7,8-pentaCDD indicated the high likelihood of obtaining reductive dechlorination (RD) products that are more toxic than the parent species. The calculated vertical electron affinities of PCDDs are consistent with the available experimental attachment energies, and the positive adiabatic electron affinities suggest that PCDDs may act as electron acceptors in living cells.

Physical origin for the nonlinear sorption of very hydrophobic organic chemicals in a membrane-like polymer film.

Bioconcentration factor (BCF) is often assumed to be linearly associated with the octanol-water partition coefficient K(ow) for hydrophobic organic chemicals (HOCs). However, a large amount of data has suggested that the correlation between the logBCF and logK(ow) is curvilinear for HOCs. Similar curvilinear relationship has also been noticed for sorption of HOCs into poly(dimethyl)siloxane (PDMS), a polymer with cross-linked interior structures. So far no satisfactory explanation has been given to account for the deviation. In this study, we acquired additional experimental data to show that the curvilinear relationship between the log-based PDMS-coated fiber-water partition coefficient (logK(f)) and logK(ow) for polychlorinated biphenyls (PCBs) was indeed a reflection of the sorption process occurring in PDMS film other than experimental defects. The physical origin of the nonlinearity was pinpointed based on the theory of phase partitioning for HOCs. The linear relationship is observed if the solute molecule is considerably smaller than the size of a monomer unit of PDMS in that the Gibbs free energy required for cavity formation in PDMS is comparable to that in octanol. Higher free energy of cavity formation is needed to create sufficient free volume if the PCB molecular size is comparable to or larger than the monomer unit of PDMS. On the other hand, the free energy of cavity formation in octanol remains almost constant when this occurs, resulting in the observed curvilinear relationship. The proposed model adequately explains the observed data, as well as sheds lights into the physical origin of the steric interactions of large molecular size solute with the PDMS polymer network.

Reaction mechanisms and kinetics for the oxidations of dimethyl sulfide, dimethyl disulfide, and methyl mercaptan by the nitrate radical.

This study examines the initial oxidation routes of the three major reduced sulfur compounds (CH(3)SH, CH(3)SCH(3), and CH(3)SSCH(3)) by the nitrate radical using density functional and ab initio methods. Stationary points along each reaction pathway are examined using different levels of theory and basis sets to ensure the convergence of the results. Kinetics calculations follow on the determined reaction pathways to obtain the rate constants. This study shows that sulfur compounds exhibit a general trend of hydrogen abstraction following the formation of an initial sulfur-nitrate complex. The results are in agreement with experimental work on CH(3)SCH(3) and CH(3)SH, while refuting a proposal of several previous studies that oxygen addition is the dominant oxidation pathway in the case of CH(3)SSCH(3). The rate constants obtained from kinetics calculations are consistent with experimental findings and exhibit negative temperature dependence. Overall, this study confirms the importance of nitrate in the oxidation of reduced sulfur compounds in the atmosphere.

Correlations and predictions of pKa values of fluorophenols and bromophenols using hydrogen-bonded complexes with ammonia.

Density functional theory calculations have been preformed on a series of the hydrogen-bonded fluorophenol-ammonia and bromophenol-ammonia complexes. Intermolecular and intramolecular properties, particularly those related to hydrogen bonding, have been carefully analyzed. Several properties, such as the bond length and stretching frequency of the hydroxyl group, the hydrogen bond length and binding energy, are shown to be highly correlated with each other and are linearly correlated with known experimental pKa values of the halogenated phenols. The linear correlations have been used to predict the pKa values of all fluorophenols and bromophenols in the series. The predicted pKa values are shown to be consistent from different molecular properties and are in good agreement with available experimental values. This study suggests that calculated molecular properties of hydrogen-bonded complexes allow the effective and systematic prediction of pKa values for a large range of organic acids using the established linear correlations.

Asymmetrical linear structures including three-electron hemibonds or other interactions in the (ABA)-type triatomic cations: Ne3+, (He-Ne-He)+, (Ar-Ne-Ar)+, (Ar-O-Ar)+, (He-O-He)+, and (Ar-He-Ar)+.

By the counterpoise geometry optimization at the level of CCSD(T)aug-cc-pVDZ, the asymmetrical linear structures with all the real frequencies were obtained for the triatomic cations of (ABA)+ type: Ne3+, (He-Ne-He)+, (Ar-Ne-Ar)+, (Ar-He-Ar)+, (He-O-He)+, and (Ar-O-Ar)+. The validity of this optimization method is confirmed by comparing with the method of the potential-energy surface for the calculations of Ne3+ and (He-Ne-He)+. Using the molecular-orbital theory, it is found that the interaction within the triatomic cations is dominated by the contribution from the first two atoms while the contribution from the third atom is small. This result is justified as a direct consequence of forming an asymmetrical linear structure. Specifically, four types of interaction within the triatomic cations are identified: three-electron sigma-type hemibond, three-electron pi-type hemibond, two-electron sigma bond, and the attraction between cation and atoms. For Ne3+, (He-Ne-He)+, and (He-O-He)+ clusters, it is shown that the electron correlation effect supports the asymmetry.

Theoretical study of hydrogen-bonded complexes of chlorophenols with water or ammonia: Correlations and predictions of pKa values.

A simple practical method for predicting the acidity constants (as pKa values) of chlorophenols is proposed based on density functional theory calculations of a series of hydrogen-bonded complexes of phenol and 19 different congeners of chlorophenol, with a single probe molecule, either water or ammonia. Relevant structural parameters and molecular properties of these complexes, primarily involving the acidic hydroxyl group, are examined and plotted against the known pKa values of 14 chlorophenols and phenol. Strong linear correlations are found for these compounds. Such correlations are used to determine the pKa values of five chlorophenols whose experimental acidities have large uncertainties. Similar predicted pKa values are obtained by using different structural parameters and molecular properties for the complexes with either probe molecule. The study may be extended to determine the acidity of other compounds with a single acidic functional group.