PubMed Health⌕ Search

Biomedical subjects

Peter J Rossky

Publications and source records attributed to Peter J Rossky.

14 recordsLinked to original sources

Quantum simulation of solution phase intramolecular electron transfer rates in betaine-30.

Mixed quantum-classical atomistic simulations have been carried out to investigate the mechanistic details of excited state intramolecular electron transfer in a betaine-30 molecule in acetonitrile. The key electronic degrees of freedom of the solute molecule are treated quantum mechanically using the semiempirical Pariser-Parr-Pople Hamiltonian, including the solvent influence on electronic structure. The intramolecular vibrational modes are also treated explicitly at a quantum level, with the remaining elements treated classically using empirical potentials. The electron-transfer rate, corresponding to S1 --> S0 relaxation, is evaluated via time-dependent perturbation theory with the explicit inclusion of the dynamics of solvation and intramolecular conformation. The calculations reveal that, while solvation dynamics is critical to the rate, the intramolecular torsional dynamics also plays an important role. The importance of the use of multiple high-frequency quantum modes is also discussed.

Betaine↗

Evaluation of quantum correlation functions from classical data: Anharmonic models.

The previously introduced method of evaluating quantum mechanical time correlation functions using as input only classical simulation data is generalized and applied to two anharmonic model systems, as a further test. The quantum correction approach utilizes the relation between a general quantum correlation function and its classical analog. For the tested models, we obtain numerical results of nonlinear correlation functions with comparable accuracy to that of the centroid molecular dynamics method, although the present method is much simpler to implement and not limited to real valued quantum correlation functions.

Journal Article↗

Quantized time correlation function approach to nonadiabatic decay rates in condensed phase: application to solvated electrons in water and methanol.

A new, alternative form of the golden rule formula defining the nonadiabatic transition rate between two quantum states in condensed phase is presented. The formula involves the quantum time correlation function of the energy gap, of the nonadiabatic coupling, and their cross terms. Those quantities can be inferred from their classical counterparts, determined via molecular dynamics simulations. The formalism is applied to the problem of the nonadiabatic p-->s relaxation of an equilibrated p-electron in water and methanol. We find that, in both solvents, the relaxation is induced by the coupling to the vibrational modes and the quantum effects modify the rate by a factor of 2-10 depending on the quantization procedure applied. The resulting p-state lifetime for a hypothetical equilibrium excited state appears extremely short, in the sub-100 fs regime. Although this result is in contrast with all previous theoretical predictions, we also illustrate that the lifetimes computed here are very sensitive to the simulated electronic quantum gap and to the strongly correlated nonadiabatic coupling.

Journal Article↗

Excess electron localization sites in neutral water clusters.

We present approximate pseudopotential quantum-mechanical calculations of the excess electron states of equilibrated neutral water clusters sampled by classical molecular dynamics simulations. The internal energy of the clusters are representative of those present at temperatures of 200 and 300 K. Correlated electronic structure calculations are used to validate the pseudopotential for this purpose. We find that the neutral clusters support localized, bound excess electron ground states in about 50% of the configurations for the smallest cluster size studied (n = 20), and in almost all configurations for larger clusters (n > 66). The state is always exterior to the molecular frame, forming typically a diffuse surface state. Both cluster size and temperature dependence of energetic and structural properties of the clusters and the electron distribution are explored. We show that the stabilization of the electron is strongly correlated with the preexisting instantaneous dipole moment of the neutral clusters, and its ground state energy is reflected in the electronic radius. The findings are consistent with electron attachment via an initial surface state. The hypothetical spectral dynamics following such attachment is also discussed.

Journal Article↗

A comparative study of imaginary time path integral based methods for quantum dynamics.

The recently introduced approximate many-body quantum simulation method, ring polymer molecular dynamics (RPMD), is compared to the centroid molecular dynamics method (CMD). Comparisons of simulation results for liquid para-hydrogen at two state points and liquid ortho-deuterium at one state point are presented. The calculated quantum correlation functions for the two methods are shown to be in good agreement with one another for a large portion of the time spectrum. However, as the quantum mechanical nature of the system increases, RPMD is less accurate in predicting the kinetic energy of the system than is CMD. A simplified and highly efficient algorithm is proposed which largely corrects this deficiency.

Journal Article↗

Effect of pressure on the phase behavior and structure of water confined between nanoscale hydrophobic and hydrophilic plates.

We perform systematic molecular dynamics simulations of water confined between two nanoscale plates at T = 300K. We investigate the effect of pressure (-0.15 GPa< or = P< or =0.2GPa) and plate separation (0.4 nm < or =d < or =1.6 nm) on the phase behavior of water when the plates are either hydrophobic or hydrophilic. When water is confined between hydrophobic plates, capillary evaporation occurs between the plates at low enough P. The threshold value of d at which this transition occurs decreases with P (e.g., 1.6 nm at P approximately equal to -0.05 GPa, 0.5 nm at P approximately equal to 0.1 GPa), until, at high P, no capillary evaporation occurs. For d approximately equal to 0.6 nm and P > or =0.1 GPa, the system crystallizes into a bilayer ice. A P-d phase diagram showing the vapor, liquid, and bilayer ice phases is proposed. When water is confined by hydrophilic (hydroxylated silica) plates, it remains in the liquid phase at all P and d studied. Interestingly, we observe for this case that even at the P at which bulk water cavitates, the confined water remains in the liquid state. We also study systematically the state of hydration at different P for both kinds of plates. For the range of conditions studied here, we find that in the presence of hydrophobic plates the effect of P is to enhance water structure and to push water molecules toward the plates. The average orientation of water molecules next to the hydrophobic plates does not change upon pressurization. In contrast, in the presence of hydrophilic plates, water structure is insensitive to P. Hence, our results suggest that upon pressurization, hydrophobic plates behave as "soft" surfaces (in the sense of accommodating pressure-dependent changes in water structure) while hydrophilic walls behave as "hard" surfaces.

Journal Article↗

Origins of nonexponential decay in single molecule measurements of rotational dynamics.

Recent reports have demonstrated that the correlation function of the fluorescence dichroism signal, measured as a probe of single molecule rotational dynamics, should not manifest a single exponential decay even for isotropic diffusion. This has called into question the attribution of observed nonexponential behavior in supercooled fluids and polymer systems to dynamical heterogeneity. We show here that, for the case of a high numerical aperture objective, the dichroism decay becomes indistinguishable from a single exponential. As a consequence, observed nonexponential decays can be associated with complex rotational dynamics. These effects are illustrated via simulated rotational trajectories for isotropic diffusion of a dipole.

Journal Article↗

Characterization of excess electrons in water-cluster anions by quantum simulations.

Water-cluster anions can serve as a bridge to understand the transition from gaseous species to the bulk hydrated electron. However, debate continues regarding how the excess electron is bound in (H2O)-n, as an interior, bulklike, or surface electronic state. To address the uncertainty, the properties of (H2O)-n clusters with 20 to 200 water molecules have been evaluated by mixed quantum-classical simulations. The theory reproduces every observed energetic, spectral, and structural trend with cluster size that is seen in experimental photoelectron and optical absorption spectra. More important, surface states and interior states each manifest a characteristic signature in the simulation data. The results strongly support assignment of surface-bound electronic states to the water-cluster anions in published experimental studies thus far.

Journal Article↗

Static and dynamic quantum effects in molecular liquids: a linearized path integral description of water.

Structure, transport properties, and IR spectra including quantum effects are calculated for a flexible simple point charge model of liquid water. A recently introduced combination of a variational local harmonic description of the liquid potential surface and the classical Wigner approximation for the dynamics is used. The potential energy and interatomic radial distribution functions are in good agreement with accurate results from the literature and are significantly closer to experiment than predictions found from classical theory. The oxygen and hydrogen velocity correlation functions are also calculated, and the corresponding molecular diffusion coefficient is in good accord with existing theoretical estimates including quantum effects. Of most interest, an ab initio quantum correction factor is obtained to correct the far IR spectrum of water. When corrected, a spectrum based on a classical simulation yields results that agree well with experiment. Combined with internal tests of consistency, these observations indicate that this quite flexible approach will be effective for a variety of molecular problems involving the dynamics of light nuclei.

Biophysics↗

An analysis of electronic dephasing in the spin-boson model.

In order to develop a more complete understanding of the limitations of mixed quantum-classical simulation methods, the origins of electronic dephasing are analyzed in a simple model of the condensed phase, namely, the spin-boson model with an ohmic spectral density. We focus on the decay of the thermally averaged nuclear overlap/phase function (NOPF). Considering the strong coupling/high temperature limit, a relationship is obtained at short time between the rate of electronic coherence loss and the electronic dephasing rate characteristic of a classical bath. Using this relationship, we clarify the origin of the decay of the NOPF. In the same limit, we also reproduce an earlier relationship between the electronic decoherence time and a solvation relaxation time. Finally, we point out that, for the spin-boson model, the exact quantum mechanical description of electronic dephasing is reproduced by mixed quantum/classical methods if a Gaussian distribution of quantum fluctuations around each classical phase space point is introduced. That spatial distribution of quantum fluctuations is functionally the same as that appearing in the Feynman-Kleinert variational local harmonic approximation, and also that implemented in existing classical trajectory-based estimates of coherence dissipation times.

Journal Article↗

Critical evaluation of approximate quantum decoherence rates for an electronic transition in methanol solution.

We present a quantum molecular dynamics calculation of a semiclassical decoherence function to evaluate the accuracy of alternative short-time approximations for coherence loss in the dynamics of condensed phase electronically nonadiabatic processes. The semiclassical function from mixed quantum-classical molecular dynamics simulations and frozen Gaussian wave packets is computed for the electronic transition of an excited state excess electron to the ground state in liquid methanol. The decoherence function decays on a 10 fs time scale that is qualitatively similar to the aqueous case. We demonstrate that it is the motion of the hydrogen atom, and, in particular, the hydrogen rotation around the oxygen-methyl bond which is predominantly responsible for destroying the quantum correlations between alternative states. Multiple time scales due to the slower diffusive nuclear modes, which dominate the solvation response of methanol, do not contribute to the coherence loss. The choice of the coordinate representation is investigated in detail and concluded to be irrelevant to the decay. Changes in both nuclear momenta and positions on the two alternative potential surfaces are found to contribute to decoherence, the former dominating at short times (t < 5 fs), the latter controlling the decay at longer times. Various short-time approximations to the full dynamics for the decoherence function are tested for the first time. The present treatment rigorously develops the short-time description and establishes its range of validity. Whereas the lowest-order short-time approximation proves to be a very good approximation up to about 5 fs, we also find that it bounds the decay of the decoherence function. After 5 fs, the coherence decay in fact becomes faster than the single Gaussian predicted in the lowest-order short-time limit. This decay is well reflected by an enhanced low-order approximation, which is also easily computed from equilibrium classical forces.

Journal Article↗