1. |
- Nørby, Morten Steen
(author)
-
Extending the capabilities of polarizable embedding through damped linear response theory
- 2016
-
Doctoral thesis (other academic/artistic)abstract
- One of the greater challenges of modern quantum chemistry is to meet the growing demands of a detailed description of macro molecules, such as proteins and RNA/DNA. The size of such systems make it an almost impossible task to, from first principle calculations, obtain any detailed information regarding electronic structure. This thesis is focused on development of new computational tools for investigation of electronic transitions within large molecular systems. Central to this is the polarizable embedding (PE) model. For many of these systems electronic transitions are localized to one or a few active sites which in turn will be described by accurate quantum chemical methods. The remaining part of the system, defining the environment, is then treated at a lower level using a classical electrostatic approach. The environment is further incorporated into the Hamiltonian describing the quantum region in an effective manner.The work presented in this thesis covers extensions to the PE model both in the quantum region and in the environmental description. The former, starts by combining the PE model with a damped linear response formalism. This coupling enables the study of electronic processes near resonance and in any frequency range desirable thus covering all from core hole transitions to valence transitions. Another important aspect of this work is the coupling of PE to a polarizable continuum model which allows for inclusion of bulk solvent effects in an effective manner. In addition to these extensions we have further considered alternative ways of acquiring parameters used in the construction of the embedding potential. This is essential for accurate embedding calculations as the embedding potential enters directly in the wave function optimization.
|
|
2. |
- Ahmadzadeh, Karan, 1991-
(author)
-
Efficient Calculation of Nonlinear Spectroscopic Properties within the Time-Dependent Density Functional Theory Approximation
- 2023
-
Doctoral thesis (other academic/artistic)abstract
- This thesis introduces a novel computational scheme tailored for efficient calculations of nonlinear spectroscopic observables. First, a derivation and implementation of an algorithm designed to harness the linearity of the Fock matrix construction in calculating two-photon absorption cross-sections within the self-consistent field approximation is presented. Subsequently, this computational scheme is extended to the density functional theory approximation for functionals belonging to the generalized gradient approximation. Lastly the derivation and implementation of the nonlinear exchange-correlation kernel for functionals belonging to the meta generalized gradient approximation are presented for the first time.Collectively, the advancements presented in this thesis contribute new methodologies and insights to the computational realm of nonlinear spectroscopic calculations, offering the potential for large-scale theoretical spectroscopy calculations at the level of density functional theory.
|
|
3. |
- Brand, Manuel
(author)
-
Implementation of complex polarization propagator theory for linear response properties of large molecular systems
- 2022
-
Doctoral thesis (other academic/artistic)abstract
- Since its beginning, the remarkable development from the first commercially available computers toward exascale supercomputers just within the span of a lifetime has been closely intertwined with the perpetual quest for the utilization of the arising computing power for the avail of theoretical chemistry. With the aim of further pushing the limits of computationally accessible molecular system sizes, this thesis includes the presentation of programming efforts, which brought forth two quantum chemical software codes, as well as a range of ab initio studies on carbon-based systems, enabled by the former.The VeloxChem and Gator programs, developed for spectroscopy simulations at the level of density-functional theory (DFT) and correlated wave function methods, respectively, employ a hybrid message passing interface (MPI)/open multiprocessing (OpenMP) parallelization scheme embedded in a modular program structure written in a Python/C++ layered fashion for the execution in contemporary high-performance computing (HPC) environments. Included numerical solvers for the evaluation of real and complex linear response functions in combination with the parallel construction of auxiliary Fock matrices enable the efficient calculation of one-photon absorption or electronic circular dichroism (CD) spectra in the ultraviolet/visible (UV/vis) or X-ray spectral region, as well as van der Waals C6 dispersion coefficients.Employing the VeloxChem program in two comprehensive investigations, the C6 dispersion coefficients of carbon fullerenes up to a system size of C540 and the CD spectra of carbohelicenes ranging from CH[5] to CH[30] have been calculated at the DFT level of theory. The revealed non-additive scaling with respect to the number of carbon atoms of N2.2 in the former and nontrivial CD band progressions, arguably linked to the number of overlapping layers in the helical conjugated systems, in the latter rectify the current conception in their respective fields.In a benchmark for the Gator program on a series of guanine oligomers, the full UV/vis spectrum for a tetrad was calculated at the level of a second-order algebraic-diagrammatic construction (ADC(2)) scheme in just under 15 hours by efficient employment of 32,768 central processing unit (CPU) cores.Exceeding the limit of 10,000 and 1,000 contracted basis functions for a treatment with the DFT and ADC(2) methods, respectively, these practical examples demonstrate the capability of VeloxChem and Gator to harness vast computational resources made available by contemporary and future HPC systems and thereby routinely address scopes of system sizes that were previously out of reach.
|
|
4. |
- Falklöf, Olle
(author)
-
Computational Studies of Photobiological Keto-Enol Reactions and Chromophores
- 2015
-
Doctoral thesis (other academic/artistic)abstract
- This thesis presents computational chemistry studies of keto-enol reactions and chromophores of photobiological signicance.The rst part of the thesis is concerned with two protein-bound chromophores that, depending on the chemical conditions, can exist in a number of dierent ketonic and enolic forms. The rst chromophore is astaxanthin, which occurs in the protein complex responsible for the deep-blue color of lobster carapace. By investigating how dierent forms of astaxanthin absorb UV-vis radiation of dierent wavelengths, a model is presented that explains the origin of the dramatic color change from deep-blue to red upon cooking of live lobsters.The second chromophore is the oxyluciferin light emitter of fireflies, which is formed in the catalytic center of the enzyme firefly luciferase. To date, there is no consensus regarding which of the possible ketonic and enolic forms is the key contributor to the light emission. In the thesis, the intrinsic tendency of oxyluciferin to prefer one particular form over other possible forms is established through calculation of keto-enol and acid-base excited-state equilibrium constants in aqueous solution.The second part of the thesis is concerned with two families of biological photoreceptors: the blue-light-absorbing LOV-domain proteins and the red-light-absorbing phytochromes. Based on the ambient light environment, these proteins regulate physiological and developmental processes by switching between inactive and active conformations. In both families, the conversion of the inactive into the active conformation is triggered by a chemical reaction of the respective chromophore.The LOV-domain proteins bind a LOV-domain proteins bidn in flavin chromophore and regulate processes such as chloroplast relocation and phototropism in plants. An important step in the activation of these photoreceptors is a singlet-triplet transition between two electronically excited states of the flavin chromophore. In the thesis, this transition is used as a prototype example for illustrating, for the rst time, the ability of rst-principles methods to calculate rate constants of inter-excited state phosphorescence events.Phytochromes, in turn, bind bilin chromophores and are active in the regulation of processes like seed germination and flowering time in plants. Following two systematic studies identifying the best way to model the UV-vis absorption and fluorescence spectra of these photoreceptors, it is demonstrated that steric interactions between the chromophore and the apoprotein play a decisive role for how phytochromes are activated by light.
|
|
5. |
- Fransson, Thomas, 1986-
(author)
-
X-ray spectroscopies through damped linear response theory
- 2016
-
Doctoral thesis (other academic/artistic)abstract
- In order to reach a fundamental understanding of interactions between electromagnetic radiation and molecular materials, experimental measurements need to be supplemented with theoretical models and simulations. With the use of this combination, it is possible to characterize materials in terms of, e.g., chemical composition and molecular structure, as well as achieve time-resolution in studies of chemical reactions. This doctoral thesis focuses on the development and evaluation of theoretical methods with which, amongst others, X-ray absorption and X-ray emission spectroscopies can be interpreted and predicted. In X-ray absorption spectroscopy the photon energy is tuned such that core electrons are targeted and excited to either bound or continuum states, and X-ray emission spectroscopy measures the subsequent decay from such an excited state. These core excitations/de-excitations exhibit strong relaxation effects, making theoretical considerations of the processes particularly challenging. While the removal of a valence electron leaves the remaining electrons relatively unaffected, removing core electrons has a substantial effect on the other electrons due to the significant change in the screening of the nucleus. Additionally, the core-excited states are embedded in a manifold of valence-excited states that needs to be considered by some computationally feasible method. In this thesis, a damped formalism of linear response theory, which is a perturbative manner of considering the interactions of (weak) external or internal fields with molecular systems, has been utilized to investigate mainly the X-ray absorption spectra of small- to medium-sized molecular systems.Amongst the standard quantum chemical methods available, coupled cluster is perhaps the most accurate, with a well-defined, hierarchical manner of approaching the correct electronic wave function. Combined with response theory, it provides a reliable theoretical method in which relaxation effects are addressed by means of an accurate treatment of electron correlation. The first part of this thesis deals with the development and evaluation of such an approach, and it is shown that the relaxation effects can be addressed by the inclusion of double excitations in the coupled cluster manifold. However, these calculations are computationally very demanding, and in order to treat larger systems the performance of the coupled cluster approach has been compared to that of the less demanding method of time-density dependent functional theory (TDDFT). Both methods have been used to investigate the X-ray absorption spectrum of water, which has been extensively debated in the scientific community following a relatively recent hypothesis concerning the underlying structure of liquid water. Water exhibits a great number of anomalous properties that stand out from those of most compounds, and the importance of reaching a fundamental understanding of this substance cannot be overstated. It has been demonstrated that TDDFT yields excellent results for liquid water, opening up possibilities of investigating the correlation between spectral features and local structures.Furthermore, recent developments in damped linear response TDDFT in the four-component relativistic regime have enabled the inclusion of spin-orbit coupling in damped linear response calculations, making black-box calculations of absorption spectra in a relativistic setting practical. With this approach, it is possible to address the spin-orbit splitting in L2,3-edge X-ray absorption spectra, and the performance of such a method has been demonstrated for a set of small molecules. Excellent agreement with experiment is obtained in terms of relative features, but an anomalous error in absolute energy has been observed for silane derivatives featuring fluorine-substitutions. This is likely a result of the strong influence of the very electronegative fluorine atoms on the electron density of the core-excited atom.Finally, the treatment of non-resonant X-ray emission spectroscopy using damped linear response theory is discussed. The expansion needed for the development of a simple method by which this spectroscopy can be treated using damped linear response theory at the TDDFT level of theory has been identified, and proof of principle calculations at the time-dependent Hartree-Fock level of theory are presented.
|
|
6. |
- Gustafsson, Camilla
(author)
-
Modeling environment effects on spectroscopic properties of biomarkers and catalytic mechanisms in enzymes
- 2020
-
Doctoral thesis (other academic/artistic)abstract
- Arguably, humans are in need of both better diagnostic tools to prevent pro- gression of diseases as well as greener catalysts for synthesis of chemicals.Neurodegenerative diseases affecting neurons in the brain leads to demen- tias, where Alzheimer’s disease (AD) is the most prevalent. It is estimated that about 50 million people worldwide suffer from AD, a number that has more than doubled during the last 30 years. Currently, there is no cure for AD, but in order to slow the progression of symptoms it is crucial to develop biomarkers for early detection and initiation of clinical interventions.With theoretical tools it is possible to better understand the optical prop- erties of fluorescent biomarkers, and thus contribute to steering the design of biomarkers for distinguishing different types of disease-associated proteins. Lu- minescent conjugated oligothiophenes (LCO) is a class of molecules that binds to aggregates of misfolded amyloid-β proteins, facilitating in vivo-detection of the pathological hallmarks of AD. By performing molecular dynamics (MD) simulations and subsequent response theory calculations of a LCO, it could be concluded that the differences in the spectroscopic fingerprints for the bound and free biomarker were predominantly due to conformational changes of the conjugated π-system in the molecular backbone. The introduction of differ- ent central units with donor properties yield donor-acceptor-donor electronic systems that increase the range of spectroscopic detection of LCO biomark- ers, without reducing the selectivity towards amyloid-β. It was also revealed that in order to capture more of the two-photon absorption (TPA) signal it would be optimal to design biomarkers with the dominant TPA signal at longer wavelenghts.The second part of this work is centered around computational enzyme design, and how single point mutations can alter the flow of water in the active site. The altered flow of water likely impacts the catalysis in the active site of the enzymes. The enzymes considered in this work belongs to two different enzyme classes, and catalyse different kinds of reactions. Squalene hopene cyclase (SHC) is a monotopic membrane enzyme that catalyses the cyclization of squalene to hopene, and ω-transaminase catalyses the transfer of an amino and keto group between an amino acid and a keto acid. Enzyme variants of both SHC and ω-transaminase, where single-point mutations have been introduced, display different experimentally observed properties compared to their corresponding wild-types (WT). By performing MD simulations, the flow of water in the active sites of both enzymes could be tracked. Distinct differences in the flow of water in the WT and enzyme variants could be detected. These changes are proposed to influence the catalysis, and help to explain the experimentally observed differences in the protein variants.
|
|
7. |
- Jakobsson, Mattias, 1981-
(author)
-
Monte Carlo Studies of Charge Transport Below the Mobility Edge
- 2012
-
Doctoral thesis (other academic/artistic)abstract
- Charge transport below the mobility edge, where the charge carriers are hopping between localized electronic states, is the dominant charge transport mechanism in a wide range of disordered materials. This type of incoherent charge transport is fundamentally different from the coherent charge transport in ordered crystalline materials. With the advent of organic electronics, where small organic molecules or polymers replace traditional inorganic semiconductors, the interest for this type of hopping charge transport has increased greatly. The work documented in this thesis has been dedicated to the understanding of this charge transport below the mobility edge.While analytical solutions exist for the transport coefficients in several simplified models of hopping charge transport, no analytical solutions yet exist that can describe these coefficients in most real systems. Due to this, Monte Carlo simulations, sometimes described as ideal experiments performed by computers, have been extensively used in this work.A particularly interesting organic system is deoxyribonucleic acid (DNA). Besides its overwhelming biological importance, DNA’s recognition and self-assembly properties have made it an interesting candidate as a molecular wire in the field of molecular electronics. In this work, it is shown that incoherent hopping and the Nobel prize-awarded Marcus theory can be used to describe the results of experimental studies on DNA. Furthermore, using this experimentally verified model, predictions of the bottlenecks in DNA conduction are made.The second part of this work concerns charge transport in conjugated polymers, the flagship of organic materials with respect to processability. It is shown that polaronic effects, accounted for by Marcus theory but not by the more commonly used Miller-Abrahams theory, can be very important for the charge transport process. A significant step is also taken in the modeling of the off-diagonal disorder in organic systems. By taking the geometry of the system from large-scale molecular dynamics simulations and calculating the electronic transfer integrals using Mulliken theory, the off-diagonal disorder is for the first time modeled directly from theory without the need for an assumed parametric random distribution.
|
|
8. |
- Sjöqvist, Jonas, 1985-
(author)
-
Luminescence properties of flexible conjugated dyes
- 2012
-
Licentiate thesis (other academic/artistic)abstract
- In this licentiate thesis the luminescence properties of two flexible conjugated dyes have been studied. The first, Pt1, is a platinum(II) acetylide chromophore used in optical power limiting materials. The second is a set of optical probes known as luminescent conjugated oligothiophenes (LCOs), which are used to detect and characterize the protein structures associated with amyloid diseases such as Alzheimer’s disease.MM3 and CHARMM force field parameters have been derived for the Pt1 chromophore and LCOs, respectively, based on potential energy surface references calculated at the density functional theory (DFT)/B3LYP level of theory. The parameters have been used to perform room temperature molecular dynamics simulations of the chromophores in solvent, where tetrahydrofuran was used for Pt1 and water for the LCOs. Conformationally averaged absorption spectra were obtained, based on response theory calculations at the time-dependent DFT(TDDFT)/CAM-B3LYP level of theory for a selection of structures from the simulations. For one of the LCOs, p-HTAA, force field parameters were also created describing the dominant first excited state, based on TDDFT/B3LYP reference potential energy surfaces. These were used for molecular dynamics simulations of the chromophore in the excited state, allowing the creation of an emission spectrum. A theoretically obtained Stokes shift of 112 nm could be computed based on the absorption and emission spectra, which is in good agreement with the experimental value of 124 nm.In addition, a quantum mechanics/molecular mechanics study of the effects of solvation on the absorption properties of the p-HTAA chromophore in water has been conducted, resulting in two models for including these effects in the averaged spectra. The first includes explicit water molecules in the form of point charges and polarizable dipole moments, and results in an absorption wavelength that is blueshifted by 2 nm from a high quality reference calculation. The second model involves the complete removal of the solvent as well as the ionic groups of the chromophore. The resulting absorption wavelength is blueshifted by an additional 4 nm as compared to the first model, but requires only one fifth of the computational resources.
|
|
9. |
- Todarwal, Yogesh
(author)
-
Accurate Force Fields for Spectroscopic Studies of Protein–Ligand Interactions and Self-Assembly Structures
- 2024
-
Doctoral thesis (other academic/artistic)abstract
- The computational prediction of complex molecular behaviors is an essen- tial component of modern chemistry, as it provides a faster and more cost- effective way to explore molecular interactions that may be difficult or even impossible to study experimentally. Molecular dynamics (MD) simulations of- ten serve as a valuable tool for such predictions; however, their accuracy is inherently dependent on the force field (FF) parameters employed. While the general amber force field (GAFF) is designed to provide reasonable results for a broad array of small molecules, it often requires further refinement when using it for a specific small organic molecule. Especially for ligands of the oligothiophene class, the dihedral potential representing the rotatable bond between the two thiophene rings (of the SCCS type) is inadequately described. An objective of this dissertation is to refine FF parameters for producing meaningful MD trajectories that capture key molecular interactions, binding modes, and thermodynamic properties, and subsequent accurate calculations of spectroscopic properties. The refined FF parameters were first tested by comparing the dihedral potential derived from the FF method to the density functional theory (DFT) based dihedral potential. They were then validated by assessing the relative energies of conformers optimized using both FF and DFT methods, and by comparing the transition wavelengths calculated based on geometries optimized with both FF and DFT approaches. Importantly, the errors in dihedral potential were kept below 1 kcal/mol, and the discrepancies in transition energies were less than 0.1 eV for molecular transitions around 5 eV. This FF parametrization methodology was used in research studies focus- ing on two classes of supramolecular systems: host-guest chemistry related to neurodegenerative diseases, and self-assembly systems for material development. Specifically, we examine host-guest interactions involving proteins such as amyloid-beta, tau, and transthyretin (TTR), which are associated with neu- rodegenerative diseases. Various fluorescent ligands are used for the detection of these proteins in pathological samples. Our results for these protein–ligand systems propose strong binding sites and modes, and include estimations of binding energies for different ligands interacting with the targeted proteins. Additionally, comparative studies among the ligands have been conducted. Interestingly, no fluorescence was observed when low binding energy ligands interacted with amyloid fibrils. In the case of bTVBT4 binding to tau associated with Alzheimer’s disease (AD), a unique binding site was identified. This site was not accessible in the tau fold found in Pick’s disease (PiD), thus explaining the specificity of bTVBT4 for AD-related tau. For self-assembly systems, our findings encompass spectral profiles altered by tyrosine substitutions in oligothiophenes, a stable self-assembly model formed by chiral sulfonimidamides that explains the involved interactions, and comparisons of experimental circular dichroism (CD) profiles to assign isomers of [4]cyclonaphthodithiophene diimides to specific spectral profiles. We also investigated the solvent effects on the spectroscopic properties of symmetric and asymmetric azaoxahelicenes. In conclusion, the methodological development of FF parameters provides a robust framework for accurately modeling the behavior of complex supramolecular systems. The improvements in the dihedral potential align closely with DFT-based calculations, thereby elevating the predictive power of MD simulations for both binding modes and subsequent spectroscopic properties. The research has direct applications in the detection of neurodegenerative diseases at an early stage by designing fluorescent ligands that specifically bind to targeted proteins. It also contributes to the creation of advanced materials with finely-tuned properties. Furthermore, the methodology employed can serve as a blueprint for future studies aiming to refine computational models for other classes of molecules.
|
|
10. |
- De Oliveira, Danilo Hirabae, et al.
(author)
-
Silk Assembly against Hydrophobic Surfaces?Modeling and Imaging of Formation of Nanofibrils
- 2023
-
In: ACS Applied Bio Materials. - : AMER CHEMICAL SOC. - 2576-6422. ; 6:3, s. 1011-1018
-
Journal article (peer-reviewed)abstract
- A detailed insight about the molecular organization behind spider silk assembly is valuable for the decoding of the unique properties of silk. The recombinant partial spider silk protein 4RepCT contains four poly-alanine/glycine-rich repeats followed by an amphiphilic C-terminal domain and has shown the capacity to self-assemble into fibrils on hydrophobic surfaces. We herein use molecular dynamic simulations to address the structure of 4RepCT and its different parts on hydrophobic versus hydrophilic surfaces. When 4RepCT is placed in a wing arrangement model and periodically repeated on a hydrophobic surface, fi-sheet structures of the poly-alanine repeats are preserved, while the CT part is settled on top, presenting a fibril with a height of similar to 7 nm and a width of similar to 11 nm. Both atomic force microscopy and cryo-electron microscopy imaging support this model as a possible fibril formation on hydrophobic surfaces. These results contribute to the understanding of silk assembly and alignment mechanism onto hydrophobic surfaces.
|
|