A new SOS-DFPT approximation for NMR shielding calculations : the Loc.3 correction applied to the catalytic mechanism of Serine Proteases

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Path Optimization in Free Energy Calculations

Free energy calculations are a computational method for determining thermodynamic quantities, such as free energies of binding, via simulation.

Currently, due to computational and algorithmic limitations, free energy calculations are limited in scope.

In this work, we propose two methods for improving the efficiency of free energy calculations.

First, we expand the state space of alchemical intermediates, and show that this expansion enables us to calculate free energies along lower variance paths.

We use Q-learning, a reinforcement learning technique, to discover and optimize paths at low computational cost.

Second, we reduce the cost of sampling along a given path by using sequential Monte Carlo samplers.

We develop a new free energy estimator, pCrooks (pairwise Crooks), a variant on the Crooks fluctuation theorem (CFT), which enables decomposition of the variance of the free energy estimate for discrete paths, while retaining beneficial characteristics of CFT.

Combining these two advancements, we show that for some test models, optimal expanded-space paths have a nearly 80% reduction in variance relative to the standard path.

Additionally, our free energy estimator converges at a more consistent rate and on average 1.8 times faster when we enable path searching, even when the cost of path discovery and refinement is considered.

A new SOS-DFPT approximation for NMR shielding calculations : the Loc.3 correction applied to the catalytic mechanism of Serine Proteases

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Alkenone stable carbon isotopes, carbonate stable carbon and oxygen isotopes, planktic foraminifera boron isotopes and atmospheric CO2 calculations and box model results from Ras il-Pellegrin section, Malta

The development of a permanent, stable ice sheet in East Antarctica happened during the middle Miocene, about 14 million years (Myr) ago. The middle Miocene therefore represents one of the distinct phases of rapid change in the transition from the "greenhouse" of the early Eocene to the "icehouse" of the present day. Carbonate carbon isotope records of the period immediately following the main stage of ice sheet development reveal a major perturbation in the carbon system, represented by the positive d13C excursion known as carbon maximum 6 ("M6"), which has traditionally been interpreted as reflecting increased burial of organic matter and atmospheric pCO2 drawdown. More recently, it has been suggested that the d13C excursion records a negative feedback resulting from the reduction of silicate weathering and an increase in atmospheric pCO2. Here we present high-resolution multi-proxy (alkenone carbon and foraminiferal boron isotope) records of atmospheric carbon dioxide and sea surface temperature across CM6. Similar to previously published records spanning this interval, our records document a world of generally low (~300 ppm) atmospheric pCO2 at a time generally accepted to be much warmer than today. Crucially, they also reveal a pCO2 decrease with associated cooling, which demonstrates that the carbon burial hypothesis for CM6 is feasible and could have acted as a positive feedback on global cooling.

Cost-effectiveness of natriuretic peptide-based screening and collaborative care: a report from the STOP-HF (St Vincent's Screening TO Prevent Heart Failure) study

AIMS: Prevention of cardiovascular disease and heart failure (HF) in a cost-effective manner is a public health goal. This work aims to assess the cost-effectiveness of the St Vincent's Screening TO Prevent Heart Failure (STOP-HF) intervention.

METHODS AND RESULTS: This is a substudy of 1054 participants with cardiovascular risk factors [median age 65.8 years, interquartile range (IQR) 57.8:72.4, with 4.3 years, IQR 3.4:5.2, follow-up]. Annual natriuretic peptide-based screening was performed, with collaborative cardiovascular care between specialist physicians and general practitioners provided to patients with BNP levels >50 pg/mL. Analysis of cost per case prevented and cost-effectiveness per quality-adjusted life year (QALY) gained was performed. The primary clinical endpoint of LV dysfunction (LVD) with or without HF was reduced in intervention patients [odds ratio (OR) 0.60; 95% confidence interval (CI) 0.38-0.94; P = 0.026]. There were 157 deaths and/or emergency hospitalizations for major adverse cardiac events (MACE) in the control group vs. 102 in the intervention group (OR 0.68; 95% CI 0.49-0.93; P = 0.01). The cost per case of LVD/HF prevented was €9683 (sensitivity range -€843 to €20 210), whereas the cost per MACE prevented was €3471 (sensitivity range -€302 to €7245). Cardiovascular hospitalization savings offset increased outpatient and primary care costs. The cost per QALY gain was €1104 and the intervention has an 88% probability of being cost-effective at a willingness to pay threshold of €30 000.

CONCLUSION: Among patients with cardiovascular risk factors, natriuretic peptide-based screening and collaborative care reduced LVD, HF, and MACE, and has a high probability of being cost-effective.

TRIAL REGISTRATION: NCT00921960.

Natriuretic peptide-based screening and collaborative care for heart failure:the STOP-HF randomized trial

IMPORTANCE: Prevention strategies for heart failure are needed.

OBJECTIVE: To determine the efficacy of a screening program using brain-type natriuretic peptide (BNP) and collaborative care in an at-risk population in reducing newly diagnosed heart failure and prevalence of significant left ventricular (LV) systolic and/or diastolic dysfunction.

DESIGN, SETTING, AND PARTICIPANTS: The St Vincent's Screening to Prevent Heart Failure Study, a parallel-group randomized trial involving 1374 participants with cardiovascular risk factors (mean age, 64.8 [SD, 10.2] years) recruited from 39 primary care practices in Ireland between January 2005 and December 2009 and followed up until December 2011 (mean follow-up, 4.2 [SD, 1.2] years).

INTERVENTION: Patients were randomly assigned to receive usual primary care (control condition; n=677) or screening with BNP testing (n=697). Intervention-group participants with BNP levels of 50 pg/mL or higher underwent echocardiography and collaborative care between their primary care physician and specialist cardiovascular service.

MAIN OUTCOMES AND MEASURES: The primary end point was prevalence of asymptomatic LV dysfunction with or without newly diagnosed heart failure. Secondary end points included emergency hospitalization for arrhythmia, transient ischemic attack, stroke, myocardial infarction, peripheral or pulmonary thrombosis/embolus, or heart failure.

RESULTS: A total of 263 patients (41.6%) in the intervention group had at least 1 BNP reading of 50 pg/mL or higher. The intervention group underwent more cardiovascular investigations (control, 496 per 1000 patient-years vs intervention, 850 per 1000 patient-years; incidence rate ratio, 1.71; 95% CI, 1.61-1.83; P<.001) and received more renin-angiotensin-aldosterone system-based therapy at follow-up (control, 49.6%; intervention, 56.5%; P=.01). The primary end point of LV dysfunction with or without heart failure was met in 59 (8.7%) of 677 in the control group and 37 (5.3%) of 697 in the intervention group (odds ratio [OR], 0.55; 95% CI, 0.37-0.82; P = .003). Asymptomatic LV dysfunction was found in 45 (6.6%) of 677 control-group patients and 30 (4.3%) of 697 intervention-group patients (OR, 0.57; 95% CI, 0.37-0.88; P = .01). Heart failure occurred in 14 (2.1%) of 677 control-group patients and 7 (1.0%) of 697 intervention-group patients (OR, 0.48; 95% CI, 0.20-1.20; P = .12). The incidence rates of emergency hospitalization for major cardiovascular events were 40.4 per 1000 patient-years in the control group vs 22.3 per 1000 patient-years in the intervention group (incidence rate ratio, 0.60; 95% CI, 0.45-0.81; P = .002).

CONCLUSION AND RELEVANCE: Among patients at risk of heart failure, BNP-based screening and collaborative care reduced the combined rates of LV systolic dysfunction, diastolic dysfunction, and heart failure.

TRIAL REGISTRATION: clinicaltrials.gov Identifier: NCT00921960.

Progression of left atrial volume index in a population at risk for heart failure:a substudy of the STOP-HF (St Vincent's Screening TO Prevent Heart Failure) trial

AIMS: Limited data are available concerning the evolution of the left atrial volume index (LAVI) in pre-heart failure (HF) patients. The aim of this study was to investigate clinical characteristics and serological biomarkers in a cohort with risk factors for HF and evidence of serial atrial dilatation.

METHODS AND RESULTS: This was a prospective substudy within the framework of the STOP-HF cohort (NCT00921960) involving 518 patients with risk factors for HF electively undergoing serial clinical, echocardiographic, and natriuretic peptide assessment. Mean follow-up time between assessments was 15 ± 6 months. 'Progressors' (n = 39) were defined as those with serial LAVI change ≥3.5 mL/m(2) (and baseline LAVI between 20 and 34 mL/m(2)). This cut-off was derived from a calculated reference change value above the biological, analytical, and observer variability of serial LAVI measurement. Multivariate analysis identified significant baseline clinical associates of LAVI progression as increased age, beta-blocker usage, and left ventricular mass index (all P < 0.05). Serological biomarkers were measured in a randomly selected subcohort of 30 'Progressors' matched to 30 'Non-progressors'. For 'Progressors', relative changes in matrix metalloproteinase 9 (MMP9), tissue inhibitor of metalloproteinase 1 (TIMP1), and the TIMP1/MMP9 ratio, markers of interstitial remodelling, tracked with changes in LAVI over time (all P < 0.05).

CONCLUSION: Accelerated LAVI increase was found to occur in up to 14% of all pre-HF patients undergoing serial echocardiograms over a relatively short follow-up period. In a randomly selected subcohort of 'Progressors', changes in LAVI were closely linked with alterations in MMP9, TIMP1, and the ratio of these enzymes, a potential aid in highlighting this at-risk group.

Benchmarking a modified version of the civ3 nonrelativistic atomic-structure code within Na-like-tungsten R -matrix calculations

In this work we explore the validity of employing a modified version of the nonrelativistic structure code civ3 for heavy, highly charged systems, using Na-like tungsten as a simple benchmark. Consequently, we present radiative and subsequent collisional atomic data compared with corresponding results from a fully relativistic structure and collisional model. Our motivation for this line of study is to benchmark civ3 against the relativistic grasp0 structure code. This is an important study as civ3 wave functions in nonrelativistic R-matrix calculations are computationally less expensive than their Dirac counterparts. There are very few existing data for the W LXIV ion in the literature with which we can compare except for an incomplete set of energy levels available from the NIST database. The overall accuracy of the present results is thus determined by the comparison between the civ3 and grasp0 structure codes alongside collisional atomic data computed by the R-matrix Breit-Pauli and Dirac codes. It is found that the electron-impact collision strengths and effective collision strengths computed by these differing methods are in good general agreement for the majority of the transitions considered, across a broad range of electron temperatures.

Practical calculations for uniform external pressure buckling in cylindrical shells with stepped walls

Metal cylindrical storage structures of significant size, such as silos and vertical-axis tanks, are almost always constructed from many short cylindrical shells of different thickness as the stress resultants on the wall progressively increase towards the base. The resulting increases in thickness are always made in step changes using metal sheets of uniform thickness because of the availability of such source materials. The result is a shell with a stepped wall with multiple discrete steps in thickness. Such shells are very susceptible to buckling under external pressure when empty or partially filled, but the buckling mode may involve only part of the shell height due to the changes in shell thickness. These changes must therefore be accounted for within the design process. A new method of determining the critical buckling resistance of such shells was recently developed, and although it has been shown to be valid, the methodology for its application in practical design has not been set out or shown. This paper therefore briefly describes the new method and demonstrates the manner in which it can be used to produce rapid, safe assessments of cylindrical shells with a wide range of patterns of wall thickness changes. The results are then suitable for direct introduction into such documents as the European standard on metal shells [1] and the ECCS Recommendations [2].

3D Ionospheric Effects on HF Propagation and Heating

The thesis uses a three-dimensional, first-principles model of the ionosphere in combination with High Frequency (HF) raytracing model to address key topics related to the physics of HF propagation and artificial ionospheric heating. In particular: 1. Explores the effect of the ubiquitous electron density gradients caused by Medium Scale Traveling Ionospheric Disturbances (MSTIDs) on high-angle of incidence HF radio wave propagation. Previous studies neglected the all-important presence of horizontal gradients in both the cross- and down-range directions, which refract the HF waves, significantly changing their path through the ionosphere. The physics-based ionosphere model SAMI3/ESF is used to generate a self-consistently evolving MSTID that allows for the examination of the spatio-temporal progression of the HF radio waves in the ionosphere. 2. Tests the potential and determines engineering requirements for ground- based high power HF heaters to trigger and control the evolution of Equatorial Spread F (ESF). Interference from ESF on radio wave propagation through the ionosphere remains a critical issue on HF systems reliability. Artificial HF heating has been shown to create plasma density cavities in the ionosphere similar to those that may trigger ESF bubbles. The work explores whether HF heating may trigger or control ESF bubbles. 3. Uses the combined ionosphere and HF raytracing models to create the first self-consistent HF Heating model. This model is utilized to simulate results from an Arecibo experiment and to provide understanding of the physical mechanism behind observed phenomena. The insights gained provide engineering guidance for new artificial heaters that are being built for use in low to middle latitude regions. In accomplishing the above topics: (i) I generated a model MSTID using the SAMI3/ESF code, and used a raytrace model to examine the effects of the MSTID gradients on radio wave propagation observables; (ii) I implemented a three- dimensional HF heating model in SAMI3/ESF and used the model to determine whether HF heating could artificially generate an ESF bubble; (iii) I created the first self-consistent model for artificial HF heating using the SAMI3/ESF ionosphere model and the MoJo raytrace model and ran a series of simulations that successfully modeled the results of early artificial heating experiments at Arecibo.

Approximate ab initio calculations and the method of molecular fragments

A two stage approach to performing ab initio calculations on medium and large sized molecules is described. The first step is to perform SCF calculations on small molecules or molecular fragments using the OPIT Program. This employs a small basis set of spherical and p-type Gaussian functions. The Gaussian functions can be identified very closely with atomic cores, bond pairs, lone pairs, etc. The position and exponent of any of the Gaussian functions can be varied by OPIT to produce a small but fully optimised basis set. The second stage is the molecular fragments method. As an example of this, Gaussian exponents and distances are taken from an OPIT calculation on ethylene and used unchanged in a single SCF calculation on benzene. Approximate ab initio calculations of this type give much useful information and are often preferable to semi-empirical approaches, since the nature of the approximations involved is much better defined.