No description available.

No description available.

No description available.

This study focused on the preparation and characterization of new indigo diimine (Nindigo) derivatives as a new atropisomeric scaffold. Trans- and cis- indigo diimines were studied and structure-property relationships were investigated regarding N-C atropisomerism using variable temperature 1H NMR studies and density functional theory calculations. Neutral trans- and protonated cis-Nindigos were prepared featuring a variety of mono ortho-substituted aryl imine groups with varying levels of steric bulk. The neutral trans-Nindigo derivatives generally have smaller N-C rotational energy barriers than their protonated cis-congeners. This finding is consistent with the latter having closer proximity of the N-aryl groups to each other, leading to steric repulsions between the two groups. The N-C rotational energy barriers are substituent dependent; the N-C rotational energy barriers of mono ortho substituted trans-Nindigos were in the range of 6.0 – 16.4 kcal/mol and can be classified as predominantly “Class 1’ atropisomers as defined by LaPlante, while the mono ortho substituted protonated cis-Nindigo analogs have N-C rotational barriers between 12.3 – 25.5 kcal/mol and are classified as “Class 1” and “Class 2” atropisomers. The introduction of additional substituents onto the other ortho position of the aryl imine subunit has significant consequences for the N-C rotational energy barriers of both the neutral trans- and protonated cis-Nindigos making them stable, or close to being, ‘Class 3’ atropisomers, having N-C rotational energy barriers between 31.5 – 276.9 kcal/mol and 29.3 – 32.6 kcal/mol respectively. Recognizing that the protonation state induced trans- to cis-isomerization process could have significant consequences regarding the potential applicability of these atropisomeric Nindigo derivatives, cis-Nindigo derivatives were synthesized that contained a tether (oxalyl or palladium (II) acetylacetonate) between the two indole type nitrogens of the Nindigo, which prevent the central -C=C- from isomerizing. The N-C rotational barriers of the tethered cis-Nindigos also displayed substituent dependent N-C rotational energy barriers. The protonation state of the N, N’-oxalyl bridged cis-Nindigos has a significant impact (higher in energy by a minimum of 5.1 kcal/mol) on the N-C rotational barriers; the neutral N, N’-oxalyl bridged cis-Nindigos have N-C rotational energy barriers ranging between 11.8 – 14.9 kcal/mol, classifying them as “Class 1” atropisomers, while their protonated congeners have N-C rotational energy barriers between 16.9 – 19.8 kcal/mol, which classifies them as “Class 1” atropisomers but are on the cusp of being “Class 2” atropisomers. The size of the tether influences the N-C rotational energy barriers of cis-Nindigos; the one-atom bridged palladium (II) acetylacetonate complexes have generally lower N-C rotational energy barriers than their protonated N, N’-oxalyl bridged cis-Nindigo congeners. The palladium acetylacetonate tethered cis-Nindigo complexes displayed substituent N-C rotational energy barrier dependence and the mono ortho substituted analogs have N-C rotational energy barriers between 12.4 – 20.2 kcal/mol and are predominantly “Class 1” atropisomers, while the bulkier analogs are “Class 2” atropisomers. The palladium (II) acetylacetonate cis-Nindigo complexes that have aryl imine groups with a 2,6-disubstitution pattern have N-C rotational energy barriers greater than 19.7 and 20.2 kcal/mol and are presumed to be stable “Class 3” atropisomers like their unbridged neutral trans- and protonated cis-Nindigo counterparts. / Graduate / 2023-09-12

he interactions of the chemical warfare agent sulfur mustard with amorphous silica were investigated using electronic structure calculations. In this thesis, the binding energies of sulfur mustard and mimic species used in the laboratory were calculated using density functional theory and fully ab initio calculations. The wB97XD and B97D functionals, which include functions to account for long-range dispersion interactions, were compared to experimental trends. The hydroxylated amorphous silica surface was approximated using a gas-phase silanol molecule and clusters containing a single hydroxyl moiety. Recent temperature programmed desorption experiments performed in UHV concluded that sulfur mustard and its less toxic mimics undergo molecular adsorption to amorphous silica. Hydrogen bonding can occur between surface silanol groups and either the sulfur or chlorine atom of the adsorbates, and the calculations indicate that the binding energies for the two hydrogen bond acceptors are similar. The adsorption of sulfur mustard and its mimics on silica also exhibits the presence of significant van der Waals interactions between alkyl of the adsorbates and the surface. These interactions, in combination with the formation of a hydrogen bond between a surface silanol group and the Cl or S atoms of the adsorbates, provide remarkably large binding energies. / Master of Science

Yes / The creation of multi-component crystals between chiral and achiral components has gained increased interest in recent years. In many cases the overall crystal structure is similar with the creation of a pseudo-inversion centre in the enantiopure case. This allows for the formation of solid solutions between the two extremes, which may have applications within chiral resolution. Utilising a combination of database mining, computational prediction and experimental screening, the frequency of formation for such materials has been investigated showing that for co-crystals this occurs more frequently than for salts, though there is a limited number of samples to draw structural conclusions. Computational modelling indicates the prediction of such systems can be challenging due to the similarities in energy of many crystal structures, so development of tools to design such systems is required to fully utilise these concepts. / The Deanship of Scientific Research at Princess Nourah bint Abdulrahman University through the Fast-track Research Funding Program.

Simulations of various hydration levels of lamellar phase 23:2 Diyne PC were performed, and subsequent, serial docking simulations of a tyrosine monomer were replicated for each system in both hydrated and dehydrated states.<br>The goal was to evaluate how hydration impacts self-assembly and crystallization on the surface, and<br>whether or not these simulations, when run sequentially, could determine the answer. It was discovered that hydrated and dehydrated surfaces behave differently, and that<br>headgroup orientation plays a role in the initial docking and self-assembly process of the tyrosine monomer. It was also determined that potential energy as a sole metric<br>for determining whether or not a specific conformation of intermolecular orientation is not entirely useful, and docking scores are likely useful metrics in discriminating between conformations with identical potential energy values. <br>

<div>Computational Chemistry has proven to be an effective means of solving chemical problems. The two main tools of Computational Chemistry - quantum mechanics and molecular mechanics, have provided viable avenues to probe such chemical problems at an electronic or molecular level, with varying levels of accuracy and speed. In this work, attempts have been made to combine the speed of molecular mechanics and the accuracy of quantum mechanics to work across multiples scales of time and length, effectively resulting in simulations of large chemical systems without compromising the accuracy.</div><div><br></div><div>The primary tool utilized for methods development and application in this work is the Effective Fragment Potential (EFP) method. The EFP method is a computational technique for studying non-covalent interactions in complex systems. EFP is an accurate \textit{ab initio} force field, with accuracy comparable to many Density Functional Theory (DFT) methods, at significantly lower computational cost. EFP decomposes intermolecular interactions into contributions from four terms: electrostatics, polarization, exchange-repulsion and dispersion.</div><div><br></div><div>In the first chapter, the possibility of applying EFP method to study large radical-water clusters is probed. An approximate theoretical model in which the transition dipole moments of excitations are computed using the information from the ground state orbitals is implemented.</div><div><br></div><div>A major challenge to broaden the scope of EFP is to overcome its limitation in describing only small and rigid molecules such as water, acetone, etc. In the second chapter, the extension of EFP method to large covalently bound biomolecules and polymers such as proteins, lipids etc., is described. Using this new method, referred to as BioEFP/mEFP, it is shown that the effect of polarization is non-negligible and must be accounted for when modeling photochemical and electron-transfer processes in photoactive proteins.</div><div><br></div><div>Another area of interest is the development of novel drug-target binding models, in which a chemically active part of the ligand is modified via functional group modification, while the rest of the system remains intact. In the third chapter, the development and application of a drug-target binding model is explained.<br></div><div><br></div><div><div>Lastly, in the fourth and final chapter, we show the derivation for working equations corresponding to the coupling gradient term describing the dispersion interactions between quantum mechanical and effective fragment potential regions.</div><div><br></div><div>The primary focus of this work is to explore and expand the boundaries of multiscale QM/MM simulations applied to chemical and biomolecular systems. We believe that the work described here leads to exciting pathways in the future in terms of modeling novel systems and processes such as heterogeneous catalysis, QSAR, crystal structure prediction, etc.</div></div>

Hydrogen storage in lightweight hydrides continues to attract significant interest as the lack of a safe and efficient storage of hydrogen remains the major technological barrier to the widespread use of hydrogen as a fuel. The metal borohydrides Ca(BH₄)₂ and LiBH₄ form the subject of this thesis; three aspects of considerable academic interest were investigated by density functional theory (DFT) and molecular dynamics (MD) modelling. (i) High-pressure crystal structures of Ca(BH₄)₂ were predicted from a structural analogy between metal borohydrides and isoelectronic metal oxides. The structural stability of hydrogen storage materials under high pressure is an important aspect, as high-pressure polymorphs may provide structures with better hydrogen desorption properties. The isoelectronic analogue of Ca(BH₄)₂ is TiO₂, and structural equivalents of Ca(BH₄)₂ in the baddeleyite, columbite and cotunnite structures of TiO₂ were found to be stable at elevated pressure. Thermodynamic stability was evaluated by computing the Gibbs energy with respect to pressure and temperature. The pressure-dependence of the Helmholtz energy was determined to described a third-order Birch-Murnaghan equation of state, and the harmonic approximation was used to compute the vibrational energy levels and the Helmholtz energy as a function of temperature. The proposed structures are consistent with reports of two hitherto unidentified high-pressure phases observed experimentally. (ii) The disordered structure of the high-temperature phase of LiBH4 was studied by ab initio molecular dynamics (MD) at temperatures ranging from 200-535 K. It was found that the model emerging from analysis of the MD simulations properly accounts for dynamical disorder and fundamentally differs from the published experimental and theoretical structures. The validity of the MD model was corroborated by comparison of calculated pair distribution functions, vibrational spectra and a crystallographic model with neutron diffraction data; good agreement was found. A reassignment of the space group from P63mc to P63/mmc is proposed based on evidence for additional symmetry from MD simulations. (iii) Finally, a new MD-based method was developed to simulate fast ionic diffusion in LiBH₄. The colour diffusion algorithm - a nonequilibrium molecular dynamics method originally developed for the study of model fluids - was adapted and applied to self-diffusion of atoms in a solid for the first time. Calculated diffusion coefficients agreed very well with published measurements, and diffusion pathways that include collective particle effects were determined directly from the simulation results, thereby opening up a promising and efficient new method for the study of phenomena such as superionic conduction.

In this thesis, mathematical modelling is used to theoretically investigate the electrochemical behaviour of surfaces which can be broadly classified as being ‘electrochemically heterogeneous’. Simulated voltammetry is used in the exploration of a number of specific systems as listed below. The cyclic voltammetry of electrodes composed of two different electroactive materials that differ in terms of their electrochemical rate constants towards any given redox couple. The effect of the distribution of the two materials was investigated and the occurrence of split peak cyclic voltammetry where two peaks are observed in the forward sweep, was studied. The technique was specifically applied to the modelling of highly-ordered pyrolytic graphite (HOPG). The steady-state voltammetry of a conducting spherical particle resting on an insulating supporting surface. An algebraic expression that completely describes the voltammetric waveform in the limit of irreversible kinetics was developed. The cyclic voltammetry of the EC′ (catalytic) mechanism at a regularly distributed array of hemispherical particles on an insulating supporting surface. Particular attention was paid to the ‘split-wave’ phenomenon, where two peaks are observed in the forward scan of a cyclic voltammogram and the conditions under which these peaks are resolvable were elucidated. The linear sweep voltammetry of micro- and nano-particle modified electrodes and other electrodes of partially covered and non-planar geometry. It was demonstrated that the apparent electrochemical rate constant of the reaction and thus the peak position of the voltammetry is dependent only on the relative electroactive surface area of the particles on the surface and not upon their shape or distribution. This has wide reaching implications as it can be used to explain some instances of a purported nano-catalytic effect without appeal to altered properties at the nanoscale. The linear sweep voltammetry of the interior of a partially electroactive cylindrical pore. Four limiting cases were observed and fully characterised. The linear sweep voltammetry of porous surfaces. It was established that if the pores are less than a certain threshold depth, then a porous surface will also display an apparent catalytic effect that is dependent on the relative electroactive surface area (including the area in the interior of the pores).