Femtosecond quantum control of molecular dynamics in the condensed phase

We review the progress in controlling quantum dynamical processes in the condensed phase with femtosecond laser pulses. Due to its high particle density the condensed phase has both high relevance and appeal for chemical synthesis. Thus, in recent years different methods have been developed to manipulate the dynamics of condensed-phase systems by changing one or multiple laser pulse parameters. Single-parameter control is often achieved by variation of the excitation pulse's wavelength, its linear chirp or its temporal subpulse separation in case of pulse sequences. Multiparameter control schemes are more flexible and provide a much larger parameter space for an optimal solution. This is realized in adaptive femtosecond quantum control, in which the optimal solution is iteratively obtained through the combination of an experimental feedback signal and an automated learning algorithm. Several experiments are presented that illustrate the different control concepts and highlight their broad applicability. These fascinating achievements show the continuous progress on the way towards the control of complex quantum reactions in the condensed phase.     

Conformation and dynamics of polymer chains in condensed systems

A Stockmayer's type-A chain such as cis-polyisoprene (PI) having dipole components parallel to the chain contour exhibits slow dielectric response, reflecting fluctuation of its end-to-end vector ( dielectric normal modes ). This article reviews our own work on dielectric spectroscopy for studying dynamics and conformational properties of probe PI chains in binary blends with PI or polybutadiene (PB) that does not exhibit normal modes. Through such studies we were able to see the crossover of the probe-chain dynamics from Rouse behaviour in low-molecular-weight (lmw) matrices to pure reptation behaviour in high-molecular-weight (hmw) matrices. We were also able to confirm experimentally Flory's earlier prediction on excluded volumes in dilute blends of hmw probe chains in lmw matrices, for which a proof has never been given before.     

Polymorphic phase behavior of cardiolipin derivatives studied by coarse-grained molecular dynamics

Cardiolipin (CL) is a negatively charged four acyl chain lipid, associated with energy production in bacterial and mitochondrial membranes. Due to the shape of CL, negative curvatures of aggregates are favorable if the charges in the head group can be reduced. The phase polymorphism of CL, and of associated derivatives with 2, 3, 4, or 5 chains, has been determined previously and offers a model system in which micellar, lamellar, and inverse hexagonal phases can be observed. We present an extension to a previously established coarse-grained molecular dynamics model with the aim of reproducing the different CL phases with two adjustable parameters: the number of acyl chains and the effective head group charge. With molecular dynamics simulations of large lipid systems, we observed transitions between different phases on the nanosecond to microsecond time scale. Charge screening by high salt or low pH was successfully modeled by a reduction of phosphate charge, which led to the adoption of aggregates with more negative curvature. Although specific ion binding at the interface and other atomistic details are sacrificed in the coarse-grained model, we found that it captures general phase features over a large range of aggregate geometries.     

On a fokker-planck type approach to molecular dynamics in condensed phase

A fast algorithm for dealing with multiplicative stochastic processes such as dielectric relaxation is developed. This technique allows us to get a straightforward evaluation of the effects of non-linearity, expecially non-Gaussian behaviour. Its use jointly to computer “experiments” on non-Gaussian behaviour is shown to be of basic importance for establishing a new approach to molecular dynamics in terms of non-linear models.     

Momentum representation of the solute/bath interaction in the dynamic theory of chemical processes in condensed phase

For the system consisting of the chemically reactive solute immersed in the oscillator bath, we consider an approach based on the solute/medium interaction expressed in terms of momenta rather than coordinates. In the adiabatic representation the medium reorganization effects are suppressed, being hidden in the solute renormalized potential and new spectral density function. The advantage proposed by the bilinear interaction in momentum representation is its spatial uniformity important for approximate dynamical treatments. The procedure of explicit transforming a standard spectral density (coordinate representation of interaction) into the spectral density in adiabatic representation (momentum representation of interaction) is the main new result of the present study. Illustrative calculations for several types of spectral functions are performed. Special discussion is devoted to clarifying the nature of the slow diffusion coordinate, to which the present approach is mainly addressed.     

Vibrational relaxation in the condensed phase

A modified Landau-Teller equation for vibrational relaxation in the condensed phase is proposed. This equation differs from previous approaches by accounting for the fluctuations of the energies of the vibrational levels that result from the interactions with the surroundings (bath). In the conventional approach the effects of the bath are only included in the coupling between the relaxing and accepting vibrational modes. It is shown that the additional inclusion of the fluctuations of the energy levels can lead to a dramatic change of the vibrational relaxation rate.     

Effective force fields for condensed phase systems from ab initio molecular dynamics simulation: a new method for force-matching

A novel least-squares fitting approach is presented to obtain classical force fields from trajectory and force databases produced by ab initio (e.g., Car-Parrinello) molecular dynamics (MD) simulations. The method was applied to derive effective nonpolarizable three-site force fields for liquid water at ambient conditions from Car-Parrinello MD simulations in the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure includes a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The various parameterizations of the water force field differ by an enforced smooth cut-off applied to the short-ranged interaction term. These were obtained by fitting to the trajectory and force data produced by Car-Parrinello MD simulations of systems of 32 and 64 H(2)O molecules. The new water force fields were developed assuming both flexible or rigid molecular geometry. The simulated structural and self-diffusion properties of liquid water using the fitted force fields are in close agreement with those observed in the underlying Car-Parrinello MD simulations. The resulting empirical models compare to experiment much better than many conventional simple point charge (SPC) models. The fitted potential is also shown to combine well with more sophisticated intramolecular potentials. Importantly, the computational cost of the new models is comparable to that for SPC-like potentials.     

Time-resolved multiphoton ionization in the organic condensed phase: Picosecond conformational dynamics of cis-stilbene and tetraphenylethylene

Picosecond time-resolved multiphoton ionization has been extended to studies of photochemical dynamics in the condensed phase. The early time dynamics     

On the direct calculation of the free energy of quantization for molecular systems in the condensed phase

Using the path integral formalism or the Feynman-Hibbs approach, various expressions for the free energy of quantization for a molecular system in the condensed phase can be derived. These lead to alternative methods to directly compute quantization free energies from molecular dynamics computer simulations, which were investigated with an eye to their practical use. For a test system of liquid neon, two methods are shown to be most efficient for a direct evaluation of the excess free energy of quantization. One of them makes use of path integral simulations in combination with a single-step free energy perturbation approach and was previously reported in the literature. The other method employs a Feynman-Hibbs effective Hamiltonian together with the thermodynamic integration formalism. However, both methods are found to give less accurate results for the excess free energy of quantization than the estimate obtained from explicit path integral calculations on the excess free energy of the neon liquid in the classical and quantum mechanical limit. Suggestions are made to make both methods more accurate.     

Kinetics of energy transfer in the condensed phase

The decay kinetics of the donor state is studied; this decay is determined by the resonance energy transfer to the surrounding donor molecules of the medium. A new element is the consideration of the strong interaction of the donor with the acceptor over short distances, thanks to which the onset of decomposition assumes an exponential character. It was shown that this circumstance does not affect the fluorescence quenching, but it can exert a considerable effect on phenomena related to radiation shielding.     

Description of quantum effects in the condensed phase

Recent experiments on light harvesting complexes have shown clear indication of coherent transport of excitations in these aggregates. We review the theoretical models that have been used to study energy transfer in molecular aggregates, beginning with the early models of Förster and Davydov. We cover the Redfield and Haken Strobl models in some detail, in order to set the nomenclature and because they are the most common and easiest models to understand and work with. We briefly discuss more complex models.     

Semiquantal time-dependent Hartree approach to condensed phase chemical dynamics: application to the system-bath model

A semiquantal analysis of condensed phase chemical dynamics, outlined recently for a double-well linearly coupled to dissipative harmonic bath, is formulated in detail to clarify its general features as well as the specifics of the linear and quadratic coupling cases. The theory may be called a "semiquantal time-dependent Hartree (SQTDH)" approach, as it assumes a factorized product of the squeezed coherent state wave packets for the variational subspace of the many-dimensional time-dependent wave function. Due to this assumption, it straightforwardly satisfies the canonicity condition introduced by Marumori et al. is described by a set of Hamilton equations of motion in an extended phase space that includes auxiliary coordinates representing the wave packet widths. The potential in the extended phase space provides a pictorial understanding of the quantum effects affected due to the bath coupling, e.g., suppression of the wave packet spreading in terms of the potential wall developing along the auxiliary coordinates. The idea is illustrated by prototypical models of quartic double-well and cubic metastable potentials linearly and quadratically coupled to the bath. Further applications and extensions, where the SQTDH method will offer a practical approach for introducing quantum effects into realistic molecular dynamics simulations, are also discussed.     

Positive muons in condensed phase ammonia

The chemical behavior of positive muons in condensed phase ammonia has been investigated in order to elucidate the phase and temperature effects on the chemical and physical behavior of the muon and muonium formation in a simple binary compound. Diamagnetic muon yield (P_D) was constant at 0.67±0.01 in both solid and liquid above 125 K. Muonium formation in solids were observed above 100 K with slow muonium spin relaxation. In liquids, the muonium yield and its spin relaxation rate showed temperature dependence. Addition of metallic sodium increased P_D in liquids.     

Radical-olefin addition reactions in condensed phase

A study has been made of the rate constant for the reaction The states of the product have quasilevels E_n — i_n, which, depending on the heat of reaction, may be in resonance with the state of the reactants after surmounting the activation barrier. Because of this, there is no monotonic variation of the propagation rate constant, even as a function of the heat of reaction; the behavior of the rate constant has a resonance character. The results obtained in this investigation are in good agreement with experimental data.Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 26, No. 1, pp. 25–32, January–February, 1990.     

A new perspective on the coarse-grained dynamics of fluids

A computational methodology is presented that is designed to model, at a coarse-grained level, the mesoscale dynamics of fluids and potentially other forms of soft matter. Within a molecular dynamics simulation, "ghost" particles of a specific size, corresponding to the fundamental length-scale of coarse-graining, are used as micro-probes designed to respond to local mesoscale fluid flows and stress gradients. A subsequent coarse-grained model is then developed that incorporates both the coarse-grained mesoscale dynamics and isothermal compressibility of the original microscopic system. The method is applied to water and methanol. A contrast with dissipative particle dynamics (DPD) is also presented.     

Problems in molecular dynamics of condensed phases

A review of the recent theoretical and computational activity at the Chemistry Department of the University of Firenze in the field of molecular simulations of condensed phases is reported. The topics considered include quantitative methods for accurate free energy calculations, molecular dynamics of liquids and ionic solutions, chemical reactions in solutions, phase transformations and polymerization reactions at high pressures.     

Intermolecular force constants of hcn in the condensed phase

A simple theory of linear lattice is applied to the hydrogen bonded linear chain system of HCN to calculate the intermolecular force constants at different temperatures in the condensed phase. The strong CN bond is assumed to remain unperturbed in the hydrogen bond formation. The sharp change in intermolecular force constant while passing from the crystalline to the liquid phase is interpreted as a characteristic of this phase transition (fusion).