Descripción: We use ring-polymer-molecular-dynamics (RPMD) techniques and the semi-empirical q-TIP4P/F water model to investigate the relationship between hydrogen bond connectivity and the characteristics of nuclear position fluctuations, including explicit incorporation of quantum effects, for the energetically low lying isomers of the prototype cluster [H2O] 8 at T = 50 K and at 150 K. Our results reveal that tunneling and zero-point energy effects lead to sensible increments in the magnitudes of the fluctuations of intra and intermolecular distances. The degree of proton spatial delocalization is found to map logically with the hydrogen-bond connectivity pattern of the cluster. Dangling hydrogen bonds exhibit the largest extent of spatial delocalization and participate in shorter intramolecular O-H bonds. Combined effects from quantum and polarization fluctuations on the resulting individual dipole moments are also examined. From the dynamical side, we analyze the characteristics of the infrared absorption spectrum. The incorporation of nuclear quantum fluctuations promotes red shifts and sensible broadening relative to the classical profile, bringing the simulation results in much more satisfactory agreement with direct experimental information in the mid and high frequency range of the stretching band. While RPMD predictions overestimate the peak position of the low frequency shoulder, the overall agreement with that reported using an accurate, parameterized, many-body potential is reasonable, and far superior to that one obtains by implementing a partially adiabatic centroid molecular dynamics approach. Quantum effects on the collective dynamics, as reported by instantaneous normal modes, are also discussed. © 2013 AIP Publishing LLC.

Descripción: Ring polymer molecular dynamics experiments have been carried out to examine effects derived from nuclear quantum fluctuations at ambient conditions on equilibrium and non-equilibrium dynamical characteristics of charge solvation by a popular simple, rigid, water model, SPC/E, and for a more recent, and flexible, q-TIP4P/F model, to examine the generality of conclusions. In particular, we have recorded the relaxation of the solvent energy gap following instantaneous, ±e charge jumps in an initially uncharged Lennard-Jones-like solute. In both charge cases, quantum effects are reflected in sharper decays at the initial stages of the relaxation, which produce up to a ∼20% reduction in the characteristic timescales describing the solvation processes. For anionic solvation, the magnitude of polarization fluctuations controlling the extent of the water proton localization in the first solvation shell is somewhat more marked than for cations, bringing the quantum solvation process closer to the classical case. Effects on the solvation response from the explicit incorporation of flexibility in the water Hamiltonian are also examined. Predictions from linear response theories for the overall relaxation profile and for the corresponding characteristic timescales are reasonably accurate for the solvation of cations, whereas we find that they are much less satisfactory for the anionic case. © 2013 AIP Publishing LLC.

Descripción: New insights into equilibrium and dynamical aspects of electron photodetachment reactions in small water clusters were given. It focuses on assessing the effects of thermal and polarization fluctuations provided by three cluster environments with different extents of spatial confinement, on the microscopic mechanisms that drive the reaction. These fluctuations, in turn, determine the characteristics of the electron localization and the subsequent detachment following photoexcitation of the probe.

Descripción: Our previously developed analytical infinite order sudden (IOS) quantum theory of triatomic photodissociation is extended to describe indirect photodissociation processes through a real or virtual intermediate state. The theory uses the IOS approximation for the dynamics in the final dissociative channels and an Airy function approximation for the continuum states. These approximations enable us to evaluate the multi-dimensional non-separable transition amplitudes analytically (as one-dimensional quadratures), despite the different natural coordinates for the initial bound, the intermediate resonant, and the final dissociative states. The fragment internal energy distributions are described as a function of the initial and final quantum states and the photon excitation energy. The theory readily permits the evaluation of rotational distributions for high values of the total angular momentum J in the initial bound molecular state, a feature that would be very difficult with close-coupled methods. In paper II we apply the theory to describe the photofragment yield spectrum of NOCl in the region of the T1(13A″)→S0(11A′) transition. © 1997 American Institute of Physics.

Descripción: The resolution of molecular electric dipole polarizabilities into atomic contributions was studied in the context of methyl and ethyl derivatives, CH3-X and CH3-CH2-X. The reliability of the partitioning schemes allowing for adhoc quantum chemical tools was also tested. Results depicted a good degree of the transferability from molecule to molecule as shown by the evaluation of gross atomic isotropic contributions for carbon, hydrogen, and heteroatoms.

Descripción: Continuum solvent models have become a standard technique in the context of electronic structure calculations, yet no implementations have been reported capable to perform molecular dynamics at solid-liquid interfaces. We propose here such a continuum approach in a density functional theory framework using plane-wave basis sets and periodic boundary conditions. Our work stems from a recent model designed for Car-Parrinello simulations of quantum solutes in a dielectric medium [D. A. Scherlis, J. Chem. Phys. 124, 074103 (2006)], for which the permittivity of the solvent is defined as a function of the electronic density of the solute. This strategy turns out to be inadequate for systems extended in two dimensions: the dependence of the dielectric function on the electronic density introduces a new term in the Kohn-Sham potential, which becomes unphysically large at the interfacial region, seriously affecting the convergence of the self-consistent calculations. If the dielectric medium is properly redefined as a function of the atomic coordinates, a good convergence is obtained and the constant of motion is conserved during the molecular dynamics simulations. The Poisson problem is solved using a multigrid method, and in this way Car-Parrinello molecular dynamics simulations of solid-liquid interfaces can be performed at a very moderate computational cost. This scheme is employed to investigate the acid-base equilibrium at the TiO2 -water interface. The aqueous behavior of titania surfaces has stimulated a large amount of experimental research, but many open questions remain concerning the molecular mechanisms determining the chemistry of the interface. Here we make an attempt to answer some of them, putting to the test our continuum model. © 2009 American Institute of Physics.

Descripción: Energetics, structural features, polarity, and melting transitions in water clusters containing up to eight molecules were studied using ab initio methods and empirical force field models. Our quantum approach was based on density functional theory performed at the generalized gradient approximation level. For the specific case of (H2O)6, we selected five conformers of similar energy with different geometries and dipolar moments. For these cases, the cyclic arrangement was found to be the only nonpolar aggregate. For (H2O)8, the most stable structures corresponded to nonpolar, cubic-like, D2d and S4 conformers. Higher energy aggregates exhibit a large spectrum in their polarities. The static polarizability was found to be proportional to the size of the aggregates and presents a weak dependence with the number of hydrogen bonds. In order to examine the influence of thermal fluctuations on the aggregates, we have performed a series of classical molecular dynamics experiments from low temperature up to the melting transition using two different effective pseudopotentials: the TIP4P and MCY models. Minimum energy structures for both classical potentials were found to reproduce reasonably well the results obtained using ab initio methods. Isomerization and phase transitions were monitored by following changes in dipole moments, number of hydrogen bonds and Lindemann's parameter. For (H2O)6 and (H2O)8, the melting transitions were found at Tm≈50 and 160 K, respectively; for both aggregates, we observed premelting transitions between well differentiated conformers as well. © 1999 American Institute of Physics.

Descripción: We present results from molecular dynamics simulations performed on reverse micelles immersed in cyclohexane. Three different inner polar phases are considered: water (W), formamide (FM), and an equimolar mixture of the two solvents. In all cases, the surfactant was sodium bis(2-ethylhexyl) sulfosuccinate (usually known as AOT). The initial radii of the micelles were R∼15 Å, while the corresponding polar solvent-to-surfactant molar ratios were intermediate between w0 =4.3 for FM and w0 =7 for W. The resulting overall shapes of the micelles resemble distorted ellipsoids, with average eccentricities of the order of ∼0.75. Moreover, the pattern of the surfactant layer separating the inner pool from the non-polar phase looks highly irregular, with a roughness characterized by length scales comparable to the micelle radii. Solvent dipole orientation polarization along radial directions exhibit steady growths as one moves from central positions toward head group locations. Local density correlations within the micelles indicate preferential solvation of sodium ionic species by water, in contrast to the behavior found in bulk equimolar mixtures. Still, a sizable fraction of ∼90% of Na+ remains associated with the head groups. Compared to bulk results, the translational and rotational modes of the confined solvents exhibit important retardations, most notably those operated in rotational motions where the characteristic time scales may be up to 50 times larger. Modifications of the intramolecular connectivity expressed in terms of the average number of hydrogen bonds and their lifetimes are also discussed. © 2008 American Institute of Physics.

Descripción: A computational scheme was developed for a fourth-rank hyprmagnetizability tensor denoted by Xαβγδ. It was shown that this instrinsic property of diamagnetic atoms and molecules was useful to rationalize their nonlinear response to intense magnetic field. The terms connected with the fourth power of the perturbing field which represent the fourth-rank hypermagnetizability of the H2, HF, H2O, NH3, and CH4 molecules were evaluated at the coupled Hartree-Fock level of accuracy. The Gaugeless basis sets of increasing size and flexibility that adopted two different coordinate systems to estimate the degree of convergence of theoretical tensor components was employed.

Descripción: Vertical excitation and electron detachment energies associated with the optical absorption of iodide ions dissolved in supercritical ammonia at 420 K have been calculated in two limiting scenarios: as a solvated free I- ion and forming a K+ I- contact ion pair (CIP). The evolution of the transition energies as a result of the gradual building up of the solvation structure was studied for each absorbing species as the solvent's density increased, i.e., changing the N H3 supercritical thermodynamic state. In both cases, if the solvent density is sufficiently high, photon absorption produces a spatially extended electron charge beyond the volume occupied by the solvated solute core; this excited state resembles a typical charge-transfer-to-solvent (CTTS) state. A combination of classical molecular dynamics simulations followed by quantum mechanical calculations for the ground, first-excited, and electron-detached electronic states have been carried out for the system consisting of one donor species (free I- ion or K+ I- CIP) surrounded by ammonia molecules. Vertical excitation and electron detachment energies were obtained by averaging 100 randomly chosen microconfigurations along the molecular dynamics trajectory computed for each thermodynamic condition (fluid density). Short- and long-range contributions of the solvent-donor interaction upon the CTTS states of I- and K+ I- were identified by performing additional electronic structure calculations where only the solvent interaction due to the first neighbor molecules was taken into account. These computations, together with previous experimental evidence that we collected for the system, have been used to analyze the solvent effects on the CTTS transition. In this paper we have established the following: (i) the CTTS electron of free I- ion or K+ I- CIP presents similar features, and it gradually localizes in close proximity of the iodine parent atom when the ammonia density is increased; (ii) for the free I- ion, the short-range solvent interaction contributes to the stabilization of the ground state more than it does for the CTTS excited state, which is evidenced experimentally as a blueshift in the maximum absorption of the CTTS transition when the density is increased; (iii) this effect is less noticeable for the K+ I- ion pair, because in this case a tight solvation structure, formed by four N H3 molecules wedged between the ions, appears at very low density and is very little affected by changes in the density; (iv) the long-range contribution to the solvent stabilization can be neglected for the K+ I- CIP, since the main features of its electronic transition can be explained on the basis of the vicinity of the cation; (v) however, the long-range solvent field contribution is essential for the free I- ion to become an efficient CTTS donor upon photoexcitation, and this establishes a difference in the CTTS behavior of I- in bulk and in clusters. © 2007 American Institute of Physics.

Descripción: CO recombination kinetics has been investigated in the type II truncated hemoglobin from Thermobifida fusca (Tf-trHb) over more than 10 time decades (from 1 ps to ~100 ms) by combining femtosecond transient absorption, nanosecond laser flash photolysis and optoacoustic spectroscopy. Photolysis is followed by a rapid geminate recombination with a time constant of ~2 ns representing almost 60% of the overall reaction. An additional, small amplitude geminate recombination was identified at ~100 ns. Finally, CO pressure dependent measurements brought out the presence of two transient species in the second order rebinding phase, with time constants ranging from ~3 to ~100 ms. The available experimental evidence suggests that the two transients are due to the presence of two conformations which do not interconvert within the time frame of the experiment. Computational studies revealed that the plasticity of protein structure is able to define a branched pathway connecting the ligand binding site and the solvent. This allowed to build a kinetic model capable of describing the complete time course of the CO rebinding kinetics to Tf-trHb. © 2012 Marcelli et al.