Michele Ceotto

Quantum instanton approximation for thermal rate constants of chemical reactions

A quantum mechanical theory for chemical reaction rates is presented which is modeled after the [semiclassical (SC)] instanton approximation. It incorporates the desirable aspects of the instanton picture, which involves only properties of the (SC approximation to the) Boltzmann operator, but corrects its quantitative deficiencies by replacing the SC approximation for the Boltzmann operator by the quantum Boltzmann operator, exp(−βĤ). Since a calculation of the quantum Boltzmann operator is feasible for quite complex molecular systems (by Monte Carlo path integral methods), having an accurate rate theory that involves only the Boltzmann operator could be quite useful. The application of this quantum instanton approximation to several one- and two-dimensional model problems illustrates its potential; e.g., it is able to describe thermal rate constants accurately (∼10–20% error) from high to low temperatures deep in the tunneling regime, and applies equally well to asymmetric and symmetric potentials.

Multiple Coherent States for First-Principles Semiclassical Initial Value Representation Molecular Dynamics

A multiple coherent states implementation of the semiclassical approximation is introduced and employed to obtain the power spectra with a few classical trajectories. The method is integrated with the time-averaging semiclassical initial value representation to successfully reproduce anharmonicity and Fermi resonance splittings at a level of accuracy comparable to semiclassical simulations of thousands of trajectories. The method is tested on two different model systems with analytical potentials and implemented in conjunction with the first-principles molecular dynamics scheme to obtain the power spectrum for the carbon dioxide molecule.

Wettability of bare and fluorinated silanes: A combined approach based on surface free energy evaluations and dipole moment calculations

The assessment of the surface free energy (SFE) of a material permits to control and predict a large number of physicochemical properties of a solid surface and its reactivity. Here, the surface energies of a series of bare and fluorinated silanes are determined by means of different semi-empirical models on the grounds of contact angle determinations for different solvents. Literature data are also considered in order to produce a series of films with increasing SFE. Wetting envelopes (WE) are obtained for the various surfaces in order to predict the wettability of the films by numerous commonly employed solvents. The polar and disperse components of the surface energies are obtained by the Owens-Wendt-Rabel-Kaelbe (OWRK) method; the values of the polar components are compared with gas-phase dipole moments obtained by theoretical calculations employing semi-empirical Hamiltonians. The sequences of the polar components of the SFE and of the calculated dipole moments for the different molecules are strictly the same. The interplay between theoretical and experimental approaches proves efficient in predicting the behavior of different systems and it could be employed in tuning the SFE of a solid surface with respect to its final applications.

Multiple coherent states semiclassical initial value representation spectra calculations of lateral interactions for CO on Cu(100)

Lateral interactions between carbon monoxide molecules adsorbed on a copper Cu(100) surface are investigated via semiclassical initial value representation (SC-IVR) molecular dynamics. A previous analytical potential is extended to include long-range dipole interactions between coadsorbed molecules and preliminary classical simulations were performed to tune the potential parameters. Then, the spectra for several coadsorbed molecules are calculated using the multiple coherent states approximation of the time-averaging representation of the SC-IVR propagator. Results show strong resonances between coadsorbed molecules as observed by past experiments. Resonances turn into dephasing when isotopical substitutions are performed.

Evaluating the Accuracy of Hessian Approximations for Direct Dynamics Simulations.

Direct dynamics simulations are a very useful and general approach for studying the atomistic properties of complex chemical systems, since an electronic structure theory representation of a systems potential energy surface is possible without the need for fitting an analytic potential energy function. In this paper, recently introduced compact finite difference (CFD) schemes for approximating the Hessian [J. Chem. Phys.2010, 133, 074101] are tested by employing the monodromy matrix equations of motion. Several systems, including carbon dioxide and benzene, are simulated, using both analytic potential energy surfaces and on-the-fly direct dynamics. The results show, depending on the molecular system, that electronic structure theory Hessian direct dynamics can be accelerated up to 2 orders of magnitude. The CFD approximation is found to be robust enough to deal with chaotic motion, concomitant with floppy and stiff mode dynamics, Fermi resonances, and other kinds of molecular couplings. Finally, the CFD approximations allow parametrical tuning of different CFD parameters to attain the best possible accuracy for different molecular systems. Thus, a direct dynamics simulation requiring the Hessian at every integration step may be replaced with an approximate Hessian updating by tuning the appropriate accuracy.

First Principles Semiclassical Calculations of Vibrational Eigenfunctions

Vibrational eigenfunctions are calculated on-the-fly using semiclassical methods in conjunction with ab initio density functional theory classical trajectories. Various semiclassical approximations based on the time-dependent representation of the eigenfunctions are tested on an analytical potential describing the chemisorption of CO on Cu(100). Then, first principles semiclassical vibrational eigenfunctions are calculated for the CO(2) molecule and its accuracy evaluated. The multiple coherent states initial value representations semiclassical method recently developed by us has shown with only six ab initio trajectories to evaluate eigenvalues and eigenfunctions at the accuracy level of thousands trajectory semiclassical initial value representation simulations.

A Close Look at the Structure of the TiO2-APTES Interface in Hybrid Nanomaterials and Its Degradation Pathway: An Experimental and Theoretical Study

The surface functionalization of TiO2-based materials with alkylsilanes is attractive in several cutting-edge applications, such as photovoltaics, sensors, and nanocarriers for the controlled release of bioactive molecules. (3-Aminopropyl)triethoxysilane (APTES) is able to self-assemble to form monolayers on TiO2 surfaces, but its adsorption geometry and solar-induced photodegradation pathways are not well understood. We here employ advanced experimental (XPS, NEXAFS, AFM, HR-TEM, and FT-IR) and theoretical (plane-wave DFT) tools to investigate the preferential interaction mode of APTES on anatase TiO2. We demonstrate that monomeric APTES chemisorption should proceed through covalent Si–O–Ti bonds. Although dimerization of the silane through Si–O–Si bonds is possible, further polymerization on the surface is scarcely probable. Terminal amino groups are expected to be partially involved in strong charge-assisted hydrogen bonds with surface hydroxyl groups of TiO2, resulting in a reduced propensity to react with other species. Solar-induced mineralization proceeds through preferential cleavage of the alkyl groups, leading to the rapid loss of the terminal NH2 moieties, whereas the Si-bearing head of APTES undergoes slower oxidation and remains bound to the surface. The suitability of employing the silane as a linker with other chemical species is discussed in the context of controlled degradation of APTES monolayers for drug release and surface patterning.

Graphics processing units accelerated semiclassical initial value representation molecular dynamics

This paper presents a Graphics Processing Units (GPUs) implementation of the Semiclassical Initial Value Representation (SC-IVR) propagator for vibrational molecular spectroscopy calculations. The time-averaging formulation of the SC-IVR for power spectrum calculations is employed. Details about the GPU implementation of the semiclassical code are provided. Four molecules with an increasing number of atoms are considered and the GPU-calculated vibrational frequencies perfectly match the benchmark values. The computational time scaling of two GPUs (NVIDIA Tesla C2075 and Kepler K20), respectively, versus two CPUs (Intel Core i5 and Intel Xeon E5-2687W) and the critical issues related to the GPU implementation are discussed. The resulting reduction in computational time and power consumption is significant and semiclassical GPU calculations are shown to be environment friendly.

Mixed semiclassical initial value representation time-averaging propagator for spectroscopic calculations

A mixed semiclassical initial value representation expression for spectroscopic calculations is derived. The formulation takes advantage of the time-averaging filtering and the hierarchical properties of different trajectory based propagation methods. A separable approximation is then introduced that greatly reduces (about an order of magnitude) the computational cost compared with a full Herman-Kluk time-averaging semiclassical calculation for the same systems. The expression is exact for the harmonic case and it is tested numerically for a Morse potential coupled to one or two additional harmonic degrees of freedom. Results are compared to full Herman-Kluk time-averaging calculations and exact quantum wavepacket propagations. We found the peak positions of the mixed semiclassical approximations to be always in very good agreement with full quantum calculations, while overtone peak intensities are lower with respect to the exact ones. Given the reduced computational effort required by this new mixed semiclassical approximation, we believe the present method to make spectroscopic calculations available for higher dimensional systems than accessible before.

Helium Isotope Enrichment by Resonant Tunneling through Nanoporous Graphene Bilayers

Graphene is impermeable to gases, but introducing subnanometer pores can allow for selective gas separation. Because graphene is only one atom thick, tunneling can play an important role, especially for low-mass gases such as helium, and this has been proposed as a means of separating (3)He from (4)He. In this paper, we consider the possibility of utilizing resonant tunneling of helium isotopes through nanoporous graphene bilayers. Using a model potential fit to previously reported DFT potential energy surfaces, we calculate the thermal rate constant as a function of interlayer separation using a recently described time-independent method for arbitrary multibarrier potentials. Resonant transmission allows for the total flux rate of (3)He to remain the same as the best-known single-barrier pores but doubles the selectivity with respect to (4)He when the optimal interlayer spacing of 4.6 Å is used. The high flux rate and selectivity are robust against variations of the interlayer spacing and asymmetries in the potential that may occur in experiment.