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.

Ability of the PM3 Quantum Mechanical Method to Model Intermolecular Hydrogen Bonding between Neutral Molecules

The PM3 semiempirical quantum-mechanical method was found to systematically describe intermolecular hydrogen bonding in small polar molecules. PM3 shows charge transfer from the donor to acceptor molecules on the order of 0.02-0.06 units of charge when strong hydrogen bonds are formed. The PM3 method is predictive; calculated hydrogen bond energies with an absolute magnitude greater than 2 kcal mol-' suggest that the global minimum is a hydrogen bonded complex; absolute energies less than 2 kcal mol-' imply that other van der Waals complexes are more stable. The geometries of the PM3 hydrogen bonded complexes agree with high-resolution spectroscopic observations, gas electron diffraction data, and high-level ab initio calculations. The main limitations in the PM3 method are the underestimation of hydrogen bond lengths by 0.1-0.2 for some systems and the underestimation of reliable experimental hydrogen bond energies by approximately 1-2 kcal mol-l. The PM3 method predicts that ammonia is a good hydrogen bond acceptor and a poor hydrogen donor when interacting with neutral molecules. Electronegativity differences between F, N, and 0 predict that donor strength follows the order F > 0 > N and acceptor strength follows the order N > 0 > F. In the calculations presented in this article, the PM3 method mirrors these electronegativity differences, predicting the F-H- - -N bond to be the strongest and the N-H- - -F bond the weakest. It appears that the PM3 Hamiltonian is able to model hydrogen bonding because of the reduction of two-center repulsive forces brought about by the parameterization of the Gaussian core-core interactions. The ability of the PM3 method to model intermolecular hydrogen bonding means reasonably accurate quantum-mechanical calculations can be applied to small biologic systems.

Atmospheric Implications for Formation of Clusters of Ammonium and 1−10 Water Molecules

A mixed molecular dynamics/quantum mechanics model has been applied to the ammonium/water clustering system. The use of the high level MP2 calculation method and correlated basis sets, such as aug-cc-pVDZ and aug-cc-pVTZ, lends confidence in the accuracy of the extrapolated energies. These calculations provide electronic and free energies for the formation of clusters of ammonium and 1−10 water molecules at two different temperatures. Structures and thermodynamic values are in good agreement with previous experimental and theoretical results. The estimated concentration of these clusters in the troposphere was calculated using atmospheric amounts of ammonium and water. Results show the favorability of forming these clusters and implications for ion-induced nucleation in the atmosphere.

Accurate Predictions of Water Cluster Formation, (H2O)n=2−10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

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.

Vibronic Spectra of Jet-Cooled 2-Aminopurine center dot H2O Clusters Studied by UV Resonant Two-Photon Ionization Spectroscopy and Quantum Chemical Calculations

For understanding the major- and minor-groove hydration patterns of DNAs and RNAs, it is important to understand the local solvation of individual nucleobases at the molecular level. We have investigated the 2-aminopurine center dot H2O. monohydrate by two-color resonant two-photon ionization and UV/UV hole-burning spectroscopies, which reveal two isomers, denoted A and B. The electronic spectral shift delta nu of the S-1 <- S-0 transition relative to bare 9H-2-aminopurine (9H-2AP) is small for isomer A (-70 cm(-1)), while that of isomer B is much larger (delta nu = 889 cm(-1)). B3LYP geometry optimizations with the TZVP basis set predict four cluster isomers, of which three are doubly H-bonded, with H2O acting as an acceptor to a N-H or -NH2 group and as a donor to either of the pyrimidine N sites. The "sugar-edge" isomer A is calculated to be the most stable form with binding energy D-e = 56.4 kJ/mol. Isomers B and C are H-bonded between the -NH2 group and pyrimidine moieties and are 2.5 and 6.9 kJ/mol less stable, respectively. Time-dependent (TD) B3LYP/TZVP calculations predict the adiabatic energies of the lowest (1)pi pi* states of A and B in excellent agreement with the observed 0(0)(0) bands; also, the relative intensities of the A and B origin bands agree well with the calculated S-0 state relative energies. This allows unequivocal identification of the isomers. The R2PI spectra of 9H-2AP and of isomer A exhibit intense low-frequency out-of-plane overtone and combination bands, which is interpreted as a coupling of the optically excited (1)pi pi* state to the lower-lying (1)n pi* dark state. In contrast, these overtone and combination bands are much weaker for isomer B, implying that the (1)pi pi* state of B is planar and decoupled from the (1)n pi* state. These observations agree with the calculations, which predict the (1)n pi* above the (1)pi pi* state for isomer B but below the (1)pi pi* for both 9H-2AP and isomer A.

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).

The impact of different nanoparticle surface chemistry and size on uptake and toxicity in a murine macrophage cell line

This study investigated the uptake, kinetics and cellular distribution of different surface coated quantum dots (QDs) before relating this to their toxicity. J774.A1 cells were treated with organic, COOH and NH2 (PEG) surface coated QDs (40 nM). Model 20 nm and 200 nm COOH-modified coated polystyrene beads (PBs) were also examined (50 microg ml(-1)). The potential for uptake of QDs was examined by both fixed and live cell confocal microscopy as well as by flow cytometry over 2 h. Both the COOH 20 nm and 200 nm PBs were clearly and rapidly taken up by the J774.A1 cells, with uptake of 20 nm PBs being relatively quicker and more extensive. Similarly, COOH QDs were clearly taken up by the macrophages. Uptake of NH2 (PEG) QDs was not detectable by live cell imaging however, was observed following 3D reconstruction of fixed cells, as well as by flow cytometry. Cells treated with organic QDs, monitored by live cell imaging, showed only a small amount of uptake in a relatively small number of cells. This uptake was insufficient to be detected by flow cytometry. Imaging of fixed cells was not possible due to a loss in cell integrity related to cytotoxicity. A significant reduction (p<0.05) in the fluorescent intensity in a cell-free environment was found with organic QDs, NH2 (PEG) QDs, 20 nm and 200 nm PBs at pH 4.0 (indicative of an endosome) after 2 h, suggesting reduced stability. No evidence of exocytosis was found over 2 h. These findings confirm that surface coating has a significant influence on the mode of NP interaction with cells, as well as the subsequent consequences of that interaction.

Quantum modeling of bioconjugated nanomaterials

In recent years, the bio-conjugated nanostructured materials have emerged as a new class of materials for the bio-sensing and medical diagnostics applications. In spite of their multi-directional applications, interfacing nanomaterials with bio-molecules has been a challenge due to somewhat limited knowledge about the underlying physics and chemistry behind these interactions and also for the complexity of biomolecules. The main objective of this dissertation is to provide such a detailed knowledge on bioconjugated nanomaterials toward their applications in designing the next generation of sensing devices. Specifically, we investigate the changes in the electronic properties of a boron nitride nanotube (BNNT) due to the adsorption of different bio-molecules, ranging from neutral (DNA/RNA nucleobases) to polar (amino acid molecules). BNNT is a typical member of III-V compounds semiconductors with morphology similar to that of carbon nanotubes (CNTs) but with its own distinct properties. More specifically, the natural affinity of BNNTs toward living cells with no apparent toxicity instigates the applications of BNNTs in drug delivery and cell therapy. Our results predict that the adsorption of DNA/RNA nucleobases on BNNTs amounts to different degrees of modulation in the band gap of BNNTs, which can be exploited for distinguishing these nucleobases from each other. Interestingly, for the polar amino acid molecules, the nature of interaction appeared to vary ranging from Coulombic, van der Waals and covalent depending on the polarity of the individual molecules, each with a different binding strength and amount of charge transfer involved in the interaction. The strong binding of amino acid molecules on the BNNTs explains the observed protein wrapping onto BNNTs without any linkers, unlike carbon nanotubes (CNTs). Additionally, the widely varying binding energies corresponding to different amino acid molecules toward BNNTs indicate to the suitability of BNNTs for the biosensing applications, as compared to the metallic CNTs. The calculated I-V characteristics in these bioconjugated nanotubes predict notable changes in the conductivity of BNNTs due to the physisorption of DNA/RNA nucleobases. This is not the case with metallic CNTs whose transport properties remained unaltered in their conjugated systems with the nucleobases. Collectively, the bioconjugated BNNTs are found to be an excellent system for the next generation sensing devices.

Spin ladders and quantum simulators for Tomonaga–Luttinger liquids

Magnetic insulators have proven to be usable as quantum simulators for itinerant interacting quantum systems. In particular the compound (C5H12N)2CuBr4 (for short: (Hpip)2CuBr4) was shown to be a remarkable realization of a Tomonaga–Luttinger liquid (TLL) and allowed us to quantitatively test the TLL theory. Substitution weakly disorders this class of compounds and thus allows us to use them to tackle questions pertaining to the effect of disorder in TLL as well, such as that of the formation of the Bose glass. In this paper we present, as a first step in this direction, a study of the properties of the related (Hpip)2CuCl4 compound. We determine the exchange couplings and compute the temperature and magnetic field dependence of the specific heat, using a finite temperature density matrix renormalization group procedure. Comparison with the measured specific heat at zero magnetic field confirms the exchange parameters and Hamiltonian for the (Hpip)2CuCl4 compound, giving the basis needed to begin studying the disorder effects.

50th anniversary of the Judd–Ofelt theory: An experimentalist's view of the formalism and its application

The theory on the intensities of 4f-4f transitions introduced by B.R. Judd and G.S. Ofelt in 1962 has become a center piece in rare-earth optical spectroscopy over the past five decades. Many fundamental studies have since explored the physical origins of the Judd–Ofelt theory and have proposed numerous extensions to the original model. A great number of studies have applied the Judd–Ofelt theory to a wide range of rare-earth doped materials, many of them with important applications in solid-state lasers, optical amplifiers, phosphors for displays and solid state lighting, upconversion and quantum-cutting materials, and fluorescent markers. This paper takes the view of the experimentalist who is interested in appreciating the basic concepts, implications, assumptions, and limitations of the Judd–Ofelt theory in order to properly apply it to practical problems. We first present the formalism for calculating the wavefunctions of 4f electronic states in a concise form and then show their application to the calculation and fitting of 4f-4f transition intensities. The potential, limitations and pitfalls of the theory are discussed, and a detailed case study of LaCl3:Er3+ is presented.

Optimizing infrared to near infrared upconversion quantum yield of β-NaYF4:Er3+ in fluoropolymer matrix for photovoltaic devices

The present study reports for the first time the optimization of the infrared (1523 nm) to near-infrared (980 nm) upconversion quantum yield (UC-QY) of hexagonal trivalent erbium doped sodium yttrium fluoride (β-NaYF4:Er3+) in a perfluorocyclobutane (PFCB) host matrix under monochromatic excitation. Maximum internal and external UC-QYs of 8.4% ± 0.8% and 6.5% ± 0.7%, respectively, have been achieved for 1523 nm excitation of 970 ± 43 Wm−2 for an optimum Er3+ concentration of 25 mol% and a phosphor concentration of 84.9 w/w% in the matrix. These results correspond to normalized internal and external efficiencies of 0.86 ± 0.12 cm2 W−1 and 0.67 ± 0.10 cm2 W−1, respectively. These are the highest values ever reported for β-NaYF4:Er3+ under monochromatic excitation. The special characteristics of both the UC phosphor β-NaYF4:Er3+ and the PFCB matrix give rise to this outstanding property. Detailed power and time dependent luminescence measurements reveal energy transfer upconversion as the dominant UC mechanism.