Absorption and emission lineshapes and solvation dynamics of NO in supercritical Ar

We perform a theoretical study of electronic spectroscopy of dilute NO in supercritical Ar fluid. Absorption and emission lineshapes for the A(2)Sigma(+)<--X(2)Pi Rydberg transition of NO in argon have been previously measured and simulated, which yielded results for the NO/Ar ground- and excited-state pair potentials [Larregaray et al., Chem. Phys. 308, 13 (2005)]. Using these potentials, we have performed molecular dynamics simulations and theoretical statistical mechanical calculations of absorption and emission lineshapes and nonequilibrium solvation correlation functions for a wide range of solvent densities and temperatures. Theory was shown to be in good agreement with simulation. Linear response treatment of solvation dynamics was shown to break down at near-critical temperature due to dramatic change in the solute-solvent microstructure upon solute excitation to the Rydberg state and the concomitant increase of the solute size.     

Absorption and emission lineshapes and ultrafast solvation dynamics of NO in parahydrogen

We perform a theoretical study on the electronic spectroscopy of dilute NO impurity embedded in parahydrogen (p-H(2)). Absorption and emission lineshapes for the A (2)Sigma(+)<--X (2)Pi Rydberg transition of NO in parahydrogen have been previously measured and simulated, which yielded results for the NO/p-H(2) ground and excited state pair potentials [L. Bonacina et al., J. Chem. Phys. 125, 054507 (2006)]. Using these potentials, we performed molecular dynamics simulation, theoretical statistical mechanical calculations of absorption and emission lineshapes, and both equilibrium and nonequilibrium solvation correlation functions for NO chromophore in parahydrogen. Theory was shown to be in good agreement with simulation. Linear response treatment of solvation dynamics was shown to break down due to a dramatic change in the solute-solvent microstructure upon solute excitation to the Rydberg state and the concomitant increase of the solute size.     

Two-color excitation of NO in a supersonic free jet. Autoionization of high rydberg states

High Rydberg states of NO above the ionization limit have been measured for the isolated molecule in a supersonic free jet by two-color multiphoton ionization. Three Rydberg series (ns, np and nf) were identified, which appeared by rotational and the vibrational autoionization. The rotational structures of the 13s(υ = 1), 13p(υ = 1) and 12f(υ = 1) states have been analyzed in detail. The fluorescence dip spectra for the intermediate A²Σ⁺(3sσ) state have been measured simultaneously and the cross sections of the one-photon absorption to the high Rydberg states from the A²Σ⁺(υ = 1) state have been determined.     

Perturbations of molecular extravalence excitations by rare-gas fluids

In this paper we report the results of an experimental study of the vacuum ultraviolet absorption spectra of molecular impurity states of methyl iodide in Ar (density range ? = 0–1.4 g cm^?3) and in Kr (? = 0–2.3 g cm^?3), of carbon disulphide in Ar (? = 0–1.4 g cm^?3) and of formaldehyde in Ar (? = 0–1.25 g cm^?3). The experimental results provide new information regarding medium perturbations of intravalenc transitions, of the lowest extravalence transitions and of transitions to mixed valence—Rydberg configurations, which serve as a diagnostic tool to distinguish between different types of electronic excitations. All the lowest extravalence molecular excitations exhibit appreciable blue spectral shifts at moderate and at high fluid densities, intravalence transitions are practically insensitive to medium effects, while excitations to mixed valence—Rydberg configurations are characterized by a moderate blue spectral shift. New information has been obtained concerning the energetics of molecular ionization processes in a dense fluid. The high n = 2–5 Rydberg states of CH₃l exhibit a large red shift at moderate (? = 0–0.5 cm^?3) Ar densities. The ionization potential E_g and the effective Rydberg constant G for CH₃I in Ar was found to decrease from G = 13.6 eV and E_g = 9.55 eV at ? = 0 and E_g = 9.08 eV and constant G for CH₃l in Ar was found to decrease from G = 13.6 eV and E_g = 9.55eV at ? = 0 and E_g = 9.08 eV and G ≈ 7.15 eV at ? = 0.5 g cm^?3. Experimental evidence was obtained for the identification of n = 2 molecular Wannier impurity states of CH₃I and of CH₂O in liquid Ar. These spectroscopic data result in E_g ≈ 8.6 eV for CH₃I in liquid Ar and E_g ≈ 10.2 eV for CH₂O in liquid Ar.     

Computational Analysis of Solute–Solvent Coupling Magnitude in the Z/E Isomerization Reaction of Nitroazobenzene and Benzylideneanilines

The dynamic solvent effect often arises in solution reactions, where coupling between chemical reaction and solvent fluctuation plays a decisive role in the reaction kinetics. In this study, the Z/E isomerization reaction of nitoroazobenzene and benzylideneanilines in the ground state was computationally studied by molecular dynamics simulations. The non-equilibrium solvation effect was analyzed using two approaches: (1) metadynamics Gibbs energy surface exploration and (2) solvation Gibbs energy evaluation using a frozen solvation droplet model. The solute–solvent coupling parameter (C_coupled) was estimated by the ratio of the solvent fluctuation Gibbs energy over the corresponding isomerization activation Gibbs energy. The results were discussed in comparison with the ones estimated by means of the analytical models based on a reaction–diffusion equation with a sink term. The second approach using a frozen solvation droplet reached qualitative agreement with the analytical models, while the first metadynamics approach failed. This is because the second approach explicitly considers the non-equilibrium solvation in the droplet, which consists of a solute at the reactant geometry immersed in the pre-organized solvents fitted with the solute at the transition state geometry.     

A 3D-analysis of the Cl?�Cbenzene dimer solvation by Ar atoms

The solvation of the Cl^??Cbenzene (Cl^??CBz) aggregate by Ar atoms has been investigated employing molecular dynamics (MD) simulations. The gradual evolution from cluster rearrangement to solvation dynamics is discussed by considering ensembles of n (n = 1,...,30) Ar atoms around the Cl^??CBz clusters. The energetic of the solvated cluster is decomposed as a sum of pairs (including both the Ar?CAr and the Cl^??CAr terms), Cl^??CBz and Ar?CBz interactions and their relative contributions are analyzed as a function of the cluster size. The geometrical distribution of Ar atoms around Cl^??CBz is investigated in terms of radial distribution functions (RDF), bidimensional (2D) angular distributions and tridimensional (3D) probability densities. The variation on the spatial distribution of the Ar atoms around Cl^??CBz when the Ar number increases is investigated from a novel prospective, employing spherical coordinates of the solvent atoms within an inertial reference frame. Isomerization processes are also studied.     

A molecular dynamics computer simulation study of room-temperature ionic liquids. II. Equilibrium and nonequilibrium solvation dynamics

The molecular dynamics (MD) simulation study of solvation structure and free energetics in 1-ethyl-3-methylimidazolium chloride and 1-ethyl-3-methylimidazolium hexafluorophosphate using a probe solute in the preceding article [Y. Shim, M. Y. Choi and H. J. Kim, J. Chem. Phys. 122, 044510 (2005)] is extended to investigate dynamic properties of these liquids. Solvent fluctuation dynamics near equilibrium are studied via MD and associated time-dependent friction is analyzed via the generalized Langevin equation. Nonequilibrium solvent relaxation following an instantaneous change in the solute charge distribution and accompanying solvent structure reorganization are also investigated. Both equilibrium and nonequilibrium solvation dynamics are characterized by at least two vastly different time scales--a subpicosecond inertial regime followed by a slow diffusive regime. Solvent regions contributing to the subpicosecond nonequilibrium relaxation are found to vary significantly with initial solvation configurations, especially near the solute. If the solvent density near the solute is sufficiently high at the outset of the relaxation, subpicosecond dynamics are mainly governed by the motions of a few ions close to the solute. By contrast, in the case of a low local density, solvent ions located not only close to but also relatively far from the solute participate in the subpicosecond relaxation. Despite this difference, linear response holds reasonably well in both ionic liquids.     

Solvation of coumarin 153 in supercritical fluoroform

We present a study of local density augmentation around an attractive solute (i.e., giving rise to more attractive interaction with the solvent than solvent-solvent interactions) in supercritical fluoroform. This work is based on molecular dynamics simulations of coumarin 153 in supercritical fluoroform at densities both above and below the critical density, ranging from dilute gas-like to liquid-like, at a reduced temperature (T/T(c)) of 1.03. We focused on studying the structure of the solvation shell and the variation of the solute electronic absorption and emission shifts with density. Quantum calculations at the density functional theory (DFT) level were run on the solute in the ground state, and time-dependent DFT calculations were performed in the solute excited state in order to determine the solute-solvent potential parameters. The results obtained for the Stokes shift are in agreement with the experimental measurements. To evaluate local density augmentation from simulations, we used two different definitions, one based on the solvation number and the other derived from solvatochromic shifts. In the former case, the agreement with experimental results is good, while, in the latter case, better agreement is achieved by perturbatively including the induced-dipole contribution to the solvation energy.     

Evaluating the Free Energies of Solvation and Electronic Structures of Lithium‐Ion Battery Electrolytes

Adaptive biasing force molecular dynamics simulations and density functional theory calculations were performed to understand the interaction of Li⁺ with pure carbonates and ethylene carbonate (EC)‐based binary mixtures. The most favorable Li carbonate cluster configurations obtained from molecular dynamics simulations were subjected to detailed structural and thermochemistry calculations on the basis of the M06‐2X/6‐311++G(d,p) level of theory. We report the ranking of these electrolytes on the basis of the free energies of Li‐ion solvation in carbonates and EC‐based mixtures. A strong local tetrahedral order involving four carbonates around the Li⁺ was seen in the first solvation shell. Thermochemistry calculations revealed that the enthalpy of solvation and the Gibbs free energy of solvation of the Li⁺ ion with carbonates are negative and suggested the ion–carbonate complexation process to be exothermic and spontaneous. Natural bond orbital analysis indicated that Li⁺ interacts with the lone pairs of electrons on the carbonyl oxygen atom in the primary solvation sphere. These interactions lead to an increase in the carbonyl (C=O) bond lengths, as evidenced by a redshift in the vibrational frequencies [ν(C=O)] and a decrease in the electron density values at the C=O bond critical points in the primary solvation sphere. Quantum theory of atoms in molecules, localized molecular orbital energy decomposition analysis (LMO‐EDA), and noncovalent interaction plots revealed the electrostatic nature of the Li⁺ ion interactions with the carbonyl oxygen atoms in these complexes. On the basis of LMO‐EDA, the strongest attractive interaction in these complexes was found to be the electrostatic interaction followed by polarization, dispersion, and exchange interactions. Overall, our calculations predicted EC and a binary mixture of EC/dimethyl carbonate to be appropriate electrolytes for Li‐ion batteries, which complies with experiments and other theoretical results.     

A MS-CASPT2 study of the low-lying electronic excited states of CH2BrCl

The low-lying singlet excited states of CH₂BrCl have been calculated using multiconfigurational CASSCF, second-order perturbation theory CASPT2 and its multistate extension MS-CASPT2. The CASSCF method shows spurious valence–Rydberg mixing and a wrong order of states. Inclusion of dynamical correlation by single root CASPT2 lowers dramatically the energy of the valences states but does not lead to a complete separation between valence and Rydberg states. This situation is improved by the MS-CASPT2 calculations, which gives two valence states for both A^′ and A″ symmetries below the lowest Rydberg state, corresponding to n(Br)→σ^*(C–Br) and n(Cl)→σ^*(C–Cl) transitions at 6.1 eV (203 nm) and 7.2 eV (173 nm), and being repulsive along C–Br and C–Cl coordinates.     

The effects of solute-solvent electrostatic interactions on solvation dynamics in supercritical CO2

We present here the results of molecular-dynamics simulation of solvation dynamics in supercritical CO(2) at a temperature of about 1.05T(c), where T(c) is the critical temperature, and at a series of densities ranging from 0.4 to 2.0 of the critical density rho(c). We focus on electrostatic solvation dynamics, representing the electronic excitation of the chromophore as a change in its charge distribution from a quadrupolar-symmetry ground state to a dipolar excited state. Two perturbations are considered, corresponding to different magnitudes of solute excited-state dipoles, denoted as d5 and d8. The d8 solute is more attractive, leading to a larger enhancement in CO(2) clustering upon solute electronic excitation. This has a large impact on solvation dynamics, especially at densities below rho(c). At these densities, solvation dynamics is much slower for the d8 than for the d5 solute. For both solutes, solvation dynamics becomes faster at densities above rho(c) at which solvent clustering diminishes. We show that the slowest solvation time scale is associated with solvent clustering and we relate it to solute-solvent mutual translational diffusion and the extent of change in effective local density resulting from solute electronic excitation.     

Solvation in supercritical water

Solvation in supercritical water under equilibrium and nonequilibrium conditions is studied via molecular dynamics simulations. The influence of solute charge distributions and solvent density on the solvation structures and dynamics is examined with a diatomic probe solute molecule. It is found that the solvation structure varies dramatically with the solute dipole moment, especially in low-density water, in accord with many previous studies on ion solvation. This electrostrictive effect has important consequences for solvation dynamics. In the case of a nonequilibrium solvent relaxation, if there are sufficiently many water molecules close to the solute at the outset of the relaxation, the solvent response measured as a dynamic Stokes shift is almost completely governed by inertial rotations of these water molecules. By contrast, in the opposite case of a low local solvent density near the solute, not only rotations but also translations of water molecules play an important role in solvent relaxation dynamics. The applicability of a linear response is found to be significantly restricted at low water densities.     

X‐Ray Emission Spectroscopy: A Spectroscopic Measure for the Determination of NO Oxidation States in Fe–NO Complexes

Extensive study of the electronic structure of Fe‐NO complexes using a variety of spectroscopic methods was attempted to understand how iron controls the binding and release of nitric oxide. The comparable energy levels of NO π* orbitals and Fe 3d orbitals complicate the bonding interaction within Fe? NO complexes and puzzle the quantitative assignment of NO oxidation state. Enemark–Feltham notation, {Fe(NO)_x}ⁿ, was devised to circumvent this puzzle. This 40‐year puzzle is revisited using valence‐to‐core X‐ray emission spectroscopy (V2C XES) in combination with computational study. DFT calculation establishes a linear relationship between ΔE_σ2s*‐σ2p of NO and its oxidation state. V2C Fe XES study of Fe? NO complexes reveals the ΔE_σ2s*‐σ2p of NO derived from NO σ_2s*/σ_2p→Fe_1s transitions and determines NO oxidation state in Fe? NO complexes. Quantitative assignment of NO oxidation state will correlate the feasible redox process of nitric oxide and Fe‐nitrosylation biology.     

Molecular solvation in water-methanol and water-sorbitol mixtures: the roles of preferential hydration, hydrophobicity, and the equation of state

Molecular dynamics simulations of aqueous mixtures of methanol and sorbitol were performed over a wide range of binary composition, density (pressure), and temperature to study the equation of state and solvation of small apolar solutes. Experimentally, methanol is a canonical solubilizing agent for apolar solutes and a protein denaturant in mixed-aqueous solvents; sorbitol represents a canonical "salting-out" or protein-stabilizing cosolvent. The results reported here show increasing sorbitol concentration under isothermal, isobaric conditions results in monotonic increases in apolar solute excess chemical potential (mu2ex) over the range of experimentally relevant temperatures. For methanol at elevated temperatures, increasing cosolvent composition results in monotonically decreasing mu2ex. However, at lower temperatures mu2ex exhibits a maximum versus cosolvent concentration, as seen experimentally for Ar in ethanol-water solutions. Both density anomalies and hydrophobic effects--characterized by temperatures of density maxima and apolar solute solubility minima, respectively--are suppressed upon addition of either sorbitol or methanol at all temperatures and compositions simulated here. Thus, the contrasting effects of sorbitol and methanol on solute chemical potential cannot be explained by qualitative differences in their ability to enhance or suppress hydrophobic effects. Rather, we find mu2ex values across a broad range of temperatures and cosolvent composition can be quantitatively explained in terms of isobaric changes in solvent density--i.e., the equation of state--along with the corresponding packing fraction of the solvent. Analysis in terms of truncated preferential interaction parameters highlights that care must be taken in interpreting cosolvent effects on solvation in terms of local preferential hydration.     

Incorporating solvent and ion screening into molecular dynamics using the finite-difference Poisson–Boltzmann method

The finite difference method for solving the Poisson–Boltzmann equation is used to calculate the reaction field acting on a macromolecular solute due to the surrounding water and ions. Comparisons with analytical test cases indicate that the solvation forces can be calculated rapidly and accurately with this method. These forces act to move charged solute atoms towards the solvent where they are better solvated, and to screen interactions between charges. A way of combining such calculations with conventional molecular dynamics force fields is proposed which requires little modification of existing molecular dynamics programs. Simulations on the alanine dipeptide show that solvent forces affect the conformational dynamics by reducing the preference for internal H-bonding forms, increasing the R-alpha helix preference and reducing transition barriers. These solvent effects are similar to previous explicit solvent simulations, but require little more computation than vacuum simulations. The method should scale up with little increase in computational cost to larger molecules such as proteins and nucleic acids.     

Electronic Structure Insights into the Solvation of Magnesium Ions with Cyclic and Acyclic Carbonates

A computational framework to rank the solvation behavior of Mg²⁺ in carbonates by using molecular dynamics simulations and density functional theory is reported. Based on the binding energies and enthalpies of solvation calculated at the M06‐2X/6‐311++G(d,p) level of theory and the free energies of solvation from ABF‐MD simulations, we find that ethylene carbonate (EC) and the ethylene carbonate:propylene carbonate (EC:PC) binary mixture are the best carbonate solvents for interacting with Mg²⁺. Natural bond orbital and quantum theory of atoms in molecules analyses support the thermochemistry calculations with the highest values of charge transfer, perturbative stabilization energies, electron densities, and Wiberg bond indices being observed in the Mg²⁺(EC) and Mg²⁺(EC:PC) complexes. The plots of the noncovalent interactions indicate that those responsible for the formation of Mg²⁺ carbonate complexes are strong‐to‐weak attractive interactions, depending on the regions that are interacting. Finally, density of state calculations indicate that the interactions between Mg²⁺ and the carbonate solvents affects the HOMO and LUMO states of all carbonate solvents and moves them to more negative energy values.     

Liquid-liquid phase separation of binary Lennard-Jones fluid in slit nanopores

The capillary phase separation of a binary mixture of two truncated and shifted Lennard-Jones (LJ) Ar liquids in slit-shaped oxygen nanopores is examined. The LJ parameters—ε(Ar(A)–Ar(A))=ε(Ar(B)–Ar(B))=0.8ε(Ar(A)–Ar(B)) and 0.5ε(Ar(A)–O)?=?ε(Ar(B)–O)—were used to distinguish the two Ar liquids. The cut off distance for Ar was 3.5σ. We employed a molecular dynamics (MD) technique in which a pore space was connected with a bulk solution to easily determine the equilibrium bulk concentration. Liquid phase isotherms were obtained for pores with widths ranging from 5.5σ to 9.5σ, and the relation between the pore width and the phase separation concentration was determined. Each simulation was run until the bulk concentration attained equilibrium (1–2 μs). The MD results show that the Patrick model overestimates the bulk concentration for a given pore size. We proposed a modified Patrick model in which the pore wall potential is considered. In our model, the Gibbs-Tolman-Koenig-Buff effect is not considered for the interfacial tension since two surfaces of tension exist on both sides of the equimolar dividing surface of the two-Ar liquid phase. The two surfaces of tension neutralized Gibbs-Tolman-Koenig-Buff effect each other. The present simple model successfully describes the relation to prove its reliability.     

Nature of the resonant states of NO−

The stabilization method of Taylor has been used to calculate the nature of the states of NO^? corresponding to the lowest resonances observed recently by Sanche and Schulz. It is confirmed that they consist of two Rydberg electrons in the core of the X ¹Σ⁺ ground state of the ion. The proposed assignments are the (Rsσ)²¹Σ⁺, (Rsσ)(Rpπ) ³Π, (Rsσ)(Rpσ) ³Σ⁺ states. The fourth resonance is attributed to an (Rpπ)²³Σ⁺ state. The Rydberg or of these states of the negative ion are more diffuse than the corresponding ones of the neutral molecule.