M. Waroquier

When adsorbing guest molecules, the porous metal-organic framework MIL-53(Cr) may vary its cell parameters drastically while retaining its crystallinity. A first approach to the thermodynamic analysis of this 'framework breathing' consists of comparing the osmotic potential in two distinct shapes only (large-pore and narrow-pore). In this paper, we propose a generic parametrized free energy model including three contributions: host free energy, guest-guest interactions, and host-guest interaction. Free energy landscapes may now be constructed scanning all shapes and any adsorbed amount of guest molecules. This allows to determine which shapes are the most stable states for arbitrary combinations of experimental control parameters, such as the adsorbing gas chemical potential, the external pressure, and the temperature. The new model correctly reproduces the structural transitions along the CO₂ and CH₄ isotherms. Moreover, our model successfully explains the adsorption versus desorption hysteresis as a consequence of the creation, stabilization, destabilization, and disappearance of a second free energy minimum under the assumptions of a first order phase transition and collective behavior. Our general thermodynamic description allows to decouple the gas chemical potential μ and mechanical pressure P as two independent thermodynamic variables and predict the complete (μ,P) phase diagram for CO₂ adsorption in MIL-53(Cr). The free energy model proposed here is an important step towards a general thermodynamics description of flexible metal-organic frameworks.

A new polarizable force field (PFF), namely atom-condensed Kohn-Sham density functional theory approximated to second order (ACKS2), is proposed for the efficient computation of atomic charges and linear response properties of extended molecular systems. It is derived from Kohn-Sham density functional theory (KS-DFT), making use of two novel ingredients in the context of PFFs: (i) constrained atomic populations and (ii) the Legendre transform of the Kohn-Sham kinetic energy. ACKS2 is essentially an extension of the Electronegativity Equalization Method (EEM) [W. J. Mortier, S. K. Ghosh, and S. Shankar, J. Am. Chem. Soc. 108, 4315 (1986)]10.1021/ja00275a013 in which two major EEM shortcomings are fixed: ACKS2 predicts a linear size-dependence of the dipole polarizability in the macroscopic limit and correctly describes the charge distribution when a molecule dissociates. All ACKS2 parameters are defined as atoms-in-molecules expectation values. The implementation of ACKS2 is very similar to that of EEM, with only a small increase in computational cost.

Radicals generated in trehalose single crystals by X radiation at room temperature were investigated by electron paramagnetic resonance (EPR), electron nuclear double resonance (ENDOR) and ENDOR-induced EPR measurements, together with periodic density functional theory calculations. In the first days after irradiation, three radical species (I1, I2 and I3) were detected, two of which (I1 and I2) dominate the EPR spectrum and could be identified as H-abstracted species centered at C3′ (I1) and C2 (I2), the latter with additional formation of a carbonyl group at C3. Annealing the sample at 40°C for 3 days or storing it in ambient conditions for three months resulted in another, more stable EPR spectrum. Two major species could be characterized in this stage (S1 and S2), only one of which was tentatively identified as an H-abstracted, C2-centered species (S1). Our findings disagree with a previous EPR study [Gräslund and Löfroth (23)] on several accounts. This work stresses the need for caution when interpreting composite EPR spectra and thermally induced spectral changes of radiation-induced species, even in these relatively simple carbohydrates. It also provides further evidence that the pathways for radiation damage critically depend on the specific conformation of a molecule and its environment, but also that carbonyl group formation is a common process in the radiation chemistry of sugars and related compounds.

Structural analysis of modified DNA with NMR is becoming ever more difficult with increasingly complex compounds under scrutiny for use in medical diagnosis, therapeutics, material science and chemical synthesis. To facilitate this process, we develop a molecular modeling approach to predict proton chemical shifts in sufficient agreement with experimental NMR measurements to guide structure elucidation. It relies on a QM/MM partitioning scheme and first principle calculations to predict the spatial structure and calculate corresponding proton chemical shifts. It is shown that molecular dynamics simulations that take into account solvent and temperature effects properly are of utmost importance to sample the conformational space sufficiently. The proposed computational procedure is universally applicable to modified oligonucleotides and DNA, attaining a mean error for the proton chemical shifts of less than 0.2 ppm. Here, it is applied on the Drew-Dickerson d(CGCGAATTCGCG)2 dodecamer as a benchmark system and the mispair-aligned N3T-ethyl-N3T cross-linked d(CGAAAT*TTTCG)2 undecamer, illustrating its universal use as computational tool to assist in structure elucidation. For the proton chemical shifts in the cross-linked system our methodology yields a strikingly superior description, surpassing the predictive power of (semi-)empirical methods. In addition, our methodology is the only one available to make an accurate prediction for the protons in the actual cross-link. To the best of our knowledge, this is the first computational study that attempts to determine the chemical shifts of oligonucleotides of this size and at this level of complexity.

An ecient protocol is presented to compensate for the basis set superposition error (BSSE) in DFT molecular dynamics (MD) simulations using localized Gaussian basis sets. We propose a classical correction term that can be added a posteriori to account for BSSE. It is tested to what extension this term will improve radial distribution functions (RDFs). The proposed term is pairwise between certain atoms in dierent molecules and was calibrated by tting reference BSSE data points computed with the counterpoise method. It is veried that the proposed exponential decaying functional form of the model is valid. This work focuses on hydrogen-bonded liquids, i.e. methanol, and more specic on the intermolecular hydrogen bond, but in principle the method is generally applicable on any type of interaction where BSSE is significant. We evaluated the relative importance of the Grimme-dispersion versus BSSE and found that they are of the same order of magnitude, but with an opposite sign. Upon introduction of the correction, the relevant RDFs, obtained from MD, have amplitudes equal to experiment.