BSc, Laurentian University
PhD, University of Toronto
My research program can be divided into three main areas:
1. Mixed quantum-classical formulation of nonlinear spectroscopy
(a) Simulating the quantum dynamics of processes occurring in classical condensed phase hosts
Starting from the mixed quantum-classical Liouville (MQCL) equation of motion, it is possible to construct practical trajectory-based algorithms for simulating the dynamics of a quantum subsystem coupled to a classical environment. Over the past few years, several methods have been successfully applied in the treatment of simple model systems, but their implementation has been less straightforward in the case of more complex systems. It is our desire to develop and implement schemes based on the MQCL equation, which overcome the pitfalls of existing algorithms when dealing with complex, many-body systems.
(b) Modeling of multidimensional nonlinear spectra
Ultrafast nonlinear spectroscopy is capable of probing molecular dynamics on the femtosecond time scale. Often the resulting spectra are complex, thereby requiring a theoretical framework for their interpretation. Moreover, with comparisons between simulated and experimentally measured spectra becoming increasingly sophisticated, the development and implementation of accurate methods for modeling spectroscopic response is timely. The MQCL approach provides a convenient way for simulating laser-driven dynamics and will thus provide a suitable platform for the development of a general framework for calculating multidimensional nonlinear spectra.
2. Applications to the multidimensional infrared spectroscopy of chemical and biological systems in nanoconfined environments
We are interested in simulating one- and two-dimensional infrared spectra to study the structure and dynamics of a wide variety of nanoconfined systems of experimental interest, some examples of which include:
(a) Nanoconfined water
In many chemical and biological systems, water molecules can be confined on nanometer length scales. Under these conditions, the molecules are in contact with different types of interfaces. Near an interface, the hydrogen bonding network of water changes considerably because it must adjust to the shape of that interface. As a result, the properties and dynamics of nanoconfined water differ substantially from those of bulk water and must therefore be studied in the presence of the interface.
(b) Nanoconfined nonaqueous polar liquid clusters
Proton transfer in nanoconfined nonaqueous polar liquid clusters represents a class of reactions that are ubiquitous in chemistry. This charge transfer reaction is strongly coupled to the polar solvent and will therefore be greatly affected by solvent confinement. For example, several experiments have shown that the proton transfer rate constant can decrease significantly upon confinement. As a result of this sensitivity, one may design materials with specific chemical purposes by simply varying properties of the cluster such as its size and shape.
(c) Hydrogen transfer in enzymatic catalysis
Hydrogen transfer reactions are ubiquitous in enzymatic catalysis. The interior of an enzyme can form a nanoconfined environment around its active site and this confinement may play an important role in its function. Studying the effects of factors such as hydrogen tunneling and enzymatic motions is crucial for a detailed understanding of the transfer mechanism.
3. Ab-initio molecular dynamics of rare chemical events in condensed phase environments
We are interested in investigating the microscopic mechanisms, energetics, and kinetics of high-barrier chemical reactions and other rare events in both liquids and solids with the aid of ab-initio molecular dynamics and metadynamics. Examples of such processes under study in our group include:
(a) Dissociation and decomposition of carbonic acid in water
The reactions of carbonic acid (H_{2}CO_{3}) in water are important in many chemical and biological processes such as those involved in the pH regulation in blood, CO_{2} transport in the lungs, and the global carbon cycle. H_{2}CO_{3}undergoes dissociation (H_{2}CO_{3} --> HCO_{3}^{-} + H^{+}) in water, which is followed by decomposition either via the hydroxide route (HCO_{3}- --> CO_{2} + OH^{-}) or the water route (HCO_{3}^{-} + H_{3}O^{+} --> CO_{2} + 2H_{2}O). Hydrogen bonding between H_{2}CO_{3} and neighbouring water molecules plays an important role in determining the mechanisms, energetics, and kinetics of these reactions, and therefore an explicit and accurate treatment of the solvation is needed.
(b) Absorption of CO_{2} in aqueous solutions of NH3 and alkanolamines
Amine scrubbing is the leading technology currently used for capturing CO_{2} from flue gases in power plants. Typically, ammonia and alkanolamines (R-NH_{2}) are used, which react with CO_{2} in aqueous solution according to the reaction CO_{2} + 2R-NH_{2} --> RNHCOO^{-} + RNH_{3}^{+}. Several mechanisms for this reaction have been proposed, but they are still speculative. Therefore, further theoretical studies are required since a detailed understanding of the mechanisms of aqueous CO_{2}/R-NH_{2} systems is important for the development and improvement of amine scrubbing technologies.
(c) Polymorphism in solid carbonic acid
Solid carbonic acid ((H_{2}CO_{3})_{n}) has been recently proposed as a potential candidate for CO_{2} capture. Experimental studies have shown that it can exist in two different polymorphs (i.e., α- and β-carbonic acid), but their exact crystal structures are still unknown. Theoretical studies can aid in the search for likely crystal structures. Moreover, in order to assess the feasibility of storing CO_{2} as solid H_{2}CO_{3}, the physical and chemical stability of the various polymorphs of solid H_{2}CO_{3} must be investigated under different conditions.An introduction to the quantum view of nature with applications to atomic and molecular structure. Methods to describe the quantum world are introduced, used to describe simple electronic, vibrational and rotational structure of model systems, and applied to the hydrogen atom, many-electron atoms, simple diatomic molecules, and the electronic structure of polyatomic molecules. The laboratory portion of the course consists of practical applications enriching and illustrating the lecture material, and incorporates the use of computers as a routine aid to processing experimental results. Prerequisites: CHEM 102 or 105; one 200-level CHEM course; MATH 115 or 136 or 146 and PHYS 124 or 144. Corequisite: PHYS 146 if PHYS 144 presented as a prerequisite instead of PHYS 124.
Winter Term 2021The fundamentals of statistical mechanics are covered to set up the theoretical framework for Molecular Dynamics (MD) simulation. The basic components of MD simulation are discussed in detail, followed by a brief foray into Monte Carlo simulation. A variety of applications are presented, including the study of structural properties of liquids, the calculation of diffusion coefficients for a solute in a solvent, and the calculation of reaction rate constants. A brief overview of methods for incorporating quantum effects into MD simulations is given. Computational exercises will be assigned to exemplify various topics encountered in the lectures. Prerequisite: CHEM 282 and CHEM 371; or consent of the instructor.
Fall Term 2020The fundamentals of statistical mechanics are covered to set up the theoretical framework for Molecular Dynamics (MD) simulation. The basic components of MD simulation are discussed in detail, followed by a brief foray into Monte Carlo simulation. A variety of applications are presented, including the study of structural properties of liquids, the calculation of diffusion coefficients for a solute in a solvent, and the calculation of reaction rate constants. A brief overview of methods for incorporating quantum effects into MD simulations is given. Computational exercises will be assigned to exemplify various topics encountered in the lectures. Not open to students with credit in CHEM 495.
Fall Term 2020