You are here

Publications

[1] Andrei Buin, Haiyan Wang, Styliani Consta, and Yining Huang. A study of conformational equilibrium of 1,1,2-trichloroethane in FAU-type zeolites. MICROPOROUS AND MESOPOROUS MATERIALS, 183:207-217, JAN 1 2014.

We study the effects of loading and nature of the zeolitic framework in the trans reversible arrow gauche equilibrium of 1,1,2-trichloroethane (TCE) in zeolites with Faujasite (FAU) structure including sodium Y (Na-Y) and siliceous Y (Si-Y). TCE is used as a prototype molecule for the study of zeolite-chlorinated hydrocarbon interactions that determine conformational equilibrium. Conformational changes of TCE were sampled using empirical force fields and quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) in combination with replica exchange molecular dynamics (REMD). The population ratio of the conformers at high loading was compared with that found by Fourier transform (FT)-Raman spectroscopy and the results are in very good agreement. Comparison between the infinitely dilute limit of loading, and higher loading shows that the population ratio of the conformers at equilibrium depends strongly on the Si/Al ratio and consequently, the presence of charge compensating cations as well as the loading. At very low loading, in both Na-Y and Si-Y the population of the trans conformer is the largest with the exception of the empirical modeling that predicts the gauche conformer to be dominant in Na-Y. At the higher loading, in Na-Y the population the gauche conformer constitutes the great majority. In Si-Y the population of the trans conformer continues to dominate but the percentage of the gauche conformer increases substantially relative to the loading at the infinitely dilute limit. Such different conformational behaviors provide insight into the factors responsible for diffusion mechanisms and adsorption energies in the FAU-type zeolitic frameworks. (C) 2013 Elsevier Inc. All rights reserved.

[2] Anatoly Malevanets and Styliani Consta. Variation of droplet acidity during evaporation. JOURNAL OF CHEMICAL PHYSICS, 138(18), MAY 14 2013.

Variation of acidity and associated chemical changes of macromolecules in evaporating droplets is of central importance in electrosprayed aerosols. We study changes in acidity during evolution of a droplet that is composed of solvent and a charge binding macromolecule. We analyze the acidity of the droplet using analytical theory and stochastic modeling. We derive a universal relation for the minimum pH of a droplet in the presence of a protein and the results are confirmed by the stochastic modeling of ubiquitin and lysozyme at varying values of pH. We establish that in acidic droplets, once the number of solvated charges reaches the macroion charge, the further droplet evaporation, counter-intuitively, reduces the number of charges on the macromolecule and increases the droplet pH. (C) 2013 AIP Publishing LLC.

[3] Styliani Consta and Anatoly Malevanets. Classification of the ejection mechanisms of charged macromolecules from liquid droplets. JOURNAL OF CHEMICAL PHYSICS, 138(4), JAN 28 2013.

The relation between the charge state of a macromolecule and its ejection mechanism from droplets is one of the important questions in electrospray ionization methods. In this article, effects of solvent-solute interaction on the manifestation of the charge induced instability in a droplet are examined. We studied the instabilities in a prototype system of a droplet comprised of charged poly(ethylene glycol) and methanol, acetonitrile, and water solvents. We observed instances of three, previously only conjectured, [S. Consta, J. Phys. Chem. B 114, 5263 (2010)] mechanisms of macroion ejection. The mechanism of ejection of charged macroion in methanol is reminiscent of “pearl†model in polymer physics. In acetonitrile droplets, the instability manifests through formation of solvent spines around the solvated macroion. In water, we find that the macroion is ejected from the droplet through contiguous extrusion of a part of the chain. The difference in the morphology of the instabilities is attributed to the interplay between forces arising from the macroion solvation energy and the surface energy of the droplet interface. For the contiguous extrusion of a charged macromolecule from a droplet, we demonstrate that the proposed mechanism leads to ejection of the macromolecule from droplets with sizes well below the Rayleigh limit. The ejected macromolecule may hold charge significantly higher than that suggested by prevailing theories. The simulations reveal new mechanisms of macroion evaporation that differ from conventional charge residue model and ion evaporation mechanisms. (C) 2013 American Institute of Physics. [http://dx.doi.org/10.1063/1.4789018]

[4] Styliani Consta and Anatoly Malevanets. Manifestations of Charge Induced Instability in Droplets Effected by Charged Macromolecules. PHYSICAL REVIEW LETTERS, 109(14), OCT 4 2012.

Ion-release processes from droplets that contain excess charge are of central importance in determining the charge-state distributions of macromolecules in electrospray ionization methods. We develop an analytical theory to describe the mechanism of contiguous extrusion of a charged macromolecule from a droplet. We find that the universal parameter determining the system behavior is the ratio of solvation energy per unit length to the square of the ion charge density per unit length. Systems with the same value of the ratio will follow the same path in the course of droplet evaporation. The analytical model is compared with molecular simulations of charged polyethylene glycol macroion in aqueous droplets, and the results are in excellent agreement.

[5] Jun Kyung Chung and Styliani Consta. Release Mechanisms of Poly(ethylene glycol) Macroions from Aqueous Charged Nanodroplets. JOURNAL OF PHYSICAL CHEMISTRY B, 116(19):5777-5785, MAY 17 2012.

Ion-release processes in nanodroplets that contain excess charge are of central importance in atmospheric aerosols as well as in determining the charge state distributions of macroions that are detected in electrospray mass spectrometry (ESMS) experiments. We performed molecular simulations of systems of a poly(ethylene glycol) (PEG) associated with various ions (Na+, Li+,Ca2+) in aqueous charged nanodroplets in order to investigate the manner that the macroion emerges from an aqueous nanodroplet as well as its final charge state. In the study we focused on a specific region of the parameter space with respect to charge and size of droplets that is close to the Rayleigh limit. We found that for sizes of droplets with linear dimensions of several nanometers and length of PEG up to 100 monomers, the PEG macroion emerges from the droplet following a three-step process: (i) phase separation, (ii) gradual extension of the macroion out of the droplet, and (iii) drying-out of the solvent or spontaneous detachment of the macroion from the droplet. The third step is determined by the ratio of charge on the macroion to the ions in the water portion of the droplet. The chemical transformation that is caused in PEG by the transfer of ions from the solvent into PEG determines its release mechanism. When the charge is carried by macroions, the charge-induced instability manifests by following one of the expected scenarios of Rayleigh instability; however, the assumptions of the Rayleigh model break down. We also examined the release of the macroion below the Rayleigh limit, and we found that the macroion emerges from the droplet by drying-out of the solvent. On balance of phenomenological evidence, we concluded that the ion-evaporation mechanism (IEM) in its most common meaning is not the followed mechanism for the parameter space of the systems that we studied. The final charge state of the macroion is in excellent agreement with the experimental data of Fenn and co-workers.

[6] M. Paliy, O. M. Braun, and S. Consta. Friction in a Thin Water Layer: Dissociative versus Nondissociative Friction. JOURNAL OF PHYSICAL CHEMISTRY C, 116(16):8932-8942, APR 26 2012.

An ultrathin water film confined between two substrates in moving contact is studied using Langevin molecular dynamics with a coordinate- and velocity-dependent damping coefficient. The water molecules are modeled with the Central-Force model that allows for the dissociation of water molecules into H+ and OH-. Two different friction scenarios are found depending on the applied pressure and the strength of the interaction of water with substrates. Under low loads, the water molecules stay intact during the frictional sliding. However, when the applied pressure increases past a value of similar to 20 GPa, the water molecules begin to dissociate and recombine immediately in the course of sliding, which results in a large increase of friction. The rate of such dissociation is found to be roughly proportional to the speed of driving. The relation of the observed phenomena to the “superionic†and “ionic fluid†states of water and its relevance in practical friction situations is discussed.

[7] Andrei Buin, Jianrong Ma, Yining Huang, Styliani Consta, and Zhang Hui. Conformational Changes of trans-1,2-Dichlorocyclohexane Adsorbed in Zeolites Studied by FT-Raman Spectroscopy and Molecular QM/MM Simulations. JOURNAL OF PHYSICAL CHEMISTRY C, 116(15):8608-8618, APR 19 2012.

The conformational behavior of trans-1,2-dichlorocyclohexane (T12D) adsorbed inside zeolites with Faujasite structure (FAU) including sodium Y (Na-Y) and siliceous Y (Si-Y) was investigated by FT-Raman spectroscopy and molecular simulations. The results have clearly shown that the conformational and dynamic properties of T12D strongly depend on the presence of charge-balancing cations as well as Si/Al ratio. Upon loading into Na-Y, the population of the diequatorial (ee) conformer increases compared with pure T12D liquid due to the strong interaction with the extra-framework Na+ ions. Molecular simulations of T12D in Na-Y and Si-Y have also been carried out. The T12D molecule was modeled by quantum and semiempirical quantum chemistry methods and embedded in a zeolite framework that was described by empirical force field. Conformational changes were sampled using quantum mechanics/molecular mechanics molecular dynamics and replica exchange molecular dynamics. Molecular simulations revealed that the ee conformer is preferable versus the diaxial (aa) conformer in both Na-Y and Si-Y frameworks. However, in the Na-Y supercage environment, the Na+ ions polarize the ee conformer stronger than in Si-Y. This leads to a larger shift of the conformational equilibrium in favor of the ee component in Na-Y relative to Si Y. Computations of the equilibrium population of aa and ee conformers using the dihedral CI-C-C-Cl angle distributions of aa and ee showed good quantitative agreement with the experimental findings with respect to the dominant ee conformer in both Na-Y and Si-Y.

[8] Doaa M. Ragab, Sohrab Rohani, and Styliani Consta. Controlled release of 5-fluorouracil and progesterone from magnetic nanoaggregates. INTERNATIONAL JOURNAL OF NANOMEDICINE, 7:3167-3189, 2012.

Background: The potential use of magnetic nanoparticles in biomedical applications has witnessed an exponential growth in recent years. Methods: In this study, we used nanoaggregates of magnetic nanoparticles as carriers for controlled drug delivery. The nanoaggregates are formed due to the presence of the block copolymer of polyethylene oxide-polypropylene oxide (Pluronic F-68) and beta-cyclodextrin that surround the magnetic core of the nanoparticles. The administration of the drug carriers occurs by inhalation, and the drug is delivered systemically via the pulmonary route. We tested the delivery of 5-fluorouracil and progesterone, which are used as models of hydrophilic and hydrophobic drugs, respectively. Results: The estimated nanoaggregates' diameters are between 293 nm +/- 14.65 nm and 90.2 nm +/- 4.51 nm, respectively. In-situ and post-synthesis techniques are two approaches for drug loading. The polymer composition of nanoaggregates and initial drug concentration showed a significant effect on both the drug entrapment efficiency and release kinetics. Average drug entrapment efficiencies ranged between 16.11% and 83.25%. In-situ loaded samples showed significantly slower release rates. The drug release mechanism is investigated by mathematical curve fitting to different drug release kinetics models. In most cases, the Peppas model has shown good correlations (coefficients of correlation, R-2, between 0.85 and 0.99) with the examined release profiles. The estimated release indices are below 0.5, which indicates the Fickian diffusion mechanism. For samples with an initial burst effect, the modified Peppas model can provide a better understanding of the drug release mechanism, both in the samples loaded with progesterone, or those high polymer concentrations. Conclusion: Our work showed prolonged delivery of drugs (5-fluorouracil and progesterone) by diffusion from nanoaggregates, with the potential to reduce dose-related adverse effects.

[9] Andrei Buin, Styliani Consta, and Tsun-Kong Sham. Mechanisms of Phase Transformations of TiO2 Nanotubes and Nanorods. JOURNAL OF PHYSICAL CHEMISTRY C, 115(45):22257-22264, NOV 17 2011.

Phase transformations of titanium dioxide (TiO2) nanotubes and nanorods at elevated temperatures are studied using molecular dynamics (MD) and replica exchange molecular dynamics (REMD) utilized here in the same way as simulated annealing. In the study, TiO2 nanotubes with amorphous (amT) and anatase structures (anT) as well as TiO2 amorphous nanorods (amR) that are amenable to experimental investigation are considered at various temperatures. It is found that amT and amR transform into a rutile rod, while anT transforms into a brookite nanotube. It is demonstrated that transformation of anT starts from TiO4 and TiO5 complexes found in the surfaces of the system, in contrast to amT and amR, where initial grains of the new phase may develop throughout the entire system starting from TiO5 and TiO6 complexes. The evolution of the number of TiOx (x = 4, 5, 6, 7) complexes indicates that the transformation of amT and amR occurs almost suddenly relative to the transformation of anT. The initial grains of transformation of amT have a structure close to rutile, while those of anT have brookite features. To our knowledge, we report the first simulations of phase transformations of TiO2 nanotubes and nanorods where simulations are performed beyond mu s.

[10] Styliani Consta and Jun Kyung Chung. Charge-Induced Conformational Changes of PEG-(Na-n(+)) in a Vacuum and Aqueous Nanodroplets. JOURNAL OF PHYSICAL CHEMISTRY B, 115(35):10447-10455, SEP 8 2011.

Charged-induced conformational changes of sodiated poly(ethylene glycol) (PEG-(Na+)(n)) in a vacuum and water droplets were studied using molecular simulations. In a vacuum, compact and partially unwound conformations were identified by analyzing occupation probabilities on reaction surfaces defined by the accessible surface area of the macroion, the distance between the centers of mass of the PEG molecule and of the cations and radius of gyration of the macromolecule. The critical charge of the macromolecule for which there is coexistence of various conformations was estimated using the Rayleigh criterion for the breakdown of highly charged droplets and compared with that observed in the simulations. The simulation findings agreed well with the Rayleigh prediction. The properties of sodiated PEG macromolecules in the presence of solvent and ions were also investigated. It was found that the macroion becomes saturated with charge. The highly charged state leads to an extended conformation that is partially expelled by the droplet. However, a portion of the chain is still in contact with the solvent. Drying-out of the solvent leaves the macroion in a stretched necklace-like conformation. The mechanism of release of sodiated PEG from an aqueous droplet was shown to be distinctly different from that of a protonated polyhistidine in a water nanodroplet, demonstrating sensitivity of the disintegration mechanism of the nanodroplet to the degree of hydrophilicity of the macroion.

[11] Yi Liu, Styliani Consta, and William A. Goddard, III. Nanoimmiscibility: Selective Absorption of Liquid Methanol-Water Mixtures in Carbon Nanotubes. JOURNAL OF NANOSCIENCE AND NANOTECHNOLOGY, 10(6, SI):3834-3843, JUN 2010.

Despite the continuing research interests in CNT-liquid systems, the microscopic structure and transport behavior of liquid mixtures in carbon nanotubes (CNTs) remain poorly understood. Methanol and water liquids are completely miscible across the entire range of concentration; however, recent research reveals that they are immiscible at a molecular level. In this work, we carried out classical molecular dynamics to study the molecular distribution, structure ordering, clustering and transport behavior of liquid methanol-water mixtures within CNT confinement. We found that CNTs preferentially absorbed methanol over water molecule even though the latter has a smaller molecular size, indicating that chemical effect such as molecular hydrophilicity plays a crucial role in the molecular absorption of CNTs. Due to the selective absorption of CNTs, methanol aqueous solution changes from microscopically immiscible to macroscopically immiscible at nanoscale. This nanoscale immiscibility may be utilized in various applications of CNTs including direct methanol fuel cells, nanosensors, molecular sieves, nanofluidic chips, and capsules for drug delivery.

[12] Styliani Consta. Manifestation of Rayleigh Instability in Droplets Containing Multiply Charged Macroions. JOURNAL OF PHYSICAL CHEMISTRY B, 114(16):5263-5268, APR 29 2010.

The Rayleigh limit and manifestations of instability in liquid droplets containing charged macroions are examined by molecular simulations. It is found that beyond the Rayleigh limit, the spherical droplets become unstable and form structures with distinct features. Regardless of the nature of the charged macroion, an assembly of spines of highly ordered polar solvent molecules form on the droplet surface. The surface charge distribution of the spiny droplet is highly nonuniform, and the macroscopic description of the droplet energy as a sum of electrostatic and surface terms is no longer valid. When the macroion is a charge-saturated polyhistidine chain, it is shown that the changes in the structure of the droplet are accompanied by the chain extension. Contrary to the conventional point of view, it is found that single ions present in droplets containing a highly charged macroion do not escape spontaneously but rather form complexes stable on the nanosecond time scale, depending on the degree of deviation from the Rayleigh limit and the nature of the ion. The effect of the instability in the disintegration mechanism of charged droplets in electrospray mass spectrometry experiments is discussed.

[13] Yi Liu, Styliani Consta, Yujun Shi, R. H. Lipson, and William A. Goddard, III. Prediction of the Size Distributions of Methanol-Ethanol Clusters Detected in VUV Laser/Time-of-Flight Mass Spectrometry. JOURNAL OF PHYSICAL CHEMISTRY A, 113(25):6865-6875, JUN 25 2009.

The size distributions and geometries of vapor clusters equilibrated with methanol-ethanol (Me-Et) liquid mixtures were recently studied by vacuum ultraviolet (VUV) laser time-of-flight (TOF) mass spectrometry and density functional theory (DFT) calculations (Liu, Y.; Consta, S.; Ogeer, F.; Shi, Y. J.; Lipson, R. H. Can. J. Client. 2007, 85, 843-852). On the basis of the mass spectra recorded, it was concluded that the formation of neutral tetramers is particularly prominent. Here we develop grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) frameworks to compute cluster size distributions in vapor mixtures that allow a direct comparison with experimental mass spectra. Using the all-atom optimized potential for liquid simulations (OPLS-AA) force field, we systematically examined the neutral cluster size distributions as functions of pressure and temperature. These neutral cluster distributions were then used to derive ionized cluster distributions to compare directly with the experiments. The simulations suggest that supersaturation at 12 to 16 times the equilibrium vapor pressure at 298 K or supercooling at temperature 240 to 260 K at the equilibrium vapor pressure can lead to the relatively abundant tetramer population observed in the experiments. Our simulations capture the most distinct features observed in the experimental TOF mass spectra: Et3H+ at m/z = 139 in the vapor corresponding to 10:90% Me-Et liquid mixture and Me3H+ at m/z = 97 in the vapors corresponding to 50:50% and 90: 10% Me-Et liquid mixtures. The hybrid GCMC scheme developed in this work extends the capability of studying the size distributions of neat clusters to mixed species and provides a useful tool for studying environmentally important systems such as atmospheric aerosols.

[14] Y. Liu, S. Consta, F. Ogeer, Y. J. Shi, and R. H. Lipson. Geometries and energetics of methanol-ethanol clusters: a VUV laser/time-of-flight mass spectrometry and density functional theory study. CANADIAN JOURNAL OF CHEMISTRY-REVUE CANADIENNE DE CHIMIE, 85(10):843-852, OCT 2007.

Hydrogen-bonded clusters, formed above liquid methanol (Me) and ethanol (Et) mixtures of various compositions, were entrained in a supersonic jet and probed using 118 nm vacuum ultraviolet (VUV) laser single-photon ionization/time-of-flight mass spectrometry. The spectra are dominated by protonated cluster ions, formed by ionizing hydrogen-bonded MemEtn neutrals, m = 0-4, n = 0-3, and m + n = 2-5. The structures and energetics of the neutral and ionic species were investigated using both the all-atom optimized potential for liquid state, OPLS-AA, and the density functional (DFT) calculations. The energetic factors affecting the observed cluster distributions were examined. Calculations indicate that the large change in binding energy going from trimer to tetramer can be attributed more to pair-wise interactions than to cooperativity effects.

[15] Kengo Ichiki and Styliani Consta. Disintegration mechanisms of charged aqueous nanodroplets studied by simulations and analytical models. JOURNAL OF PHYSICAL CHEMISTRY B, 110(39):19168-19175, OCT 5 2006.

The mechanism of fragmentation processes in aqueous nanodroplets charged with ions is studied by molecular dynamics (MD) simulations. By using constant-temperature MD, the evaporation of the water is naturally taken into account and sequences of ion fragmentation events are observed. The size of the critical radius of the charged droplet just before the fragmentation and the distribution of the sizes of the fragments are estimated. Comparison of the Rayleigh critical radius for fragmentation and simulation data is within 0.23 nm. This seemingly small difference arises from a large difference in the number of water molecules that makes fragmentation an activated process as in the ion evaporation mechanism (IEM). This finding is in agreement with the predictions of Labowsky et al. [Anal. Chim. Acta 2000, 406, 105-118] for charged aqueous drops. The size of the daughter droplets is larger than the prediction of Born's theory by 0.1 to 0.15 nm. The nature and the dynamics of the intermediate states of the fragmentation process characterized by a bridge formed between the mother droplet and the evaporating ion or thorned structures where the ion sits on the tip are important for the outcome of the size-distribution of the fragments, while they are is missing in Born's theory.

[16] Styliani Consta. Detecting reaction pathways and computing reaction rates in condensed phase. THEORETICAL CHEMISTRY ACCOUNTS, 116(1-3):373-382, AUG 2006.

Methods for the computation of rate constants that characterize classical reactions occurring in the condensed phase are discussed. While microscopic expressions for these transport properties are well known, their computation presents challenges for simulation since reactive events often occur rarely, and the long time scales that are typical for reactive processes are not accessible using simple molecular dynamics methods. Furthermore, the underlying free energy surface is very complex with many saddle points that prevent sampling of possible reaction pathways. As a result, the reaction coordinate may be a complex many-body function of the system's degrees of freedom. Since there is not an a priori way to define a “good†reaction coordinate, methods are being developed to assist in a systematic construction of a reaction coordinate. These methods are reviewed and examples of non-trivial reaction coordinates are presented.

[17] M. Paliy, O. M. Braun, and S. Consta. The friction properties of an ultrathin confined water film. TRIBOLOGY LETTERS, 23(1):7-14, JUL 2006.

An ultrathin water film confined between two substrates in moving contact is studied using Langevin molecular dynamics with coordinate- and velocity-dependent damping coefficient. The correlations between the structure of the water lubricant film and its frictional properties are found in a wide range of applied loads and for various strengths of interactions between water and surfaces. A self-organization of the film into a low-friction state under driving is observed to occur. Surprisingly, it is found that the “hydrophilic†surfaces exhibit a lower friction than the “hydrophobic†ones. The viscosity of the confined water estimated from the simulations is in a good agreement with the experimental value.

[18] S. Consta. Detecting reaction pathways and computing reaction rates in condensed phase. In Simos, T and Maroulis, G, editor, ADVANCES IN COMPUTATIONAL METHODS IN SCIENCES AND ENGINEERING 2005, VOLS 4 A & 4 B, volume 4A-4B of LECTURE SERIES ON COMPUTER AND COMPUTATIONAL SCIENCES, pages 709-712. Amer Chem Soc; Amer Phys Soc, 2005. International Conference on Computational Methods in Sciences and Engineering (ICCMSE 2005), Corinth, GREECE, OCT 21-26, 2005.

In this contribution, methods for the computation of rate constants that characterize classical reactions occurring in the condensed phase axe discussed. Methods of transition path sampling and minimal action are reviewed and their adaption in the study on physical fragmentation of charged aqueous clusters is discussed.

[19] OM Braun, M Paliy, and S Consta. Ordering of a thin lubricant film due to sliding. PHYSICAL REVIEW LETTERS, 92(25), JUN 25 2004.

A thin lubricant film confined between two substrates in moving contact is studied using Langevin molecular dynamics with the coordinate- and velocity-dependent damping coefficient. It is shown that an optimal choice of the interaction within the lubricant can lead to minimal kinetic friction as well as to low critical velocity of the stick-slip to smooth-sliding transition. The strength of this interaction should be high enough (relative to the strength of the interaction of lubricant atoms with the substrates) so that the lubricant remains in a solid state during sliding. At the same time, the strength of the interaction should not be too high, in order to allow annealing of defects in the lubricant at slips.

[20] S Consta, KR Mainer, and W Novak. Fragmentation mechanisms of aqueous clusters charged with ions. JOURNAL OF CHEMICAL PHYSICS, 119(19):10125-10132, NOV 15 2003.

Fragmentation processes of mesoscopic aqueous clusters charged with ions of similar sign are studied by computer simulations. In order to examine differences in the fragmentation that depend on the nature of the ions and the charge distribution, the clusters contain positive (Na+ and Ca2+) or negative (Cl-) ions. Insight into the fragmentation mechanism is obtained by theories of activated processes. Critical to this approach is the use of a new reaction coordinate that captures the shape fluctuations of the droplet that are responsible for the reaction. Reversible work profiles for the reaction are constructed along the reaction coordinate, and dynamics is performed. The dynamics validates the use of the reaction coordinate, and shows diffusive barrier crossing. It is found that clusters with even number of charges fragment unevenly in contrast to analytical theories that predict even fission by considering only the energetic factors that determine the stability of charged droplets. (C) 2003 American Institute of Physics.

[21] GS Fanourgakis, YJ Shi, S Consta, and RH Lipson. A spectroscopic and computer simulation study of butanol vapors. JOURNAL OF CHEMICAL PHYSICS, 119(13):6597-6608, OCT 1 2003.

Clusters of butanol formed above neat liquid samples were entrained in a supersonic jet and probed using 10.5 eV vacuum ultraviolet laser single-photon ionization/time-of-flight mass spectrometry. The four different isomers of butanol (n-butanol, sec-butanol, iso-butanol, and tert-butanol) were studied separately to assess the influence of the structure of the alkyl chain on the formation and stability of the hydrogen bonded clusters. Most of the higher mass features observed in the mass spectra could be assigned to protonated alcohol clusters, H(ROH)(n)(+), nless than or equal to3; R=C4H9, that arise from facile proton-alkoxy radical/alkoxide anion dissociation. Signals due to protonated trimers were only evident in the spectra of tert- and sec-butanol. Empirical force fields, density functional theory and ab initio methods were used to identify the geometries of all clusters up to the pentamers for the different isomers. Monte Carlo simulations established vapor-phase cluster distributions, while molecular dynamics was used to assess the relative stability of the isomeric tetramers. Together, these experimental and theoretical results suggest that butanol tetramers are “magic-number†structures, and that the protonated ion signals of size n could be correlated with the neutral cluster of size n+1, provided the vapor pressures sampled in the supersonic jet exceeded equilibrium values. (C) 2003 American Institute of Physics.

[22] S Consta. Fragmentation reactions of charged aqueous clusters. JOURNAL OF MOLECULAR STRUCTURE-THEOCHEM, 591:131-140, AUG 30 2002.

Fragmentation reaction of mesoscopic aqueous clusters loaded with several sodium ions is investigated by computer simulations. The ratio of the number of ions over that of the solvent molecules is such that the clusters are stable for several nanoseconds. Fragmentation is caused by infrequent shape fluctuations. A new reaction coordinate is introduced that captures the shape changes that determine the reaction mechanism and allows for the study of the reaction using theories of activated processes. The nature of the shape fluctuations responsible for the fragmentation is characterized by free energy profiles computed as function of the new reaction coordinate. Dynamics of the fragmentation shows that the barrier crossing is diffusive so that the dynamical corrections to the transition state are large. The structure of the decay of the time-dependent rate constant reflects the diffusive character of the recrossing dynamics so that a plateau is established after a long transient time of 5 ps. The free energy and dynamics of the reaction demonstrate that clusters fragment unevenly in contrast to predictions of analytical theories. (C) 2002 Elsevier Science B.V. All rights reserved.

[23] YJ Shi, S Consta, AK Das, B Mallik, D Lacey, and RH Lipson. A 118 nm vacuum ultraviolet laser/time-of-flight mass spectroscopic study of methanol and ethanol clusters in the vapor phase. JOURNAL OF CHEMICAL PHYSICS, 116(16):6990-6999, APR 22 2002.

Clusters of methanol and ethanol formed above neat liquid samples were entrained in a supersonic jet of helium and probed in the expansion using 118 nm vacuum ultraviolet laser single-photon ionization/time-of-flight (TOF) mass spectrometry. Almost every cluster ion observed in the TOF mass spectra could be represented by the formula H(ROH)(n)(+), where R=CH3 or C2H5, and n=1-5. Formation of these species is attributed to a well-established ionization pathway where each protonated (n-1)-mer originates from its n-mer neutral parent. Signals in the TOF mass spectra due to the protonated trimers H(CH3OH)(3)(+) and H(CH3CH2OH)(3)(+) were found to be the most intense and provides direct evidence that these particular cluster ions are “magic-number†structures. The possible relationships between the observed ion data and the neutral cluster vapor phase distributions are discussed. In this context, methanol and ethanol vapor cluster distributions at 298.15 K and at several pressuresgreater than or equal tothe equilibrium vapor pressure were computed using the grand canonical Monte Carlo and molecular dynamics techniques. Lastly, differences between these experiments and the results of bimolecular reaction studies are discussed. (C) 2002 American Institute of Physics.

[24] S Consta, TJH Vlugt, JW Hoeth, B Smit, and D Frenkel. Recoil growth algorithm for chain molecules with continuous interactions. MOLECULAR PHYSICS, 97(12):1243-1254, DEC 20 1999.

The recoil growth (RG) scheme is a dynamic Monte Carlo algorithm that has been suggested as an improvement over the configurational bias Monte Carlo (CBMC) method (Consta, S., Wilding, N. B., Frenkel, D. and Alexandrowicz, Z., 1999, J. chem. Phys., 110, 3220). The RG method had originally been tested for hard core polymers on a lattice, and it was found that RG outperforms CBMC for dense systems and long chain molecules. In the present paper, the RG scheme is extended to the practically more relevant case of off-lattice chain molecules with continuous interactions. It is found that for longer chain molecules RG becomes over an order of magnitude more efficient than CBMC. However, other schemes are better suited to the computation of the excess chemical potential. Moreover, it is more difficult to parallelize RG than CBMC.

[25] S Consta and R Kapral. Ionization reactions of ion complexes in mesoscopic water clusters. JOURNAL OF CHEMICAL PHYSICS, 111(22):10183-10191, DEC 8 1999.

The free energy and dynamics of the dissociation reactions of the [Na+(Cl-)(2)] ion complex in mesoscopic water clusters are examined. The free energy surface shows the existence of stable single and double solvent-separated complex species formed from ionization of the stable double-contact ion complex. The reaction occurs on the cluster surface for mesoscale clusters composed of tens of water molecules. Passage between stable species is an activated process but barrier crossing has a large diffusive component so that dynamical corrections to transition state theory are large. The structure of the decay of the time-dependent rate constant reflects the diffusive character of the recrossing dynamics so that a plateau is not established on a 10 ps time scale in contrast to ionization dynamics in bulk fluids. (C) 1999 American Institute of Physics. [S0021-9606(99)51941-7].

[26] S Consta, NB Wilding, D Frenkel, and Z Alexandrowicz. Recoil growth: An efficient simulation method for multi-polymer systems. JOURNAL OF CHEMICAL PHYSICS, 110(6):3220-3228, FEB 8 1999.

We present a new Monte Carlo scheme for the efficient simulation of multi-polymer systems. The method permits chains to be inserted into the system using a biased growth technique. The growth proceeds via the use of a retractable feeler, which probes possible pathways ahead of the growing chain. By recoiling from traps and excessively dense regions, the growth process yields high success rates for both chain construction and acceptance. Extensive tests of the method using self-avoiding walks on a cubic lattice show that for long chains and at high densities it is considerably more efficient than configurational bias Monte Carlo, of which it may be considered a generalization. (C) 1999 American Institute of Physics. [S0021-9606(99)51406-2].

[27] R Kapral, S Consta, and D Laria. 1996 Polanyi Award Lecture Proton reactions in clusters. CANADIAN JOURNAL OF CHEMISTRY-REVUE CANADIENNE DE CHIMIE, 75(1):1-8, JAN 1997. 79th Annual Conference of the Canadian-Society-for-Chemistry, ST JOHNS, CANADA, JUN 23-26, 1996.

Reactions in mesoscopic, molecular clusters may proceed by mechanisms and with rates that differ from those in bulk solvents. Two examples of reactions in large, liquid-state, molecular clusters are described to illustrate the distinctive features of these reactions: acid dissociation and proton transfer in aprotic, polar solvents. Both of these reactions involve proton dynamics so methods for dealing with mixed quantum-classical systems must be utilized to investigate the reaction dynamics. Surface versus bulk solvation effects play an important role in determining the reaction mechanisms as do the strong cluster fluctuations. Mechanisms for proton transfer within clusters that have no bulk analogs will be described.

[28] S Consta and R Kapral. Dynamics of proton transfer in mesoscopic clusters. JOURNAL OF CHEMICAL PHYSICS, 104(12):4581-4590, MAR 22 1996.

Proton transfer rates and mechanisms are studied in mesoscopic, liquid-state, molecular clusters. The proton transfer occurs in a proton-ion complex solvated by polar molecules comprising the cluster environment. The rates and mechanisms of the reaction are studied using both adiabatic and nonadiabatic molecular dynamics. For large molecular clusters, the proton-ion complex resides primarily on the surface of the cluster or one layer of solvent molecules inside the surface. The proton transfer occurs as the complex undergoes orientational fluctuations on the cluster surface or penetrates one solvent layer into the cluster leading to solvent configurations that favor the transfer. For smaller clusters the complex resides mostly on the surface of the cluster and proton transfer is observed only when the complex penetrates the cluster and solvent configurations that favor the proton transfer are achieved. Quantitative information on the cluster reaction rate constants is also presented. (C) 1996 American Institute of Physics.

[29] S CONSTA and R KAPRAL. PROTON-TRANSFER IN MESOSCOPIC, MOLECULAR CLUSTERS. JOURNAL OF CHEMICAL PHYSICS, 101(12):10908-10914, DEC 15 1994.
[30] ID PETSALAKIS, G THEODORAKOPOULOS, and S CONSTA. THEORETICAL CALCULATIONS ON THE ELECTRONIC STATES OF XEH. MOLECULAR PHYSICS, 75(4):805-810, MAR 1992.

Configuration interaction calculations have been carried out on the ground and ten excited states of XeH. The results are in agreement with the existing experimental data on the Rydberg spectra of this system. In addition, secondary minima are found in the first two excited states (A 2-SIGMA+ and B'2 PI) at internuclear distance of 5.0 bohr, where the two states have charge-transfer character, in agreement with recent experimental findings. The present calculations are not of sufficient accuracy to reproduce the Van der Waals minimum of the ground state.