Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SOLVATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68

Many-body effects on the structures and stability of Ba{sup 2+}Xe{sub n} (n = 1–39, 54) clusters

The structures and relative stabilities of mixed Ba{sup 2+}Xe{sub n} (n = 1–39, 54) clusters have been theoretically studied using basin-hopping global optimization. Analytical potential energy surfaces were constructed from ab initio or experimental data, assuming either purely additive interactions or including many-body polarization effects and the mutual contribution of self-consistent induced dipoles. For both models the stable structures are characterized by the barium cation being coated by a shell of xenon atoms, as expected from simple energetic arguments. Icosahedral packing is dominantly found, the exceptional stability of the icosahedral motif at n = 12 being further manifested at the size n = 32 where the basic icosahedron is surrounded by a dodecahedral cage, and at n = 54 where the transition to multilayer Mackay icosahedra has occurred. Interactions between induced dipoles generally tend to decrease the Xe-Xe binding, leading to different solvation patterns at small sizes but also favoring polyicosahedral growth. Besides attenuating relative energetic stability, many-body effects affect the structures by expanding the clusters by a few percents and allowing them to deform more.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How far do electrons delocalize?

Electron delocalization is central to chemical bonding, but it is also a fundamentally nonclassical and nonintuitive quantum mechanical phenomenon. Tools to quantify and visualize electron delocalization help to understand, teach, and predict chemical reactivity. We develop a new approach to quantify and visualize electron delocalization in real space. Our electron delocalization range function EDR(r{sup -vector};u) quantifies the degree to which electrons at point r{sup -vector} in a calculated wavefunction delocalize over length scale u. Its predictions are physically reasonable. For example, EDR(r{sup -vector};u=0.25 bohr) is close to one at points r{sup -vector} in the cores of first-row atoms, consistent with the localization of core electrons to ∼0.25 bohr. EDR(r{sup -vector};u=1 bohr) is close to one at points r{sup -vector} in typical covalent bonds, consistent with electrons delocalizing over the length of the bond. Our approach provides a rich representation of atomic shell structure; covalent and ionic bonding; the delocalization of excited states, defects, and solvated electrons; metallic and insulating systems; and bond stretching and strong correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical structural and vibrational study of 5-trifluoromethyluracil. A comparison with uracil

In the present work, a comparative study on the structural and vibrational properties of the 5-trifluoromethyluracil (TFMU) derivative with those corresponding to uracil in gas and aqueous solution phases was performed combining the available H{sup 1}-NMR, C{sup 13}-NMR, F{sup 19}-NMR and FTIR spectra with Density Functional Theory (DFT) calculations. Three stable conformers were theoretically determined in both media by using the hybrid B3LYP/6-31G* method. The solvent effects were simulated by means of the self-consistent reaction field (SCRF) method employing the integral equation formalism variant (IEFPCM). Complete assignments of the vibrational spectra in both phases were performed combining the internal coordinates analysis and the DFT calculations with the Scaled Quantum Mechanics Force Field (SQMFF) methodology. The atomic charges, bond orders, solvation energies, dipole moments, molecular electrostatic potentials and force constants parameters were calculated for the three conformers of TFMU in gas phase and aqueous solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simplified continuum solvent model with a smooth cavity based on volumetric data

We present a continuum solvent model (CSM) with a smooth cavity for the application in grid-based electronic structure methods. The cavity is identified with the inherently smooth distribution function of a binary mixture at infinite dilution. We obtain a cavity model based on atomic van der Waals radii and one free parameter controlling the overall size. This single parameter is sufficient to adequately reproduce experimental partial molar volumes. The CSM based on this cavity is of similar accuracy in the prediction of aqueous solvation Gibbs energies of small neutral molecules and ions as other CSMs with a smooth cavity. We apply the model to systems in non-aqueous solution, i.e., spiropyran/merocyanin energetics, a proton transfer reaction in dimethyl sulfoxide, and the electrostatic screening of charged gold clusters in an ionic liquid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical investigation of the relative stability of Na{sup +}He{sub n} (n = 2–24) clusters: Many-body versus delocalization effects

The solvation of the Na{sup +} ion in helium clusters has been studied theoretically using optimization methods. A many-body empirical potential was developed to account for Na{sup +}–He and polarization interactions, and the most stable structures of Na{sup +}He{sub n} clusters were determined using the basin-hopping method. Vibrational delocalization was accounted for using zero-point energy corrections at the harmonic or anharmonic levels, the latter being evaluated from quantum Monte Carlo simulations for spinless particles. From the static perspective, many-body effects are found to play a minor role, and the structures obtained reflect homogeneous covering up to n = 10, followed by polyicosahedral packing above this size, the cluster obtained at n = 12 appearing particularly stable. The cationic impurity binds the closest helium atoms sufficiently to negate vibrational delocalization at small sizes. However, this snowball effect is obliterated earlier than shell completion, the nuclear wavefunctions of {sup 4}He{sub n}Na{sup +} with n = 5–7, and n > 10 already exhibiting multiple inherent structures. The decrease in the snowball size due to many-body effects is consistent with recent mass spectrometry measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum, classical, and hybrid QM/MM calculations in solution: General implementation of the ddCOSMO linear scaling strategy

We present the general theory and implementation of the Conductor-like Screening Model according to the recently developed ddCOSMO paradigm. The various quantities needed to apply ddCOSMO at different levels of theory, including quantum mechanical descriptions, are discussed in detail, with a particular focus on how to compute the integrals needed to evaluate the ddCOSMO solvation energy and its derivatives. The overall computational cost of a ddCOSMO computation is then analyzed and decomposed in the various steps: the different relative weights of such contributions are then discussed for both ddCOSMO and the fastest available alternative discretization to the COSMO equations. Finally, the scaling of the cost of the various steps with respect to the size of the solute is analyzed and discussed, showing how ddCOSMO opens significantly new possibilities when cheap or hybrid molecular mechanics/quantum mechanics methods are used to describe the solute.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energy deposition dynamics of femtosecond pulses in water

We exploit inverse Raman scattering and solvated electron absorption to perform a quantitative characterization of the energy loss and ionization dynamics in water with tightly focused near-infrared femtosecond pulses. A comparison between experimental data and numerical simulations suggests that the ionization energy of water is 8 eV, rather than the commonly used value of 6.5 eV. We also introduce an equation for the Raman gain valid for ultra-short pulses that validates our experimental procedure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Solvatochromism of 9,10-phenanthrenequinone: An electronic and resonance Raman spectroscopic study

Solvent effects play a vital role in various chemical, physical, and biological processes. To gain a fundamental understanding of the solute-solvent interactions and their implications on the energy level re-ordering and structure, UV-VIS absorption, resonance Raman spectroscopic, and density functional theory calculation studies on 9,10-phenanthrenequinone (PQ) in different solvents of diverse solvent polarity has been carried out. The solvatochromic analysis of the absorption spectra of PQ in protic dipolar solvents suggests that the longest (1n-π{sup 1}*; S{sub 1} state) and the shorter (1π-π{sup 1}*; S{sub 2} state) wavelength band undergoes a hypsochromic and bathochromic shift due to intermolecular hydrogen bond weakening and strengthening, respectively. It also indicates that hydrogen bonding plays a major role in the differential solvation of the S{sub 2} state relative to the ground state. Raman excitation profiles of PQ (400–1800 cm{sup −1}) in various solvents followed their corresponding absorption spectra therefore the enhancements on resonant excitation are from single-state rather than mixed states. The hyperchromism of the longer wavelength band is attributed to intensity borrowing from the nearby allowed electronic transition through vibronic coupling. Computational calculation with C{sub 2ν} symmetry constraint on the S{sub 2} state resulted in an imaginary frequency along the low-frequency out-of-plane torsional modes involving the C=O site and therefore, we hypothesize that this mode could be involved in the vibronic coupling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A first principles approach to the electronic properties of liquid and supercritical CO{sub 2}

The electronic absorption spectra of liquid and supercritical CO{sub 2} (scCO{sub 2}) are investigated by coupling a many-body energy decomposition scheme to configurations generated by Born-Oppenheimer molecular dynamics. A Frenkel exciton Hamiltonian formalism was adopted and the excitation energies were calculated with time dependent density functional theory. A red-shift of ∼ 0.2 eV relative to the gas-phase monomer is observed for the first electronic absorption maximum in liquid and scCO{sub 2}. The origin of this shift, which is not very dependent on deviations from the linearity of the CO{sub 2} molecule, is mainly related to polarization effects. However, the geometry changes of the CO{sub 2} monomer induced by thermal effects and intermolecular interactions in condensed phase lead to the appearance of an average monomeric electric dipole moment 〈μ〉 = 0.26 ± 0.04 D that is practically the same at liquid and supercritical conditions. The predicted average quadrupole moment for both liquid and scCO{sub 2} is 〈Θ〉 = − 5.5 D Å, which is increased by ∼ −0.9 D Å relative to its gas-phase value. The importance of investigating the electronic properties for a better understanding of the role played by CO{sub 2} in supercritical solvation is stressed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interactions of model biomolecules. Benchmark CC calculations within MOLCAS

We present results using the OVOS approach (Optimized Virtual Orbitals Space) aimed at enhancing the effectiveness of the Coupled Cluster calculations. This approach allows to reduce the total computer time required for large-scale CCSD(T) calculations about ten times when the original full virtual space is reduced to about 50% of its original size without affecting the accuracy. The method is implemented in the MOLCAS computer program. When combined with the Cholesky decomposition of the two-electron integrals and suitable parallelization it allows calculations which were formerly prohibitively too demanding. We focused ourselves to accurate calculations of the hydrogen bonded and the stacking interactions of the model biomolecules. Interaction energies of the formaldehyde, formamide, benzene, and uracil dimers and the three-body contributions in the cytosine – guanine tetramer are presented. Other applications, as the electron affinity of the uracil affected by solvation are also shortly mentioned.

74 ATOMIC AND MOLECULAR PHYSICS↗

Design of molecular switches and sensors based on proton transfer - Theory vs. experiment

Substituent and solvent effects on tautomeric equilibria of isomeric phenylazo naphthols and the analogous Schiff bases are calculated by various quantum chemical methods and solvation models. Special emphasis is put on the reliability of the computational procedures used. The results obtained should aid in resolving discrepancies of experimental data resulting from UV/Vis and NMR spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Maximum probability domains for the analysis of the microscopic structure of liquids

We introduce the concept of maximum probability domains (MPDs), developed in the context of the analysis of electronic densities, in the study of the microscopic spatial structures of liquids. The idea of locating a particle in a three dimensional region, by determining the domain where the probability of finding that, and only that, particle is maximum, gives an interesting characterization of the local structure of the liquid. The optimization procedure, required for the search of the domain of maximum probability, is carried out by the implementation of the level set method. Results for a couple of case studies are presented, to illustrate the structure of liquid water at ambient conditions and upon increasing pressure from the point of view of MPDs and to compare the information encoded in the solvation shells of sodium in water with, once again, that extracted from the MPDs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electrolyte pore/solution partitioning by expanded grand canonical ensemble Monte Carlo simulation

Using a newly developed grand canonical Monte Carlo approach based on fractional exchanges of dissolved ions and water molecules, we studied equilibrium partitioning of both components between laterally extended apolar confinements and surrounding electrolyte solution. Accurate calculations of the Hamiltonian and tensorial pressure components at anisotropic conditions in the pore required the development of a novel algorithm for a self-consistent correction of nonelectrostatic cut-off effects. At pore widths above the kinetic threshold to capillary evaporation, the molality of the salt inside the confinement grows in parallel with that of the bulk phase, but presents a nonuniform width-dependence, being depleted at some and elevated at other separations. The presence of the salt enhances the layered structure in the slit and lengthens the range of inter-wall pressure exerted by the metastable liquid. Solvation pressure becomes increasingly repulsive with growing salt molality in the surrounding bath. Depending on the sign of the excess molality in the pore, the wetting free energy of pore walls is either increased or decreased by the presence of the salt. Because of simultaneous rise in the solution surface tension, which increases the free-energy cost of vapor nucleation, the rise in the apparent hydrophobicity of the walls has not been shown to enhance the volatility of the metastable liquid in the pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultrafast vibrational spectroscopy (2D-IR) of CO{sub 2} in ionic liquids: Carbon capture from carbon dioxide’s point of view

The CO{sub 2}ν{sub 3} asymmetric stretching mode is established as a vibrational chromophore for ultrafast two-dimensional infrared (2D-IR) spectroscopic studies of local structure and dynamics in ionic liquids, which are of interest for carbon capture applications. CO{sub 2} is dissolved in a series of 1-butyl-3-methylimidazolium-based ionic liquids ([C{sub 4}C{sub 1}im][X], where [X]{sup −} is the anion from the series hexafluorophosphate (PF{sub 6}{sup −}), tetrafluoroborate (BF{sub 4}{sup −}), bis-(trifluoromethyl)sulfonylimide (Tf{sub 2}N{sup −}), triflate (TfO{sup −}), trifluoroacetate (TFA{sup −}), dicyanamide (DCA{sup −}), and thiocyanate (SCN{sup −})). In the ionic liquids studied, the ν{sub 3} center frequency is sensitive to the local solvation environment and reports on the timescales for local structural relaxation. Density functional theory calculations predict charge transfer from the anion to the CO{sub 2} and from CO{sub 2} to the cation. The charge transfer drives geometrical distortion of CO{sub 2}, which in turn changes the ν{sub 3} frequency. The observed structural relaxation timescales vary by up to an order of magnitude between ionic liquids. Shoulders in the 2D-IR spectra arise from anharmonic coupling of the ν{sub 2} and ν{sub 3} normal modes of CO{sub 2}. Thermal fluctuations in the ν{sub 2} population stochastically modulate the ν{sub 3} frequency and generate dynamic cross-peaks. These timescales are attributed to the breakup of ion cages that create a well-defined local environment for CO{sub 2}. The results suggest that the picosecond dynamics of CO{sub 2} are gated by local diffusion of anions and cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measuring correlated electronic and vibrational spectral dynamics using line shapes in two-dimensional electronic-vibrational spectroscopy

Two-dimensional electronic-vibrational (2DEV) spectroscopy is an experimental technique that shows great promise in its ability to provide detailed information concerning the interactions between the electronic and vibrational degrees of freedom in molecular systems. The physical quantities 2DEV is particularly suited for measuring have not yet been fully determined, nor how these effects manifest in the spectra. In this work, we investigate the use of the center line slope of a peak in a 2DEV spectrum as a measure of both the dynamic and static correlations between the electronic and vibrational states of a dye molecule in solution. We show how this center line slope is directly related to the solvation correlation function for the vibrational degrees of freedom. We also demonstrate how the strength with which the vibration on the electronic excited state couples to its bath can be extracted from a set of 2DEV spectra. These analytical techniques are then applied to experimental data from the laser dye 3,3′-diethylthiatricarbocyanine iodide in deuterated chloroform, where we determine the lifetime of the correlation between the electronic transition frequency and the transition frequency for the backbone C = C stretch mode to be ∼1.7 ps. Furthermore, we find that on the electronic excited state, this mode couples to the bath ∼1.5 times more strongly than on the electronic ground state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical analysis of co-solvent effect on the proton transfer reaction of glycine in a water–acetonitrile mixture

The co-solvent effect on the proton transfer reaction of glycine in a water–acetonitrile mixture was examined using the reference interaction-site model self-consistent field theory. The free energy profiles of the proton transfer reaction of glycine between the carboxyl oxygen and amino nitrogen were computed in a water–acetonitrile mixture solvent at various molar fractions. Two types of reactions, the intramolecular proton transfer and water-mediated proton transfer, were considered. In both types of the reactions, a similar tendency was observed. In the pure water solvent, the zwitterionic form, where the carboxyl oxygen is deprotonated while the amino nitrogen is protonated, is more stable than the neutral form. The reaction free energy is −10.6 kcal mol{sup −1}. On the other hand, in the pure acetonitrile solvent, glycine takes only the neutral form. The reaction free energy from the neutral to zwitterionic form gradually increases with increasing acetonitrile concentration, and in an equally mixed solvent, the zwitterionic and neutral forms are almost isoenergetic, with a difference of only 0.3 kcal mol{sup −1}. The free energy component analysis based on the thermodynamic cycle of the reaction also revealed that the free energy change of the neutral form is insensitive to the change of solvent environment but the zwitterionic form shows drastic changes. In particular, the excess chemical potential, one of the components of the solvation free energy, is dominant and contributes to the stabilization of the zwitterionic form.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An EQT-cDFT approach to determine thermodynamic properties of confined fluids

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗