Finite difference calculations of permeability in large domains in a wide porosity range.

Determining effective hydraulic, thermal, mechanical and electrical properties of porous materials by means of classical physical experiments is often time-consuming and expensive. Thus, accurate numerical calculations of material properties are of increasing interest in geophysical, manufacturing, bio-mechanical and environmental applications, among other fields. Characteristic material properties (e.g. intrinsic permeability, thermal conductivity and elastic moduli) depend on morphological details on the porescale such as shape and size of pores and pore throats or cracks. To obtain reliable predictions of these properties it is necessary to perform numerical analyses of sufficiently large unit cells. Such representative volume elements require optimized numerical simulation techniques. Current state-of-the-art simulation tools to calculate effective permeabilities of porous materials are based on various methods, e.g. lattice Boltzmann, finite volumes or explicit jump Stokes methods. All approaches still have limitations in the maximum size of the simulation domain. In response to these deficits of the well-established methods we propose an efficient and reliable numerical method which allows to calculate intrinsic permeabilities directly from voxel-based data obtained from 3D imaging techniques like X-ray microtomography. We present a modelling framework based on a parallel finite differences solver, allowing the calculation of large domains with relative low computing requirements (i.e. desktop computers). The presented method is validated in a diverse selection of materials, obtaining accurate results for a large range of porosities, wider than the ranges previously reported. Ongoing work includes the estimation of other effective properties of porous media.

Ligand effects on the structure and vibrational properties of the thiolated Au18 cluster

Most of the studies devoted to thiolated gold clusters suppose that their core and Au-S framework do not suffer from distortion independently of the protecting ligands (-SR) and it is assumed as correct to simplify the ligand as SCH3. In this work is delivered a systematic study of the structure and vibrational properties (IR and Raman) of the Au18(SR)14 cluster. The pursued goal is to understand the dependency of the displayed vibrational properties of the thiolated Au18 cluster with the ligands type. A set of six ligands was considered during calculations of the vibrational properties based on density functional theory (DFT) and in its dispersioncorrected approach (DFT-D)

Thermodynamic calculations for the salt crystallisation damage in porous built heritage using PHREEQC

This work considers the crystallisation mechanisms of the most common and aggressive salts that generate stress in porous building stones as a result of changing ambient conditions. These mechanisms include the salt crystallisation that result from decreasing relative humidity and changes in temperature and, in hydrated salts, the dissolution of the lower hydrated form and the subsequent precipitation of the hydrated salt. We propose a new methodology for thermodynamic calculations using PHREEQC that includes these crystallisation mechanisms. This approach permits the calculation of the equilibrium relative humidity and the parameterization of the critical relative humidity and crystallisation pressures for the dissolution–precipitation transitions. The influence of other salts on the effectives of salt crystallisation and chemical weathering is also assessed. We review the sodium and magnesium sulphate and sodium chloride systems, in both single and multicomponent solutions, and they are compared to the sodium carbonate and calcium carbonate systems. The variation of crystallisation pressure, the formation of new minerals and the chemical dissolution by the presence of other salts is also evaluated. Results for hydrated salt systems show that high crystallisation pressures are possible as lower hydrated salts dissolve and more hydrated salts precipitate. High stresses may be also produced by decreasing temperature, although it requires that porous materials are wet for long periods of time. The presence of other salts changes the temperature and relative humidity of salt transitions that generates stress rather than reducing the pressure of crystallisation, if any salt has previously precipitated. Several practical conclusions derive from proposed methodology and provide conservators and architects with information on the potential weathering activity of soluble salts. Furthermore, the model calculations might be coupled with projections of future climate to give as improved understanding of the likely changes in the frequency of phase transitions in salts within porous stone.

Double Hybrid Functionals and the Π-System Bond Length Alternation Challenge: Rivaling Accuracy of Post-HF Methods

Predicting accurate bond length alternations (BLAs) in long conjugated oligomers has been a significant challenge for electronic-structure methods for many decades, made particularly important by the close relationships between BLA and the rich optoelectronic properties of π-delocalized systems. Here, we test the accuracy of recently developed, and increasingly popular, double hybrid (DH) functionals, positioned at the top of Jacobs Ladder of DFT methods of increasing sophistication, computational cost, and accuracy, due to incorporation of MP2 correlation energy. Our test systems comprise oligomeric series of polyacetylene, polymethineimine, and polysilaacetylene up to six units long. MP2 calculations reveal a pronounced shift in BLAs between the 6-31G(d) basis set used in many studies of BLA to date and the larger cc-pVTZ basis set, but only modest shifts between cc-pVTZ and aug-cc-pVQZ results. We hence perform new reference CCSD(T)/cc-pVTZ calculations for all three series of oligomers against which we assess the performance of several families of DH functionals based on BLYP, PBE, and TPSS, along with lower-rung relatives including global- and range-separated hybrids. Our results show that DH functionals systematically improve the accuracy of BLAs relative to single hybrid functionals. xDH-PBE0 (N4 scaling using SOS-MP2) emerges as a DH functional rivaling the BLA accuracy of SCS-MP2 (N5 scaling), which was found to offer the best compromise between computational cost and accuracy the last time the BLA accuracy of DFT- and wave function-based methods was systematically investigated. Interestingly, xDH-PBE0 (XYG3), which differs to other DHs in that its MP2 term uses PBE0 (B3LYP) orbitals that are not self-consistent with the DH functional, is an outlier of trends of decreasing average BLA errors with increasing fractions of MP2 correlation and HF exchange.

Theoretical Rationalization of the Singlet–Triplet Gap in OLEDs Materials: Impact of Charge-Transfer Character

New materials for OLED applications with low singlet–triplet energy splitting have been recently synthesized in order to allow for the conversion of triplet into singlet excitons (emitting light) via a Thermally Activated Delayed Fluorescence (TADF) process, which involves excited-states with a non-negligible amount of Charge-Transfer (CT). The accurate modeling of these states with Time-Dependent Density Functional Theory (TD-DFT), the most used method so far because of the favorable trade-off between accuracy and computational cost, is however particularly challenging. We carefully address this issue here by considering materials with small (high) singlet–triplet gap acting as emitter (host) in OLEDs and by comparing the accuracy of TD-DFT and the corresponding Tamm-Dancoff Approximation (TDA), which is found to greatly reduce error bars with respect to experiments thanks to better estimates for the lowest singlet–triplet transition. Finally, we quantitatively correlate the singlet–triplet splitting values with the extent of CT, using for it a simple metric extracted from calculations with double-hybrid functionals, that might be applied in further molecular engineering studies.

A THEORETICAL STUDY OF INTERACTION OF NANOPARTICLES WITH BIOMOLECULE

Many types of materials at nanoscale are currently being used in everyday life. The production and use of such products based on engineered nanomaterials have raised concerns of the possible risks and hazards associated with these nanomaterials. In order to evaluate and gain a better understanding of their effects on living organisms, we have performed first-principles quantum mechanical calculations and molecular dynamics simulations. Specifically, we will investigate the interaction of nanomaterials including semiconducting quantum dots and metallic nanoparticles with various biological molecules, such as dopamine, DNA nucleobases and lipid membranes. Firstly, interactions of semiconducting CdSe/CdS quantum dots (QDs) with the dopamine and the DNA nucleobase molecules are investigated using similar quantum mechanical approach to the one used for the metallic nanoparticles. A variety of interaction sites are explored. Our results show that small-sized Cd4Se4 and Cd4S4 QDs interact strongly with the DNA nucleobase if a DNA nucleobase has the amide or hydroxyl chemical group. These results indicate that these QDs are suitable for detecting subcellular structures, as also reported by experiments. The next two chapters describe a preparation required for the simulation of nanoparticles interacting with membranes leading to accurate structure models for the membranes. We develop a method for the molecular crystalline structure prediction of 1,2-Dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC), 1,2-Dimyristoyl-sn-glycero-3-phosphorylethanolamine (DMPE) and cyclic di-amino acid peptide using first-principles methods. Since an accurate determination of the structure of an organic crystal is usually an extremely difficult task due to availability of the large number of its conformers, we propose a new computational scheme by applying knowledge of symmetry, structural chemistry and chemical bonding to reduce the sampling size of the conformation space. The interaction of metal nanoparticles with cell membranes is finally carried out by molecular dynamics simulations, and the results are reported in the last chapter. A new force field is developed which accurately describes the interaction forces between the clusters representing small-sized metal nanoparticles and the lipid bilayer molecules. The permeation of nanoparticles into the cell membrane is analyzed together with the RMSD values of the membrane modeled by a lipid bilayer. The simulation results suggest that the AgNPs could cause the same amount of deformation as the AuNPs for the dysfunction of the membrane.

Enhanced oxidation activity from modified ceria: MnOx-ceria, CrOx-ceria and Mg doped VOx-ceria

Ceria is an important component of catalysts for oxidation reactions that proceed through the Mars-van Krevelen mechanism, promoting activity. A paradigm example of this is the VOx–CeO2 system for oxidative dehydrogenation reactions, where vanadium oxide species are supported on ceria and a special synergy between them is behind the enhanced activity: reduction of the catalyst is promoted by ceria undergoing reduction. This leads to favourable oxygen vacancy formation and hydrogen adsorption energies—useful descriptors for the oxidation activity of VOx–CeO2 catalysts. In this paper, we examine if this promoting effect on ceria-based catalysts holds for other metal oxide modifiers and we investigate MnOn– and CrOn–CeO2(111) (n = 0 − 4) as examples. We show, combining density functional theory calculations and statistical thermodynamics that similarly to the vanadia modifier, the stable species in each case is MnO2– and CrO2–CeO2. Both show favourable energetics for oxygen vacancy formation and hydrogen adsorption, indicating that VO2–CeO2 is not the only system of this type that can have an enhanced activity for oxidation reactions. However, the mechanism involved in each case is different: CrO2–CeO2 shows similar properties to VO2–CeO2 with ceria reduction upon oxygen removal stabilising the 5+ oxidation state of Cr. In contrast, with MnO2–CeO2, Mn is preferentially reduced. Finally, a model system of VO2–Mg:CeO2 is explored that shows a synergy between VO2 modification and Mg doping. These results shed light on the factors involved in active oxidation catalysts based on supported metal oxides on ceria that should be taken into consideration in a rational design of such catalysts.

First-Principles Calculations of the Electronic and Structural Properties of GaSb

In this paper, we carried out first-principles calculations in order to investigate the structural and electronic properties of the binary compound gallium antimonide (GaSb). This theoretical study was carried out using the Density Functional Theory within the plane-wave pseudopotential method. The effects ofexchange and correlation (XC) were treated using the functional Local Density Approximation (LDA), generalized gradient approximation (GGA): Perdew–Burke–Ernzerhof (PBE), Perdew-Burke-Ernzerhof revised for solids (PBEsol), Perdew-Wang91 (PW91), revised Perdew–Burke–Ernzerhof (rPBE), Armiento–Mattson 2005 (AM05) and meta-generalized gradient approximation (meta-GGA): Tao–Perdew– Staroverov–Scuseria (TPSS) and revised Tao–Perdew–Staroverov–Scuseria (RTPSS) and modified Becke-Johnson (MBJ). We calculated the densities of state (DOS) and band structure with different XC potentials identified and compared them with the theoretical and experimental results reported in the literature. It was discovered that functional: LDA, PBEsol, AM05 and RTPSS provide the best results to calculate the lattice parameters (a) and bulk modulus (B0); while for the cohesive energy (Ecoh), functional: AM05, RTPSS and PW91 are closer to the values obtained experimentally. The MBJ, Rtpss and AM05 values found for the band gap energy is slightly underestimated with those values reported experimentally.

DFT STUDIES ON THE MECHANISM OF RING-OPENING POLYMERIZATION OF LACTIDE WITH ADENINE

The use of organic molecules as catalysts for the ring-opening polymerization (ROP) of cyclic esters has gained much interest last years.[1] The use of a molecule of biological interest, able to initiate ROP of cyclic esters without any cocatalyst is even more interesting, as the resulting material will not contain any catalytic residue. Nucleobase-polymer conjugates development is thus an emerging area envisaging biomedical applications.[2] However, they are usually synthesized by tedious multistep procedures. Recently, adenine was used as organoinitiator for the ROP of L-lactide.[3] Reaction conditions involving short reaction times and relatively low temperatures enable the access to adenine-polylactide(Adn-PLA)conjugates in a simple one-step procedure, without additional catalyst and in the absence of solvent. In this study, computational investigations with density functional theory (DFT) were performed in order to clarify the reaction mechanism leading to the desired Adn-PLA. The results show that a hydrogen bond catalytic mechanism, involving a nucleophilic attack of the activated amine group of adenine onto the carbonyl group of lactide, seem to be plausible.

A DFT study on chemodivergent preparation of piperidinic and morpholinic heterocycles

The topic of this thesis is the DFT computational study of the mechanisms for the synthesis of chiral 3,4,5-trisubstituted piperidines and 2,6-disubstituted morpholines. The goal of this synthesis is to use, the same substrate containing two electrophilic sites: an α,β-unsaturated ester and a ketone, which evolve according to the nucleophile used (cyanide, phenyl sulfide) through different addition and cyclization reactions. A quaternary ammonium salt is used as a catalyst for these reactions, which leads to a diastereoisomeric excess both for the reactions of morpholine and piperidine products. Studies in silico of the pathways of these reactions explain the chemoselection and diasteroselection deriving from the two nucleophiles used. In this case of piperidine products, it was also possible to validate the hypothesis of a concerted nucleophilic addition mechanism on the α,β-unsaturated site and cyclization due to an intramolecular Michael addition.

Studio computazionale DFT e TD-DFT dei dimeri della Rodamina-B

L’oggetto di questo elaborato è lo studio computazionale, a livello della teoria del funzionale della densità (DFT) e della sua formulazione dipendente dal tempo (TD-DFT), dei dimeri della molecola di rodamina-B, parallelo allo sviluppo di una procedura di tuning ottimale del funzionale CAM-B3LYP. Questa molecola, che assume notevole rilevanza nei sistemi light harvesting grazie alle sue proprietà fotochimiche di emissione nel visibile, è impiegata nella sintesi di nanoparticelle (NPs) fluorescenti in ambito di diagnostica medica e bio-imaging, che sfruttano il fenomeno di trasferimento di energia per risonanza (FRET). Per via della notevole importanza che questa molecola riveste nell’ambito della fotochimica, essa è stata oggetto di esperimenti del gruppo di ricerca del laboratorio di biofotonica e farmacologia “Nanochemistry and Bioimaging”, che collabora con il gruppo di chimica computazionale dell’area chimico/fisica del Dipartimento. La dimerizzazione della rodamina all’interno delle NPs può innescare canali di self-quenching che abbassano la resa quantica di fluorescenza, pregiudicando l’efficienza dei dispositivi: l’obiettivo dello studio è la caratterizzazione dei dimeri, in solventi e con controioni diversi, impiegando dei modelli molecolari, per identificarne le specie più stabili e descrivere la fotofisica degli stati elettronici eccitati. Il carattere generalmente charge-transfer (CT) di questi stati elettronici richiede un “tuning ottimale” della metodologia computazionale DFT/TD-DFT per una descrizione quantitativa accurata.

Theoretical models and ab initio calculations of phonons

Questa tesi intende approfondire da un punto di vista, sia teorico sia computazionale, le proprietà fondamentali dei fononi. A tal fine, sono presentati i modelli quantistici di Einstein e di Debye che permettono la derivazione analitica degli osservabili macroscopici principali di un solido, come l’energia media e la capacità termica. Ciò è possibile tramite una trattazione meccano-statistica basata sull’approssimazione armonica dei modi normali di vibrazione degli ioni reticolari. Quindi, all’inizio si mostrano brevemente i risultati principali riguardanti l’oscillatore armonico quantistico. Successivamente, si approfondiscono i temi della dispersione fononica e della densità degli stati vibrazionali per reticoli cristallini 1D e 3D. Si ottiene che la prima non può essere considerata lineare se non nel limite di alte lunghezze d’onda, e che la seconda può presentare punti di singolarità correlati alla forma della relazione di dispersione. Infine, sono state svolte alcune analisi computazionali ab initio relative alla dispersione fononica, la densità degli stati vibrazionali e la frequenza di Debye del Carbonio (diamante) tramite i programmi VASP e Phonopy, confrontando i risultati con dati sperimentali presenti in letteratura.

Compliance verification of virtual Phasor Measuring Unit using Ip-DFT algorithm according to the IEEE-60255-118-1(2018).

With the increase in load demand for various sectors, protection and safety of the network are key factors that have to be taken into consideration over the electric grid and distribution network. A phasor Measuring unit is an Intelligent electronics device that collects the data in the form of a real-time synchrophasor with a precise time tag using GPS (Global positioning system) and transfers the data to the grid command to monitor and assess the data. The measurements made by PMU have to be very precise to protect the relays and measuring equipment according to the IEEE 60255-118-1(2018). As a device PMU is very expensive to research and develop new functionalities there is a need to find an alternative to working with. Hence many open source virtual libraries are available to replicate the exact function of PMU in the virtual environment(Software) to continue the research on multiple objectives, providing the very least error results when verified. In this thesis, I executed performance and compliance verification of the virtual PMU which was developed using the I-DFT (Interpolated Discrete Fourier transforms) C-class algorithm in MATLAB. In this thesis, a test environment has been developed in MATLAB and tested the virtually developed PMU on both steady state and dynamic state for verifying the latest standard compliance(IEEE-60255-118-1).

Desorption/ionization Efficiencies Of Triacylglycerols And Phospholipids Via Easi-ms.

Knowledge of the major effects governing desorption/ionization efficiency is required for the development and application of ambient mass spectrometry. Although all triacylglycerols (TAG) have the same favorable protonation and cationization sites, their desorption/ionization efficiencies can vary dramatically during easy ambient sonic-spray ionization because of structural differences in the carbon chain. To quantify this somewhat surprising and drastic effect, we have performed a systematic investigation of desorption/ionization efficiencies as a function of unsaturation and length for TAG as well as for diacylglycerols, monoacylglycerols and several phospholipids (PL). Affinities for Na(+) as a function of unsaturation level have also been assayed via comprehensive metadynamics calculations to understand the influence of this phenomenon on the ionization efficiency. The results suggest that dipole-dipole interactions within a carbon chain tuned by unsaturation sites govern ionization efficiency of TAG and PL.