A Computationally Efficient Procedure for Modeling the First Step in the Alkaline Hydrolysis of Esters

A computationally efficient procedure for modeling the alkaline hydrolysis of esters is proposed based on calculations performed on methyl acetate and methyl benzoate systems. Extensive geometry and energy comparisons were performed on the simple ester methyl acetate. The effectiveness of performing high level single point ab initio energy calculations on the geometries obtained from semiempirical and ab initio methods was determined. The AM1 and PM3 semiempirical methods are evaluated for their ability to model the transition states and intermediates for ester hydrolysis. The Cramer/Truhlar SM3 solvation method was used to determine activation energies. The most computationally efficient way to model the transition states of large esters is to use the PM3 method. The PM3 transition structure can then be used as a template for the design of haptens capable of inducing catalytic antibodies.

Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1 clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Hydration of the Bisulfate Ion: Atmospheric Implications

Using molecular dynamics configurational sampling combined with ab initio energy calculations, we determined the low energy isomers of the bisulfate hydrates. We calculated the CCSD(T) complete basis set (CBS) binding electronic and Gibbs free energies for 53 low energy isomers of HSO4–(H2O)n=1–6 and derived the thermodynamics of adding waters sequentially to the bisulfate ion and its hydrates. Comparing the HSO4–/H2O system to the neutral H2SO4/H2O cluster, water binds more strongly to the anion than it does to the neutral molecules. The difference in the binding thermodynamics of HSO4–/H2O and H2SO4/H2O systems decreases with increasing number of waters. The thermodynamics for the formation of HSO4–(H2O)n=1–5 is favorable at 298.15 K, and that of HSO4–(H2O)n=1–6 is favorable for T < 273.15 K. The HSO4– ion is almost always hydrated at temperatures and relative humidity values encountered in the troposphere. Because the bisulfate ion binds more strongly to sulfuric acid than it does to water, it is expected to play a role in ion-induced nucleation by forming a strong complex with sulfuric acid and water, thus facilitating the formation of a critical nucleus.

Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Hydration of the Sulfuric Acid−Methylamine Complex and Implications for Aerosol Formation

The binary H2SO4−H2O nucleation is one of the most important pathways by which aerosols form in the atmosphere, and the presence of ternary species like amines increases aerosol formation rates. In this study, we focus on the hydration of a ternary system of sulfuric acid (H2SO4), methylamine (NH2CH3), and up to six waters to evaluate its implications for aerosol formation. By combining molecular dynamics (MD) sampling with high-level ab initio calculations, we determine the thermodynamics of forming H2SO4(NH2CH3)(H2O)n, where n = 0−6. Because it is a strong acid−base system, H2SO4−NH2CH3 quickly forms a tightly bound HSO4−−NH3CH3+ complex that condenses water more readily than H2SO4 alone. The electronic binding energy of H2SO4−NH2CH3 is −21.8 kcal mol−1 compared with −16.8 kcal mol−1 for H2SO4−NH3 and −12.8 kcal mol−1 for H2SO4−H2O. Adding one to two water molecules to the H2SO4−NH2CH3 complex is more favorable than adding to H2SO4 alone, yet there is no systematic difference for n ≥ 3. However, the average number of water molecules around H2SO4−NH2CH3 is consistently higher than that of H2SO4, and it is fairly independent of temperature and relative humidity.

Hydration of the Sulfuric Acid-Methylamine Complex and Implications for Aerosol Formation

The binary H2SO4-H2O nucleation is one of the most important pathways by which aerosols form in the atmosphere, and the presence of ternary species like amines increases aerosol formation rates. In this study, we focus on the hydration of a ternary system of sulfuric acid (H2SO4), methylamine (NH2CH3), and up to six waters to evaluate its implications for aerosol formation. By combining molecular dynamics (MD) sampling with high-level ab initio calculations, we determine the thermodynamics of forming H2SO4(NH2CH3)(H2O)n, where n = 0-6. Because it is a strong acid-base system, H2SO4-NH2CH3 quickly forms a tightly bound HSO4(-)-NH3CH3(+) complex that condenses water more readily than H2SO4 alone. The electronic binding energy of H2SO4-NH2CH3 is -21.8 kcal mol(-1) compared with -16.8 kcal mol(-1) for H2SO4-NH3 and -12.8 kcal mol(-1) for H2SO4-H2O. Adding one to two water molecules to the H2SO4-NH2CH3 complex is more favorable than adding to H2SO4 alone, yet there is no systematic difference for n ≥ 3. However, the average number of water molecules around H2SO4-NH2CH3 is consistently higher than that of H2SO4, and it is fairly independent of temperature and relative humidity.

Photodissociation dynamics of tert-butylnitrite following excitation to the S(1) and S(2) states. A study by velocity-map ion-imaging and 3D-REMPI spectroscopy

Excitation of tert-butylnitrite into the first and second UV absorption bands leads to efficient dissociation into the fragment radicals NO and tert-butoxy in their electronic ground states (2)Π and (2)E, respectively. Velocity distributions and angular anisotropies for the NO fragment in several hundred rotational and vibrational quantum states were obtained by velocity-map imaging and the recently developed 3D-REMPI method. Excitation into the well resolved vibronic progression bands (k = 0, 1, 2) of the NO stretch mode in the S(1) ← S(0) transition produces NO fragments mostly in the vibrational state with v = k, with smaller fractions in v = k - 1 and v = k - 2. It is concluded that dissociation occurs on the purely repulsive PES of S(1) without barrier. All velocity distributions from photolysis via the S(1)(nπ*) state are monomodal and show high negative anisotropy (β ≈ -1). The rotational distributions peak near j = 30.5 irrespective of the vibronic state S(1)(k) excited and the vibrational state v of the NO fragment. On average 46% of the excess energy is converted to kinetic energy, 23% and 31% remain as internal energy in the NO fragment and the t-BuO radical, respectively. Photolysis via excitation into the S(2) ← S(0) transition at 227 nm yields NO fragments with about equal populations in v = 0 and v = 1. The rotational distributions have a single maximum near j = 59.5. The velocity distributions are monomodal with positive anisotropy β ≈ 0.8. The average fractions of the excess energy distributed into translation, internal energy of NO, and internal energy of t-BuO are 39%, 23%, and 38%, respectively. In all cases ∼8500 cm(-1) of energy remain in the internal degrees of freedom of the t-BuO fragment. This is mostly assigned to rotational energy. An ab initio calculation of the dynamic reaction path shows that not only the NO fragment but also the t-BuO fragment gain large angular momentum during dissociation on the purely repulsive potential energy surface of S(2).

Charge and spin transport in nanoscale junction from first principles

Recently nanoscale junctions consisting of 0-D nanostructures (single molecule) or 1-D nanostructures (semiconducting nanowire) sandwiched between two metal electrodes are successfully fabricated and characterized. What lacks in the recent developments is the understanding of the mechanism behind the observed phenomena at the level of atoms and electrons. For example, the origin of observed switching effect in a semiconducting nanowire due to the influence of an external gate bias is not yet understood at the electronic structure level. On the same context, different experimental groups have reported different signs in tunneling magneto-resistance for the same organic spin valve structure, which has baffled researchers working in this field. In this thesis, we present the answers to some of these subtle questions by investigating the charge and spin transport in different nanoscale junctions. A parameter-free, single particle Green’s function approach in conjunction with a posteriori density functional theory (DFT) involving a hybrid orbital dependent functional is used to calculate the tunneling current in the coherent transport limit. The effect of spin polarization is explicitly incorporated to investigate spin transport in a nanoscale junction. Through the electron transport studies in PbS nanowire junction, a new orbital controlled mechanism behind the switching of the current is proposed. It can explain the switching behavior, not only in PbS nanowire, but in other lead-chalcogenide nanowires as well. Beside this, the electronic structure properties of this nanowire are studied using periodic DFT. The quantum confinement effect was investigated by calculating the bandgap of PbS nanowires with different diameters. Subsequently, we explain an observed semiconducting to metallic phase transition of this nanowire by calculating the bandgap of the nanowire under uniform radial strain. The compressive radial strain on the nanowire was found to be responsible for the metallic to semiconducting phase transition. Apart from studying one dimensional nanostructure, we also present transport properties in zero dimensional single molecular junctions. We proposed a new codoping approach in a single molecular carborane junction, where a cation and an anion are simultaneously doped to find the role of a single atom in the device. The main purpose was to build a molecular junction where a single atom can dictate the flow of electrons in a circuit. Recent observations of both positive and negative sign in tunneling magnetoresistance (TMR) the using same organic spin-valve structure hasmystified researchers. From our spin dependent transport studies in a prototypical organic molecular tunneling device, we found that a 3% change in metal-molecule interfacial distance can alter the sign of TMR. Changing the interfacial distance by 3%, the number of participating eigenstates as well as their orbital characteristic changes for anti-parallel configuration of the magnetization at the two electrodes, leading to the sign reversal of the TMR. Apart from this, the magnetic proximity effect under applied bias is investigated quantitatively, which can be used to understand the observed unexpectedmagnetismin carbon basedmaterials when they are in close proximity with magnetic substrates.

Supramolecular Organization of Heptapyrenotide Oligomers—An in Depth Investigation by Molecular Dynamics Simulations

We present a molecular modeling study based on ab initio and classical molecular dynamics calculations, for the investigation of the tridimensional structure and supramolecular assembly formation of heptapyrenotide oligomers in water solution. Our calculations show that free oligomers self-assemble in helical structures characterized by an inner core formed by π- stacked pyrene units, and external grooves formed by the linker moieties. The coiling of the linkers has high ordering, dominated by hydrogen-bond interactions among the phosphate and amide groups. Our models support a mechanism of longitudinal supramolecular oligomerization based on interstrand pyrene intercalation. Only a minimal number of pyrene units intercalate at one end, favoring formation of very extended longitudinal chains, as also detected by AFM experiment. Our results provide a structural explanation of the mechanism of chirality amplification in 1:1 mixtures of standard heptapyrenotides and modified oligomers with covalently linked deoxycytidine, based on selective molecular recognition and binding of the nucleotide to the groove of the left-wound helix.

The origin of the low-spin character of the resting state of cytochrome P450cam investigated by means of active site analogues

Crown-capped iron(S−) porphyrins 1·H2O and 2·H2O and their corresponding Ba2+ complexes have been prepared as active site analogues of the resting state of cytochrome P450cam. cw-EPR studies and electronic structure calculations at the density functional theory (DFT) level of model systems suggest a functional role of the water cluster of P450cam.

Liquid–Vapor Interface of Formic Acid Solutions in Salt Water: A Comparison of Macroscopic Surface Tension and Microscopic in Situ X-ray Photoelectron Spectroscopy Measurements

The liquid–vapor interface is difficult to access experimentally but is of interest from a theoretical and applied point of view and has particular importance in atmospheric aerosol chemistry. Here we examine the liquid–vapor interface for mixtures of water, sodium chloride, and formic acid, an abundant chemical in the atmosphere. We compare the results of surface tension and X-ray photoelectron spectroscopy (XPS) measurements over a wide range of formic acid concentrations. Surface tension measurements provide a macroscopic characterization of solutions ranging from 0 to 3 M sodium chloride and from 0 to over 0.5 mole fraction formic acid. Sodium chloride was found to be a weak salting out agent for formic acid with surface excess depending only slightly on salt concentration. In situ XPS provides a complementary molecular level description about the liquid–vapor interface. XPS measurements over an experimental probe depth of 51 Å gave the C 1s to O 1s ratio for both total oxygen and oxygen from water. XPS also provides detailed electronic structure information that is inaccessible by surface tension. Density functional theory calculations were performed to understand the observed shift in C 1s binding energies to lower values with increasing formic acid concentration. Part of the experimental −0.2 eV shift can be assigned to the solution composition changing from predominantly monomers of formic acid to a combination of monomers and dimers; however, the lack of an appropriate reference to calibrate the absolute BE scale at high formic acid mole fraction complicates the interpretation. Our data are consistent with surface tension measurements yielding a significantly more surface sensitive measurement than XPS due to the relatively weak propensity of formic acid for the interface. A simple model allowed us to replicate the XPS results under the assumption that the surface excess was contained in the top four angstroms of solution.

The calculation of Afρ and mass loss rate for comets

Ab initio calculations of Afρ are presented using Mie scattering theory and a Direct Simulation Monte Carlo (DSMC) dust outflow model in support of the Rosetta mission and its target 67P/Churyumov-Gerasimenko (CG). These calculations are performed for particle sizes ranging from 0.010 μm to 1.0 cm. The present status of our knowledge of various differential particle size distributions is reviewed and a variety of particle size distributions is used to explore their effect on Afρ , and the dust mass production View the MathML sourcem˙. A new simple two parameter particle size distribution that curtails the effect of particles below 1 μm is developed. The contributions of all particle sizes are summed to get a resulting overall Afρ. The resultant Afρ could not easily be predicted a priori and turned out to be considerably more constraining regarding the mass loss rate than expected. It is found that a proper calculation of Afρ combined with a good Afρ measurement can constrain the dust/gas ratio in the coma of comets as well as other methods presently available. Phase curves of Afρ versus scattering angle are calculated and produce good agreement with observational data. The major conclusions of our calculations are: – The original definition of A in Afρ is problematical and Afρ should be: qsca(n,λ)×p(g)×f×ρqsca(n,λ)×p(g)×f×ρ. Nevertheless, we keep the present nomenclature of Afρ as a measured quantity for an ensemble of coma particles.– The ratio between Afρ and the dust mass loss rate View the MathML sourcem˙ is dominated by the particle size distribution. – For most particle size distributions presently in use, small particles in the range from 0.10 to 1.0 μm contribute a large fraction to Afρ. – Simplifying the calculation of Afρ by considering only large particles and approximating qsca does not represent a realistic model. Mie scattering theory or if necessary, more complex scattering calculations must be used. – For the commonly used particle size distribution, dn/da ∼ a−3.5 to a−4, there is a natural cut off in Afρ contribution for both small and large particles. – The scattering phase function must be taken into account for each particle size; otherwise the contribution of large particles can be over-estimated by a factor of 10. – Using an imaginary index of refraction of i = 0.10 does not produce sufficient backscattering to match observational data. – A mixture of dark particles with i ⩾ 0.10 and brighter silicate particles with i ⩽ 0.04 matches the observed phase curves quite well. – Using current observational constraints, we find the dust/gas mass-production ratio of CG at 1.3 AU is confined to a range of 0.03–0.5 with a reasonably likely value around 0.1.

Microwave Spectroscopy and Structures of Several Substituted Acetylenes and Fluorocarbons

The structures of two substituted acetylene compounds have been characterized from their microwave rotational spectra. In the first study, two structures of 6-methyl-3-heptyne have been determined. This compound can be thought of as an ethyl group separated from an isobutyl group by a C≡C spacer. Both structures have the ethyl and isobutyl groups eclipsed, consistent with the dominant interaction determining the orientation about the acetylene axis being the weak dispersion attraction between the end groups. One structure is with the isobutyl group in a symmetric conformation and the other with the isobutyl group asymmetric. In addition, the microwave spectrum of the butane analogue 3,5-octadiyne has been observed. This compound consists of two ethyl groups separated by two C≡C spacers. The study is still in progress, but it appears that the ethyl end groups are freely rotating. Therefore, it seems that the dispersion attractions between the end groups are too weak at this longer distance of about 7 Å. The structures of several fluorocarbons have also been studied by microwave spectroscopy. The structures of perfluoropentane and perfluorohexane have been shown to be helical, like the polymer polytetrafluoroethylene (Teflon©). The structure of perfluoropropane and two conformers of 1H-heptafluoropropane have been determined to be non-helical. It is apparent that the steric and dipole repulsions between fluorine atoms that have been attributed to the helical structure of longer fluorocarbon chains are not sufficient in a three carbon chain to cause a twist in the structures.

Bank Performance: Market Power or Efficient Structure?

Regulatory change not seen since the Great Depression swept the U.S. banking industry beginning in the early 1980s, culminating with the Interstate Banking and Branching Efficiency Act of 1994. Significant consolidations have occurred in the banking industry. This paper considers the market-power versus the efficient-structure theories of the positive correlation between banking concentration and performance on a state-by-state basis. Temporal causality tests imply that bank concentration leads bank profitability, supporting the market-power, rather than the efficient-structure, theory of that positive correlation. Our finding suggests that bank regulators, by focusing on local banking markets, missed the initial stages of an important structural change at the state level.

Estudio a escala laboratorio y planta piloto, de la adsorción de NO y SO2, sobre CR2O3/AL2O3 a altas temperaturas, provenientes de fuentes fijas

Resumen: Se planificaron las experiencias con el objeto de analizar el comportamiento del catalizador en la columna metálica de mayor diámetro. Se modificaron las masas usadas para verificar la eficiencia de retención respecto de la masa. Se realizaron ciclos de adsorción, desorción y readsorción sobre una misma muestra para determinar variaciones en la eficiencia del catalizador. En otra fase, en colaboración con el Dr. V. A. Ranea y el Prof. E. E. Mola (INIFTA, UNLP), se desarrolló el estudio teórico de la adsorción de moléculas de SO2, CH4, CO2, O2 y CO sobre Cr2O3(0001) mediante Teoría del Funcional Densidad (programa VASP, Vienna Ab-initio Simulation Package), y el estudio de la cinética de la reacción entre CH4, SO2 y el O2 junto con la presencia de especies sulfito y sulfato. Este estudio permitió hallar los sitios preferenciales de adsorción de Sº y la posible competencia con SO2 experimentalmente y por cálculos teóricos. Dentro del marco de la presente línea de investigación, la Ing. Sabrina Hernández Guiance continúa realizando experiencias en el marco del proyecto conjunto con el INIFTA, las cuales forman parte del desarrollo de su tesis doctoral. Experimentalmente, se observa que la eficiencia de adsorción del catalizador respecto al SO2 es cercana al 100%. Se observa un pico de termodesorción a 1120 K. Luego, se estudió la oxidación de CH4 con SO2. Se observa que hay producción de CO2 desde temperatura inicial, seguida de un aumento significativo en la formación de CO2 hasta 330-340 K. Luego, la producción de CO2 se mantiene aproximadamente constante. Mediante el empleo de la ecuación de Arrhenius y resultados experimentales, se obtuvo la energía de activación de la reacción global, de 7 Kcal/mol. Mediante estudios teóricos, se determinó que la energía de quimisorción del SO2 sobre el Cr2O3 es de -3.09 eV para la configuración más estable, una energía de adsorción de O2 en estado disociativo de -1.567 eV, una energía para CH4 sobre O2 adsorbido previamente de -0.335 eV, y - 0.812 eV para la configuración más estable de CO2 sobre el sustrato.