Non-adiabatic direct quantum dynamics using force fields: Toward solvation.

Quantum dynamics simulations are becoming a powerful tool for understanding photo-excited molecules. Their poor scaling, however, means that it is hard to study molecules with more than a few atoms accurately, and a major challenge at the moment is the inclusion of the molecular environment. Here, we present a proof of principle for a way to break the two bottlenecks preventing large but accurate simulations. First, the problem of providing the potential energy surfaces for a general system is addressed by parameterizing a standard force field to reproduce the potential surfaces of the molecule's excited-states, including the all-important vibronic coupling. While not shown here, this would trivially enable the use of an explicit solvent. Second, to help the scaling of the nuclear dynamics propagation, a hierarchy of approximations is introduced to the variational multi-configurational Gaussian method that retains the variational quantum wavepacket description of the key quantum degrees of freedom and uses classical trajectories for the remaining in a quantum mechanics/molecular mechanics like approach. The method is referred to as force field quantum dynamics (FF-QD), and a two-state ππ*/nπ* model of uracil, excited to its lowest bright ππ* state, is used as a test case.

Structural, dynamic and photophysical properties of a fluorescent dye incorporated in an amorphous hydrophobic polymer bundle.

The properties of a low molecular weight organic dye, namely 4-naphthyloxy-1-methoxy-2,2,6,6-tetramethylpiperidine, covalently bound to an apolar polyolefin were investigated by means of a multi-level approach, combining classical molecular dynamics simulations, based on purposely parameterized force fields, and quantum mechanical calculations based on density functional theory (DFT) and its time-dependent extension (TD-DFT). The structure and dynamics of the dye in its embedding medium were analyzed and discussed taking the entangling effect of the surrounding polymer into account, and also by comparing the results to those obtained for a different environment, i.e. toluene solution. Finally, the influence was investigated of long lived cages found in the polymeric embedding on photophysical properties, in terms of the slow and fast dye's internal dynamics, by comparing computed IR and UV spectra with their experimental counterparts.

Modeling benzene with single-site potentials from ab initio calculations: a step toward hybrid models of complex molecules.

Extensive ab initio calculations at the MP2/6-31G* level have been carried out to sample the energy surface for the interactions of the benzene dimers. This database has been used to parameterize two anisotropic single-site models, meant to be used as building blocks in hybrid models of complex, liquid crystal forming molecules. A quadrupolar Gay-Berne (GBQIII) and an S-function (SF) Corner potentials have been obtained in this way. Their ability to reproduce, qualitatively at least, the phase diagram as well as energetic and structural properties of benzene has been tested with Monte Carlo simulations and compared with previous literature potentials, GBQI [S. Gupta et al., Mol. Phys. 65, 961 (1988)] and GBQII [T. R. Walsh, Mol. Phys. 100, 2867 (2002)]. It turned out that GBQI showed no melting transition in the temperature range explored (100-400 K), while GBQII underwent a phase transition from solid to gas, with no liquid phase. Conversely, both models parameterized on our database of ab initio interaction energies (GBQIII and SF) gave rise to a stable liquid phase. Melting has been observed between 100 and 150 K (GBQIII) and in the range 300-350 K (SF), i.e., substantially below and slightly above the experimental value at ambient pressure, 278 K. The description of the crystal structure of benzene at atmospheric pressure is also in better agreement with experimental data if the SF model is used, while positional correlations in the liquid are better described by the GBQIII potential. The S-function potential is also computationally more convenient. These results could be useful in the semirealistic modeling of more complex molecules.