No description available.

A molecule can crystallise in more than one crystal structure, a common phenomenon in organic compounds known as polymorphism. Different polymorphic forms may have significantly different physical properties, and a reliable prediction would be beneficial to the pharmaceutical industry. However, crystal structure prediction (CSP) based on the knowledge of the chemical structure had long been considered impossible. Previous failures of some CSP attempts led to speculation that the thermodynamic calculations in CSP methodologies failed to predict the kinetically favoured structures. Similarly, regarding the stabilities of co-crystals relative to their pure components, the results from lattice energy calculations and full CSP studies were inconclusive. In this thesis, these problems are addressed using the state-of-the-art CSP methodology implemented in the GRACE software. Firstly, it is shown that the low-energy predicted structures of four organic molecules, which have previously been considered difficult for CSP, correspond to their experimental structures. The possible outcomes of crystallisation can be reliably predicted by sufficiently accurate thermodynamic calculations. Then, the polymorphism of 5- chloroaspirin is investigated theoretically. The order of polymorph stability is predicted correctly and the isostructural relationships between a number of predicted structures and the experimental structures of other aspirin derivatives are established. Regarding the stabilities of co-crystals, 99 out of 102 co-crystals and salts of nicotinamide, isonicotinamide and picolinamide reported in the Cambridge Structural Database (CSD) are found to be more stable than their corresponding co-formers. Finally, full CSP studies of two co-crystal systems are conducted to explain why the co-crystals are not easily obtained experimentally. / University of Bradford

Covariant density functional theory (CDFT) is one of the modern theoretical tools for describing the nuclear structure physics of finite nuclei. Its performance is defined by underlying covariant energy density functionals (CEDFs). In this dissertation and within the framework of the CDFT, different physical properties of the ground and the excited states of rotating and non-rotating nuclei have been investigated. A systematic global investigation of pairing properties based on all available experimental data on pairing indicators has been performed for the first time in the framework of covariant density functional theory. It is based on separable pairing interaction of Ref.\ \cite{TMR.09}. The optimization of the scaling factors of this interaction to experimental data clearly reveals its isospin dependence in the neutron subsystem. However, the situation is less certain in the proton subsystem since similar accuracy of the description of pairing indicators can be achieved both with isospin-dependent and mass-dependent scaling factors. The differences in the functional dependencies of scaling factors lead to the uncertainties in the prediction of proton and neutron pairing properties which are especially pronounced at high isospin and could have a significant impact on some physical observables. Although the present investigation is based on the NL5(E) covariant energy density functional (CEDF), its general conclusions are expected to be valid also for other CEDFs built at the Hartree level. It is shown for the first time that rotational bands which are proton unbound at zero or low spins can be transformed into proton bound ones at high spin by collective rotation of nuclear systems. This is due to strong Coriolis interaction, which acts on high-$N$ or strongly mixed M orbitals and drives the highest in energy occupied single-particle states of nucleonic configurations into the negative energy domain. Proton emission from such proton bound rotational states is suppressed by the disappearance of static pairing correlations at high spins of interest. These physical mechanisms lead to a substantial extension of the nuclear landscape beyond the spin zero proton drip line. In addition, a new phenomenon of the formation of giant proton halos in rotating nuclei emerges: it is triggered by the occupation of strongly mixed M intruder orbitals. Possible experimental fingerprints of the transition from particle bound to particle unbound part of rotational bands are discussed and compared for proton and neutron rich nuclei near and beyond respective drip lines.

<p dir="ltr">Nanomaterials exhibit many unique properties compared to traditional bulk materials, interfaces play a more important role in nanoscale systems by significantly influencing the mechanical performance. In this thesis, we focus on an intricate exploration of various interfaces, ranging from simple GBs in bicrystal models to intricate GB networks within polycrystalline structures and interfaces within nanocomposite materials. Various computational methodologies, including MD, DFT, and advanced machine learning algorithms<del>,</del> were employed to simulate and predict the mechanical properties of interfaces with microstructural complexity.</p><p dir="ltr">Firstly, utilizing MC/MD simulations, we established a distinct correlation between GB motion in the Cantor alloy and the Cr concentration within the GBs. A formulation is calculated to link the GB mobility with the Cr concentration. Subsequently, DFT simulations highlight that vacancies in Tungsten GBs prefer to appear in the layer adjacent to the GB plane rather than the GB plane itself. These vacancies, as the findings suggest, cause the strength to decrease under tensile loading. Then, to expedite the prediction of interfacial properties, a cGAN model was developed to predict GB network evolution in polycrystalline samples based on the training data of MD simulation results. Finally, two modified deep learning models are introduced including the CNN-Prob and FNN-Prob, to predict the yield stress of a composite material, Cu-Cu/Zr. These models encompassed dual components for predicting both mean values and associated standard deviations.</p>

Utilizing density functional theory (DFT) the Vienna Ab initio Simulation Package (VASP) was used to calculate the Raman spectra for five single-walled carbon nanotubes (SWCNTs) with chiralities (4,4), (6,6), (8,0), (12,0) and (7,1). The radial breathing mode (RBM), when compared with experimental frequencies, shows good correlation. When compared to RBM:s calculated with tight binding the frequencies calculated with DFT displayed higher accuracy. The precision of G-band frequencies were inconclusive due to lack of experimental data. The frequencies did not agree well with the results from tight-binding theory. The correctness of the Raman activity estimations using results from DFT calculations was found to be questionable. An unknown mode, which was found to be highly Raman active in the calculated spectra of (4,4), (6,6), and possibly (8,0), and (12,0), is also discussed. It was concluded that further calculations on larger tubes, especially armchair tubes are relevant for future studies. Further verification of the determination of Raman activity is also needed. / F7042T - Project in Engineering Physics

Computational modeling of organic interface formation on semiconductors poses a challenge to a density functional theory-based description due to structural and chemical complexity. A hierarchical approach is presented, where parts of the interface are successively removed in order to increase computational efficiency while maintaining the necessary accuracy. First, a benchmark is performed to probe the validity of this approach for three model reactions and five dispersion corrected density functionals. Reaction energies are generally well reproduced by generalized gradient approximation-type functionals but accurate reaction barriers require the use of hybrid functionals. Best performance is found for the model system that does not explicitly consider the substrate but includes its templating effects. Finally, this efficient model is used to provide coverage dependent reaction energies and suggest synthetic principles for the prevention of unwanted growth termination reactions for organic layers on semiconductor surfaces.

In much of the literature about organic photovoltaics, the topic is framed within the current landscape of energy production and the research on these materials is cited as a possible solution to the energy crisis looming ahead. Despite being the most frequent, this is by no means the only perspective that can be offered. Indeed, the same research may also be set within the larger perspective offered by the field of functional materials. These materials are usually exploited for their particular responses to electrical, magnetic and chemical stimuli and are at the basis of many technologies fundamental to our society. The prominent position of functional materials in modern science is due to the emergence of novel technological needs that such materials have been able to satisfy thanks to their peculiar properties. These properties have been rationalised and mastered by expanding the theoretical description of the underlying physical mechanisms. This theoretical body, combined with the growth and diffusion of computational capabilities has fostered a change in the scientific paradigm underpinning the research effort. More and more, the predictive power of numerical approaches is exploited to lead the way in the exploration of the immense chemical space. The ultimate promise is to achieve the purpose-driven design of compounds thanks to which the molecular structure can be engineered before the actual synthesis to meet the demands dictated by a specific application. To fulfil this role, computational approaches need to be able to simulate the solid state properties at the most relevant time and length scales. If this can be accomplished then a reliable prediction of the performance can be achieved. The current work deals with the development and application of one such protocol, for the particular case of organic photovoltaic semiconductors. Given the specific application, the properties targeted are light absorption and charge transport. Particular effort is put in the simulation of local morphologies at scales above the molecular one to describe supramolecular organisation with sufficient resolution. In this thesis, the protocol is applied to two molecular systems employed in solar devices. Both systems have been selected on the basis of data suggesting that a detailed microscopic description of their behaviour could be highly informative about the aspects responsible for their photovoltaic performance. In particular, chapter 3 details the investigation of a small-molecule donor that has been shown in the literature to have a remarkable behaviour in absorption. While chapter 4 reports the study of a donor-acceptor dyad used as active layer in single-component solar devices with relatively high conversion efficiency. In both cases, the computational protocol has proven capable of providing a detailed microscopic description of the systems. The picture drawn has allowed to clarify the plausible mechanisms behind the observations and to rationalise these behaviours in a broader and more general theoretical framework.

The work in this thesis is basically divided into two related but separate investigations. The first part treats simple chemical reactions of adsorbate molecules on metallic surfaces, induced by means of a scanning tunneling probe (STM). The investigation serves as a parameter free extension to existing theories. The theoretical framework is based on a combination of density functional theory (DFT) and non-equilibrium Green's functions (NEGF). Tunneling electrons that pass the adsorbate molecule are assumed to heat up the molecule, and excite vibrations that directly correspond to the reaction coordinate. The theory is demonstrated for an OD molecule adsorbed on a bridge site on a Cu(110) surface, and critically compared to the corresponding experimental results. Both reaction rates and pathways are deduced, opening up the understanding of energy transfer between different configurational geometries, and suggests a deeper insight, and ultimately a higher control of the behaviour of adsorbate molecules on surfaces. The second part describes a method to calculate STM images in the low bias regime in order to overcome the limitations of localized orbital DFT in the weak coupling limit, i.e., for large vacuum gaps between a tip and the adsorbate molecule. The theory is based on Bardeen's approach to tunneling, where the orbitals computed by DFT are used together with the single-particle Green's function formalism, to accurately describe the orbitals far away from the surface/tip. In particular, the theory successfully reproduces the experimentally well-observed characteristic dip in the tunneling current for a carbon monoxide (CO) molecule adsorbed on a Cu(111) surface. Constant height/current STM images provide direct comparisons to experiments, and from the developed method further insights into elastic tunneling are gained.

This dissertation focuses on computational study of the geometry and energetics small molecules and nanoclusters involving transition metals (TM). These clusters may be used for various industrial applications including catalysis and photonics. Specifically, in this work we have studied hydrides and carbides of 3d-transition metal systems (Sc through Cu), small nickel and gold clusters. Qualitatively correct description of the bond dissociation is ensured by allowing the spatial and spin symmetry to break. We have tested applicability of new exchange-correlation functional and alternative theoretical descriptions (spin-contamination correction in broken symmetry DFT and ensemble Kohn-Sham (EKS)) as well. We studies TM hydrides and carbides systems to understand the importance of underlying phenomenon of bond breaking in catalytic processes. We have tested several exchange-correlation functionals including explicit dependence on kinetic energy density for the description of hydrides (both neutral and cationic) and carbides formed by 3d-transition metals. We find M05-2x and BMK dissociation energies are in better agreement with experiment (where available) than those obtained with high level wavefunction theory methods, published previously. This agreement with experiment deteriorates quickly for other functionals when the fraction of the Hartree-Fock exchange in DFT functional is decreased. Higher fraction of HF exchange is also essential in EKS formalism, but it does not help when spin-adapted unrestricted approach is employed. We analyze the electron spin densities using Natural Bond Orbital population analysis and find that simple description of 3d electrons as non-bonding in character is rarely correct. Unrestricted formalism results in appreciable spin-contamination for some of the systems at equilibrium, which motivated us to investigate it further in details. In order to correct the spin contamination effect on the energies, we propose a new scheme to correct for spin contamination arising in broken-symmetry DFT approach. Unlike conventional schemes, our spin correction is introduced for each spin-polarized electron pair individually and therefore is expected to yield more accurate energy values. We derive an expression to extract the energy of the pure singlet state from the energy of the broken-symmetry DFT description of the low spin state and the energies of the high spin states (pentuplet and two spin-contaminated triplets in the case of two spin-polarized electron pairs). We validate our spin-contamination correction approach by a simple example of H2 and applied to more complex MnH system. Ensemble KS formalism is also applied to investigate the dissociation of C2 molecule. We find that high fraction of HF exchange is essential to reproduce the results of EKS treatment with exact exchange-correlation functional. We analyze the geometry and energetics of small nickel clusters (Ni2-Ni5) for several lowest energy isomers. We also study all possible spin states of small nickel cluster isomers and report observed trends in energetics. Finally we determine the geometry and energetics of ten lowest energy isomers of four small gold clusters (Au2, Au4, Au6, and Au8). We have also investigated the influence of cluster geometry, ligation, solvation and relativistic effects on electronic structure of these gold clusters. The effect of one-by-one ligand attachment in vacuum and solvent environment is also studied. Performance of five DFT functionals are tested as well; Local Spin Density Approximation (SVWN5), Generalized Gradient Approximation (PBE), kinetic energy density-dependent functional (TPSS), hybrid DFT (B3LYP), and CAM-B3LYP which accounts for long-range exchange effects believed to be important in the analysis of metal bonding in gold complexes and clusters. Our results exhibit the ligand induced stability enhancement of otherwise less stable isomers of Au4, Au6 and Au8. Ligands are found to play a crucial role in determining the 2D to 3D transition realized in small gold clusters. In order to select an appropriate theory level to use in this study, we investigate the effect of attachment of four different ligands (NH3, NMe3, PH3, PMe3) on cluster geometry and energetics of Au2 and Au4 in vacuum and in solution. Our results benchmark the applicability of DFT functional model and polarization functions in the basis set for calculations of ligated gold cluster systems. We employ five different basis sets with increasing amount of polarization and diffuse functions; LANL2DZ, LANL2DZ-P, def2-SVP, def2-TZVP, and def2-QZVP. We obtain NMe3 = NH3 > PH3 > PMe3 order of ligand binding energies and observe shallow potential energy surfaces in all molecules. Our results suggest appropriate quantum-chemical methodologies to model small noble metal clusters in realistic ligand environment to provide reliable theoretical analysis in order to complement experiments.

Covariant density functional theory (CDFT) is a modern theoretical tool for the description of nuclear structure phenomena. Different physical observables of the ground and excited states in even-even nuclei have been studied within the CDFT framework employing three major classes of the state-of-the-art covariant energy density functionals. The global assessment of the accuracy of the description of the ground state properties and systematic theoretical uncertainties of atomic nuclei have been investigated. Large-scale axial relativistic Hartree-Bogoliubov (RHB) calculations are performed for all Z < 106 even-even nuclei between the two-proton and two-neutron drip lines. The sources of theoretical uncertainties in the prediction of the two-neutron drip line are analyzed in the framework of CDFT. We concentrate on single-particle and pairing properties as potential sources of these uncertainties. The major source of these uncertainties can be traced back to the differences in the underlying single-particle structure of the various CEDFs. A systematic search for axial octupole deformation in the actinides and superheavy nuclei with proton numbers Z = 88 - 126 and neutron numbers from two-proton drip line up to N = 210 has been performed in CDFT. The nuclei in the Z ~ 96, N ~ 196 region of octupole deformation have been investigated in detail and their systematic uncertainties have been quantified. The structure of superheavy nuclei has been reanalyzed with inclusion of quadrupole deformation. Theoretical uncertainties in the predictions of inner fission barrier heights in superheavy elements have been investigated in a systematic way. The correlations between global description of the ground state properties and nuclear matter properties have been studied. It was concluded that the strict enforcement of the constraints on the nuclear matter properties (NMP) defined in Ref. [1] will not necessary lead to the functionals with good description of ground state properties. The different aspects of the existence and stability of hyperheavy nuclei have been investigated. For the first time, we demonstrate the existence of three regions of spherical hyperheavy nuclei centered around (Z ~ 138, N ~ 230), (Z ~ 156, N ~ 310) and (Z ~ 174, N ~ 410) which are expected to be reasonably stable against spontaneous fission.