Complete thermodynamically consistent kinetic model of particle nucleation and growth: numerical study of the applicability of the classical theory of homogeneous nucleation

A complete thermodynamically consistent elementary reaction kinetic model of particle nucleation and growth from supersaturated vapor was developed and numerically evaluated to determine the conditions for the steady-state regime. The model treats all processes recognized in the aerosol science (such as nucleation, condensation, evaporation, agglomerationcoagulation, etc.) as reversible elementary reactions. It includes all possible forward reactions (i.e., of monomers, dimers, trimers, etc.) together with the thermodynamically consistent reverse processes. The model is built based on the Kelvin approximation, and has two dimensionless parameters: S0-the initial supersaturation and Theta-the dimensionless surface tension. The time evolution of the size distribution function was obtained over the ranges of parameters S0 and Theta. At low initial supersaturations, S0, the steady state is established after a delay, and the steady-state distribution function corresponds to the predictions of the classical nucleation theory. At high initial supersaturations, the depletion of monomers due to condensation on large clusters starts before the establishing of the steady state. The steady state is never reached, and the classical nucleation theory is not applicable. The boundary that separates these two regimes in the two dimensionless parameter space, S0 and Theta, was determined. The model was applied to several experiments on water nucleation in an expansion chamber [J. Wolk and R. Strey, J. Phys. Chem. B 105, 11683 (2001)] and in Laval nozzle [Y. J. Kim et al., J. Phys. Chem. A 108, 4365 (2004)]. The conditions of the experiments performed using Laval nozzle (S0=40-120) were found to be close to the boundary of the non-steady-state regime. Additional calculations have shown that in the non-steady-state regime the nucleation rate is sensitive to the rate constants of the initial steps of the nucleation process, such as the monomer-monomer, monomer-dimer, etc., reactions. This conclusion is particularly important for nucleation from supersaturated water vapor, since these processes for water molecules at and below the atmospheric pressure are in the low pressure limit, and the rate constants can be several orders of magnitude lower than the gas kinetic. In addition, the impact of the thermodynamic inconsistency of the previously developed partially reversible kinetic numerical models was assessed. At typical experimental conditions for water nucleation, S0=10 and Theta=10 (T=250 K), the error in the particle nucleation rate introduced by the thermodynamic inconsistency exceeds one order of magnitude.     

Diversity through isosterism: the case of boron-substituted 1,2-dihydro-1,2-azaborines

The first general synthesis of boron-substituted 1,2-dihydro-1,2-azaborines is described. The versatile 1,2-dihydro-1,2-azaborine precursor 4 is synthesized through a ring-closing metathesis-oxidation sequence. Treatment of 4 with a wide range of anionic nucleophiles furnishes the desired adducts 5 in good yields. The scope includes hydrogen- and a variety of carbon- and heteroatom-based nucleophiles. Furthermore, the boron-containing isostere (7) of the potent hypolipidemic agent, methyl 2-ethylphenoxyacetate (8), is readily prepared through our method.     

Does alpha-fluorination affect the structural trans-influence and kinetic trans-effect of an alkyl ligand? Molecular structures of Pd(TMEDA)(CH(3))(R(F)) and a kinetic study of the trans to cis isomerization of Pt(TMEDA)(CH(3))(2)I(R(F)) [R(F) = CF(2)CF(3), CFHCF(3), CH(2)CF(3)

Reaction of Pd(TMEDA)(CH(3))(2) [TMEDA = tetramethylethylenediamine] with fluoroalkyl iodides R(F)I affords a series of square planar Pd(II) complexes Pd(TMEDA)(CH(3))(R(F)) [R(F) = CF(2)CF(3) (9), CFHCF(3) (10), CH(2)CF(3) (11)], presumably by oxidative addition followed by reductive elimination of CH(3)I. The solid-state structures of each compound have been determined by single crystal X-ray diffraction studies, allowing the effect of increasing alpha-fluorination on the structural trans-influence of alkyl ligands to be examined. In these compounds there is no significant difference observed in the trans-influence of the three fluorinated alkyl ligands toward the trans-N atom, although a significant cis-influence on the neighboring methyl ligand is apparent. Oxidative addition of the same series of fluoroalkyl ligands to the corresponding Pt(TMEDA)(CH(3))(2) affords octahedral Pt(IV) complexes trans-Pt(TMEDA)(CH(3))(2)(R(F))I [R(F) = CF(2)CF(3) (12), CFHCF(3) (13), CH(2)CF(3) (14)] as the kinetic products. In each case, subsequent isomerization to the corresponding all cis-isomers is observed; in the case of 13, the stereocenter at the alpha-carbon results in two diastereomeric cis-isomers, which are formed at different rates. The molecular structures of 13 and its more stable all cis-isomer 16b have been crystallographically determined. Kinetic studies of the trans-cis isomerization reactions show the mechanism to involve a polar transition state, presumably involving iodide dissociation, followed by rearrangement of the cation, and iodide recombination. High dielectric solvents increase the rate, but solvent coordinating ability has no effect. Dissolved salts (LiI, LiOTf) show normal accelerative salt effects, with no inhibition in the case of added iodide, consistent with the formation of an intimate ion pair intermediate. The kinetic parameters show that the trans-effects of fluoroalkyl ligands in these compounds follow the order expected from the relative sigma-donor properties of the ligands, with CF(2)CF(3) < CFHCF(3) < CH(2)CF(3).     

Investigation of Functionalized α‐Chloroalkyllithiums for a Stereospecific Reagent‐Controlled Homologation Approach to the Analgesic Alkaloid (−)‐Epibatidine

Four putative functionalized α‐chloroakyllithiums RCH₂CHLiCl, where R=CHCH₂ ( 18 a ), CCH ( 18 b ), CH₂OBn ( 18 c ), and CH[O(CH₂)₂O] ( 18 d ), were generated in situ by sulfoxide–lithium exchange from α‐chlorosulfoxides, and investigated for the stereospecific reagent‐controlled homologation (StReCH) of phenethyl and 2‐chloropyrid‐5‐yl ( 17 ) pinacol boronic esters. Deuterium labeling experiments revealed that α‐chloroalkyllithiums are quenched by proton transfer from their α‐chlorosulfoxide precursors and it was established that this effect compromises the yield of StReCH reactions. Use of α‐deuterated α‐chlorosulfoxides was discovered to ameliorate the problem by retarding the rate of acid‐base chemistry between the carbenoid and its precursor. Carbenoids 18 a and 18 b showed poor StReCH efficacy, particularly the propargyl group bearing carbenoid 18 b , the instability of which was attributed to a facile 1,2‐hydride shift. By contrast, 18 d , a carbenoid that benefits from a stabilizing interaction between O and Li atoms gave good StReCH yields. Boronate 17 was chain extended by carbenoids 18 a , 18 b , and 18 d in 16, 0, and 68 % yield, respectively; α‐deuterated isotopomers D ‐ 18 a and D ‐ 18 d gave yields of 33 and 79 % for the same reaction. Double StReCH of 17 was pursued to target contiguous stereodiads appropriate for the total synthesis of (?)‐epibatidine ( 15 ). One‐pot double StReCH of boronate 17 by two exposures to (S)‐D ‐ 18 a (≤66 % ee), followed by work‐up with KOOH, gave the expected stereodiad product in 16 % yield (d.r.～67:33). The comparable reaction using two exposures to (S)‐D ‐ 18 d (≤90 % ee) delivered the expected bisacetal containing stereodiad (R,R)‐DD ‐ 48 in 40 % yield (≥98 % ee, d.r.=85:15). Double StReCH of 17 using (S)‐D ‐ 18 d (≤90 % ee) followed by (R)‐D ‐ 18 d (≤90 % ee) likewise gave (R,S)‐DD ‐ 48 in 49 % yield (≥97 % ee, d.r.=79:21). (R,S)‐DD ‐ 48 was converted to a dideuterated isotopomer of a synthetic intermediate in Corey’s synthesis of 15 .     

Rh(II) and Rh(I) two-legged piano-stool complexes: structure,reactivity, and electronic properties

The ligand 1,4-bis[4-(diphenylphosphino)butyl]-2,3,5,6-tetramethylbenzene, 3, was used to synthesize a mononuclear Rh(II) complex [(eta(1):eta(6):eta(1)-1,4-bis[4-(diphenylphosphino)butyl]-2,3,5,6-tetramethylbenzene)Rh][PF(6)](2), 6+, in a two-legged piano-stool geometry. The structural and electronic properties of this novel complex including a single-crystal EPR analysis are reported. The complex can be cleanly interconverted with its Rh(I) form, allowing for a comparison of the structural properties and reactivity of both oxidation states. The Rh(I) form 6 reacts with CO, tert-butyl isocyanide, and acetonitrile to form a series of 15-membered mononuclear cyclophanes [(eta(1):eta(1)-1,4-bis[4-(diphenylphosphino)butyl]-2,3,5,6-tetramethylbenzene)Rh(CO)(3)][PF(6)] (8), [(eta(1):eta(1)-1,4-bis[4-(diphenylphosphino)butyl]-2,3,5,6-tetramethylbenzene)Rh(CNC(CH(3))(3))(2)][PF(6)] (10), and [(eta(1):eta(1)-1,4-bis[4-(diphenylphosphino)butyl]-2,3,5,6-tetramethylbenzene)Rh(CO)(CH(3)CN)][PF(6)] (11). The Rh(II) complex 6+ reacts with the same small molecules, but over shorter periods of time, to form the same Rh(I) products. In addition, a model two-legged piano-stool complex [(eta(1):eta(6):eta(1)-1,4-bis[3-(diphenylphosphino)propoxy]-2,3,5,6-tetramethylbenzene)Rh][B(C(6)F(5))(4)], 5, has been synthesized and characterized for comparison purposes. The solid-state structures of complexes 5, 6, 6+, and 11 are reported. Structure data for 5: triclinic; P(-)1; a = 10.1587(7) A; b = 11.5228(8) A; c = 17.2381(12) A; alpha = 96.4379(13) degrees; beta = 91.1870(12) degrees; gamma = 106.1470(13) degrees; Z = 2. 6: triclinic; P(-)1; a = 11.1934(5) A; b = 12.4807(6) A; c = 16.1771(7) A; alpha = 81.935(7) degrees; beta = 89.943(1) degrees; gamma = 78.292(1) degrees; Z = 2. 6+: monoclinic; P2(1)/n; a = 11.9371(18) A; b = 32.401(5) A; c = 12.782(2) A; beta = 102.890(3) degrees; Z = 4. 11: triclinic; P(-)1; a = 13.5476(7) A; b = 13.8306(7) A; c = 14.9948(8) A; alpha = 74.551(1) degrees; beta = 73.895(1) degrees; gamma = 66.046(1) degrees; Z = 2.     

Simulation of the radiation-chemical transformations of acetic acid in aqueous solutions

A general model for the radiolysis of acetic acid and its aqueous solutions is proposed. The model adequately describes experimental data on the degradation of the acid and the formation of gases (H₂, CO₂, and CH₄) in aqueous solutions at various pH values.     

Solution-phase structural characterization of supramolecular assemblies by molecular diffraction

Structures of four molecular squares based on rhenium coordination chemistry have been characterized in the solution phase using pair distribution function (PDF) analyses of wide-angle X-ray scattering measured to better than 1 A spatial resolution. In this report we have focused, in particular, on a comparison of structures for pyrazine- and bipyridine-edged squares measured in solution with structures determined for these molecules in the solid state using X-ray crystallography and models derived from geometry optimization and molecular dynamics simulations using a classical force field. The wide-angle scattering for these assemblies is dominated by pair correlations involving one or more rhenium atoms, with both edge and diagonal Re-Re interactions appearing prominently in PDF plots. The pyrazine square is characterized by a relatively rigid structure in solution, with PDF peak positions and linewidths corresponding closely to those calculated from crystal structure data. For the bipyridine-edged square, the experimental PDF peaks measured along the molecular sides match the positions and linewidths of the PDF peaks calculated from static models. In contrast, PDF peaks measured across the diagonal distances of the molecular square deviate significantly from those calculated from the static crystallographic and energy minimized models. The experimental data are instead indicative of configurational broadening of the diagonal distances. In this respect, molecular dynamics simulations point to the importance of butterfly type motions that modulate the Re-Re diagonal distance. Indeed, the experimental data are reasonably well fit by assuming a bimodal distribution of butterfly conformers differing by approximately 25 degrees in the Re-Re-Re-Re torsion angle. Additionally, the measurements provide evidence for solvent ordering by the supramolecular assemblies detected as regions of solvent association and exclusion.     

Molecular Structure and Conformational Composition of 1-[(Methylthio)methyl]-2-nitrobenzene (MTMNB): A Theoretical and Experimental Study

The conformational composition of gaseous MTMNB and the molecular structures of the rotational forms have been studied by electron diffraction at 130^∘C aided by results from ab initio and density functional theory calculations. The conformational potential energy surface has been investigated by using the B3LYP/6-31G(d,p) method. As a result, six minimum-energy conformers have been identified. Geometries of all conformers were optimized using MP2/6-31G(d,p), B3LYP/6-31G(d,p), and B3LYP/cc-pVTZ methods. These calculations resulted in accurate geometries, relative energies, and harmonic vibrational frequencies for all conformers. The B3LYP/cc-pVTZ energies were then used to calculate the Boltzmann distribution of conformers. The best fit of the electron diffraction data to calculated values was obtained for the six conformer model, in agreement with the theoretical predictions. Average parameter values (r_a in angstroms, angle α in degrees, and estimated total errors given in parentheses) weighted for the mixture of six conformers are r(C–C) = 1.507(5), r(C–C)_{ring, av} = 1.397(3), r(C–S)_av = 1.814(4), r(C–N) = 1.495(4), r(N–O)_av = 1.223(3), ∠(C–C–C)_ring = 116.0–122.5, ∠ C6–C4–C7 = 118.2(4), ∠ C–C–S = 113.6(6), ∠ C–S–C = 98.5(12), ∠ N–C–C4 = 121.9(3), ∠(O–N–C)_av = 116.8(3), ∠ O–N–O = 127.0(4). Torsional angles could not be refined. Theoretical B3LYP/cc-pVTZ torsional angles for the rotation about C–N bond, φ_C−N, were found to be 30.5–36.5^∘ for different conformers. As to internal rotation about C–C and C–S bonds, values of φ_C−C = 68–118^∘ and φ_C−S = 66–71^∘ were obtained for the three most stable conformers with gauche orientation with respect to these bonds. Some conclusions of this work were presented in a short communication in Russ. J. Phys. Chem. 2005, 79, 1701.     

Synthesis of aza-, oxa-, and thiabicyclo[3.1.0]hexane heterocycles from a common synthetic intermediate

[reaction: see text] An efficient and stereospecific approach to the synthesis of structurally constrained aza-, oxa-, and thiabicyclo[3.1.0]hexane heterocycles has been achieved through application of the intramolecular cyclopropanation reaction of diazoacetates. The various constrained heterocycles (X = N, O, or S) are conveniently prepared from a common diol intermediate accessible from readily available cinnamyl alcohols. Application of the methodology to the synthesis of conformationally constrained oxazolidinone antibacterials is also discussed.     

Bis(m-phenylene)-32-crown-10-based cryptands, powerful hosts for paraquat derivatives

Four new bis(m-phenylene)-32-crown-10-based cryptands with different third bridges were prepared. Their complexes with paraquat derivatives were studied by proton NMR spectroscopy, mass spectrometry, and X-ray analysis. It was found that these cryptands bind paraquat derivatives very strongly. Specifically, a diester cryptand with a pyridyl nitrogen atom located at a site occupied by either water or a PF(6) anion in analogous complexes exhibited the highest association constant K(a) = 5.0 x 10(6) M(-1) in acetone with paraquat, 9000 times greater than the crown ether system. X-ray structures of this and analogous complexes demonstrate that improved complexation with this host is a consequence of preorganization, adequate ring size for occupation by the guest, and the proper location of the pyridyl N-atom for binding to the beta-pyridinium hydrogens of the paraquat guests. This readily accessible cryptand is one of the most powerful hosts reported for paraquats.