Search

The solubility of a nonpolar solute in water is changed upon addition of a salt or cosolvent. Hereafter, "solvent" is formed by water molecules for pure water, by water molecules, cations, and anions for water-salt solution, and by water and cosolvent molecules for water-cosolvent solution. Decrease and increase in the solubility, respectively, are ascribed to enhancement and reduction of the hydrophobic effect. Plenty of experimental data are available for the change in solubility of argon or methane arising from the addition. We show that the integral equation theory combined with a rigid-body model, in which the solute and solvent particles are modeled as hard spheres with different diameters, can reproduce the data for the following items: salting out by an alkali halide and salting in by tetramethylammonium bromide, increase in solubility by a monohydric alcohol, and decrease in solubility by sucrose or urea. The orders of cation or anion species in terms of the power of decreasing the solubility can also be reproduced for alkali halides. With the rigid-body model, the analyses are focused on the roles of entropy originating from the translational displacement of solvent particles. It is argued by decomposing the solvation entropy of a nonpolar solute into physically insightful constituents that the solvent crowding in the bulk is a pivotal factor of the hydrophobic effect: When the solvent crowding in the bulk becomes more serious, the effect is strengthened, and when it becomes less serious, the effect is weakened. It is experimentally known that the thermal stability of a protein is also influenced by the salt or cosolvent addition. The additions which decrease and increase the solubility of a nonpolar solute, respectively, usually enhance and lower the thermal stability. This suggests that the enhanced or reduced hydrophobic effect is also a principal factor governing the stability change. However, urea decreases the solubility but lowers the stability. Bromide and iodide ions decrease the solubility but lower the stability of a protein with a large, positive total charge. In these cases, the urea- or ion-protein van der Waals interaction energy as well as the hydrophobic effect needs to be taken into account in arguing the stability change. We also present a new view on the so-called Hofmeister series: We show how it is expressed when the change in hydrophobic effect dominates and how it is modified when other factors are also influential.

At iThemba Laboratory for Accelerator Based Sciences (iThemba LABS) an electron cyclotron ion source was installed and commissioned. This source is a copy of the Grenoble Test Source (GTS) for the production of highly charged ions. The source is similar to the GTS-LHC at CERN and named GTS2. A collaboration between the Accelerators and Beam Physics Group of CERN and the Accelerator and Engineering Department of iThemba LABS was proposed in which the development of high intensity argon and xenon beams is envisaged. In this paper, we present beam experiments with the GTS2 at iThemba LABS, in which the results of continuous wave and afterglow operation of xenon ion beams with oxygen as supporting gases are presented.

4-propylaniline complexes with one and two argon atoms formed in the molecular beam were studied in the first excited electronic state, S1, using resonance enhanced two-photon ionization spectroscopy and in the cation ground state, D0, using mass analyzed threshold ionization spectroscopy. The combination of electronic and cationic spectra of the clusters allows two conformations to be identified in both aniline-Ar1 and aniline-Ar2, which are assigned to either the gauche configuration or anti-configuration of 4-propylaniline. The gauche isomer exhibits complex bands shifted 29 cm(-1) and 89 cm(-1) from the S1 origin bands and 83 cm(-1) and 148 cm(-1) from the ionization potential assigned to the Ar1 and Ar2 complexes, respectively. For the anti-rotamer, the corresponding shifts actually become nearly additive, 53 cm(-1) and 109 cm(-1) for the S1 origin bands, and 61 cm(-1) and 125 cm(-1) for the ionization potentials. Ab initio calculations provide insights into the influences of the propyl and amino groups on the positions of the argon atoms within the clusters. In addition, the binding energy of one argon with the gauche isomer of 4-propylaniline has been measured to be 550 ± 5 cm(-1) in the D0 state, 496 ± 5 cm(-1) in the S1 state, and 467 ± 5 cm(-1) in the neutral ground state, S0.

We present a continuum-based approach to predict the structure and thermodynamic properties of confined fluids at multiple length-scales, ranging from a few angstroms to macro-meters. The continuum approach is based on the empirical potential-based quasi-continuum theory (EQT) and classical density functional theory (cDFT). EQT is a simple and fast approach to predict inhomogeneous density and potential profiles of confined fluids. We use EQT potentials to construct a grand potential functional for cDFT. The EQT-cDFT-based grand potential can be used to predict various thermodynamic properties of confined fluids. In this work, we demonstrate the EQT-cDFT approach by simulating Lennard-Jones fluids, namely, methane and argon, confined inside slit-like channels of graphene. We show that the EQT-cDFT can accurately predict the structure and thermodynamic properties, such as density profiles, adsorption, local pressure tensor, surface tension, and solvation force, of confined fluids as compared to the molecular dynamics simulation results.

Tunneling corrections to Landau-Zener rate coefficients for the vibrational relaxation NO(X(2)Π, v = 1) + Ar → NO(X(2)Π, v = 0) + Ar between 300 and 2000 K are determined employing ab initio potential energy surfaces calculated by the code provided by Alexander [J. Chem. Phys. 111, 7426 (1999)]. The calculations use a reaction coordinate approach and lead to vibronically nonadiabatic transition probabilities within the generalized Airy approximation as extended to the WKB underbarrier Landau-Lifshitz limit. The calculations confirm experimental evidence for an onset of major tunneling contributions to the relaxation rate at temperatures below about 900 K and rationalize large tunneling contributions at 300 K. These effects increase the rate coefficients by several orders of magnitude over the uncorrected Landau-Zener values and remove the large gap between the latter and experimental results.

Momentum-transfer collision cross-sections and integral collision cross-sections for the collision-induced dissociation are calculated for collisions of ionized argon dimers with argon atoms using a nonadiabatic semiclassical method with the electronic Hamiltonian calculated on the fly via a diatomics-in-molecules semiempirical model as well as inverse-method modeling based on simple isotropic rigid-core potential. The collision cross-sections are then used in an optimized Monte Carlo code for evaluations of the Ar 2 (+) mobility in argon gas, longitudinal diffusion coefficient, and collision-induced dissociation rates. A thorough comparison of various theoretical calculations as well as with available experimental data on the Ar 2 (+) mobility and collision cross-sections is performed. Good agreement is found between both theoretical approaches and the experiment. Analysis of the role of inelastic processes in Ar 2 (+)/Ar collisions is also provided.

We report a new three-dimensional potential energy surface for Ar-CS2 involving the Q3 normal mode for the υ3 antisymmetric stretching vibration of the CS2 molecule. The potential energies were calculated using the supermolecular method at the coupled-cluster singles and doubles level with noniterative inclusion of connected triples, using augmented correlation-consistent quadruple-zeta basis set plus midpoint bond functions. Two vibrationally averaged potentials with CS2 at both the ground (υ = 0) and the first excited (υ = 1)υ3 vibrational states were generated from the integration of the three-dimensional potential over the Q3 coordinate. Each potential was found to have a T-shaped global minimum and two equivalent linear local minima. The radial discrete variable representation /angular finite basis representation method and the Lanczos algorithm were applied to calculate the rovibrational energy levels. The calculated band origin shift of the complex (0.0622 cm(-1)) is very close to the observed one (0.0671 cm(-1)). The predicted infrared spectra and spectroscopic parameters based on the two averaged potentials are in excellent agreement with the available experimental data.

Theoretical studies of the potential energy surfaces (PESs) and bound states are performed for rare gas-carbon disulfide complexes, He-CS2, Ne-CS2, and Ar-CS2. Three two-dimensional intermolecular PESs are constructed from ab initio data points which are calculated at the CCSD(T) level with aug-cc-pVTZ basis set supplemented with bond functions. We find that the three PESs have very similar features and each PES can be characterized by a global T-shaped minimum, two equivalent local linear minima, and the saddle points between them. The T-shaped isomer is energetically more stable than the linear isomer for each complex. The linear isomers, which have not been observed in experiment so far, are predicted from our PESs and further identified by bound state calculations. Moreover, we assign several intermolecular vibrational states for both the T-shaped and linear isomers of the three complexes via the analysis of wavefunctions. The corresponding vibrational frequencies are calculated from the bound state energies for these assigned states. These frequencies could be helpful for further experimental studies, especially for the linear isomers. We also calculate the rovibrational transition frequencies for the three T-shaped isomers and the pure rotational transition frequencies for the linear isomers, respectively. The accuracy of the PESs is validated by the good agreement between theoretical and experimental results for the rovibrational transition frequencies and spectroscopic parameters.

Combining rules, such as the Lorentz-Berthelot rules, are routinely used to calculate the thermodynamic properties of mixtures using molecular simulations. Here we extend the expanded Wang-Landau simulation approach to determine the impact of the combining rules on the value of the partition function of binary systems, and, in turn, on the phase coexistence and thermodynamics of these mixtures. We study various types of mixtures, ranging from systems of rare gases to biologically and technologically relevant mixtures, such as water-urea and water-carbon dioxide. Comparing the simulation results to the experimental data on mixtures of rare gases allows us to rank the performance of combining rules. We find that the widely used Lorentz-Berthelot rules exhibit the largest deviations from the experimental data, both for the bulk and at coexistence, while the Kong and Waldman-Hagler provide much better alternatives. In particular, in the case of aqueous solutions of urea, we show that the use of the Lorentz-Berthelot rules has a strong impact on the Gibbs free energy of the solute, overshooting the value predicted by the Waldman-Hagler rules by 7%. This result emphasizes the importance of the combining rule for the determination of hydration free energies using molecular simulations.

A matrix isolation study of the infrared spectra and structure of anethole (1-methoxy-4-(1-propenyl)benzene) has been carried out, showing the presence of two E conformers (AE1, AE2) of the molecule in the as-deposited matrices. Irradiation using ultraviolet-tunable laser light at 308-307 nm induced conformationally selective phototransformations of these forms into two less stable Z conformers (AZ1, AZ2). The back reactions were also detected upon irradiation at 301 nm. On the whole, the obtained results allow for full assignment of the infrared spectra of all the four experimentally observed anethole isomers and showed that the narrowband UV-induced E-Z photoisomerization is an efficient and selective way to interconvert the two isomers of anethole into each other, with conformational discrimination. Photolysis of anethole was observed as well, with initial methoxyl O-C bond cleavage and formation of CH3 and p-propenylphenoxy (AR) radicals, followed by radical recombination to form 2-methyl-4-propenyl-2,4-cyclohexadienone, which subsequently undergoes ring-opening generating several conformers of long-chain conjugated ketenes. Interpretation of the experimental observations was supported by density functional theory (B3LYP and B2PLYD) calculations.

Rotational diffusion processes are correlated with nanoparticle visualization and manipulation techniques, widely used in nanocomposites, nanofluids, bioscience, and so on. However, a systematical methodology of deriving this diffusivity is still lacking. In the current work, three molecular dynamics (MD) schemes, including equilibrium (Green-Kubo formula and Einstein relation) and nonequilibrium (Einstein-Smoluchowski relation) methods, are developed to calculate the rotational diffusion coefficient, taking a single rigid carbon nanotube in fluid argon as a case. We can conclude that the three methods produce same results on the basis of plenty of data with variation of the calculation parameters (tube length, diameter, fluid temperature, density, and viscosity), indicative of the validity and accuracy of the MD simulations. However, these results have a non-negligible deviation from the theoretical predictions of Tirado et al. [J. Chem. Phys. 81, 2047 (1984)], which may come from several unrevealed factors of the theory. The three MD methods proposed in this paper can also be applied to other situations of calculating rotational diffusion coefficient.

The two photon resonant, three photon ionization spectra of ArXe were recorded in the spectral region of 88,500-90,100 cm(-1). Seven new molecular band progressions dissociating to ArXe* → Ar(1)S0 + Xe* 7p[1/2]0, Xe* 7p[3/2]2, Xe* 6p'[3/2]2, Xe* 6p'[1/2]1, Xe* 6p'[1/2]0 have been selected and analyzed. The molecular constants for the excited states of ArXe* of these vibrational progressions were determined in the approximation of the anharmonic oscillator, the Morse potential and the Franck-Condon principle. The photoelectron spectra were recorded by several excited electronic-vibrational transitions of ArXe, the dissociation channels of the excited molecules were determined and extra information about the electron structure of the excited molecular states was obtained.

A high fidelity computational fluid dynamic model is used to simulate the flow, pressure, and density fields generated in a cylindrical and a conical resonator by a vibrating end wall/piston producing high-amplitude standing waves. The waves in the conical resonator are found to be shock-less and can generate peak acoustic overpressures that exceed the initial undisturbed pressure by two to three times. A cylindrical (consonant) acoustic resonator has limitations to the output response observed at one end when the opposite end is acoustically excited. In the conical geometry (dissonant acoustic resonator) the linear acoustic input is converted to high energy un-shocked nonlinear acoustic output. The model is validated using past numerical results of standing waves in cylindrical resonators. The nonlinear nature of the harmonic response in the conical resonator system is further investigated for two different working fluids (carbon dioxide and argon) operating at various values of piston amplitude. The high amplitude nonlinear oscillations observed in the conical resonator can potentially enhance the performance of pulse tube thermoacoustic refrigerators and these conical resonators can be used as efficient mixers.

We present the experimental and simulated (2+1) REMPI spectrum of the C(2)Π state of the NO-Ar complex, in the vicinity of the 3p Rydberg state of NO. Two Rydberg states of NO are expected in this energy region: the C(2)Π (3pπ) and D(2)Σ(+) (3pσ) states, and we concentrate on the former here. When the C(2)Π (3pπ) state interacts with Ar at nonlinear orientations, the symmetry is lowered to C(s), splitting the degeneracy of the (2)Π state to yield C((2)A") and C((2)A') states. For these two states of NO-Ar, we calculate potential energy surfaces using second order Møller-Plesset perturbation theory, exploiting a procedure to converge the reference Hartree-Fock wavefunction to describe the excited states, the maximum overlap method. The bound rovibrational states obtained from the surfaces are used to simulate the electronic spectrum, which is in excellent agreement with experiment, providing assignments for the observed spectral lines from the calculated rovibrational wavefunctions.

Spin-adaptation of virtual functions in state-specific multireference perturbation theory is examined. Redundancy occurring among virtual functions generated by unitary group based excitation operators on a model-space function is handled by canonical orthogonalization. The treatment is found to remove non-physical kinks observed earlier on potential energy surfaces. Sensitivity analysis of the new approach confirms the elimination of the drastic increase in singular values of sensitivity matrices, reported earlier.

Pure acetylene and mixed Ar-acetylene clusters are formed in supersonic expansions of acetylene/argon mixtures and analysed using reflectron time-of-flight mass spectrometer with variable electron energy ionization source. Acetylene clusters composed of more than a hundred acetylene molecules are generated at the acetylene concentration of ≈8%, while mixed species are produced at low concentrations (≈0.7%). The electron energy dependence of the mass spectra revealed the ionization process mechanisms in clusters. The ionization above the threshold for acetylene molecule of 11.5 eV results in the main ionic fragment progression (C2H2)n(+). At the electron energies ≥21.5 eV above the CH+CH(+) dissociative ionization limit of acetylene the fragment ions nominally labelled as (C2H2)nCH(+), n ≥ 2, are observed. For n ≤ 7 these fragments correspond to covalently bound ionic structures as suggested by the observed strong dehydrogenation [(C2H2)n - k × H](+) and [(C2H2)nCH - k × H](+). The dehydrogenation is significantly reduced in the mixed clusters where evaporation of Ar instead of hydrogen can stabilize the nascent molecular ion. The C3H3(+) ion was previously assigned to originate from the benzene molecular ion; however, the low appearance energy of ≈13.7 eV indicates that a less rigid covalently bound structure of C6H6(+) ion must also be formed upon the acetylene cluster electron ionization. The appearance energy of Arn(C2H2)(+) fragments above ≈15.1 eV indicates that the argon ionization is the first step in the fragment ion production, and the appearance energy of Arn≥2(C2H2)m≥2(+) at ≈13.7 eV is discussed in terms of an exciton transfer mechanism.

The investigation of the chemical complexity of hot environments, ranging from combustion flames to circumstellar envelopes of evolved stars, relies on the determination of the reaction kinetics and product branching ratio. We have designed a chemical reactor for the exploration of high temperature chemistry. This apparatus is employed in the present study to measure the reaction kinetics of the CN radical with C2H4 and C2H6 over the 300-1200 K temperature range. In our setup and in some environments, the CN radical is partially produced in a vibrationally excited state, before relaxing by collision with the surrounding gas. We complement the experimental kinetic studies of hydrocarbons reactions with CN(v = 0) with a theoretical study of vibrational relaxation of CN(v = 1) by He and Ar atoms, the main collisional partners in our apparatus. Calculations are carried out to determine the collisional elastic and inelastic cross sections versus the kinetic energy as well as the corresponding vibrationally elastic and inelastic rate coefficients. The results are compared with empirical calculations and with a few experimental observations. The range of validity of the empirical model is discussed and potential applications sketched.

Three rovibrational bands of Ar-D2O and two rovibrational bands of Ar-HDO were observed in the 1.5 μm range by continuous wave cavity ringdown spectroscopy. Their analyses led to the determination of rotational constants for the upper states and vibrational shifts indicating that the potential energy surface is only slightly affected by the vibrational excitation. Some Coriolis couplings were identified. The observed lines were fitted to retrieve a 3.5 ns lifetime of the upper state, showing that even with a triple or double excitation quanta in the water subunit, the Ar-D2O and Ar-HDO complexes are long-lived species.

An explicit polarizable potential for the naphthalene-argon complex has been derived assuming only atomic contributions, aiming at large scale simulations of naphthalene under argon environment. The potential was parametrized from dedicated quantum chemical calculations at the CCSD(T) level, and satisfactorily reproduces available structural and energetic properties. Combining this potential with a tight-binding model for naphthalene, collisional energy transfer is studied by means of dedicated molecular dynamics simulations, nuclear quantum effects being accounted for in the path-integral framework. Except at low target temperature, nuclear quantum effects do not alter the average energies transferred by the collision or the collision duration. However, the distribution of energy transferred is much broader in the quantum case due to the significant zero-point energy and the higher density of states. Using an ab initio potential for the Ar-Ar interaction, the IR absorption spectrum of naphthalene solvated by argon clusters or an entire Ar matrix is computed via classical and centroid molecular dynamics. The classical spectra exhibit variations with growing argon environment that are absent from quantum spectra. This is interpreted by the greater fluxional character experienced by the argon atoms due to vibrational delocalization.

Exact relation for contributions to heat capacity of liquids is obtained from hydrodynamic theory. It is shown from analysis of the long-wavelength limit of heat density autocorrelation functions that the heat capacity of simple liquids is represented as a sum of two contributions due to "phonon-like" collective excitations and heat relaxation. The ratio of both contributions being the analogy of Landau-Placzek ratio for heat processes depends on the specific heats ratio. The theory of heat density autocorrelation functions in liquids is verified by computer simulations. Molecular dynamics simulations for six liquids having the ratio of specific heats γ in the range 1.1-2.3, were used for evaluation of the heat density autocorrelation functions and predicted Landau-Placzek ratio for heat processes. The dependence of contributions from collective excitations and heat relaxation process to specific heat on γ is shown to be in excellent agreement with the theory.

Injection of photoelectrons into gaseous or liquid dielectrics is a widely used technique to produce cold plasmas in weakly ionized systems for investigating the transport properties of electrons. We report measurements of the collection efficiency of photoelectrons injected into dense argon gas for T = 152.7 K, close to the critical temperature T(c) ≈ 150.9 K, and for T = 200.0 K. The high-field data agree with the Young-Bradbury model and with previous measurements below T(c) and at an intermediate temperature above T(c). The effective, density-dependent electron-atom momentum transfer scattering cross section can be deduced. However, the weak-field data near T(c) show large deviations from the theoretical model. We show that the electron behavior at weak field is influenced by electrostriction effects that are only important near the critical point.

The Ca(2) → Ca(4s4p(1)P) + Ca(4s(2)(1)S) photodissociation was investigated in a He droplet isolation experiment where the droplets are doped by Ar atoms. Fluorescence spectra associated with the Ca(4s4p(1)P → 4s(2)(1)S) emission were recorded as a function of the average number of Ar atoms per droplet. Three contributions were observed depending on whether the emitting Ca atoms are free, bound to helium atoms or bound to argon atoms. Moreover, the full Ca(4s4p(1)P → 4s(2)(1)S) fluorescence emission was recorded as a function of the wavelength of the photodissociation laser, hence providing the action spectrum of the Ca(2) → Ca(4s4p(1)P) + Ca(4s(2)(1)S) process. The latter spectrum suggests that in He droplets doped by argon, Ca atoms are attracted inside the droplet where they associate as Ca(2). Full analysis of the spectra indicate that the emission of Ca bound to a single Ar atom is redshifted by 94 cm(-1) with respect to the emission of free Ca.

The stability of the ground or excited state calcium atom in an argon-doped helium droplet has been investigated using an extension of the helium density functional method to treat clusters. This work was motivated by the experimental study presented in a companion paper, hereafter called Paper I [A. Masson, M. Briant, J. M. Mestdagh, M. A. Gaveau, A. Hernando, and N. Halberstadt, J. Chem. Phys. 137, 184310 (2012)], which investigated Ca(2) photodissociation in an argon-doped helium droplet and the nature of the fluorescent species. It is found that one single argon atom is sufficient to bring the calcium atom inside the droplet, for droplets of over 200 helium atoms. The absorption and emission spectra of CaAr(M) (M = 0-7) clusters have been simulated using the recently developed density sampling method to describe the influence of the helium environment. Absorption spectra exhibit broad, double bands that are significantly blueshifted with respect to the calcium atomic line. The emission spectra are less broad and redshifted with respect to the calcium resonance line. The shifts are found to be additive only for M ≤ 2, because only the first two argon atoms are located in equivalent positions around the calcium p orbital. This finding gives a justification for the fit presented in the companion paper, which uses the observed shifts in the emission spectra as a function of argon pressure to deduce the shifts as a function of the number of argon atoms present in the cluster. An analysis of this fit is presented here, based on the calculated shifts. It is concluded that the emitting species following Ca(2) photodissociation in an argon-doped droplet in Paper I could be Ca∗Ar(M) in a partly evaporated droplet where less than 200 helium atoms remain.

We report on the focusing and guiding of the van der Waals complex formed between benzonitrile molecules (C(6)H(5)CN) and argon atoms in a cold molecular beam using an ac electric quadrupole guide. The distribution of quantum states in the guided beam is non-thermal, because the transmission efficiency depends on the state-dependent effective dipole moment in the applied electric fields. At a specific ac frequency, however, the excitation spectrum can be described by a thermal distribution at a rotational temperature of 0.8 K. From the observed transmission characteristics and a combination of trajectory and Stark-energy calculations we conclude that the permanent electric dipole moment of benzonitrile remains unchanged upon the attachment of the argon atom to within ±5%. By exploiting the different dipole-moment-to-mass ([micro sign]/m) ratios of the complex and the benzonitrile monomer, transmission can be selectively suppressed for or, in the limit of 0 K rotational temperature, restricted to the complex.

We present infrared spectra of mass-selected C(10)H(8)(-)·(H(2)O)(n)·Ar(m) cluster anions (n = 1-6) obtained by Ar predissociation spectroscopy. The experimental spectra are compared with predicted spectra from density functional theory calculations. The OH groups of the water ligands are involved in H-bonds to other water molecules or to the π system of the naphthalene anion, which accommodates the excess electron. The interactions in the water network are generally found to be more important than those between water molecules and the ion. For 2 ≤ n ≤ 4 the water molecules form single layer water networks on one side of the naphthalene anion, while for n = 5 and 6, cage and multilayer structures become more energetically favorable. For cluster sizes with more than 3 water molecules, multiple conformers are likely to be responsible for the experimental spectra.

Sixteen intermolecular vibrational levels of the S(0) state of the fluorobenzene-Ar van der Waals complex have been observed using dispersed fluorescence. The levels range up to ~130 cm(-1) in vibrational energy. The vibrational energies have been modelled using a complete set of harmonic and quartic anharmonic constants and a cubic anharmonic coupling between the stretch and long axis bend overtone that becomes near ubiquitous at higher energies. The constants predict the observed band positions with a root mean square deviation of 0.04 cm(-1). The set of vibrational levels predicted by the constants, which includes unobserved bands, has been compared with the predictions of ab initio calculations, which include all vibrational levels up to 70-75 cm(-1). There are small differences in energy, particularly above 60 cm(-1), however, the main differences are in the assignments and are largely due to the limitations of assigning the ab initio wavefunctions to a simple stretch, bend, or combination when the states are mixed by the cubic anharmonic coupling. The availability of these experimental data presents an opportunity to extend ab initio calculations to higher vibrational energies to provide an assessment of the accuracy of the calculated potential surface away from the minimum. The intermolecular modes of the fluorobenzene-Ar(2) trimer complex have also been investigated by dispersed fluorescence. The dominant structure is a pair of bands with a ~35 cm(-1) displacement from the origin band. Based on the set of vibrational modes calculated from the fluorobenzene-Ar frequencies, they are assigned to a Fermi resonance between the symmetric stretch and symmetric short axis bend overtone. The analysis of this resonance provides a measurement of the coupling strength between the stretch and short axis bend overtone in the dimer, an interaction that is not directly observed. The coupling matrix elements determined for the fluorobenzene-Ar stretch-long axis bend overtone and stretch-short axis bend overtone couplings are remarkably similar (3.8 cm(-1) cf. 3.2 cm(-1)). Several weak features seen in the fluorobenzene-Ar(2) spectrum have also been assigned.

We report finite temperature quantum mechanical simulations of structural and dynamical properties of Ar(N)-CO(2) clusters using a path integral Monte Carlo algorithm. The simulations are based on a newly developed analytical Ar-CO(2) interaction potential obtained by fitting ab initio results to an anisotropic two-dimensional Morse/Long-range function. The calculated distributions of argon atoms around the CO(2) molecule in Ar(N)-CO(2) clusters with different sizes are consistent to the previous studies of the configurations of the clusters. A first-order perturbation theory is used to quantitatively predict the CO(2) vibrational frequency shift in different clusters. The first-solvation shell is completed at N = 17. Interestingly, our simulations for larger Ar(N)-CO(2) clusters showed several different structures of the argon shell around the doped CO(2) molecule. The observed two distinct peaks (2338.8 and 2344.5 cm(-1)) in the υ(3) band of CO(2) may be due to the different arrangements of argon atoms around the dopant molecule.

We evaluate the phenylacetylene-argon intermolecular potential energy surface by fitting a representative number of ab initio interaction energies to an analytic function. These energies are calculated at a grid of intermolecular geometries, using the CCSD(T) method and the aug-cc-pVDZ basis set extended with a series of 3s3p2d1f1g midbond functions. The potential is characterized by two equivalent global minima where the Ar atom is located above and below the phenylacetylene plane at a distance of 3.5781 Å from the molecular center of mass and at an angle of 9.08° with respect to the axis perpendicular to the phenylacetylene plane and containing the center of mass. The calculated interaction energy is -418.9 cm(-1). To check further the potential, we obtain the rovibrational spectrum of the complex and the results are compared to the available experimental data.

The two-electrode extractor system made of 0.1 mm diameter tungsten wires separated by 0.7 mm has formed an argon ion beam with 50 V extraction potential. Energy spreads of the extracted beams were typically less than 2 eV when the beam current density was low. The beam intensity rapidly decreased as the distance between the extractor and the beam detector increased, indicating space charge limited transport of the beam. Problems associated with the emittance measurements are also discussed.

We investigate the hydrodynamic boundary condition for simple nanofluidic systems such as argon and methane flowing in graphene nanochannels using equilibrium molecular dynamics simulations (EMD) in conjunction with our recently proposed method [J. S. Hansen, B. D. Todd, and P. J. Daivis, Phys. Rev. E 84, 016313 (2011)]. We first calculate the fluid-graphene interfacial friction coefficient, from which we can predict the slip length and the average velocity of the first fluid layer close to the wall (referred to as the slip velocity). Using direct nonequilibrium molecular dynamics simulations (NEMD) we then calculate the slip length and slip velocity from the streaming velocity profiles in Poiseuille and Couette flows. The slip lengths and slip velocities from the NEMD simulations are found to be in excellent agreement with our EMD predictions. Our EMD method therefore enables one to directly calculate this intrinsic friction coefficient between fluid and solid and the slip length for a given fluid and solid, which is otherwise tedious to calculate using direct NEMD simulations at low pressure gradients or shear rates. The advantages of the EMD method over the NEMD method to calculate the slip lengths/flow rates for nanofluidic systems are discussed, and we finally examine the dynamic behaviour of slip due to an externally applied field and shear rate., (© 2011 American Institute of Physics)

We introduce a multiscale framework to simulate inhomogeneous fluids by coarse-graining an all-atom molecular dynamics (MD) trajectory onto sequential snapshots of hydrodynamic fields. We show that the field representation of an atomistic trajectory is quantitatively described by a dynamic field-theoretic model that couples hydrodynamic fluctuations with a Ginzburg-Landau free energy. For liquid-vapor interfaces of argon and water, the parameters of the field model can be adjusted to reproduce the bulk compressibility and surface tension calculated from the positions and forces of atoms in an MD simulation. These optimized parameters also enable the field model to reproduce the static and dynamic capillary wave spectra calculated from atomistic coordinates at the liquid-vapor interface. In addition, we show that a density-dependent gradient coefficient in the Ginzburg-Landau free energy enables bulk and interfacial fluctuations to be controlled separately. For water, this additional degree of freedom is necessary to capture both the bulk compressibility and surface tension emergent from the atomistic trajectory. The proposed multiscale framework illustrates that bottom-up coarse-graining and top-down phenomenology can be integrated with quantitative consistency to simulate the interfacial fluctuations in nanoscale transport processes., (© 2011 American Institute of Physics)

High-resolution photoelectron spectra of RgS(-) (Rg = Ne, Ar, and Kr) were obtained using slow electron velocity-map imaging (SEVI). The SEVI spectra reveal well-resolved vibrational transitions between multiple spin-orbit states of RgS(-) and RgS, both of which are open-shell species. Detailed assignments are made by comparison with theoretical simulations based on high level ab initio calculations and a atoms-in-molecule model that accounts for spin-orbit coupling in the anion and neutral. Several RgS(-) and RgS vibrational frequencies and excited-state term energies are accurately determined from the analysis of the experimental spectra and are found to be in excellent agreement with the calculated values.

Short intense laser pulses of visible and infrared light can dramatically accelerate crystal nucleation from transparent solutions; previous studies invoke mechanisms that are only applicable for nucleation of ordered phases or high dielectric phases. However, we show that similar laser pulses induce CO(2) bubble nucleation in carbonated water. Additionally, in water that is cosupersaturated with argon and glycine, argon bubbles escaping from the water can induce crystal nucleation without a laser. Our findings suggest a possible link between laser-induced nucleation of bubbles and crystals.

This work derives a quantitative description of the conformational distribution in self-guided Langevin dynamics (SGLD) simulations. SGLD simulations employ guiding forces calculated from local average momentums to enhance low-frequency motion. This enhancement in low-frequency motion dramatically accelerates conformational search efficiency, but also induces certain perturbations in conformational distribution. Through the local averaging, we separate properties of molecular systems into low-frequency and high-frequency portions. The guiding force effect on the conformational distribution is quantitatively described using these low-frequency and high-frequency properties. This quantitative relation provides a way to convert between a canonical ensemble and a self-guided ensemble. Using example systems, we demonstrated how to utilize the relation to obtain canonical ensemble properties and conformational distributions from SGLD simulations. This development makes SGLD not only an efficient approach for conformational searching, but also an accurate means for conformational sampling.

The hydrophobicity of surfaces has a strong influence on their interactions with biomolecules such as proteins. Therefore, for in vitro studies of bio-surface interactions model surfaces with tailored hydrophobicity are of utmost importance. Here, we present a method for tuning the hydrophobicity of atomically flat mica surfaces by hyperthermal Ar ion irradiation. Due to the sub-100 eV energies, only negligible roughening of the surface is observed at low ion fluences and also the chemical composition of the mica crystal remains almost undisturbed. However, the ion irradiation induces the preferential removal of the outermost layer of K(+) ions from the surface, leading to the exposure of the underlying aluminosilicate sheets which feature a large number of centers for C adsorption. The irradiated surface thus exhibits an enhanced chemical reactivity toward hydrocarbons, resulting in the adsorption of a thin hydrocarbon film from the environment. Aging these surfaces under ambient conditions leads to a continuous increase of their contact angle until a fully hydrophobic surface with a contact angle >80° is obtained after a period of about 3 months. This method thus enables the fabrication of ultrasmooth biological model surfaces with precisely tailored hydrophobicity.

Electronically excited states of atoms and molecules in an environment may decay in interatomic processes by transferring excess energy to neighboring species and ionizing them. The corresponding interatomic decay width is the most important characteristic of the decay allowing to calculate its efficiency and the final states' distribution. In this paper we present calculations of interatomic widths by the Fano-Stieltjes method applied to Lanczos pseudospectra, which has been previously shown to provide accurate autoionization widths in atoms and molecules. The use of Lanczos pseudospectra allows one to avoid the full diagonalization bottleneck and makes the method applicable to larger systems. We apply the present method to the calculation of interatomic decay widths in NeMg, NeAr and HCN[middle dot]Mg(n), n = 1, 2 clusters. The results are compared with widths obtained analytically and by other ab initio methods where available., (© 2011 American Institute of Physics.)

We investigate the Ar(+)-Fe asymmetric charge transfer (ACT) reaction using a combination of plasma diagnostics methods and a kinetic model of the afterglow plasma, which allow monitoring of the temporal evolution of the densities of different species. The iron vapor is created inside a discharge cell by cathode sputtering; its density is measured by atomic absorption spectroscopy. The rate coefficient of the reaction is evaluated from the emission intensity decay of Fe(+)∗ lines pumped by the ACT process in the He-Ar-Fe and Ar-Fe afterglow plasmas. The measurements yield a rate coefficient k = 7.6( ± 3.0) × 10(-9) cm(3) s(-1) at T = 300 K.

In this work, we have tested the efficiency of two scaling approaches aiming at relating shear viscosity to a single thermodynamic quantity in dense fluids, namely the excess entropy and the thermodynamic scaling methods. Using accurate databases, we have applied these approaches first to a model fluid, the flexible Lennard-Jones chain fluid (from the monomer to the hexadecamer), then to real fluids, such as argon and normal alkanes. To enlarge noticeably the range of thermodynamics conditions for which these scaling methods are applicable, we have shown that the use of the residual viscosity instead of the total viscosity is preferable in the scaling procedures. It has been found that both approaches, using the adequate scaling, are suitable for the Lennard-Jones chain fluid model for a wide range of thermodynamic conditions whatever the chain length when scaling law exponents and prefactors are adjusted for each chain length. Furthermore, these results were found to be well respected by the corresponding real fluids.

The intermolecular potential energy surface (PES) of the naphthalene-argon (NpAr) complex is constructed using an ab initio method. The molecule-argon interaction energy is computed at the level of the second-order Møller-Plesset (MP2) theory combined with the augmented correlation consistent polarized valence double-ζ basis set. The analytical PES fitted to a large set of single energy values is further improved with the help of correction functions determined by calculations of the interaction energy at the coupled cluster level including single and double excitations supplemented by triple excitations performed for a limited set of intermolecular configurations. The PES determined is very flat near its four equivalent global minima of -493 cm(-1) located from both sides of the Np plane at a distance of 3.435 Å and shifted from the center of Np by ±0.43 Å along its long symmetry axis. The large-amplitude motion of Ar in the complex is investigated, and dynamical consequence of a strong intermode coupling is discovered in the excited vibrational states. The theoretical results obtained allow for the reassignment of the spectral bands observed in the electronic transition S(1) ← S(0) of the NpAr complex.

Hyperthermal Ar atoms were scattered under grazing incidence (θ(i) = 60°) from a CO-saturated Ru(0001) surface held at 180 K. Collision-induced desorption involving the ejection of fast CO (∼1 eV) occurs. The angularly resolved in-plane CO desorption distribution has a peak along the surface normal. However, the angular distribution varies with the fractional coverage of the surface. As the total CO coverage decreases, the instantaneous desorption maximum shifts to larger outgoing angles. The results are consistent with a CO desorption process that involves lateral interaction with neighboring molecules. Furthermore, the data indicate that the incident Ar cannot readily penetrate the saturated CO overlayer. Time-of-flight measurements of scattered Ar exhibit two components-fast and slow. The slow component is most evident when scattering from the fully covered surface. The ratio and origin of these components vary with the CO coverage.

Hyperthermal Ar atoms were scattered from a Ru(0001) surface held at temperatures of 180, 400 and 600 K, and from a Ru(0001)-(1×1)D surface held at 114 and 180 K. The resultant angular intensity and energy distributions are complex. The in-plane angular distributions have narrow (FWHM ≤ 10°) near-specular peaks and additional off-specular features. The energy distributions show an oscillatory behavior as a function of outgoing angle. In comparison, scattered Ar atoms from a Ag(111) surface exhibit a broad angular intensity distribution and an energy distribution that qualitatively tracks the binary collision model. The features observed for Ru, which are most evident when scattering from the clean surface at 180 K and from the Ru(0001)-(1×1)D surface, are consistent with rainbow scattering. The measured TOF profiles cannot be adequately described with a single shifted Maxwell-Boltzmann distribution. They can be fitted by two components that exhibit complex variations as a function of outgoing angle. This suggests at least two significantly different site and∕or trajectory dependent energy loss processes at the surface. The results are interpreted in terms of the stiffness of the surface and highlight the anomalous nature of the apparently simple hcp(0001) ruthenium surface.

We have optimized the lowest energy structures and calculated interaction energies for the CO(2)-Ar, CO(2)-N(2), CO(2)-CO, CO(2)-H(2)O, and CO(2)-NH(3) dimers with the recently developed explicitly correlated coupled cluster singles doubles and perturbative triples [CCSD(T)]-F12 methods and the associated VXZ-F12 (where X = D,T,Q) basis sets. For a given cardinal number, we find that results obtained with the CCSD(T)-F12 methods are much closer to the CCSD(T) complete basis set limit than the conventional CCSD(T) results. The relatively modest increase in the computational cost between explicit and conventional CCSD(T) is more than compensated for by the impressive accuracy of the CCSD(T)-F12 method. We recommend use of the CCSD(T)-F12 methods in combination with the VXZ-F12 basis sets for the accurate determination of equilibrium geometries and interaction energies of weakly bound electron donor acceptor complexes.

Singly and doubly charged scandium-water ion-molecule complexes are produced in a supersonic molecular beam by laser vaporization. These ions are mass analyzed and size selected in a specially designed reflectron time-of-flight spectrometer. To probe their structure, vibrational spectroscopy is measured for these complexes in the O-H stretching region using infrared laser photodissociation and the method of rare gas atom predissociation, also known as "tagging." The O-H stretches in these systems are shifted to lower frequency than those for the free water molecule, and the intensity of the symmetric stretch band is strongly enhanced relative to the asymmetric stretch. These effects are more prominent for the doubly charged ions. Partially resolved rotational structure for the Sc(+)(H(2)O)Ar complex shows that the H-O-H bond angle is larger than it is in the free water molecule. Fragmentation and spectral patterns indicate that the coordination of the Sc(2+) ion is filled with six ligands (one water and five argons).

The supermolecule approach has been used to model molecules embedded in solid argon matrix, wherein interaction between the guest and the host atoms in the first solvation shell is evaluated with the use of density functional calculations. Structural stability and simulated spectra have been obtained for formic acid dimer (FAD)-Ar(n) (n = 21-26) clusters. The calculations at the B971∕6-31++G(3df,3pd) level have shown that the tetrasubstitutional site on Ar(111) plane is likely to incorporate FAD most stably, in view of consistency with the matrix shifts available experimentally.

Photodissociation of the ionized argon pentamer, Ar(5)(+), is studied using an extended diatomics-in-molecules interaction model with the inclusion of the spin-orbit coupling and various dynamical approaches. A thorough comparison with the experimental data available in the literature is presented, including photofragment abundances and their kinetic and internal energy distributions. New predictions are reported for ultraviolet photoexcitation energies, a range that has not been studied before either experimentally or theoretically.

The dissociation energetics in the phenol(+)⋯Ar(2)(2π) cluster ion have been investigated using photoionization efficiency and mass analyzed threshold ionization spectroscopy. The appearance energies for the loss of one and two Ar atoms are determined as ∼210 and ∼1115 cm(-1), respectively. The difference between the appearance energy for the first Ar ligand in phenol(+)⋯Ar(2)(2π) and the dissociation energy of the phenol(+)⋯Ar(π) dimer (535cm(-1)) is explained by the isomerization of one π-bound Ar ligand to the OH binding site (H-bond) upon ionization. The energy difference between phenol(+)⋯Ar(2)(2π) and phenol(+)⋯Ar(2)(H/π) could also be estimated to be around 325cm(-1), which corresponds roughly to the difference of the binding energy of a π-bound and H-bound Ar ligands. The binding energy of the H-bound Ar atom in phenol(+)⋯Ar(2)(H/π) is derived to be ∼905cm(-1).

Using high-speed visualization we demonstrate that ultrasound irradiation of pressurized carbon dioxide (CO(2)) induces phenomena that do not occur in ordinary liquids at ambient conditions. For a near-critical mixture of CO(2) and argon, sonication leads to extremely fast local phase separation, in which the system enters and leaves the two-phase region with the frequency of the imposed sound field. This phase transition can propagate with the speed of sound, but can also be located at fixed positions in the case of a standing sound wave. Sonication of a vapor-liquid interface creates a fine dispersion of liquid and vapor, irrespective whether the ultrasound horn is placed in the liquid or the vapor phase. In the absence of an interface, sonication of the liquid leads to ejection of a macroscopic vapor phase from the ultrasound horn with a velocity of several meters per second in the direction of wave propagation. The findings reported here potentially provide a tunable and noninvasive means for enhancing mass and heat transfer in high-pressure fluids.

The Axilrod-Teller equation is incorporated into nonlocal density functional theory (NLDFT) to improve the description of low-temperature argon adsorption isotherm on graphitized carbon black. Using the NLDFT based on the Tarazona smoothed density approximation and the test-particle method, the pair and triplet distribution functions for liquid argon were analyzed. The latter allowed us to quantitatively account for the three-body Axilrod-Teller nonadditive term. An analogous scheme to that of Weeks, Chandler, and Andersen is proposed to incorporate the Axilrod-Teller equation into the NLDFT in the framework of mean field approximation. The three-body nonadditive term appeared to decrease the internal energy of liquid argon by about 5%, which is in agreement with molecular simulation results. The new approach has been applied to the vapor-liquid coexistence and to argon adsorption on graphitized carbon black at 87.3 K. The approach has been shown to fit the experimental adsorption isotherm within the relative error of about 4% as opposed to 37% in the case of the standard NLDFT. This is due to weakening of the intermolecular interaction potential by about 20% in the molecular layer nearest to the graphite surface resulted mainly from the positive Axilrod-Teller contribution of triplets comprising one carbon atom.

A novel threshold photoelectron-photoion coincidence (TPEPICO) imaging spectrometer at the U14-A beamline of the Hefei National Synchrotron Radiation Laboratory is presented. A set of open electron and ion lenses are utilized to map velocity imaging of photoelectrons and photoions simultaneously, in which a repelling electric field using an extra lens is applied to magnify images of photoelectrons instead of traditional accelerating electric field in order to suppress the contribution of energetic electrons in the threshold photoelectron spectroscopy (TPES) and the mass-selected TPEPICO spectroscopy. The typical energy resolution of TPES is measured to be 9 meV (full width at half maximum), as shown on the (2)P(1/2) ionization of argon. The measured mass resolving power for the present TPEPICO imaging spectrometer is above 900 of M/DeltaM. Subsequently as a benchmark, oxygen molecule is photoionized by monochromatic synchrotron radiation at 20.298 eV and dissociates to an oxygen atomic ion and a neutral oxygen atom, and the translation energy distribution of oxygen atomic ion is measured by the time-sliced imaging based on mass-selected TPEPICO experiment. The kinetic energy resolution of the present ion velocity imaging is better than 3% of DeltaE/E.

An investigation of the behavior of O(2) molecules in and on O(2)-doped large (N approximately 8000) Ar host clusters has been performed by means of core and valence photoelectron spectroscopy. Data from pure O(2) and Ar clusters, as well as from O(2)-doped Ar clusters, are presented. The experimental data together with calculations of the binding energy shifts of oxygen molecular ions in and on the surface of a large host Ar cluster show that the diffusion behavior has a strong dependence on the doping pressure. We conclude that the oxygen molecules in the doped Ar host do not partake in band formation, since there is clear vibrational resolution in the spectral features stemming from screened O(2) (+) ions. This implies that valence photoelectron spectroscopy can be used to determine the geometrical structure of this and certain, similar, cluster systems.