% Template for Thesis, Universität Bielefeld, Center for Mathematical Economics. % % Dieses Template kann als Grundlage für die Erstellung einer Masterarbeit oder Bachelorarbeit verwendet werden. % 18.02.2021, Annika Kemper, Bielefeld. % \documentclass[12pt, a4paper, english]{report} %--- Layout \usepackage[a4paper, left=2.5cm, right=2.5cm, top=1cm, bottom=1cm, marginparsep=5mm, marginparwidth=10mm, headheight=7mm, headsep=1.5cm, footskip=1.5cm, includeheadfoot]{geometry} %--- CONFIG \usepackage[ngerman]{babel} \usepackage{listings} \usepackage{lipsum} \usepackage[utf8]{inputenc} \usepackage[T1]{fontenc} \usepackage{fancyhdr} \usepackage[onehalfspacing]{setspace} \usepackage{graphicx} \usepackage{svg} \usepackage{amsmath} \usepackage{quiver} \usepackage{tikz-cd} \usepackage{tikzscale} \usepackage[style=numeric, sorting=none,backend=biber]{biblatex} \addbibresource{refs.bib} \pagestyle{fancy} \renewcommand{\headrulewidth}{1pt} %Linie oben \fancyhf{} \fancyhead[R]{\rightmark} %Kopfzeile rechts bzw. außen \renewcommand{\footrulewidth}{0pt} \fancyfoot[C]{\thepage} \lfoot{} \fancypagestyle{plain}{ \pagestyle{fancy} \renewcommand{\headrulewidth}{1pt} %Linie oben \fancyhf{} \fancyhead[R]{\rightmark} \renewcommand{\footrulewidth}{0pt} \fancyfoot[C]{\thepage} \lfoot{}} \usepackage[bottom]{footmisc} \renewcommand{\chaptermark}[1]{\markboth{\textsc{\thechapter.\ ~~#1}}{\textsc{\thechapter.\ ~~#1}}} \renewcommand{\sectionmark}[1]{\markboth{\textsc{\thesection.\ ~~#1}}{\textsc{\thesection.\ ~~#1}}} \usepackage{xcolor} \makeatletter \definecolor{codegreen}{rgb}{0,0.6,0} \definecolor{codegray}{rgb}{0.5,0.5,0.5} \definecolor{codepurple}{rgb}{0.58,0,0.82} \definecolor{backcolour}{rgb}{0.95,0.95,0.92} \lstdefinestyle{mystyle}{ backgroundcolor=\color{backcolour}, commentstyle=\color{codegreen}, keywordstyle=\color{magenta}, numberstyle=\tiny\color{codegray}, stringstyle=\color{codepurple}, basicstyle=\ttfamily\footnotesize, breakatwhitespace=false, breaklines=true, captionpos=b, keepspaces=true, numbers=left, numbersep=5pt, showspaces=false, showstringspaces=false, showtabs=false, tabsize=2 } \lstset{style=mystyle} \def\@makechapterhead#1{% \vspace*{10\p@}% % {\parindent \z@ \centering \reset@font {\parindent \z@ \reset@font \par\nobreak \vspace*{2\p@}% {\Huge \bfseries \thechapter\quad #1\par\nobreak} \par\nobreak \vspace*{2\p@}% \vskip 40\p@ %\vskip 100\p@ }} \def\@makeschapterhead#1{% \vspace*{10\p@}% {\parindent \z@ \reset@font \par\nobreak \vspace*{2\p@}% {\Huge \bfseries #1\par\nobreak} \par\nobreak \vspace*{2\p@}% \vskip 40\p@ }} \makeatother %%--- Mathematischer Formelsatz \usepackage {mathrsfs, amssymb} \usepackage {amsmath} \usepackage {amsthm} \usepackage{amsfonts} \usepackage{latexsym} \usepackage{amssymb} \usepackage[section]{placeins} %--- add hyperref %\usepackage[hyphens]{url} \usepackage{hyperref} \urlstyle{same} %--- add Theorem-Style \usepackage{amsthm} \newtheorem{definition}{Definition} \newtheorem{prop}{Proposition} \newtheorem{theorem}{Theorem} \newtheorem*{remark}{Remark} %--- change indent of table of content \setcounter{secnumdepth}{5} \setcounter{tocdepth}{5} %\renewcommand*\thesubsubsection{\Alph{subsubsection}} \usepackage[titles]{tocloft} \cftsetindents{chapter}{0cm}{.5cm} \cftsetindents{section}{0.5 cm}{0.8cm} \cftsetindents{subsection}{0.5cm}{1.1cm} \cftsetindents{subsubsection}{0.5 cm}{1.4cm} %\renewcommand{\listfigurename}{\textsc{List of Figures}} %--- Quotes \usepackage[babel,english=american]{csquotes} %\usepackage[square,comma,numbers,authoryear]{natbib} %%================================================================================== \begin{document} \thispagestyle{empty} \begin{center} \large{Bielefeld University}\\ \large{Faculty of Technology }\\ \vspace{2cm} \LARGE{{\bf \textcolor{black}{Graph Rewriting and NMR Spectroscopy in Chemical Reaction Network Modelling}} } \vspace{0.4 cm} \LARGE{{\bf \textcolor{black}{[Graph Transformationen und NMR-Spektroskopie in Chemischen Reaktions Netzwerken]}} } \vspace{1.5cm} \Large{Masterarbeit} % \Large{Bachelorarbeit} \vspace{0.2cm} \Large{\today} \vspace{1 cm} \Large{Kilian Seng} \vspace{0.1cm} \large{4330700} \vspace{3cm} \end{center} \begin{flushleft} \vspace{3cm} \begin{tabular}{ll} \large{Erstprüfer:}& \large{Prof. Dr. Daniel Merkle} \end{tabular} \vspace{0.2 cm} \begin{tabular}{ll} \large{Zweitprüfer:} & \large{Dr. Marcel Friedrichs} \end{tabular} \end{flushleft} %%================================================================================== \newpage \pagenumbering{Roman} \onehalfspacing \phantomsection %\addcontentsline{toc}{chapter}{ Inhaltsverzeichnis} \raggedbottom \pagebreak \thispagestyle{empty} \tableofcontents \raggedbottom \pagebreak \thispagestyle{empty} \clearpage \phantomsection \addcontentsline{toc}{chapter}{Abbildungsverzeichnis} \raggedbottom \pagebreak \thispagestyle{empty} \listoffigures \raggedbottom \pagebreak \thispagestyle{empty} \clearpage \phantomsection %\addcontentsline{toc}{chapter}{Tabellenverzeichnisverzeichnis} \listoftables \raggedbottom \pagebreak \thispagestyle{empty} \newpage \chapter*{Abkürzungsverzeichnis} \textbf{ILP}: Integer-Linear-Programming \newline \textbf{DPO}: Double-Push-Out \newline \textbf{NMR}: Nucleus-Magnetic-Resonance \newline \textbf{MS}: Massenspektroskopie \newline \textbf{IR}: Infrarot \newline \textbf{UV/Vis}: Ultraviolet und Visuelles Licht \newline \textbf{DFT}: Dichte-Funktional-Theorie \newline \textbf{TMS}: Tetramethylsilan \newline \textbf{ppm}: Parts per million \newline \textbf{M\o{}d}: Med\O{}lDatschgerl \newline \textbf{GML}: Game maker language \newline \textbf{SMILES}: Simplified Molecular Input Line Entry System \newline \textbf{Exp}: Experimentell \newpage \pagenumbering{arabic} %%================================================================================== \chapter*{Abstract} Cheminformatic has been established as a valuable tool in the development and analysis of large chemical networks. Opening the possibility to tackle system sizes that would not be feasible for a human researcher. The possibilities both for the development of reactions as well as the the selection and verification of large reaction networks are presented in this thesis. This thesis aims to find sensible synthesis in chemical reaction networks using NMR-spectra and their characteristics. \chapter{Einleitung} Die Chemieinformatik ist ein aufstrebendes Forschungsfeld, das zum Ziel hat die Beschränkungen der chemischen Industrie mit der Rechenstärke der Informatik anzugehen. So sollen neue Reaktion und Moleküle entwickelt, bestehende Prozesse optimiert und Eigenschaften analysiert werden \cite{cheminformatics}. Die Chemieinformatik eignet sich dabei besonders in der Untersuchung großer Systeme , wie der Entwicklung neuer Medikamente \cite{cheminformatics2}. Diese Arbeit thematisiert eine Möglichkeit chemische Reaktionen näher zu analysieren und auf Basis dessen einen Prozess zu optimieren. Dabei werden zueinander ähnliche Moleküle über ihre zugehörigen NMR-Spektren verglichen. Basierend auf diesen Gemeinsamkeiten und Unterschieden wird eine Synthese bestimmt. \section{ILP und Graphtransformation} Das Ziel bei der Synthese chemischer Moleküle ist die Maximierung der Ausbeute und die Minimierung des nötigen Aufwands. Bekannte Reaktionen können dabei auf eine Vielzahl von Startmolekülen angewandt werden. Aber auch für ein gewähltes Startmolekül ist sowohl in der Auswahl der verwendeten Reaktionen, als auch der gewählten Reihenfolge, große Varianz möglich. Um aus diesem Netzwerk an Reaktionen eine optimale Syntheseplanung zu entwickeln eignet sich das Integer-Linear-Programming (ILP). Dabei werden alle über die einzelnen Reaktionsschritte bekannten Informationen in einem Reaktionsnetzwerk zusammengefasst. Zu den geeigneten Informationen zählen bevorzugte Konfigurationen, Energiebarrieren oder Kosten für Katalysatoren und Reaktanten. All diese Informationen können mittels ILP gewichtet zu einem geeigneten Syntheseplan optimiert werden \cite{ilpmot}. Es ist dabei möglich das Reaktionsnetzwerk ausschließlich aus bereits einmal durchgeführten Reaktionen, mit bekannten Edukten und Produkten zusammenzusetzen. In der Chemie werden Reaktionen aber meistens für Kategorien von Molekülen mit bestimmten Eigenschaften angegeben. Dabei kann das Vorliegen einer funktionellen Gruppe oder spezieller Bindungs- und Atomkonfigurationen entscheidend sein, ob eine Reaktion auf die Edukte anwendbar ist und welche Produkte zu erwarten sind. Diese Entscheidung kann bei kleinen Reaktionsnetzwerken händisch durchgeführt werden, steigt in ihrem Aufwand aber schnell an. Eine zeitsparende Alternative ist die Entwicklung des Reaktionsnetzwerks durch Graphtransformationen. Dabei werden Moleküle als Graphen und chemische Reaktionen als Graphentransformation abstrahiert \cite{andersenhyperflow}. So sind Reaktionsnetzwerke mit mehr Reaktionstypen, sowie mehr und größere Moleküle möglich. \section{NMR-Spektroskopie} Die NMR-Spektroskopie ist eine der wichtigsten chemischen Werkzeuge zur Strukturaufklärung unbekannter Moleküle. Die Spektren sind dabei charakteristisch für eine Struktur und ermöglichen so eine Bestätigung des Erfolgs einer geplanten Synthese. Über die vorliegenden Zusammenhänge ist zudem auch eine weitgehende Analyse der Nachbarschaftsbeziehungen in unbekannten Molekülen möglich. Diese Analyse wird meistens mit weiteren Methoden kombiniert, wobei die verschiedenen NMR-Typen einen wesentlichen Anteil ausmachen \cite{ocwiley}. Insbesondere in Artikeln über die Synthese neuer Moleküle wird so meistens das NMR-Spektrum angegeben, um eine externe Verifikation der Synthese durch unabhängige Chemiker zu ermöglichen. Auch bei der eigenen Synthese ist der Vergleich des NMR-Spektrums mit bekannten Spektren der Literatur eine zentrale Möglichkeit der Verifikation. Dieser Vergleich erfolgt meistens empirisch und ist zeitaufwändig, was die Automatisierung dieser Analyse nahelegt. \section{Zielsetzung} Das Ziel dieser Masterarbeit ist die Kombination des mathematischen Integer-Linear-Programming, der informatischen Struktur des Hypergraphen und der chemischen Methode des NMR-Spektrums. Dabei soll für bekannte Reaktionen und Moleküle ein chemisches Reaktionsnetzwerk in der Form eines Hypergraphen entwickelt werden. Für die so gefundenen Moleküle sollen NMR-Spektren bestimmt werden. Diese NMR-Spektren können dabei sowohl durch quantenchemische Programme berechnet werden, als auch aus den Literaturangaben zu experimentell bestimmten Spektren stammen. Diese NMR-Spektren sollen miteinander verglichen werden und ihre Ähnlichkeit, sowie ihre charakteristischen Unterschiede quantifiziert werden. Dieser Grad der Ähnlichkeit soll anschließend genutzt werden, um einen Reaktionspfad zu bestimmen, der diese Ähnlichkeit maximiert. Das angestrebte Ziel ist dabei die automatisierte der Verifizierung von Molekülen. \chapter{Theoretische Grundlagen} Für diese Arbeit werden sowohl Werkzeuge der Graphentheorie (Kapitel \ref{chgraph}), als auch ein Verständnis der NMR-Spektroskopie (Kapitel \ref{chnmr}) und der quantenchemischen Berechnung selbiger benötigt (Kapitel \ref{chdft}), welche im Folgenden näher erläutert werden. \section{Graphen Theorie}\label{chgraph} Die Graphen Theorie befasst sich mit dem Abbilden von Netzwerken und Entitätenbeziehungen. Dabei werden die Netzwerkteilnehmer und Entitäten als Menge von Knoten ($V$) und ihre Verbindungen oder Beziehungen als Menge von Kanten ($E$) abstrahiert. Ein Knoten $v$ ist dabei ein Element der Menge $V$ und eine Kante ein Element der Teilmenge $E \subseteq P_{2}(V)$. Das Tupel $(u,v)$ aus zwei Knoten u und v, mit $u,v \in V$ ist ein Element der Menge $(u,v) \in E$. Ein einfacher Graph ($G$) besteht aus einer Menge von Knoten und Kanten $G = (V, E)$ \cite{bender2006foundations}. Bei einem gerichteten Graph gilt dabei, dass die Kante $(u,v)$ von Knoten $u$ zu Knoten $v$ ein geordnetes Paar darstellt. Es gilt $(u,v) \neq (v,u)$, weshalb beide in der Menge von Kanten vorkommen können. In einem ungerichteten Graphen gilt, dass die beiden Tupel dieselbe Kante beschreiben. Es gilt also $(u,v) = (v,u)$. Ein Multigraph ist darüber definiert, dass anders als in einem einfachen Graphen, eine identische Kante mehrmals enthalten ist. Die Kanten sind also keine Menge. Die beiden Kanten haben dann die selbe Repräsentation $(u,v)\stackrel{r}{=}(u,v)$, aber sind nicht die gleiche Kante $(u,v) \neq (u,v)$. Es besteht auch die Möglichkeit, dass eine Kante von einem Knoten zu sich selber führt, also $(v,v)$, was als Schleife bezeichnet wird. Sowohl gerichtete, als auch ungerichtete Graphen können gelabelt sein. Sie werden dann als gelabelter Graph bezeichnet. Dabei lässt sich zwischen Kanten und Knotenlabeln unterscheiden. Für Kantenlabel (\ref{edgelabel}) und Knotenlabel (\ref{verticelabel}) in einem klassischen gelabelten Graphen gilt, dass jeder Knoten oder jede Kante ein Label aus der Menge der Label $\Omega$ zugeordnet haben muss. \begin{equation} \label{edgelabel} \forall e \in E: \exists label(e) \in \Omega_{E} \end{equation} \begin{equation} \label{verticelabel} \forall v \in V: \exists label(v) \in \Omega_{V} \end{equation} \subsection{Das chemische Molekül als nicht einzigartig gelabelter Graph} Ein chemisches Molekül bestehend aus kovalenten Bindungen kann als ein ungerichteter Graph $G = (V_{G}, E_{G})$ betrachtet werden. Dabei können die Kanten nach den Bindungstypen: einfach, doppelt, dreifach, vierfach und aromatisch gelabelt werden. Die Knoten können nach dem Elementtyp und mit der Elementarladung gelabelt werden \cite{andersenphd}. Die Label für Knoten und Kanten in einem Molekül müssen nicht einzigartig sein. Es ist möglich das Kantenlabel beispielsweise einer Doppelbindung durch zwei Einfachbindungen zwischen den selben Knoten zu ersetzen. Diese Darstellung kann aber keine aromatischen Bindungen repräsentieren. Etwaige sterische Informationen sind dabei mit einem Graphen nicht darstellbar. Eine chemische Reaktion kann dabei als eine Graphenoperation beschrieben werden, welche in Kapitel \ref{chdoublepushout} näher beschrieben wird. Um zu überprüfen, ob die nötigen funktionellen Gruppen für eine Reaktion vorhanden sind, muss ein Monomorphismus gefunden werden. Um zu bestimmen, ob ein erzeugtes Molekül bereits gefunden wurde, muss ein Isomorphismus gefunden werden. Beide werden in Kapitel \ref{chmophisms} definiert. Die Standard Darstellung von chemischen Molekülen als Graphen folgt dabei der Lewis-Repräsentation \cite{lewisrep, ocwiley}. \begin{itemize} \item Oktettregel: In organischen Molekülen hat Wasserstoff 2 und alle anderen Hauptgruppenelemente 8 Elektronen. \item In kovalenten Bindungen werden 2 mal die Bindungsstärke Elektronen verwendet. Diese zählen für beide an der Bindung beteiligten Atomen. \item Einfachbindungen werden als Kanten zwischen den benachbarten Atomen gezeichnet. Eine Einfachbindung repräsentiert 2 Elektronen. \item Ein Atom kann maximal einen Bindungsgrad von der Anzahl seiner Valenzelektronen haben. Zum Beispiel hat Wasserstoff maximale eine Einfachbindung und Chlor bis zu 7. \item Aromatische Bindungen zählen als 1.5 Elektronen, treten aber meistens in gerader Anzahl auf. \end{itemize} \subsection{Morphismen} \label{chmophisms} Die Definition der Morphismen wurde aus der Doktorarbeit von Jakob Andersen entnommen \cite{andersenphd}. Es seien $G = (V_{G}, E_{G})$ und $H = (V_{H}, E_{H})$, zwei einfache, ungerichtete Graphen. \subsubsection{Morphismus} Ein Graphmorphismus $m: G \rightarrow H$ ist dann eine Abbildung von $G$ nach $H$. In dieser Abbildung müssen alle Knoten und Kanten in $G$ auch eine Abbildung in $H$ haben. Also gilt: Wenn $e = (u,v) \in E_{G} \implies \exists m(e) = (m(u),m(v)) \in E_{H}$. Die Abbildung muss nicht injektiv sein. Ein Beispiel für einen Morphismus, der kein Monomorphismus ist, ist wenn der Graph $G$ aus zwei verbundenen Knoten und der Graph $H$ aus einem Knoten mit einer Schleife besteht (Abbildung \ref{morphisms}). Die Kante kann zwischen beiden Graphen injektiv abgebildet werden, aber beide Knoten aus Graph $G$ werden auf den selben Knoten in Graph $H$ abgebildet. \subsubsection{Monomorphismus} Der Monomorphismus ist ein Morphismus, bei dem alle Knoten in $H$ nur von einem Knoten in $G$ abgebildet werden. Die Abbildung muss Injektiv sein. Es gilt: $\forall u,v \in V_{G}, u \neq v \implies m(u) \neq m(v)$. Ein Beispiel für einen Monomorphismus, der aber kein Subgraph Isomorphismus ist, ist wenn der Graph $G$ aus drei in Reihe verbundenen Knoten besteht und der Graph $H$ aus drei Knoten besteht, die alle untereinander verbunden sind. Dadurch würde zwischen Knoten, die Teil der Abbildung sind eine Kante bestehen, die wiederum nicht abgebildet wird. \subsubsection{Subgraph Isomorphismus} Der Subgraph Isomorphismus ist ein Monomorphismus, bei dem zwischen allen abgebildeten Knoten keine nicht abgebildeten Kanten bestehen. Es gilt: $(u,v) \in E_{G} \iff (m(u), m(v)) \in E_{H}$. Ein Beispiel für einen Subgraph Isomorphismus, der kein Isomorphismus ist, ist wenn ein Graph $G$ aus drei linear verbundenen Knoten besteht und Graph $H$ aus einem Graphen mit vier Knoten besteht,wobei drei linear verbunden sind und der letzte beliebig mit den drei anderen verbunden ist. \subsubsection{Isomorphismus} Der Isomorphismus ist ein Subgraph Isomorphismus, bei dem alle Knoten und Kanten jeweils genau einen Knoten und eine Kante haben auf den sie abgebildet werden können. Die Abbildung ist also bijektiv. \begin{figure} [h!] \includegraphics[width=\textwidth]{morphismen.png} \caption{Beispiele der Morphismen dargestellt als rote Pfeile. (a) Ein Morphismus, der kein Monomorphismus ist, (b) ein Monomorphismus, der kein Subgraph Isomorphismus ist, (c) ein Subgraph Isomorphismus, der kein Isomorphismus ist und (d) ein Isomorphismus \cite{andersenphd}. } \label{morphisms} \end{figure} \newline Um Morphismen in einem chemisch-gelabelten Graphen zu finden, müssen sie zu der richtigen Kantenzahl und Konformation, auch das richtige Label aufweisen. Wir definieren einen chemisch-gelabelten Graphen $G = (V_{G}, E_{G}, l_{G}^{V}, l_{G}^{E})$ mit $l_{G}^{V}: V_{G} \rightarrow \Omega_{V}$ als eine Abbildung $l_{G}^{V}$ der Knoten $V_{G}$ auf die Element Label $\Omega_{V}$ und $l_{G}^{E}: E_{G} \rightarrow \Omega_{E}$ als eine Abbildung $l_{G}^{E}$ der Kanten $E_{G}$ auf die Bindungslabel $\Omega_{V} = \{-, =, :, \# \}$ als Einfach, Doppel, Aromatische und Dreifachbindung respektive. Die Abbildungen sind dabei surjektiv, da jeder Knoten und jede Kante gelabelt sind. Bei einem Morphismus $m: G \rightarrow H$ kommt die Label-Assoziation hinzu, welche einen Knoten und eine Kante jeweils einem Gegenstück in der Abbildung zuordnet (Gleichung \ref{verticelabelmorphismen}, \ref{edgelabelmorphismen}) \cite{andersenphd}. \begin{equation} \label{verticelabelmorphismen} A_{V} (m) = {\langle l_{G}^{V}(v) , l_{H}^{V}(m(v))\rangle | v \in V_{G}} \end{equation} \begin{equation} \label{edgelabelmorphismen} A_{E} (m) = {\langle l_{G}^{E}(e) , l_{H}^{E}(m(e))\rangle | e \in E_{G}} \end{equation} Da die Label als einfache Strings definiert sind gilt: $s_{1} = s_{2} , \forall(s_{1} , s_{2} ) \in A_{V} (m) \cup A_{E} (m)$ \cite{andersenphd}. Für jedes assoziierte Paar in der Knoten und Kantenmenge müssen beide Label gleich sein. \subsection{Double-Push-Out (DPO)} \label{chdoublepushout} Bisher wurde der chemische Graph statisch betrachtet. Um chemische Reaktionen zu definieren, müssen Regeln für das Umschreiben des Graphen definieren werden. Die verwendete Methode ist der Double-Push-Out-Ansatz. Zuerst wird für die Graphen eine Kategorie definiert. Die Kategorie \textbf{C} besteht aus: \begin{itemize} \item einer Klasse von Objekten $Obj(\textbf{C})$ \item einer Klasse von Morphismen $Mor(\textbf{C})$, der zwei Elemente $A,B \in Obj(\textbf{C})$ aufeinander Abbildet $m: A \rightarrow B$ \item einem Morphismus-Kompositions-Operator $\circ : Mor(\textbf{C}) \times Mor(\textbf{C}) \rightarrow Mor(\textbf{C})$ für den gilt: \begin{itemize} \item die Komposition ist abgeschlossen also: $\forall m,n \in Mor(\textbf{C})$ mit \newline $m: A \rightarrow B, n: B \rightarrow C, \exists n \circ m \in Mor(\textbf{C})$ mit $n \circ m: A \rightarrow C$ \item es existiert ein Identitätsmorphismus $\text{id}_{A} \in Mor(\textbf{C}$ mit $\text{id}_{A}: A \rightarrow A$. \item die Komposition ist assoziativ also: $\forall m,n,g \in Mor(\textbf{C})$ mit \newline $m: A \rightarrow B, n: B \rightarrow C$, $g: C \rightarrow D$ gilt $(n \circ m) \circ g = n \circ (m \circ g)$ \end{itemize} \item zwei Morphismen können den selben Start und Endgraphen haben \end{itemize} Auf dieser Kategorie von Graphen lassen sich drei unterliegende Graphentypen definieren. \subsubsection{Pushout} Bei zwei gegebenen Morphismen $g: A \rightarrow C$ und $f: A \rightarrow B$ erzeugen die Morphismen $g': B \rightarrow D$ und $f': C \rightarrow D$ genau dann ein Pushout D von f und g (Abbildung \ref{pushout}), wenn: \begin{itemize} \item die Morphismen kommutativ sind: $g' \circ f = f' \circ g$ \item für alle Morphismen Paare $f'': C \rightarrow D''$ und $g'': B \rightarrow D''$ mit $f'' \circ g = g'' \circ f$ genau ein einzigartiger Morphismus $d'': D \rightarrow D''$ existiert, sodass $f'' = d'' \circ f'$ und $g'' = d'' \circ g'$. \end{itemize} \begin{figure} [h] \includegraphics[width=\textwidth]{pushout.png} \caption{Beispiel für ein kommutatives Diagramm mit einem Pushout D. Die Knoten $A,B,C,D,D'' \in Obj(\textbf{C})$ stellen selber Graphen dar. Die Kanten $g,g',g'',f,f',f'',d'' \in Mor(\textbf{C})$ stellen Morphismen dieser Graphen dar. \cite{andersenphd}. } \label{pushout} \end{figure} \subsubsection{Pullback} Bei zwei gegebenen Morphismen $g': B \rightarrow D$ und $f': C \rightarrow D$ erzeugen die Morphismen $g: A \rightarrow C$ und $f: A \rightarrow B$ genau dann ein Pullback A von f' und g' (Abbildung \ref{pullback}), wenn: \begin{itemize} \item die Morphismen kommutativ sind: $g' \circ f = f' \circ g$ \item für alle Morphismen Paare $f'': A'' \rightarrow B$ und $g'': A'' \rightarrow C$ mit $f' \circ g'' = g'' \circ f''$ genau ein einzigartiger Morphismus $a'': A'' \rightarrow A$ existiert, sodass $f'' = f \circ a''$ und $g'' = g \circ a''$. \end{itemize} \begin{figure} [h] \includegraphics[width=\textwidth]{pullback.png} \caption{Beispiel für ein kommutatives Diagramm mit einem Pullback A. Die Knoten $A,B,C,D,A'' \in Obj(\textbf{C})$ stellen selber Graphen dar. Die Kanten $g, g', g'',f,f',f'', a'' \in Mor(\textbf{C})$ stellen Morphismen dieser Graphen dar. \cite{andersenphd}. } \label{pullback} \end{figure} \begin{figure} \begin{center} \begin{tikzcd}[column sep=small, row sep=small] \bullet && \bullet & \bullet &&&& \bullet && \bullet & \bullet && \\ \\ \bullet && \bullet & \bullet &&&& \bullet && \bullet \\ \bullet && \bullet &&&&& \bullet &&& \bullet \\ &&&&&&&&&&& \bullet & \bullet \\ &&&&&&&&&&& \bullet \\ &&&&&& \bullet \\ \bullet && \bullet & \bullet &&&& \textcolor{white}{\bullet} && \bullet & \bullet \\ \\ \bullet && \bullet & \bullet &&&& \bullet && \bullet & \bullet \\ \bullet && \bullet & \bullet && \bullet && \bullet && \bullet \\ &&&& \bullet & \bullet \\ &&& \bullet & \bullet \arrow[color={rgb,255:red,214;green,92;blue,92}, dashed, from=1-1, to=1-3] \arrow[color={rgb,255:red,214;green,92;blue,92}, dashed, from=1-1, to=3-1] \arrow[no head, from=1-3, to=1-4] \arrow[color={rgb,255:red,214;green,92;blue,92}, dashed, from=1-3, to=3-3] \arrow[color={rgb,255:red,214;green,92;blue,92}, dashed, from=1-4, to=3-4] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=1-8, to=1-10] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=1-8, to=3-8] \arrow[no head, from=1-10, to=1-11] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=1-10, to=3-10] \arrow[color={rgb,255:red,214;green,92;blue,92}, dashed, from=1-11, to=4-11] \arrow[color={rgb,255:red,214;green,92;blue,92}, dashed, from=3-1, to=3-3] \arrow[no head, from=3-3, to=3-4] \arrow[no head, from=3-3, to=4-3] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=3-8, to=3-10] \arrow[no head, from=3-10, to=4-11] \arrow[color={rgb,255:red,214;green,92;blue,92}, curve={height=-6pt}, dashed, from=3-10, to=5-12] \arrow[color={rgb,255:red,214;green,92;blue,92}, curve={height=-12pt}, dashed, from=3-10, to=5-13] \arrow[color={rgb,255:red,214;green,92;blue,92}, curve={height=12pt}, dashed, from=3-10, to=6-12] \arrow[no head, from=4-1, to=3-1] \arrow[color={rgb,255:red,214;green,92;blue,92}, dashed, from=4-1, to=4-3] \arrow[no head, from=4-8, to=3-8] \arrow[color={rgb,255:red,214;green,92;blue,92}, dashed, from=4-8, to=4-11] \arrow[no head, from=5-12, to=5-13] \arrow[no head, from=5-12, to=6-12] \arrow[color={rgb,255:red,214;green,92;blue,92}, curve={height=-6pt}, dashed, from=7-7, to=8-10] \arrow[color={rgb,255:red,214;green,92;blue,92}, curve={height=6pt}, dashed, from=7-7, to=10-8] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=8-1, to=8-3] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=8-1, to=10-1] \arrow[no head, from=8-3, to=8-4] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=8-3, to=10-3] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=8-4, to=10-4] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=8-8, to=8-10] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=8-8, to=10-8] \arrow[no head, from=8-10, to=8-11] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=8-10, to=10-10] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=8-11, to=10-11] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=10-1, to=10-3] \arrow[no head, from=10-3, to=10-4] \arrow[no head, from=10-3, to=11-3] \arrow[no head, from=10-3, to=11-4] \arrow[color={rgb,255:red,214;green,92;blue,92}, curve={height=-6pt}, dashed, from=10-3, to=12-5] \arrow[color={rgb,255:red,214;green,92;blue,92}, dashed, from=10-4, to=11-6] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=10-8, to=10-10] \arrow[no head, from=10-10, to=10-11] \arrow[no head, from=10-10, to=11-10] \arrow[no head, from=11-1, to=10-1] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=11-1, to=11-3] \arrow[color={rgb,255:red,214;green,92;blue,92}, dashed, from=11-3, to=13-4] \arrow[color={rgb,255:red,92;green,214;blue,92}, curve={height=-6pt}, dashed, from=11-4, to=12-6] \arrow[color={rgb,255:red,92;green,92;blue,214}, curve={height=6pt}, dashed, from=11-4, to=13-5] \arrow[no head, from=11-8, to=10-8] \arrow[draw={rgb,255:red,214;green,92;blue,92}, dashed, from=11-8, to=11-10] \arrow[no head, from=12-5, to=11-6] \arrow[no head, from=12-5, to=12-6] \arrow[no head, from=12-5, to=13-5] \arrow[no head, from=13-4, to=12-5] \end{tikzcd} \end{center} \caption{Beispiel für ein tatsächliches Pushout-Pullback-Konstrukt (oben-links), ein zu kleines Pushout (oben-rechts), ein zu großes Pushout (unten-links) und ein zu kleines Pullback (unten-rechts).} \label{pushpullex} \end{figure} Im Folgenden werden zum Pushout und Pullback vier Beispiele näher untersucht (Abbildung \ref{pushpullex}). Das oben-linke Beispiel stellt ein korrektes Pushout-Pullback-Konstrukt dar. Bei jedem Morphismus wird genau eine Kante und ein zu dieser Kante gehöriger Knoten eingefügt. Im oben-rechten Beispiel wird für der Pushout nur insgesamt eine Kante hinzugefügt. Die beiden zuvor erzeugten Kanten werden also kombiniert, weshalb zwei unterschiedliche Morphismusmen $d'':D \rightarrow D''$ mit D'' als tatsächliches Pushout existieren würden. Bei dem unten-linken Beispiel wird für den Pushout-Kandidat eine Kante zu viel eingefügt, wodurch der Pushout-Kandidat ein Dreistern-Graphen ist. Dadurch existieren zwei verschiedene Morphismen (grün und blau) zu einem weiteren Vierstern-Graphen als Pushout-Kandidat. Ebenso existiert kein Morphismus von dem Dreistern-Graphen zu dem tatsächlichen Pushout. Der Pushout ist also zu groß. In dem Beispiel unten.rechts wird für den Pullback-Kandidat angenommen, dass wenn die Kante entfernt wird, auch die zugehörigen Knoten entfernt werden. Dadurch ist kein Morphismus von dem tatsächlichen Pullback zu diesem Pullback-Kandidat möglich. Ebenso würden, wenn jeweils zwei Knoten mit einer Kante mit $g$ und $f$ hinzugefügt würden, diese im Pushout-Kandidat kombiniert werden. Beides widerspricht den etablierten Regeln. \subsubsection{Pushout-Komplement} Bei zwei gegebenen Morphismen $g': B \rightarrow D$ und $f: A \rightarrow B$ erzeugen die Morphismen $g: A \rightarrow C$ und $f': C \rightarrow D$ genau dann ein Pushout-Komplement C von f und g', wenn f' und g' der Pushout von f und g ist oder f und g der Pullback von f' und g' ist. \FloatBarrier \subsubsection{Double-Push-Out} Eine Graphtransformation lässt sich mit dem Double-Push-Out (DPO) Ansatz realisieren. Der DPO ist dabei eine Kombination aus den Morphismen l, r, l', r', m, m' (Abbildung \ref{dpo}) und einer Transformationsregel $p = (L \stackrel{l}{\leftarrow} K \stackrel{r}{\rightarrow} R)$. Die Regel besteht aus dem Kontext $K$, den zwei Morphismen $l$ und $r$ und den Bedingungen $L$ und $R$. Bei einer angenommen Transformation von links nach rechts stellt $L$ die Vorbedingung, $R$ die Folgebedingung, $G$ den Vorzustand und $H$ den Folgezustand dar. Die Graphransformation eines Graphen $G$ zu einem Graphen $H$ mit der Regel $p$ und dem Morphismus $m$ wird über den Ausdruck $G \stackrel{p,m}{\implies} H$ dargestellt. Dafür sind die folgenden Schritte nötig. \begin{itemize} \item Finde einen Morphismus $m: L \rightarrow G$ von der Vorbedingung $L$ zum Startgraphen $G$. \item Bestimme aus den Morphismen $m: L \rightarrow G$ und $l: K \rightarrow L$ das Pushout-Komplement $D$ mit den Morphismen $l': D \rightarrow G$ und $d: K \rightarrow D$. \item Bestimme aus den Morphismen $r: K \rightarrow R$ und $d: K \rightarrow D$ den Pushout $H$ mit den Morphismen $r': D \rightarrow H$ und $m': R \rightarrow H$. \end{itemize} \begin{figure} [h] \includegraphics[width=\textwidth]{dpo.png} \caption{Ein Beispiel für den Double-Push-Out. Die Knoten K, L, R, D, G und H stellen Graphen dar. Die Kanten l, r, l', r', m, m' und d stellen Morphismen auf diesen Graphen dar. Bei einer Transformation von $L$ zu $R$ stellt $L$ die Vorbedingung, $R$ die Folgebedingung, $G$ den Vorzustand und $H$ den Folgezustand dar. Die Knoten K und D werden als Kontext bezeichnet \cite{andersenphd}. } \label{dpo} \end{figure} Die rechte Seite $R$ und die linke Seite $L$ sind dabei austauschbar. Wenn die Graphentransformation $G \stackrel{p,m}{\implies} H$ in eine Richtung abzubilden ist, existiert auch eine inverse Graphtransformation $H \stackrel{p^{-1},m'}{\implies} G$, welche die inverse Regel $p^{-1} = (R \stackrel{r}{\leftarrow} K \stackrel{l}{\rightarrow} L)$ verwendet. Die Existenz des Morphismus $m'$ ist schon für die ursprüngliche Graphtransformation gegeben. Im chemischen Kontext lassen sich die meisten Reaktionen als Gleichgewichtsreaktionen definieren, bei welchen sowohl Hin- als auch Rückreaktion simultan stattfinden \cite{equilibrium}. Für einen gegebenen Morphismus $m: L \rightarrow G$ kann also $D$ bestimmt werden, wenn: \begin{itemize} \item es keine Kanten in $E_{G} \setminus m(E_{L})$ die an einen Knoten in $m(V_{L} \setminus l(V_{K}))$ gebunden sind, also es keine Kanten in $E_{H} \setminus m(E_{R})$ gibt, die an einen Knoten in $m(V_{R} \setminus r(V_{K}))$ gebunden sind. Also gilt für alle Knoten die nicht in einer der beiden Regelseiten, aber nicht im Kontext sind, also nach der Transformation nicht mehr vorhanden sind, nicht mit Kanten verbunden sein können, die Teil der Vorbedingung waren. Wenn $p$ also einen Knoten entfernt, müssen auch von ihm ausgehenden Kanten entfernt werden. \item es keine distinkten Knoten $u,v \in V_{L}$ gibt, mit $m(u) = m(v)$ und $u \in l(V_{k}), v \notin l(V_{k})$ und es keine distinkten Kanten $e,f \in E_{L}$ gibt, mit $m(e) = m(f)$ und $e \in l(E_{k}), f \notin l(E_{k})$. Wenn also eine Kante oder ein Knoten in der Vorbedingung ist, aber durch die Transformation entfernt wird, so kann dieser nicht durch den Morphismus $m$ dem gleichen Knoten zugeordnet werden, wie ein Knoten oder eine Kante die nicht entfernt wurde. \end{itemize} Im Folgenden werden die Morphismen $l$ und $r$ als Monomorphismen festgelegt, da so nur eine Möglichkeit für $D$ besteht. Dafür würde es reichen, dass einer der beiden Morphismen ein Monomorphismus ist. Die Beschränkung beider Morphismen ist nötig, damit die Invertierbarkeit gewahrt bleibt. Ebenso werden die Morphismen $m$ und $m'$ als Monomorphismen festgelegt, damit keine Knoten fusioniert werden können \cite{andersenphd}. Für den DPO-Ansatz bei einer Transformation an einem chemisch-gelabelten Graphen wird eine Menge von Labeln $\Omega$ angenommen. Die Graphen $L$ und $R$ sind dabei mit $\Omega$ gelabelt, während der Kontext $K$ beide Label vereint und somit mit $\Omega \times \Omega$ gelabelt ist. Es wäre möglich das Umbenennen von Labeln zu verhindern wie in Abbildung \ref{labeldpo} Beispiel (a) zusehen, aber diese Definition erlaubt das Umbenennen nach Beispiel (b) für eine effizientere Darstellung. \begin{figure} [h] \includegraphics[width=\textwidth]{labeldpo.png} \caption{Wenn in einer Regel $r$ ein Kantenlabel verändert wird, kann die Kante im Kontext nicht auftreten (a) oder als Kombination beider Kanten gelabelt werden (b). In Beispiel (a) wird die Kante mit Label a in der Transformation entfernt und eine neue Kante mit Label b hinzugefügt. In Beispiel (b) wird das Label einer Kante a zum Label b geändert. In Beispiel (a) wäre es möglich $K$ über den Labelraum $\Omega$ zu definieren. In Beispiel (b) wird der Labelraum $\Omega \times \Omega$ nötig \cite{andersenphd}. } \label{labeldpo} \end{figure} Wird dies auf den chemisch-gelabelten Graphen angewendet, so stellen die Knotenlabel weiterhin die chemischen Elemente mit etwaigen Ladungen dar, jedoch ist für die Kantenlabel neben den chemischen Bindungen noch ein \textbf{nil} Label nötig, um das Entfernen einer Bindung zu beschreiben. Bei einer Graphtransformation müssen weiterhin die Regeln für einen chemisch-gelabelten Graphen gewahrt werden. \subsection{Gerichteter Multi-Hypergraph}\label{hypergraph} Ein gerichteter Multi-Hypergraph $\mathcal{H} = (V,E)$ ist ein Tupel. Dabei ist $V$ eine Menge von Knoten und $E$ eine Menge von gerichteten Hyperkanten. Eine Hyperkante $e \in E$ ist selber ein Tupel $ e = (U, W)$ aus einer Multimenge von Startknoten (tail) $U$ und einer Multimenge von Zielknoten (head) $W$. Beide Multimengen enthalten ausschließlich Knoten und sind daher Teilmengen der Menge von Knoten $U,W \subseteq V$ \cite{andersenhyperflow}. Da es sich bei den Startkonten und Zielknoten um Multimengen handelt kann eine Kante den gleichen Knoten mehrmals in head und tail haben (Abbildung \ref{hypergraphex}). \begin{figure} [h] \includegraphics[width=\textwidth]{hypergraphex.png} \caption{Drei Beispiele für Hyperkanten in Hypergraphen. Die linke Hyperkante beschreibt den Übergang von den Graphen $g_{1}$ und $g_{2}$ zu $g_{3}$. Die mittlere Kante beschreibt den Übergang von dem Graph $g_{4}$ zu $g_{5}$. Bei einem einfachen Übergang wird eine einfache Kante verwendet. Die rechte Hyperkante beschreibt den Übergang von dem Graphen $g_{6}$ zu $g_{7}$. Es werden zwei Exemplare von $g_{6}$ benötigt, was durch den doppelten Pfeil dargestellt wird \cite{andersenphd}.} \label{hypergraphex} \end{figure} Ein $s$,$t$-Superpfad in einem Hypergraphen $\mathcal{H} = (V,E)$ für die unterscheidbaren Knoten $s,t \in V$ ist definiert als eine Kantenuntermenge $F \subseteq E$, sodass für die Kanten geordnet als $e_{1}$ bis $e_{k}$ gilt \cite{hypergraphpath}: \begin{itemize} \item $tail(e_{1}) = {s}$ \item $\forall 1 < i \leq k: tail(e_{i}) \subseteq {s} \bigcup\limits_{1 \leq j < i} head(e_{j})$ \item $head(e_{k}) = {t}$ \end{itemize} Ein Hyperpfad $P$ ist ein minimaler Superpfad. Also gilt: $\forall e \in P: P \setminus {e}$ ist kein Superpfad. \subsection{Zyklen} Für einen endlichen ungerichteten Graphen $G =(V, E)$ ist ein einfacher Pfad definiert als endliche Sequenz $V_{P} = \langle v_{1}, ... , v_{k} \rangle$ mit $v_{1}, ... , v_{k} \in V$ und $E_{P} = \{(v_{1}, v_{2}), ..., (v_{k-1}, v_{k})\}$ mit $(v_{i}, v_{i+1}) \in E$ mit $i \in {1, ..., k-1}$. Ein solcher Pfad ist ein Zyklus, wenn gilt $(v_{k}, v_{1}) \in E$. Eine Sehne in einem Pfad ist eine Kante $e_{c} \in E$ zwischen zwei Knoten des Pfades $e_{i}, e_{j} \in V_{P}$, welche wiederum nicht Teil des Pfades ist $(e_{i}, e_{j}) \notin E_{P}$. Ein sehnenloser Zyklus ist ein sehnenloser Pfad, bei dem der erste und letzte Knoten über eine Kante verbunden sind \cite{chordlesscycle}. \subsection{Integer-Linear-Programming (ILP)} Ein Lineares-Programmier (LP) Problem ist ein mathematisches Problem, bei dem für eine Menge an Variablen $\{x_{1},..., x_{n}\}$ Werte gefunden werden, so dass eine Zielfunktion maximiert oder minimiert wird. Ein LP kann mathematisch als eine Zielfunktion und eine Menge von Bedingungen beschrieben werden. \begin{itemize} \item $max \: z = \sum_{j} c_{j} \times x_{j}$ \item Bedingung 1: $\sum_{j} a_{i,j} \times x_{j} \geq b_{i}$ \item Bedingung 2: $0 \leq x_{j}$ \end{itemize} Dabei sind alle Werte außer einer Variable $x_{j}$ Konstanten. Bei einem Integer-LP Problem ist mindestens eine der Komponenten, meistens die Variable auf Integer also ganzzahlige Werte beschränkt. Es ist möglich, ein Maximierungsproblem in ein Minimierungsproblem zu überführen, indem die Zielfunktion mit einem negativen Wert multipliziert wird. Ebenso können mehrere Zielfunktionen kombiniert werden, indem sie addiert oder subtrahiert werden. Die Kombination aus Multiplikation mit einem großen Faktor und Addition erlaubt es eine primäre Zielfunktion mit einer sekundären oder tertiären Zielfunktion zu kombinieren. Um das ILP zu lösen, werden alle Variablen und Konstanten auf die gleiche Seite der Gleichung gebracht und Ungleichungen durch positive Konstanten $s$ ausgedrückt. \begin{itemize} \item $max \: z - \sum_{j} c_{j} \times x_{j}$ = 0 \item Bedingung 1: $\sum_{j} a_{i,j} \times x_{j} + s_{i} = b_{i}$ \item Bedingung 2: $0 \leq x_{j}$ \item Bedingung 3: $0 \leq s_{i}$ \end{itemize} Diese Menge an Gleichungen kann in eine Matrix $\textbf{M}$ überführt werden, bei der die Zeilen die Gleichungen und die Spalten die Variablen $x_{j}, s{i}$ darstellen (Gleichung \ref{matrixilp}). \begin{equation} \label{matrixilp} \begin{bmatrix} z & -c_{1}& ... & -c_{j} & 0 & ... & 0 &| 0\\ 0 & -a_{1,1} & ... & -a_{1,j} & s_{1} & ... & 0 &| b_{1}\\ 0 & ...& ... & ... & ... & ... & 0 &| ...\\ 0 & -a_{i,1}& ... & -a_{i,j} & 0 & ... & s_{i} &| b_{i} \end{bmatrix} \times \begin{bmatrix} x_{1}\\ ... \\ x_{j}\\ s_{1}\\ ... \\ s_{i} \end{bmatrix} \end{equation} Das lineare Gleichungssystem kann durch das Invertieren der Matrix gelöst werden (Gleichung \ref{matinv}): \begin{equation}\label{matinv} \textbf{A} \times x = b \iff \textbf{A}^{-1} \times b = x \end{equation} Die daraus resultierende Lösung ist nicht unbedingt ganzzahlig. Um eine ganzzahlige, optimale Lösung zu finden, wird ein Baum mit der über die Matrixinversion erhaltenen Lösung als Wurzel aufgespannt. Durch eine Abfolge von Schritten kann eine ganzzahlige Lösung ermittelt werden. \begin{itemize} \item Bestimme \textbf{A}$^{-1}$ aus \textbf{A} und $b$. \item Bestimme $x$ aus \textbf{A}$^{-1}$ und $b$. \item Sind alle Variablen in $x$ nur Integer? \item Wähle eine Variable $x_{i}$ oder $s_{i}$ aus der Menge von $x$ und $s$ aus. \item Teile in zwei neue lineare Gleichungssysteme auf, für welche erneut eine Lösung gesucht wird. \begin{itemize} \item Größer: Füge $x_{i} \geq ceiling(x_{i})$ als neue Bedingung hinzu. \item Kleiner: Füge $x_{i} \leq floor(x_{i})$ als neue Bedingung hinzu. \end{itemize} \item Ein Pfad muss nicht weiter verfolgt werden, wenn: \begin{itemize} \item Unmöglich: Es gibt keine optimale Lösung mit der Bedingung. \item Optimal: Der neue Knoten ist ganzzahlig. \item Begrenzung: Die Parameterbegrenzung liegt außerhalb der Anfangsbedingungen. \end{itemize} \end{itemize} Die Entscheidung, welche Variable als nächstes eingeschränkt wird, geht nach der folgenden Reihenfolge. \begin{itemize} \item Am nächsten an .5 \item Stärkster Faktor in der Zielfunktion \item niedrigster Index \end{itemize} Um zu entscheiden, welcher Knoten als nächstes untersucht wird, können verschiedene Baumsuchalgorithmen verwendet werden \cite{ilp}. \subsection{Petrinetze} Ein Petrinetz ist ein Werkzeug zur Untersuchung der Erfüllbarkeit parallel ablaufender Prozesse \cite{petri1962kommunikation}. Es besteht aus Knoten und Kanten. Die Knoten lassen sich dabei in Zustände $p$ und Übergänge $t$ unterteilen. Die Zustände werden dabei als Kreise und die Übergänge als Querbalken dargestellt (Abbildung \ref{petrinetbase}). \begin{figure} [bh!] \includegraphics[width=\textwidth]{Petrinet1.png} \caption{Ein einfaches Petrinetz bestehend aus den Zuständen $p_{1}$, $p_{2}$, $p_{3}$, $p_{4}$, $p_{5}$, $p_{6}$, $p_{7}$, dargestellt als Kreise und den Übergängen $t_{1}$, $t_{2}$, $t_{3}$, $t_{4}$, $t_{5}$, $t_{6}$, dargestellt als Querbalken \cite{petrinets}.} \label{petrinetbase} \end{figure} \newline In einem Petrinetz kann auf jedem Zustand ein Token platziert werden, um einen Gesamtzustand des Systems zu repräsentieren. Ein Übergang kann ausgelöst werden, wenn alle Zustände, die auf ihn zeigen, mindestens einen Token beinhalten. In dem Beispielgraphen (Abbildung \ref{petrinetmarked}) kann nur der Übergang $t_{2}$ auslösen, da ein Token in $p_{1}$ ist. Bei den Übergängen $t_{3}$ und $t_{5}$ liegt nur ein Token in $p_{5}$ und $p_{7}$ vor. Damit die Übergänge auslösen können, muss ein Token in $p_{3}$ vorliegen. \begin{figure} [h!] \includegraphics[width=\textwidth]{Petrinetmarked.png} \caption{Ein einfaches Petrinetz bestehend aus den Zuständen $p_{1}$, $p_{2}$, $p_{3}$, $p_{4}$, $p_{5}$, $p_{6}$, $p_{7}$, dargestellt als Kreise und den Übergängen $t_{1}$, $t_{2}$, $t_{3}$, $t_{4}$, $t_{5}$, $t_{6}$, dargestellt als Querbalken. Die Zustände $p_{1}$, $p_{5}$ und $p_{7}$ haben ein Token, weshalb nur Übergang $t_{2}$ ausgelöst werden kann \cite{petrinets}.} \label{petrinetmarked} \end{figure} \newline Nachdem Übergang $t_{2}$ ausgelöst hat, liegt ein neuer Token in $p_{2}$ und $p_{3}$ (Abbildung \ref{petrinetmarkedfree}). Der Token in $p_{3}$ ist nötig, damit die Übergänge $t_{3}$ und $t_{5}$ ausgelöst werden können. Ebenso kann ein erneutes sequenzielles auslösen von $t_{1}$ und $t_{2}$ eine unbegrenzte Anzahl weiterer Token in $p_{3}$ erzeugen. Dieser Teil des Petrinetzes ist also ein Erzeuger. Die erzeugten Token in $p_{3}$ können entweder in $t_{3}$ und $t_{4}$ oder $t_{5}$ und $t_{6}$ verbraucht werden. Die Übergangsequenzen sind also im Konflikt zueinander. \begin{figure} [h!] \includegraphics[width=\textwidth]{Petrinetmarkedfree.png} \caption{Ein einfaches Petrinetz bestehend aus den Zuständen $p_{1}$, $p_{2}$, $p_{3}$, $p_{4}$, $p_{5}$, $p_{6}$, $p_{7}$, dargestellt als Kreise und den Übergängen $t_{1}$, $t_{2}$, $t_{3}$, $t_{4}$, $t_{5}$, $t_{6}$, dargestellt als Querbalken. Die Zustände $p_{2}$, $p_{5}$, $p_{7}$ und $p_{3}$ haben ein Token, weshalb die Übergänge $t_{1}$, $t_{3}$ und $t_{5}$ ausgelöst werden können \cite{petrinets}.} \label{petrinetmarkedfree} \end{figure} Die Feststellung, welche sequenziellen Übergänge und Gesamtzustände aus einem initialen Marking oder der initialen Tokenverteilung $\mu$ möglich sind, wird \textit{Reachability} genannt. Ein Zustand $\mu'$ ist direkt erreichbar, wenn er durch das Auslösen eines möglichen Übergangs erzeugt wird. Ein Zustand $\mu$ ist allgemein erreichbar, wenn er durch einen Zustand direkt erreichbar ist, welcher selber allgemein oder direkt erreichbar ist. Eine \textit{Reachability} Menge $R(M)$ kann für ein Petrinetz $M = (P, T, I, O, \mu)$ mit der Zustandsmenge $P$, der Übergangsmenge $T$, der Menge von Input und Output Funktionen für Übergänge $I$ und $O$ und dem initialen Marking $\mu$ definiert werden, als die Menge aller Markings, die durch das markierte Petrinetz erreicht werden können \cite{petrinets}. Mit Blick auf den zuvor definierten Hypergraphen (Kapitel \ref{hypergraph}) können die Zustände in einem Petrinetz als Knoten in einem Hypergraphen und die Übergänge als Hyperkanten verstanden werden. Dann erlaubt es der Einsatz von Petrinetzten die Erfüllbarkeit einer Abfolge von Hyperkantenübergängen zu untersuchen. Insbesondere bei Reaktionsnetzwerken kann so durch ein initiales Marking der Startmoleküle oder eine erzeugende Kante ohne Input gewährleisten, dass nur Reaktionen stattfinden, für die alle nötigen Edukte vorliegen. \FloatBarrier \section{NMR-Spektroskopie}\label{chnmr} Um die im Hypergraphen entwickelten organischen Moleküle zu untersuchen eignet sich NMR-Spektroskopie. Für die Erklärung der Grundlagen der NMR-Spektroskopie wurde das Lehrbuch Organische Chemie von Peter Vollhardt \cite{ocwiley} verwendet. \subsection{Spektren} Die Spektroskopie (lat. Erscheinung, Sichtbarwerden) befasst sich mit der Darstellung realer Eigenschaften durch eine Auftragung einer kontinuierlichen oder kleinschrittig gestaffelten Größe gegen die Intensität. In der organischen Chemie werden vor allem vier Spektroskopie Typen verwendet. Die Massenspektroskopie (MS) hat in der Auftragung der x-Achse kleinschrittig gestaffelte Werte. Bei der MS wird ein Molekül in kleinere, geladene Fragmente geteilt. Diese Fragmente können über ein meist ganzzahliges Masse-Ladungs-Verhältnis $m/z$ und ihre Häufigkeit charakterisiert werden. Zerstörungsfrei hingegen sind die kontinuierlichen, auf elektromagnetischen Wellen basierenden Absorptions- oder Resonanzspektren. Die verbreitetsten Absorptionsspektren sind dabei UV/Vis und IR. Die verwendete elektromagnetische Strahlung liegt dabei im UV/Vis- (200 - 800 nm) und IR- (800 nm - 1 mm) Bereich. Die meist verbreitetste Resonanzspektroskopie hingegen ist die Kern-Magnet-Resonanz-Spektroskopie (Nucleus-Magnetic-Resonance, kurz NMR-Spektroskopie). UV-Resonanzspektroskopie ist möglich, wobei die Menge von Molekülen, die dafür verwendbar sind, sehr gering ist. Die elektromagnetischen Wellen für NMR sind im Radiowellenbereich (über 100 mm) \cite{ocwiley}. \subsection{Kernspin} Die im Atomkern befindlichen Neutronen und Protonen können sich in einem von zwei für NMR relevanten Zuständen befinden. Diese Zustände des Kernspins sind die $\alpha$- und $\beta$-Konfiguration. Dieser Kernspin ist keine visuelle Rotation, hat aber die Eigenschaft, wie Ladung in Bewegung, ein Magnetfeld zu erzeugen. Die verschiedenen Zustände unterscheiden sich dabei in der theoretischen Drehrichtung und damit auch in der Richtung des induzierten Magnetfeldes. Unter normalen Umständen treten beide Konfigurationen im Verhältnis 1:1 auf. Ein entgegen gerichtetes, externes Magnetfeld führt dazu, dass diese Konfiguration energetisch gehindert wird. Dadurch befinden sich deutlich weniger Kerne in diesem Zustand. Die kinetische Energie der Teilchen erlaubt jedoch einen Konfigurationswechsel, wenn der Energieunterschied klein genug ist. Es ist ebenfalls möglich diesen Konfigurationswechsel durch Radiowellen einer zu dem Kern gehörigen, charakteristischen Wellenlänge zu erzwingen, wobei die Rückänderung der Konfiguration, auch Relaxation genannt, nach kurzer Zeit eintritt. Diese Relaxation emittiert dabei eine für den Kern charakteristische Wellenlänge. Sowohl die absorbierte, als auch die emittierte Wellenlänge hängen von der Energiedifferenz der beiden Konfigurationen ab, welche wiederum von der Umgebung des Kerns, aber auch von dem angelegten Magnetfeld abhängen. Die Frequenz der Strahlung ist dabei linear-proportional zur Stärke des Magnetfeldes. \subsection{NMR aktive Kerne} Da für NMR-Spektren der Kernspin nötig ist, kann diese Methode nur für Spektren verwendet werden, die Kernresonanz aufweisen. Das sind alle Kerne, die keine gerade Anzahl von Protonen und gerade Anzahl von Neutronen haben. Für die organische Chemie besonders wichtige Kerne sind Wasserstoff, Kohlenstoff, Sauerstoff und Stickstoff (Tabelle \ref{tabelleisotope}). Für Wasserstoff ist dabei $^{1}$H besonders gut geeignet. Es tritt mit hoher Häufigkeit von 99.985 $\%$ auf und hat mit einer ungefähren Resonanzfrequenz von 300 MHz eine deutlich verschiedene Frequenz zu dem mit 0.015 $\%$ sehr seltenen $^{2}$H (oder D Deuterium), welches bei etwa 45.7 MHz resoniert. Bei Stickstoff bietet sich das Isotop $^{14}$N ebenso gut wie Deuterium an. Für Kohlenstoff ist $^{12}$C mit relativer Häufigkeit von 98.89 $\%$ nicht geeignet, da dieses Isotop 6 Neutronen und 6 Protonen hat und somit nicht NMR aktiv ist. Darum wird das deutlich seltenere $^{13}$C für NMR-Spektren verwendet. Für Sauerstoff ist $^{17}$O die einzige Option, wobei die relative Häufigkeit mit 0.037 $\%$ noch geringer ist. \begin{table}[h] \caption{Tabelle der chemischen Isotope, der in der organischen Chemie relevanten Elemente, mit natürlicher relativer Häufigkeit und ob sie NMR aktiv sind \cite{ocwiley}.} \begin{center} \begin{tabular}{ |c|c|c|c| } \label{tabelleisotope} Isotop N & Protonen/Neutronen & NMR-Aktivität & Natürliche Häufigkeit in $\%$ \\ \hline $^{1}$H & 1/0 & Ja & 99.985 \\ $^{2}$H (D) & 1/1 & Ja & 0.0015 \\ $^{3}$H (T) & 1/2 & Ja & 0 \\ $^{12}$C & 6/6 & Nein & 98.89 \\ $^{13}$C & 6/7 & Ja & 1.11 \\ $^{14}$N & 7/7 & Ja & 99.63 \\ $^{15}$N & 7/8 & Ja & 0.37 \\ $^{16}$O & 8/8 & Nein & 99.759 \\ $^{17}$O & 8/9 & Ja & 0.037 \\ $^{18}$O & 8/10 & Nein & 0.037 \\ \hline \end{tabular} \end{center} \end{table} Die verschiedenen Elemente und Isotope haben meistens unterschiedliche Frequenzbereiche, mit denen sie angeregt werden können. \subsection{NMR-Spektrometer} Um ein NMR-Spektrum aufzunehmen, werden wenige Milligramm einer zu untersuchenden Probe in 0.3 bis 0.5 mL Lösungsmittel aufgelöst. Das Lösungsmittel selbst sollte dabei möglichst wenig der zu untersuchenden Kerne enthalten. Bei Wasserstoff wird zum Beispiel ein deuteriertes Lösungsmittel verwendet. Die Lösung wird in ein elektrisch induziertes Magnetfeld eingebracht (Abbildung \ref{nmrdevice}). Ein Radiofrequenzpulsgenerator sendet einen kurzen Impuls an Radiowellen in die Probe, der das gesamte zu untersuchende Spektrum umfasst. Diese Radiowellen treffen auf die Röhre, wobei ein überwiegender Teil nicht absorbiert wird. Die absorbierten Radiowellen regen die Kernspins an, welche nacheinander wieder relaxieren. Die durch die Relaxation emittierten Radiowellen werden in alle Richtungen abgestrahlt und von einem Radiofrequenzdetektor gemessen. Dieser ist orthogonal zum Radiofrequenzpulsgenerator positioniert, damit möglichst keine der durchgehenden Strahlung direkt auf den Detektor trifft. Der Radiofrequenzdetektor misst eine gedämpfte Schwingung mit verschiedenen Frequenzen. Diese Frequenzen können mit der Fouriertransformation von einem zeitabhängigen Oszillationsspektrum zu einem frequenzabhängigen Peakspektrum umgerechnet werden. \begin{figure} [h] \includegraphics[width=\textwidth]{nmrdevice.png} \caption{Schematischer Aufbau eines NMR-Spektrometers \cite{ocwiley}. } \label{nmrdevice} \end{figure} \subsection{Chemische Verschiebung} Da die Rotation des geladenen Atomkerns ein Magnetfeld induziert liegt nahe, dass die Bewegung von Elektronen ebenfalls zu einer solchen Induktion führt. Das von den Elektronen induzierte Magnetfeld ist nach der Lenz'schen Regel dem induzierenden Magnetfeld entgegen gerichtet. Es lässt sich sagen, dass bei Kernen mit hoher umliegender Elektronendichte ein schwächeres effektives Feld anliegt. Um solche Kerne anzuregen, wird eine geringere Frequenz benötigt. Da diese Frequenz zwar für den Kern charakteristisch ist, aber mit dem angelegten Magnetfeld variiert, wird ein interner Standard verwendet. Für Wasserstoff- und Kohlenstoff-Spektren ist dieser Tetramethylsilan (TMS (CH$_{3}$)$_{4}$Si). Die Frequenz $\nu$ wird daher als Differenz zu der Frequenz dieses internen Standards $\nu_{TMS}$ betrachtet. Da aber die Differenz eines Kerns zum internen Standard proportional zur Feldstärke zunimmt, wird die Differenz durch die Frequenz des Spektrometers $\nu_{0}$ geteilt, um die chemische Verschiebung $\delta$ in ppm zu erhalten (Gleichung \ref{eqppm}). \begin{equation} \label{eqppm} \delta[ppm] = \frac{\nu[Hz] - \nu_{TMS}[Hz]}{\nu_{0}[MHz]} \end{equation} Kerne, bei denen die gleiche Elektronendichte vorliegt, haben ein gleichstarkes effektive Magnetfeld, absorbieren die gleiche Wellenlänge und ihre Peaks liegen somit im entstehenden Spektrum bei der gleichen ppm Zahl. Sie werden darum als chemisch-äquivalent bezeichnet. Um zu bestimmen, ob zwei Kerne in einem Molekül chemisch-äquivalent zueinander sind, müssen beide Kerne durch Symmetrieoperationen, wie Rotationen um Einfachbindungen oder Spiegelungen aufeinander abzubilden sein. Als Beispiel wird Chlorethan (H$_{3}$C-CH$_{2}$Cl) betrachtet (Abbildung \ref{symmetrie}). Die Wasserstoffe in der Methylgruppe (H$_{b}$) sind, projiziert auf eine Ebene, alle 120° auseinander. Wenn die Gruppe also sequenziell 120° um die C-C-Einfachbindung rotiert wird, können alle Wasserstoffe aufeinander abgebildet werden. Wenn die Methylgruppe so rotiert wird, dass ein Wasserstoff parallel zum Chlor steht, kann durch die H-C-Bindung und die C-C-Bindung eine Spiegelebene aufgespannt werden. Diese Spiegelebene kann die beiden chlornahen Wasserstoffe H$_{a}$ aufeinander abbilden. Ein komplizierteres Beispiel ist Cyclohexan (Abbildung \ref{symmetrie}). Die energetisch-günstigste Konfiguration für Cyclohexan ist die Sesselkonfiguration, bei welcher aus einem Sechseck zwei der Ecken in entgegengesetzte Richtungen geklappt wurden. In dieser Form sind alle Kohlenstoffe äquivalent zueinander, da die Winkel und benachbarten Atome gleich sind. Jeder Kohlenstoff hat dabei einen Wasserstoff, welcher planar oder äquatorial zum Sechsring liegt (rot) und einen Wasserstoff, welcher orthogonal oder axial zum Sechsring steht. Ohne Rotation lassen sich diese beiden Wasserstoffe voneinander unterscheiden. Durch Rotation um die zu den Spitzen führenden C-C-Einfachbindungen lässt sich die Richtung der beiden Spitzen umkehren, wodurch alle äquatorialen zu axialen Wasserstoffen werden und umgekehrt. Diese Rotation findet bei Raumtemperatur schnell genug statt, um im NMR-Spektrum nur ein Signal zu erzeugen. \begin{figure} [h] \includegraphics[width=\textwidth]{symmetrie.png} \caption{Symmetrieoperationen am Beispiel Chlorethan (A) und Cyclohexan (B). Bei Chlorethan kann eine Spiegelachse vertikal von Chlor (Cl) zum unteren Wasserstoff (grün) gezogen werden. Ebenso ist eine Rotationsachse zwischen den beiden Kohlenstoffen (Kreismitte) möglich. In Cyclohexan stehen sechs Wasserstoff orthogonal (axial) und sechs waagerecht (äquatorial) zum 6-Ring. Durch ein Ändern der Sesselkonfiguration, können beide ineinander überführt werden \cite{ocwiley}. } \label{symmetrie} \end{figure} Eine Restgruppe an einem der Kohlenstoffe oder eine Temperatur von -90 °C würde es aber erlauben, zwischen äquatorialen und axialen Positionen zu unterscheiden. \subsection{Signalintensität} Nachdem die chemische Verschiebung der Peaks diskutiert wurde, wird nun thematisiert, was an den ppm-Verschiebungswerten zu sehen ist. Das NMR-Spektrum besteht aus einer berauschten Grundlinie, aus der an manchen Stellen hohe Intensitätswerte hervorstechen. Diese Werte werden Peaks genannt. Die ppm-Werte können dabei den einzelnen Kernen zugeordnet werden. Ein Peak kann aber auch in einem Peakcluster angeordnet sein, welche in Kapitel \ref{chmultiplet} thematisiert werden. Die Peaks selber haben in den meisten Fällen eine breite Basis, die mit steigender Intensität immer spitzer zuläuft. Besonders bei sauren Wasserstoffen ist es jedoch möglich, dass der Peak auseinander gezogen wird und nur eine deutlich breitere und niedrigere Gaußglocke zu erkennen ist. In manchen Messungen ist dieser Peak nur als ein Anheben der Grundlinie zu erkennen. Die Peakhöhe und -breite kann dabei durch äußere Umstände variieren. Um für das Signal trotzdem vergleichbare Werte zu bekommen, wird die Fläche unter den Peaks integriert. Da die Flächen nur relativ zueinander konstante Werte ergeben, wird die kleinste Fläche als 1 definiert. Es wird davon ausgegangen, dass dieser Peak genau einem Wasserstoff im zu untersuchenden Molekül entspricht. Dabei gilt auch, dass die Peaks von chemisch-äquivalenten Kernen an der gleichen Verschiebung auftreten und dadurch überlagern. Die Peakfläche entspricht somit der Anzahl an Kernen, die diese Verschiebung aufweisen. Sollte bei dieser Analyse herauskommen, dass größere Peaks eindeutig Brüche darstellen, wird die kleinste Fläche in ihrem Wert erhöht, bis alle Flächen ganzzahlig sind. In dem $^{1}$H-Spektrum vorherigen Beispiel Chlorethan, gibt es zwei Peaks. Der kleinere Peak entspricht den chlornahen Wasserstoffen und der größere der Methylgruppe. Wenn für diesen Peak nun eine relative Fläche von 1 angenommen werden würde, so würde der Peak der Methylgruppe eine Fläche von 1.5 aufweisen. Die Fläche des kleinsten Peaks wird also auf die nächstgrößere, natürliche Zahl erhöht. Dadurch sind beide Peakflächen mit 2 und 3 ganzzahlig, was auch der tatsächlichen Häufigkeit in Chlorethan entspricht. \subsection{Multiplettaufspaltung} \label{chmultiplet} Da das Anlegen eines Magnetfeldes nur eine Kernspinkonfiguration bevorzugt, jedoch nicht alle Kerne in diese überführt ist immer noch ein Großteil der Kernspins zufällig ausgerichtet. Liegt ein angeregter Kern nun neben nicht angeregten Kernen, so beeinflusst deren Magnetfeld den angeregten Kern. Bei der Ausrichtung der nicht angeregten Kerne mit dem Magnetfeld, wird das Signal zu höherer Frequenz und bei entgegengesetzter Ausrichtung zu niedrigere Frequenz verschoben. Beide Signale treten nun im Spektrum auf, wobei sie zu dem selben Kern gehören. Sie werden deshalb als ein Peak zusammengefasst. Chemisch-äquivalente Kerne spalten das Multiplett mit gleicher Stärke auf. Dabei kann binomialverteilt jede Kombination an mit und entgegen gerichteten Kernen auftreten (Tabelle \ref{tablepeakheight}). Liegt in der näheren Umgebung ein weiterer, nicht zu den anderen benachbarten Kernen chemisch-äquivalenter Kern, so kann dieser das Spektrum mit einer anderen Frequenz aufspalten. Da die Frequenz der Aufspaltung verschieden ist, kommt es zu einem kombinierten Multiplett. Zum Beispiel gib es bei 1,1-Dichlorpropan zu einer Aufspaltung der zentralen 2-Wasserstoffe zu einer Dublett in Quartett Aufspaltung. Die Methylgruppe spaltet dabei in ein Quartett auf und der chlornahe Wasserstoff spaltet dieses Quartett wiederum in ein Dublett auf. Je nach Auspaltungsfrequenz können die beiden Multipletts verschieden ineinander verschränkt sein. Es können zwei nebeneinanderliegende Quartetts (Quartett in Dublett) oder jedes Signal des Quartetts aufgespalten (Dublett in Quartett) sein. Bei Wasserstoffkernen folgt die Multiplett Zahl der $N+1$-Regel. Bei anderen Kernen kann es aber auch zu einer $2N+1$-Verteilung kommen. Da die Aufspaltung durch Interaktion der Kerne miteinander entsteht und der Interaktionsabstand der Kerne zueinander gleich ist, ist die Aufspaltungsfrequenz symmetrisch. Es ist also möglich zwei benachbarte Kerne darüber zu identifizieren, dass der Abstand innerhalb des Multipletts gleich ist. Diese Interaktion ist dabei über bis zu 3 Bindungen gut erkennbar. Bei geminalen, also an den gleichen Kohlenstoff gebundenen Wasserstoffen, können Aufspaltungen von 0 - 18 Hz auftreten und bei vicinalen, also an benachbarten Kohlenstoff gebundenen Wasserstoffen, kann diese 6 - 8 Hz betragen. \begin{table}[h] \caption{Relative Intensitäten der einzelnen Peaks in den verschiedenen Multipletts, sowie die dafür nötige Anzahl chemisch-äquivalenter benachbarter Wasserstoffe \cite{ocwiley}.} \begin{center} \begin{tabular}{ |c|c|c| } \label{tablepeakheight} Anzahl chemisch-äquivalenter & & \\ Wasserstoffe N & Multiplett Name & relative Peakhöhe \\ \hline 0 & Singulett (s) & 1 \\ 1 & Dublett (d) & 1 : 1 \\ 2 & Triplett (t) & 1 : 2 : 1 \\ 3 & Quartett (q) & 1 : 3 : 3 : 1 \\ 4 & Quintett (quin) & 1 : 4 : 6 : 4 : 1 \\ 5 & Sextett & 1 : 5 : 10 : 10 : 5 : 1 \\ 6 & Septett & 1 : 6 : 15 : 20 : 15 : 6 : 1 \\ \hline \end{tabular} \end{center} \end{table} \subsection{Spektrumanalyse} Im Folgenden wird die Analyse eines NMR-Spektrums am Beispiel von 1,1,2-Trichlorpropan erläutert (Abbildung \ref{nmrex}). Nach der Rotationssymmetrie lassen sich alle Wasserstoffe der Methylgruppe aufeinander abbilden. Durch die an der 1-C-Position gebundenen Chloratome gibt es aber keine Symmetrie von den Wasserstoffen der Endpositionen aufeinander. Es gibt also 3 verschiedene Wasserstofftypen. Ebenso lässt sich sagen, dass das Signal für die Methylgruppe mit 3 Wasserstoffen eine Gesamtintensität von 3 aufweist. Da die beiden anderen Wasserstofftypen nur jeweils ein Proton umfassen, haben sie jeweils eine Intensität von 1. Die H$_{c}$-Methylwassestoffe haben zu dem an den 2-Kohlenstoff gebundenen H$_{b}$-Wasserstoff einen Abstand von 3 Bindungen. Die beiden Signale spalten einander also jeweils mit der Kopplungskonstante $J_{bc} = 6.8 Hz$ auf. Das H$_{c}$-Signal wird dabei als ein Dublett aufgespalten, da nur ein chemisch-äquivalenter Kern benachbart ist. Der Wasserstoff der 1-C-Position H$_{a}$ und der Wasserstoff der 2-C-Position sind zueinander ebenfalls 3 Bindungen entfernt und spalten einander mit der Kopplungskonstante $J_{ab} = 3.6 Hz$ auf. Da zu dem H$_{a}$ nur ein chemisch-äquivalentes H$_{b}$ benachbart ist, wird das H$_{a}$-Signal zu einem Dublett aufgespalten. Der H$_{b}$ hat sowohl H$_{a}$, als auch drei H$_{c}$ in seiner Nachbarschaft. Darum wird es von den drei H$_{c}$ zu einem Quartett und von dem H$_{a}$ zu einem Dublett aufgespalten. Es entsteht ein Dublett in Quartett. Da Chlor elektronenziehend ist, wird mit zunehmender Chlorsubstitution das Proton immer weiter entschirmt und das Signal erscheint bei höherer Verschiebung. Der H$_{a}$ hat mit zwei benachbarten Chloratomen die größte, H$_{b}$ mit einem benachbarten Chloratom die mittlere und H$_{c}$ mit keinen benachbarten Chloratomen die niedrigste chemische Verschiebung. \begin{figure} [h] \includegraphics[width=\textwidth]{nmrex.png} \caption{$^{1}$H-NMR von 1,1,2-Trichlorpropan. Der Methylwasserstoff (H$_{c}$ rot) hat ein Dublett bei 1.68 ppm mit Intensität 3, der zentrale Wasserstoff (H$_{b}$ grün) hat bei 4.35 ppm ein Dublett in Quartett mit Signalstärke 1 und der Endwasserstoff (H$_{a}$ blau) erzeugt ein Dublett bei 5.96 ppm mit Signalstärke 1 \cite{ocwiley}. } \label{nmrex} \end{figure} \subsection{Schreibweise} In einem wissenschaftlichen Kontext wird nicht immer das gesamte Spektrum als Bild angegeben. Um trotzdem eine Angabe zu machen, die das Spektrum ausreichend beschreibt, muss ein aussagekräftiger String verwendet werden. Für Chlorethan wäre dieser String zum Beispiel $^{1}$H-NMR (300 MHz): $\delta$ 1.488 (-CH$_{3}$, 3, t, $^{3}J = 7.232 Hz$), 1.488 (-CH$_{2}$-, 2, q, $^{3}J = 7.232 Hz$). Für jeden Kern wird dabei die Verschiebung in ppm, eine kurze Zuordnung zum Kern, die Intensität, die Multiplettaufspaltung und die Stärke der Aufspaltung genannt. Für die Stärke der Aufspaltung wird zudem angegeben, über wie viele Bindungen diese Interaktion stattfindet. Für das ganze Spektrum wird auch angegeben, auf welche Frequenz oder Feldstärke das Spektrometer eingestellt wurde. \subsection{Kohlenstoff NMR}\label{chdft} Der neben Wasserstoff am häufigsten untersuchte Kern ist Kohlenstoff. Für sein NMR wird $^{13}$C verwendet und ein Lösungsmittel mit $^{12}$C eingesetzt. Das Spektrum verhält sich dabei größtenteils genauso wie das von Wasserstoff, mit der Ausnahme dass $^{13}$C fast nie nebeneinander vorliegt, also keine $^{13}$C-$^{13}$C-Multiplett-Interaktion stattfindet. $^{13}$C kann aber durchaus mit $^{1}$H interagieren. Diese Kopplung wird häufig gezielt neutralisiert, da der $^{13}$C Kern weniger leicht angeregt wird und so ein schlechteres Signal/Rausch-Verhältnis vorliegt. Kleinere Peaks können so leichter im Rauschen untergehen. Darum wird oft ein $^{1}$H Breitband entkoppeltes $^{13}$C Spektrum aufgenommen. Die Notation dafür ist $^{13}$C\{$^{1}$H\}. Um dies zu erreichen, wird ein breites Wasserstofffrequenzsignal während der gesamten Messung eingestrahlt. Das resultierende Spektrum weist für jeden chemisch-äquivalenten Kohlenstoff nur einen Peak auf. Die Peak Fläche wird dabei je nach ausgeblendeten Kern unterschiedlich beeinflusst, weshalb die Intensität nur noch als Richtwert und nicht mehr als absolute Größe betrachtet werden kann. \section{Dichte Funktional Theorie (DFT)} Die Dichte Funktional Theorie (DFT) vereinfach die quantenchemischen Atomorbitale als eine über den Raum verteilte Elektronendichtefunktion $\rho[r]$. Für ein gegebens Molekül lässt sich daraus die kinetische Energie $T$ (Gleichung \ref{dftkin}), die Energie der Elektronen-Protonen-Wechselwirkung $V_{eN}$ (Gleichung \ref{dftpotn}), sowie die Energie der Elektronen-Elektronen-Wechselwirkung $V_{ee}$ (Gleichung \ref{dftpote}) als Integral dieser Dichtefunktion über den gesamten Raum $r$ bestimmen. \begin{equation} \label{dftkin} T[\rho] = C \int \rho(r)^{\frac{5}{5}} dr \end{equation} \begin{equation} \label{dftpotn} V_{eN}[\rho] = -j_{0}Z_{0} \int \frac{\rho(r)}{r} dr \end{equation} \begin{equation} \label{dftpote} V_{ee}[\rho] = \frac{1}{2}j_{0} \int \frac{\rho(r_{1})\rho(r_{1})}{|r_{1}-r_{2}|} dr_{1}dr_{2} \end{equation} Die Gesamtenergie des Moleküls lässt sich als Summe aller Teilenergien ausdrücken (Gleichung \ref{dftEtot}). \begin{equation} \label{dftEtot} E_{TF}[\rho] = T[\rho] + V_{eN}[\rho] + V_{ee}[\rho] \end{equation} Nach dem Variationsprinzip ist die tatsächliche Elektronenkonfiguration genau dann erreicht, wenn die Gesamtenergie für die feste Elektronenzahl und Kernzusammensetzung den niedrigsten Wert erreicht. Während die Elektronendichte variiert wird, ist auch die Position der Kerne festgelegt. Die beste gefundene Konfiguration ist dabei immer energetisch ungünstiger, als die tatsächliche Konfiguration, wobei der gefundene Minimalwert das System am besten beschreibt \cite{atkinsdft}. Um die Bestimmung der optimalen Dichteverteilung zu vereinfachen, wird diese Funktion $\rho(r)$ nach dem Fourierprinzip als Linearkombination von Basisfunktionen $\eta_{K}$, mit einem Faktor $d_{K}^{\mu\nu}$ ausgedrückt (Gleichung \ref{gaussian}). \begin{equation} \label{gaussian} \rho_{\mu\nu}(r) = \sum_{K} d_{K}^{\mu\nu} \eta_{K} \end{equation} Die hierfür verwendeten Basisfunktionen können variieren, sind aber meistens Slater oder Gauss-Funktionen \cite{DFTORCA}. \chapter{Methoden und Tools} Für die Entwicklung der Reaktionsnetzwerke (Kapitel \ref{chmod}), die Beschränkung der entstehenden Graphen (Kapitel \ref{chnetx}), die Berechnung der NMR-Spektren (Kapitel \ref{chorca}) und die Lösung (Kapitel \ref{chgurobi}) und Verifikation (Kapitel \ref{chmodpet}) des ILP wurden verschiedene Programme verwendet, welche im Folgenden näher erläutert werden. Ebenso wird die Methode definiert, mit der die NMR-Spektren untereinander verglichen werden (Kapitel \ref{chnmrsi}). \section{Graptransformation Derivation Graph}\label{chmod} M\o{}d kurz für Med\O{}lDatschgerl ist ein pythonbasiertes Graphtransformtions Tool, welches mit dem DPO-Ansatz aus Graphtransformationsregeln und aus Startgraphen neue Graphen bestimmt. Aus den durch Graphtransformation gefundenen neuen Graphen wird ein Hypergraph definiert. Die Graphen in M\o{}d können dabei erzeugt werden, indem sie über einen SMILES-String geladen oder über einen GML-String definiert werden. Im Folgenden wird das selbe Molekül Cyclohexan über seinen SMILES- und einen GML-String definiert. \begin{lstlisting}[language=Python] cyclohexan = smiles ("C1CCCCC1") cyclohexan = Graph.fromGMLString( """ graph [ node [ id 0 label "C" ] node [ id 1 label "H" ] node [ id 2 label "C" ] node [ id 3 label "H" ] node [ id 4 label "C" ] node [ id 5 label "H" ] node [ id 6 label "C" ] node [ id 7 label "H" ] node [ id 8 label "C" ] node [ id 9 label "H" ] node [ id 10 label "C" ] node [ id 11 label "H" ] edge [ source 0 target 1 label "-" ] edge [ source 0 target 2 label "-" ] edge [ source 2 target 3 label "-" ] edge [ source 2 target 4 label "-" ] edge [ source 4 target 5 label "-" ] edge [ source 4 target 6 label "-" ] edge [ source 6 target 7 label "-" ] edge [ source 6 target 8 label "-" ] edge [ source 8 target 9 label "-" ] edge [ source 8 target 10 label "-" ] edge [ source 10 target 11 label "-" ] edge [ source 10 target 0 label "-" ] ] """ , name = "Cyclohexan) \end{lstlisting} Wenngleich die Definition des Moleküls als String deutlich einfacher ist, erlaubt der GML-String die Definition nicht chemischer Graphen und auch von Molekülzuständen, die keinen SMILES-String zugeordnet haben. Eine Regel kann nur über einen GML-String definiert werden. Für eine Regel müssen dabei linke und rechte Seite sowie der Kontext angegeben werden. Um zum Beispiel den Ring des Cyclohexans aufzuspalten um n-Hexan zu erhalten, würde eine Regel benötigt, welche eine C-C-Einfachbindung aufspaltet und die Kohlenstoffe mit Wasserstoff sättigt. \begin{lstlisting}[language=Python] cracking = Rule.fromGMLString( """ rule [ left [ edge [ source 1 target 2 label "-" ] ] context [ node [ id 1 label "C" ] node [ id 2 label "C"] ] right [ node [ id 3 label "H" ] node [ id 4 label "H" ] edge [ source 1 target 3 label "-" ] edge [ source 2 target 4 label "-" ] ] ] ) \end{lstlisting} Alle definierten Graphen werden in der Variable inputGraphs und alle definierten Regeln in der Variable inputRules gespeichert. Mit der Variable inputGraphs kann ein Derivationgraph DG erzeugt werden. Auf die in dem Derivationgraph enthaltenen Moleküle können die Regeln mit einer benutzerdefinierten Abfolge angewendet werden. Die in der Variable inputRules zusammengefassten Regeln können dabei allgemein oder eine spezifische Regel mit ihrem Variablennamen angewendet werden. Es ist ebenfalls möglich nach jeder Regelanwendung eine restriction-Funktion aufzurufen, in welcher unerwünschte Eigenschaften der durch die Regeln erzeugten Graphen abgelehnt werden können. \begin{lstlisting}[language=Python] dg = DG(graphDatabase=inputGraphs) dg.build().execute( addSubset(inputGraphs) >> rightPredicate[ restriction ]( repeat(inputRules)) ) ) \end{lstlisting} Der durch die Regelanwendung erzeugte Hypergraph kann über verschiedene print-Befehle in einer latex-generierten PDF-Datei dargestellt werden. Dafür muss ein FlowPrinter, DGPrinter oder GraphPrinter Objekt erzeugt werden \cite{mod}. \begin{lstlisting} flowPrinter = FlowPrinter() flowPrinter.printUnfiltered = False postSection("Loaded Graphs") for a in inputGraphs: a.print() postSection("Loaded Rules") for a in inputRules: a.print() postSection("Product Graphs") for a in dg.vertices: a.graph.print() \end{lstlisting} M\o{}d wurde verwendet, um die Einwirkung von natürlich auftretenden Reaktionen auf natürliche oder synthetische Moleküle zu untersuchen. \section{networkx}\label{chnetx} Networkx ist ein pythonbasiertes Netzwerk- und Graphenanalysetool \cite{networkx}. Eine Funktion ist dabei das Identifizieren aller sehnenlosen Zyklen. Während der Laufzeit des Graphtransformationstools wurden verschiedene Analysen dieser sehnenlosen Zyklen als Beschränkung verwendet. \section{ORCA}\label{chorca} ORCA ist ein quantenchemisches Rechnungsprogramm auf DFT Basis. Es erlaubt unter anderem Geometrieoptimierung und NMR-Spektrenberechnung von chemischen Molekülen \cite{ORCA}. Für diese Arbeit wurde ORCA 6.1.1 verwendet. Für die durch das Graphtransformationstool generierten Moleküle wurden die Graphen und SMILES-Strings bestimmt. Die zugehörigen XYZ-Koordinatendateien wurden aus der RDKit-Bibliothek über ihren SMILES-String bezogen \cite{rdkit} oder mit dem chemischen Zeichenprogramm Avogadro2 \cite{avogadro} gezeichnet. In beiden Fällen wurde die Geometrie von ORCA mit den der Dichtefunktion PBE, D4 Dispersionskorrektur und dem Basissatz DEF2-SVP optimiert. Anschließend wurden für die optimierte Struktur mit der Dichtefunktion TPSS, dem Basissatz pcSseg-3 und mit autoaux und tightscf die chemischen Verschiebungen berechnet. \begin{lstlisting} !PBE D4 DEF2-SVP OPT !TPSS pcSseg-3 autoaux tightscf NMR \end{lstlisting} Die selben Berechnungen für TMS durchgeführt und in der finalen Bestimmung der NMR-Spektren angegeben. Die Kohlenstoff Verschiebung von TMS lag dabei bei 179.73 ppm. Die Bestimmung chemisch-äquivalenter Kerne wurde von Hand durchgeführt. Für jedes mit dem Graptransformationstool erzeugte Molekül wurde ein $^{13}$C-NMR berechnet. \section{NMR-Vergleich}\label{chnmrsi} Um die von dem quantenchemischen Rechnungsprogramm berechneten NMR-Spektren zu Vergleichen und ihre Ähnlichkeit zu quantifizieren wurde die Binning-Methode mit anschließender Crosscorrelation angewandt \cite{nmrbin}. Jeder Peak $j$ des Spektrums setzt sich aus seiner Intensität $I_{j}$ und seiner Verschiebung $\nu_{j}$ zusammen. Zuerst wird für ein gegebenes NMR-Spektrum die beiden Peaks mit der größten chemischen Verschiebung $\nu_{max}$ und der kleinsten chemischen Verschiebung $\nu_{min}$ bestimmt und eine Bingröße $b$ festgelegt. Der größte Wert $\nu_{max}$ wird dabei aufgerundet und um eine Bingröße $b$ erhöht (Gleichung \ref{binupper}). Der kleinste Wert $\nu_{min}$ wird abgerundet und um eine Bingröße $b$ verringert (Gleichung \ref{binlower}). \begin{equation} \label{binupper} \nu_{max} = ceiling(max) + b \end{equation} \begin{equation} \label{binlower} \nu_{min} = floor(min) - b \end{equation} Aus der Obergrenze $\nu_{max}$ und Untergrenze $\nu_{min}$ wird mit der Bingröße die nötige Anzahl an Bins $n$ bestimmt (Gleichung \ref{binnum}). \begin{equation} \label{binnum} n = ceiling(\frac{\nu_{max} - \nu_{min}}{b}) \end{equation} In welchem Bin $i$ ein Peak $j$ liegt wird bestimmt, indem der kleinste Bin bestimmt wird, dessen unter Grenze kleiner als die chemische Verschiebung des Peaks $\nu_{j}$ ist (Gleichung \ref{binzuor}) \begin{equation} \label{binzuor} i = floor((\nu_{j} - \nu{min}) / b) \end{equation} Die Intensität des Bins $I_{i}$ ist dabei die Summe aller in ihm liegenden Peaks (Gleichung \ref{binning}). \begin{equation} \label{binning} I_{i} = \sum_{b \times i \leq \nu_{j} - \nu_{min} \leq b \times (i + 1)} I_{j} \end{equation} Die Intensität wird über alle Bins normalisiert. \newline Um zwei Spektren A und B zu vergleichen wird anstelle der Randwerte des einen Spektrums gemeinsame Randwerte für beide Spektren nach dem gleichen System bestimmt. Diese gemeinsamen Randwerte werden verwendet um die Peaks auf die Bins aufzuteilen. Anschließend werden die Bins für jedes Spektrum einzeln normiert. Der Überlapp eines Bins des ersten Spektrums $i$ mit einem Bin des zweiten Spektrums $j$ ist dabei das Produkt der Werte beider Bins. Der Teilüberüberlapp für einen Bin $O_{i}$ ist die Summe aus dem Überlapp mit dem Bin mit gleichem Index und dem mit 0.5 gewichteten Überlapp der direkt benachbarten Bins (Gleichung \ref{overlap}). \begin{equation} \label{overlap} O_{i} = 0.5 \times (I_{i} \times I_{j-1}) + (I_{i} \times I_{j}) + 0.5 \times (I_{i} \times I_{j+1}) : i = j \end{equation} Der Gesamtüberlapp $O$ ist die Summe aller Teilüberlappe $O_{i}$ (Gleichung \ref{fulloverlap}) \begin{equation} \label{fulloverlap} O_{AB} = \sum_{i} O_{i} \end{equation} Aus dieser Crosskorrelation beider Spektren $O_{AB}$ und der Selbstkorrelation beider Spektren $O_{AA}$ und $O_{BB}$, wird der Ähnlichkeitsindex $SI$ berechnet (Gleichung \ref{similarity}) \begin{equation} \label{similarity} SI_{ij} = \frac{O_{AB}}{\sqrt{O_{AA} \times O_{BB}}} \end{equation} Für alle Reaktionen eines Derivationgraphes wurden Referenzspektren bestimmt. Jedes mit dem quantenchemisches Rechnungsprogramm berechnete $^{13}$C-NMR, wurde mit den Referenzspektren verglichen. Die für ein Referenzspektrum normierten Ähnlichkeitsindizes wurden als Wahrscheinlichkeit angenommen, dass beide NMR-Spektren die Gleiche Substanz abbilden. \section{Gurobipy}\label{chgurobi} Gurobipy \cite{gurobi} ist ein pythonbasierter ILP-Solver. Zur Bestimmung der optimalen Parameter wird ein ILP-Baum verwendet. Die bestimmten Wahrscheinlichkeiten wurden mit dem Derivationgraph kombiniert. Darauf basierend wurde Gruobipy verwendet, um einen diese Wahrscheinlichkeiten maximierenden Reaktionpfad in dem Derivationgraph zu identifizieren. \section{Petrinetzzertifikat}\label{chmodpet} Eine weitere Funktion von M\o{}d ist das Zertifizieren eines bestimmten Reaktionpfades des Hypergraphen mittels Petrinetzen. Dafür müssen die mittels Gurobi bestimmten Parameter für den Derivationgraph in M\o{}d festgelegt werden. Ein mittels ILP bestimmter Reaktionspfad kann so nähergehend untersucht werden, indem Reaktionen auf Erfüllbarkeit basierend auf den gegebenen Startmolekülen überprüft werden. \chapter{Ergebnisse} Mit den erläuterten Werkzeugen und Methoden werden drei Reaktionsnetzwerke analysiert. Dabei werden zwei natürlich auftretende Ein-Molekül-Reaktionen (Kapitel \ref{chonmolre}) und eine Multi-Molekül-Reaktion (Kapitel \ref{chmultmore}) aus der Organischensynthesechemie betrachtet. Die Beispiele sind dabei nicht erschöpfend und dienen lediglich als Demonstration der Methode und als ein mögliches Anwendungsschema. \section{Ein-Molekül-Reaktionen}\label{chonmolre} Für eine einfache Demonstration der Kombination von ILP wurden simple Reaktionen mit hoher natürlicher Selektivität gesucht. Dabei wurden die natürlichen Synthesen der Stoffe Koffein und Vanillin näher untersucht. Diese Synthesen erfolgen im Wesentlichen durch Reaktionen an einem zentralen Molekül, wobei die Reaktionspartner nicht weiter relevant bleiben. Sie werden im Folgenden als Ein-Molekül-Reaktionen bezeichnet. \subsection{Koffein} Koffein ist ein dreifachmethyliertes Xanthin, welches vor allem in Schwarztee und Kaffeebohnen vorkommt. Es wurde erstmals in den 1820 Jahren isoliert. Heute sind koffeinhaltige Lebensmittel aus dem Alltag nicht mehr wegzudenken. Der Wirkstoff erfüllt aber auch in zahlreichen weiteren Anwendungsfeldern eine signifikante Rolle. In der Natur wird Koffein aus den Nukleotidbasen Adenin und Guanin über die Oxidierung zu Xanthosin und anschließende Methylierung der Stickstoffe des Rings hergestellt (Abbildung \ref{coffsynth}). Die Methylierung erfolgt erst in der 7, dann in der 3 und zuletzt in der 1 Position. Diese geordnete Methylierung erfolgt mit hoher Selektivität und eignet sich so für eine nähere Untersuchung \cite{caffeinesynth}. \begin{figure} [h] \includegraphics[width=\textwidth]{coffeinesynthesis.png} \caption{Synthese von Koffein in Pflanzen aus den Nukleotiden Adenin und Guanin. Erst werden die Purinringstrukturen oxidiert und nach der Abspaltung der Ribose die Stickstoffe methyliert. Die Methylierung erfolgt in der 7, 3 und zuletzt der 1 Position \cite{caffeinesynth}. } \label{coffsynth} \end{figure} \FloatBarrier \subsection{Vanillin und 3,4-Dihydroxyaldehyd} Vanillin ist eines der beliebtesten aromatischen Moleküle mit hoher Verwendung in Lebensmitteln Der überwiegende Teil wird synthetisch hergestellt. In der Natur wird Vanillin in Vanille Orchideen aus Phenylalanin synthetisiert. Zuerst wird dabei der Aminrest in einer Eliminierungsreaktion entfernt, wodurch Zimtsäure entsteht. Die Zimtsäure wird dann in der 3- und 4-Position des Phenylrests hydroxiliert und die Doppelbingung zu einem Aldehyd oxidiert (Abbildung \ref{vansynth}). Die Hydroxylierung und die Oxidation zu einem Aldehyd sind in ihrer Reihenfolge selektiv und eignet sich für eine nähere Untersuchung \cite{vanillinesynth}. \begin{figure} [h] \includegraphics[width=0.9\textwidth]{vanillsynth.png} \caption{Synthese von Vanilline in der Vanille Orchidee. Die Synthese erfolgt aus der Aminosäure Phenylalanin. Dabei wird zuerst die Aminogruppe eliminiert, anschließend der Phenylrest hydroxyliert und die aus der Eliminierung entstandene Doppelbindung oxidiert (rot). Der letzte Schritt ist die Methylierung des 4-Hydroxyrests \cite{vanillinesynth}. } \label{vansynth} \end{figure} \FloatBarrier \section{Graphtransformation: Ein-Molekül-Reaktionen} Für die Ein-Molekül-Reaktionen wurden zuerst in dem Graphtransformationstool Transformationsregeln definiert. Diese Regeln wurden dann auf ein Startmolekül angewandt, welches Teil des natürlichen Syntheseprozesses ist. Der betrachtete Syntheseausschnitt wurde dabei mit möglichst wenigen und möglichst einfachen Regeln dargestellt. \subsection{Koffein} Für die Synthese von Koffein wurde mit dem Purinring Xanthin gestartet (Abbildung \ref{xanthine}). \begin{figure} [h] \begin{center} \includegraphics[width=0.3\textwidth]{xanthine.png} \caption{Darstellung von Xanthin in dem Graphtransformationstool.} \label{xanthine} \end{center} \end{figure} Ebenso wurde für die Transformationsregel zum Austausch eines an einen Stickstoff gebundenen Wasserstoff durch eine Methylgruppe definiert. (Abbildung \ref{nmeth}). \begin{figure} [h] \includegraphics[width=\textwidth]{nmeth.png} \caption{Darstellung der Regel zum Methylieren eines protonierten Stickstoffs in dem Graphtransformationstool.} \label{nmeth} \end{figure} Durch das wiederholte anwenden der Regel auf Xanthin konnte ein Reaktionsnetzwerk entwickelt werden, an dessen Ende die Synthese von Koffein steht (Abbildung \ref{coffsynthmod}). Dabei wurden die folgenden Zwischenprodukte bestimmt. \begin{itemize} \item $p_{0,0}$: 1-Methylxanthin \item $p_{0,1}$: 3-Methylxanthin \item $p_{0,2}$: 7-Methylxanthin \item $p_{0,3}$: 1,3-Dimethylxanthin oder Theophyllin \item $p_{0,4}$: 1,7-Dimethylxanthin oder Paraxanthin \item $p_{0,5}$: 3,7-Dimethylxanthin oder Theobromin \end{itemize} \begin{figure} [p] \includegraphics[width=\textwidth]{coffsynthmod.png} \caption{Darstellung aller Synthesepfade von Koffein aus Xanthin in dem Graphtransformationstool. Die Graphen $p_{0,0}$ bis $p_{0,5}$ stellen dabei die Zwischenprodukte 1-Methylxanthin, 3-Methylxanthin, 7-Methylxanthin, 1,3-Dimethylxanthin oder Theophylline, 1,7-Dimethylxanthin oder Paraxanthin und 3,7-Dimethylxanthin oder Theobromin dar.} \label{coffsynthmod} \end{figure} \subsection{Vanillin und 3,4-Dihydroxyaldehyd} Für die Synthese von Vanillin wurde mit Zimtsäure gestartet (Abbildung \ref{zimt}). Das natürliche Produkt ist dabei E-Zimtsäure. \begin{figure} [h] \begin{center} \includegraphics[width=0.3\textwidth]{zimt.png} \caption{Darstellung von Zimtsäure in dem Graphtransformationstool.} \label{zimt} \end{center} \end{figure} Zur Umsetzung der Zimtsäure zu dem Vanillin Zwischenprodukts, 3,4-Dihydroxybenzaldehyd, wurde eine Regel zur Hydroxylierung der aromatischen Kohlenstoffe (Abbildung \ref{hydrox}) und zur Oxidierung einer Doppelbindung mit Carboxylgruppe zu einem Aldehyd (Abbildung \ref{aldhyd}) definiert. \begin{figure} [h] \includegraphics[width=\textwidth]{hydrox.png} \caption{Darstellung der Regel zum Hydroxylieren des Phenylrests in dem Graphtransformationstool.} \label{hydrox} \end{figure} \begin{figure} [h] \begin{center} \includegraphics[width=0.5\textwidth]{aldehyd.png} \caption{Darstellung der Regel zum oxidieren der Zimtsäure Doppelbindung zu einem Aldehyd in dem Graphtransformationstool.} \label{aldhyd} \end{center} \end{figure} Durch das wiederholte Anwenden beider Regeln auf Zimtsäure konnte ein Reaktionsnetzwerk entwickelt werden, an dessen Ende die Synthese von 3,4-Dihydroxybenzaldehyd steht (Abbildung \ref{vanisynthmod}). Die breite Regelanwendung bedeutet die Ermittlung vieler Moleküle, welche nicht mehr zu dem Zielmolekül umgesetzt werden können. Darum wurden die Zwischenprodukte auf die Moleküle beschränkt, die tatsächlich Zwischenprodukte darstellen. Dabei wurden die folgenden Zwischenprodukte bestimmt. \begin{itemize} \item $p_{0,4}$: 3-Hydroxzimtsäure oder m-Coumarsäure \item $p_{0,5}$: 4-Hydroxzimtsäure oder p-Coumarsäure \item $p_{0,6}$: Benzaldehyd \item $p_{0,13}$: 3,4-Dihydroxzimtsäure oder Kaffeesäure \item $p_{0,16}$: 3-Hydroxybenzaldehyd \item $p_{0,17}$: 4-Hydroxybenzaldehyd \end{itemize} \begin{figure} [p] \includegraphics[width=\textwidth]{vanisynthmod.png} \caption{Darstellung aller Synthesepfade von 3,4-Dihydroxybenzaldehyd aus Zimtsäure in dem Graphtransformationstool. Nur ein Teil aller Zwischenprodukte wird dabei tatsächlich zu dem gewünschten Produkt.} \label{vanisynthmod} \end{figure} \section{NMR: Ein-Molekül-Reaktionen} Für alle Zwischenprodukte wurde mit dem quantenchemischen Rechnungsprogramm ein $^{13}$C-NMR berechnet und ein dazugehöriges Literaturspektrum recherchiert. \FloatBarrier \subsection{Koffein} Für alle Xanthinringe lässt sich sagen, dass ausschließlich die Protonen der Methygruppen chemisch-äquivalent zueinander sind. In einem $^{13}$C-NMR weisen also alle Kohlenstoffe ein eigenes Signal mit Intensität 1 auf. Die mit dem quantenchemischen Rechnungsprogramm berechneten NMR-Spektren wurden den recherchierten Literaturspektren gegenübergestellt (Tabellen \ref{tabelxanthin}, \ref{tabel1xanthin}, \ref{tabel3xanthin}, \ref{tabel7xanthin}, \ref{tabel13xanthin}, \ref{tabel17xanthin}, \ref{tabel37xanthin}, \ref{tabel137xanthin}) und ihre Differenz zueinander bestimmt. \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für Xanthin mit einem Literatur Spektrum \cite{xanthinnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 3.14 ppm.} \begin{center} \begin{tabular}{ |c|c|c| } \label{tabelxanthin} Berechnet & Exp & Differenz: Berechnet - Exp\\ \hline 112.57 & 108.8 & 3.77 \\ 138.61 & 137.2 & 1.413 \\ 151.03 & 140.5 & 10.53 \\ 152.87 & 152.4 & 0.47 \\ 155.64 & 156.1 & -0.46 \\ & & 3.14 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für 1-Methylxanthin mit einem Literatur Spektrum \cite{manyxanthinesnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 2.39 ppm.} \begin{center} \begin{tabular}{ |c|c|c| } \label{tabel1xanthin} Berechnet & Exp & Differenz: Berechnet - Exp\\ \hline 30.18 & 27 & 3.18 \\ 112.28 & 106.4 & 5.88 \\ 137.88 & 140.6 & -2.72 \\ 150.37 & 147 & 3.37 \\ 153.13 & 151.2 & 1.93 \\ 157.91 & 155.2 & 2.71 \\ & & 2.39 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für 3-Methylxanthin mit einem Literatur Spektrum \cite{manyxanthinesnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 2.35 ppm.} \begin{center} \begin{tabular}{ |c|c|c| } \label{tabel3xanthin} Berechnet & Exp & Differenz: Berechnet - Exp\\ \hline 31.34 & 28.8 & 2.54 \\ 112.87 & 106.9 & 5.97 \\ 138.12 & 140.5 & -2.38 \\ 153.46 & 149.4 & 4.06 \\ 154.78 & 151.2 & 3.58 \\ 155 & 154.7 & 0.30 \\ & & 2.35 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für 7-Methylxanthin mit einem Literatur Spektrum \cite{7xanthinnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 2.73 ppm.} \begin{center} \begin{tabular}{ |c|c|c| } \label{tabel7xanthin} Berechnet & Exp & Differenz: Berechnet - Exp\\ \hline 36.69 & 33.03 & 3.66 \\ 113.78 & 106.9 & 6.88 \\ 143.18 & 143.01 & 0.17 \\ 151.13 & 149.3 & 1.83 \\ 153.94 & 151.35 & 2.59 \\ 157.1 & 155.85 & 1.25 \\ & & 2.73 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für Theophyllin mit einem Literatur Spektrum \cite{manyxanthinesnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 2.61 ppm.} \begin{center} \begin{tabular}{ |c|c|c| } \label{tabel13xanthin} Berechnet & Exp & Differenz: Berechnet - Exp\\ \hline 29.23 & 27.7 & 1.53 \\ 32.05 & 29.7 & 2.35 \\ 112.36 & 106.4 & 5.96 \\ 137.43 & 140.5 & -3.07 \\ 152.62 & 147.8 & 4.82 \\ 155.4 & 151.2 & 4.20 \\ 156.89 & 154.4 & 2.49 \\ & & 2.61 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für Paraxanthin mit einem Literatur Spektrum \cite{manyxanthinesnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 3.19 ppm.} \begin{center} \begin{tabular}{ |c|c|c| } \label{tabel17xanthin} Berechnet & Exp & Differenz: Berechnet - Exp\\ \hline 28.59 & 26.7 & 1.89 \\ 36.8 & 32.9 & 3.9 \\ 113.44 & 106.5 & 6.94 \\ 142.72 & 143 & -0.28 \\ 151.81 & 147.4 & 4.41 \\ 153.04 & 151.1 & 1.94 \\ 158.8 & 155.3 & 3.5 \\ & & 3.19 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für Theobromin mit einem Literatur Spektrum \cite{theobrominnoMenmr}. Die Methylsignale (etwa 30 ppm) wurden aus einem anderen Spektrum genommen, da diese Angabe fehlte \cite{theobrominnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 2.63 ppm.} \begin{center} \begin{tabular}{ |c|c|c| } \label{tabel37xanthin} Berechnet & Exp & Differenz: Berechnet - Exp\\ \hline 30.96 & 29.3 & 1.66 \\ 36.7 & 33.9 & 2.8 \\ 114.14 & 108.1 & 6.04 \\ 142.73 & 142.6 & 0.13 \\ 153.39 & 150 & 3.39 \\ 155.92 & 153.1 & 2.82 \\ 156.49 & 154.9 & 1.59 \\ & & 2.63 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für Koffein mit einem Literatur Spektrum \cite{manyxanthinesnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 3.77 ppm.} \begin{center} \begin{tabular}{ |c|c|c| } \label{tabel137xanthin} Berechnet & Exp & Differenz: Berechnet - Exp\\ \hline 28.8 & 27.5 & 1.3 \\ 31.72 & 29.3 & 2.42 \\ 36.86 & 33.1 & 3.76 \\ 113.66 & 106.6 & 7.06 \\ 142.13 & 142.8 & -0.67 \\ 153.86 & 148.1 & 5.76 \\ 155.29 & 151 & 4.29 \\ 158.29 & 154.5 & 3.79 \\ & & 3.77 \\ \hline \end{tabular} \end{center} \end{table} \newline Die Methylwasserstoffe lassen sich dabei den Signalen um etwa 30 ppm zuordnen. Das 1-Methylsignal weist mit 28.6 bis 30.18 von dem quantenchemischen Rechnungsprogramm berechnet und 26.7 bis 27.7 in den experimentellen Spektren den niedrigsten Wert auf. Die mittlere Verschiebung weist das 3-Methylsignal mit 31.0 bis 32.1 von dem quantenchemischen Rechnungsprogramm berechnet und 28.8 bis 29.7 in den experimentellen Spektren auf. Die größte Verschiebung ist dem 7-Methylsignal mit 36.7 bis 36.9 von dem quantenchemischen Rechnungsprogramm berechnet und 32.9 bis 33.9 in den experimentellen Spektren zuzuordnen. Die Signale des Purinrings sind für den größten Teil konstant über alle Moleküle, wobei Xanthin mit 151.03 von dem quantenchemischen Rechnungsprogramm berechnet, zu 140.5 im Experiment eine besonders große Differenz aufweist. Die größten Unterschiede zwischen den berechneten und experimentellen Spektren sind über alle Xanthine bei dem Kohlenstoff neben dem 7-Stickstoff. Über alle Moleküle gemittelt lässt sich bestimmen, dass die von dem quantenchemischen Rechnungsprogramm berechneten Spektren im Durchschnitt um 2.5 ppm nach oben verschoben sind. \FloatBarrier \subsection{Vanillin und 3,4-Dihydroxyaldehyd} Für 3-hydroxylierte Zwischenprodukte gilt, dass keine Kerne chemisch-äquivalent zueinander sind. Wenn nicht oder nur in der 4-Position hydroxyliert ist, sind jeweils die 2 mit 6 und 3 mit 5 Kohlenstoffe im Phenylrest, sowie die jeweils an diese gebundenen Wasserstoffe zueinander chemisch-äquivalent. In einem $^{13}$C-NMR weisen also alle Kohlenstoffe in 3-hydroxilierten Moleküle ein eigenes Signal mit Intensität 1 auf. In nicht oder 4-hydroxylierten Molekülen verschmelzen die Signale der Kohlenstoffe in der 2 und 6 Position sowie die in der 3 und 5 Position. Es resultieren zwei Signalen mit Intensität 2. Die mit dem quantenchemischen Rechnungsprogramm berechneten NMR-Spektren wurden den recherchierten Literaturspektren gegenübergestellt (Tabellen \ref{tabelzimt}, \ref{tabelmcoum}, \ref{tabelpcoum}, \ref{tabelald}, \ref{tabelcafac}, \ref{tabel3ald}, \ref{tabel4ald}, \ref{tabel34ald}) und ihre Differenz zueinander bestimmt. \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für Zimtsäure mit einem Literatur Spektrum \cite{cinamnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 4.02 ppm.} \begin{center} \begin{tabular}{ |c|c|c|c| } \label{tabelzimt} Height & Exp & Berechnet & Differenz: Berechnet - Exp\\ \hline 1 & 117.33 & 120.33 & 3 \\ 2 & 128.39 & 133.48 & 5.09 \\ 2 & 128.98 & 134.2 & 5.22 \\ 1 & 130.77 & 135.91 & 5.14 \\ 1 & 134.04 & 140.47 & 6.43 \\ 1 & 147.14 & 151.41 & 4.27 \\ 1 & 172.6 & 171.56 & -1.04 \\ & & & 4.02 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für m-Coumarsäure mit einem Literatur Spektrum \cite{mcoumnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 3.53 ppm.} \begin{center} \begin{tabular}{ |c|c|c|c| } \label{tabelmcoum} Height & Exp & Berechnet & Differenz: Berechnet - Exp\\ \hline 1 & 115.3 & 118.6 & 3.3 \\ 1 & 118.6 & 121.09 & 2.49 \\ 1 & 119.2 & 121.69 & 2.49 \\ 1 & 120.8 & 122.36 & 1.56 \\ 1 & 131 & 134.46 & 3.46 \\ 1 & 137.1 & 141.29 & 4.19 \\ 1 & 146.5 & 153.69 & 7.19 \\ 1 & 159.1 & 164.15 & 5.05 \\ 1 & 170.5 & 172.55 & 2.05 \\ & & & 3.53 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für p-Coumarsäure mit einem Literatur Spektrum \cite{pcoumnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 4.71 ppm.} \begin{center} \begin{tabular}{ |c|c|c|c| } \label{tabelpcoum} Height & Exp & Berechnet & Differenz: Berechnet - Exp\\ \hline 1 & 115.41 & 117.18 & 1.77 \\ 2 & 115.83 & 118.65 & 2.82 \\ 1 & 125.36 & 132.76 & 7.40 \\ 2 & 130.17 & 135.75 & 5.58 \\ 1 & 144.27 & 149.67 & 5.40 \\ 1 & 159.67 & 166.09 & 6.42 \\ 1 & 168.05 & 171.61 & 3.56 \\ & & & 4.71 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für Benzaldehyd mit einem Literatur Spektrum \cite{aldhnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 4.80 ppm.} \begin{center} \begin{tabular}{ |c|c|c|c| } \label{tabelald} Height & Exp & Berechnet & Differenz: Berechnet - Exp\\ \hline 2 & 128.9 & 133.45 & 4.55 \\ 2 & 129.5 & 135.86 & 6.36 \\ 1 & 134.2 & 139.15 & 4.95 \\ 1 & 136.4 & 141.47 & 5.07 \\ 1 & 192 & 195.09 & 3.09 \\ & & & 4.80 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für Kaffeesäure mit einem Literatur Spektrum \cite{caffacdnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 2.18 ppm.} \begin{center} \begin{tabular}{ |c|c|c|c| } \label{tabelcafac} Height & Exp & Berechnet & Differenz: Berechnet - Exp\\ \hline 1 & 147.9 & 155.79 & 7.89 \\ 1 & 146.5 & 149.88 & 3.38 \\ 1 & 129.3 & 132.99 & 3.69 \\ 1 & 123.1 & 131.94 & 8.84 \\ 1 & 121.7 & 118.15 & -3.55 \\ 1 & 116.4 & 116.16 & -0.24 \\ 1 & 141.5 & 149.28 & 7.78 \\ 1 & 114.6 & 110.95 & -3.65 \\ 1 & 176.2 & 171.66 & -4.55 \\ & & & 2.18 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für 3-Hydroxybenzaldehyd mit einem Literatur Spektrum \cite{mhydaldhnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 3.83 ppm.} \begin{center} \begin{tabular}{ |c|c|c|c| } \label{tabel3ald} Height & Exp & Berechnet & Differenz: Berechnet - Exp\\ \hline 1 & 115 & 122.95 & 7.95 \\ 1 & 122 & 123.33 & 1.33 \\ 1 & 123 & 124.61 & 1.61 \\ 1 & 131 & 133.96 & 2.96 \\ 1 & 138 & 143.41 & 5.41 \\ 1 & 158 & 164.39 & 6.39 \\ 1 & 194 & 195.16 & 1.16 \\ & & & 3.83 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für 4-Hydroxybenzaldehyd mit einem Literatur Spektrum \cite{phydaldhnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt 3.35 ppm.} \begin{center} \begin{tabular}{ |c|c|c|c| } \label{tabel4ald} Height & Exp & Berechnet & Differenz: Berechnet - Exp\\ \hline 2 & 116.5 & 118.46 & 1.96 \\ 1 & 130.3 & 135.08 & 4.78 \\ 2 & 132.7 & 137.69 & 4.99 \\ 1 & 163.8 & 168.28 & 4.48 \\ 1 & 191 & 191.53 & 0.53 \\ & & & 3.35 \\ \hline \end{tabular} \end{center} \end{table} \begin{table}[p] \caption{Vergleich des mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum für 3,4-Dihydroxybenzaldehyd mit einem Literatur Spektrum \cite{mphydaldhnmr}. Die Werte sind in ppm angegeben. Die durchschnittliche Differenz beträgt -3.81 ppm.} \begin{center} \begin{tabular}{ |c|c|c|c| } \label{tabel34ald} Height & Exp & Berechnet & Differenz: Berechnet - Exp\\ \hline 1 & 113.1 & 114.22 & 1.12 \\ 1 & 117.79 & 115.21 & -2.58 \\ 1 & 133.4 & 124.59 & -8.81 \\ 1 & 134.66 & 130.13 & -4.53 \\ 1 & 149.32 & 145.44 & -3.88 \\ 1 & 158.08 & 151.27 & -6.81 \\ 1 & 191.41 & 190.26 & -1.15 \\ & & & -3.81 \\ \hline \end{tabular} \end{center} \end{table} \newline Bei den Zwischenprodukten der 3,4-Dihydroxybenzaldehyd-Synthese lässt sich im Vergleich zu den Zwischenprodukten der Koffein-Synthese kein so konstantes Muster erkennen. Das ist unter anderem darauf zurückzuführen, dass der deutliche Elektronegativitäten-Unterschied zwischen Sauerstoff und Wasserstoff einen großen Einfluss auf das Ringsystem erwirkt. Dadurch sind die Signale des Kerngerüsts deutlich weniger konstant. Das Endprodukt 3,4-Dihydroxybenzaldehyd sticht dabei trotzdem heraus, da es im Gegensatz zu allen anderen Spektren durchschnittlich eine größere berechnete chemische Verschiebung als experimentell bestimmte chemische Verschiebung aufweist. Über alle Spektren gemittelt haben die mit dem quantenchemischen Rechnungsprogramm berechneten Spektren trotzdem eine um 2.8 ppm größere chemische Verschiebung. \section{NMR $SI$: Ein-Molekül-Reaktionen} In bisherigen Analysen zu der Genauigkeit von berechneten NMR-Spektren im Vergleich zu experimentellen Spektren wurde festgestellt, dass einfache DFT Methoden die Effekte von schweren Atomen wie dem in TMS enthaltenden Silizium nicht akkurat bestimmen können. Insbesondere kann davon ausgegangen werden, dass die Frequenz des TMS-Kohlenstoff zu gering ausfällt \cite{heavyateff1,heavyatef2}. Das führt wiederum dazu, dass die chemische Verschiebung der Signale im mit dem quantenchemischen Rechnungsprogramm berechneten Spektrum zu hoch ausfällt. Dieser Wert ist systematisch, weshalb er als ungefähres Mittel der Unterschiede zwischen den mit dem quantenchemischen Rechnungsprogramm berechneten und experimentellen Spektren quantifiziert wurde. Auf alle Peaks der experimentellen Spektren wurden 2.5 ppm hinzuaddiert. Die Bingröße $b$ variiert in Literaturangaben zwischen 1 und 0.1 ppm \cite{nmrbin,nmrbinsize} mit einer Genauigkeit von 0.1 ppm. Deshalb wurde die Bingröße $b$ (Gleichung \ref{binsize}) variiert und die resultierenden Ähnlichkeitsindizes $SI_{b}$ über alle sinnvollen Werte gemittelt. \begin{equation} \label{binsize} b \in \{0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0\} \end{equation} \FloatBarrier \subsection{Koffein} Den Erwartungen entsprechend weist das mit dem quantenchemischen Rechnungsprogramm berechnete Spektrum von 7-Methylxanthin eine hohe Ähnlichkeit zu dem experimentellen Spektrum von 7-Methylxanthin auf. Das Spektrum von Paraxanthin ist jedoch noch etwas ähnlicher. Die Spektren von Theobromin haben im Vergleich dazu die größten Übereinstimmungen, wobei Koffein am zweit ähnlichsten ist. Allgemein scheint eine Identifikation des 7-Methylrest gut zu funktionieren, jedoch keine einfache Unterscheidung der restlichen Methylgruppen möglich zu sein. Da in allen Molekülen das Xanthin-Grundgerüst vorliegt und die Addition einer Methylgruppe dieses nicht zu sehr verändert, haben alle Moleküle eine signifikante Ähnlichkeit zueinander (Tabelle \ref{nmrsimcoff}). \begin{table}[h!] \caption{Resultierende normalisierte Wahrscheinlichkeiten beim Vergleich aller Zwischenprodukte der Koffeinsynthese, sowie Xanthin und Koffein mit den natürlichen Zwischenprodukten 7-Methyxanthin und Theobromin.} \begin{center} \begin{tabular}{ |c|c|c| } \label{nmrsimcoff} Molekül & 7-Methylxanthin & Theobromin\\ \hline Xanthin & 0.08 & 0.14\\ 1-Methylxanthin & 0.11 & 0.09\\ 3-Methylxanthin & 0.08 & 0.11\\ 7-Methylxanthin & 0.16 & 0.12\\ Theophylline & 0.11 & 0.09\\ Paraxanthin & 0.19 & 0.1\\ Theobromin & 0.07 & 0.18\\ Koffein & 0.19 & 0.17\\ \hline \end{tabular} \end{center} \end{table} \FloatBarrier \subsection{Vanillin und 3,4-Dihydroxyaldehyd} Bei der Vanillin Synthese liegen die chemischen Verschiebungen des Grundgerüsts weiter auseinander, weshalb manche Moleküle, wie 3-Hydroxybenzaldehyd nur eine sehr geringe Ähnlichkeit zu den natürlichen Zwischenprodukten aufweisen. Die p-Coumarsäure hat die größte Ähnlichkeit zu sowohl seinem eigenen experimentellen Spektrum, sowie zu dem von 4-Hydoxybenzaldehyd. Die zweitgrößte Ähnlichkeit zu dem experimentellen Spektrum von p-Coumarsäure haben Kaffeesäure und 4-Hydroxybenzaldehyd. Bei dem experimentellen Spektrum von 4-Hydroxybenzaldehyd hat das richtige Molekül die zweit größte Ähnlichkeit (Tabelle \ref{nmrsimvani}). \newline Im Vergleich zu den Xanthinen sind die Abweichungen der berechneten zu den experimentellen Spektren bei den Vanillevorprodukten deutlich weniger konstant. So hat Kaffeesäure mit 2.35 ppm die geringste Abweichung zwischen beiden zugehörigen Spektren. Trotzdem ist die absolute Verschiebung der einzelnen Peaks zueinander relativ groß. Darum ist die Varianz der Werte ist deutlich größer als bei den Spektren der Xanthine, weshalb die nicht normalisierte Crosskorrelation allgemein geringer ausfällt. Das berechnete Spektrum von Kaffeesäure hat dabei zu vielen Peaks in dem Literaturspektrum von p-Coumarsäure eine Distanz von unter 1.0 ppm, weshalb der Überlapp größer als erwartet ausfällt. Diese überlappenden Peaks weichen besonders stark von den Peaks im eigenen Referenzspektrum ab. Die berechneten Spektren von p-Coumarsäure und 4-Hydroxybenzaldehyd haben nur wenige überlappende Peaks mit beiden Literaturspektren. Diese überlappenden Peaks sind haben jedoch eine Intensität von 2, welche durch die Multiplikation der Intensitäten stärker gewichtet werden. Allgemein haben alle para-hydroxylierten Phenylgruppen größeren Überlapp mit der p-Coumarsäure. Bei 4-Hydroxybenzaldehyd scheint der Überlapp größer bei Molekülen ohne 3-Hydroxygruppe zu sein. \begin{table}[h!] \caption{Resultierende normalisierte Wahrscheinlichkeiten beim Vergleich aller Zwischenprodukte der Vanillesynthese mit den natürlichen Zwischenprodukten p-Coumarsäure und 4-Hydroxybenzaldehyd.} \begin{center} \begin{tabular}{ |c|c|c| } \label{nmrsimvani} Molekül & p-Coumarsäure & 4-Hydroxybenzaldehyd\\ \hline Zimtsäure & 0.06 & 0.11\\ p-Coumarsäure & 0.25 & 0.29\\ m-Coumarsäure & 0.08 & 0.08\\ Benzaldehyd & 0.06 & 0.12\\ Kaffeesäure & 0.22 & 0.1\\ 3-Hydroxybenzaldehyd & 0.01 & 0.02\\ 4-Hydroxybenzaldehyd & 0.22 & 0.17\\ 3,4-Dihydroxybenzaldehyd & 0.11 & 0.11\\ \hline \end{tabular} \end{center} \end{table} \FloatBarrier \section{ILP: Ein-Molekül-Reaktionen} Der Hypergraph der komplett explorierten Synthese wurde an Gurobi übergeben und um eine Kante, die in das Startmolekül hineinfließt und aus dem Zielmolekül herausfließt, erweitert. Die Wahrscheinlichkeiten wurden als konstante Werte $p_{i,e}$ den Kanten $e$ des Hypergraphen zugeordnet. Die Wahrscheinlichkeit des Moleküls wurde den Kanten zugeordnet, welche in den Knoten fließen. Die Kante 4 von Zimtsäure zu Benzaldehyd hat also die Wahrscheinlichkeit $p_{1,4} = 0.09$, dass diese Reaktion dem ersten Referenzspektrum mit p-Coumarsäure entspricht. Für jede Kante werden vier ganzzahlige Variablen $x_{e}$, $b_{e}$, $n_{1,e}$ und $n_{2,e}$ definiert. Die Variable $b_{e}$ ist dabei eine Wahrheitsvariable die definiert, ob eine Kante verwendet wird. Die Variable $x_{e}$ beschreibt wie oft die selbe Kante verwendet wird. Die Variablen $n_{1,e}$ und $n_{2,e}$ sind ebenfalls Wahrheitsvariablen und legen fest, welche zwei Reaktionen den Referenz NMR-Spektren entsprechen und deswegen zu der Zielfunktion dazuzählen. Dabei ist $n_{1,e}$ die Zuordnungsvariable zu dem ersten Referenzspektrum mit 7-Methylxanthin oder p-Coumarsäure und $n_{2,e}$ die Zuordnungsvariable zu dem zweiten Referenzspektrum mit 3,7-Methylxanthin oder 4-Hydroxybenzaldehyd. Die Variablen unterliegen den Beschränkungen, dass: \begin{itemize} \item $x_{e} \geq b_{e}$ \item $b_{e} \geq n_{1,e}$ \item $b_{e} \geq n_{2,e}$ \item s: ([], [Startmolekül]), $x_{s} = 1$ \item f: ([Zielmolekül], []), $x_{f} = 1$ \item $\forall v \in V$ inflow = outflow \item $n_{1,e} + n_{2,e} \leq 1$ \item $\sum_{e} n_{1,e} = 1$ \item $\sum_{e} n_{2,e} = 1$ \end{itemize} Unter diesen Beschränkungen soll ein Pfad durch den Hypergraphen bestimmt werden, welcher ein geordnetes Paar von Molekülen auswählt. Dieses Paar von Molekülen soll alle Zwischenprodukte der Synthese beinhalten und addiert die größte Wahrscheinlichkeit aufweisen. Die ausgehenden Kanten sind nicht in die Zielfunktion eingeflossen (Gleichung \ref{objective}). \begin{equation} \label{objective} max \: z = \sum_{e} 1000 \times (p_{1,e} \times n_{1,e} + p_{2,e} \times n_{2,e}) \newline - \frac{x_{e}}{ max(p_{1,e}, p_{2,e})} \end{equation} \FloatBarrier \subsubsection{Koffein} Der Hypergraph für die Koffeinsynthese wird in Tabelle \ref{hypergraphcaffsynth} dargestellt. \begin{table}[h!] \caption{Der für die Koffeinsynthese definierte Hypergraph in Gurobi mit den für jede Kante definierten Wahrscheinlichkeiten.} \begin{center} \scriptsize \begin{tabular}{ |c|c|c|c|c| } \label{hypergraphcaffsynth} Index & Edukt & Produkt & $n_{\text{7-Methylxanthin},e}$ & $n_{\text{3,7-Dimethylxanthin},e}$\\ \hline 1 & & Xanthin & 0.12 & 0.25\\ 2 &Xanthin & 1-Methylxanthin & 0.18 & 0.16\\ 3 &Xanthin & 3-Methylxanthin & 0.12 & 0.2\\ 4 &Xanthin & 7-Methylxanthin & 0.25 & 0.21\\ 5 &1-Methylxanthin & Theophyllin & 0.18 & 0.16\\ 6 &1-Methylxanthin & Paraxanthin & 0.3 & 0.18\\ 7 &3-Methylxanthin & Theophyllin & 0.18 & 0.16\\ 8 &3-Methylxanthin & Theobromin & 0.1 & 0.31\\ 9 &7-Methylxanthin & Paraxanthin & 0.3 & 0.18\\ 10 &7-Methylxanthin & Theobromin & 0.1 & 0.31\\ 11 &Theophyllin & Koffein & 0.3 & 0.29\\ 12 &Paraxanthin & Koffein & 0.3 & 0.29\\ 13 &Theobromin & Koffein & 0.3 & 0.29\\ 14 &Koffein & & 0.0 & 0.0\\ \hline \end{tabular} \end{center} \end{table} Die beste gefundene Lösung ist für die Kanten 1, 4, 10, 13 und 14 jeweils ein flow $x_{e} = 1$ und alle anderen Kanten ein flow $x_{e} = 0$ (Tabelle \ref{ilpsolcaff}). Dies entspricht der natürlich auftretenden Lösung. \begin{table}[h!] \caption{Lösung von Gurobi für die Koffeinsynthese.} \begin{center} \scriptsize \begin{tabular}{ |c|c|c|c| } \label{ilpsolcaff} Index & Edukt & Produkt & flow $x_{e}$\\ \hline 1 & & Xanthin & 1\\ 2 &Xanthin & 1-Methylxanthin & 0\\ 3 &Xanthin & 3-Methylxanthin & 0\\ 4 &Xanthin & 7-Methylxanthin & 1\\ 5 &1-Methylxanthin & Theophyllin & 0\\ 6 &1-Methylxanthin & Paraxanthin & 0\\ 7 &3-Methylxanthin & Theophyllin & 0\\ 8 &3-Methylxanthin & Theobromin & 0\\ 9 &7-Methylxanthin & Paraxanthin & 0\\ 10 &7-Methylxanthin & Theobromin & 1\\ 11 &Theophyllin & Koffein & 0\\ 12 &Paraxanthin & Koffein & 0\\ 13 &Theobromin & Koffein & 1\\ 14 &Koffein & & 1\\ \hline \end{tabular} \end{center} \end{table} \FloatBarrier \subsection{Vanillin und 3,4-Dihydroxyaldehyd} Der Hypergraph für die Vanillinsynthese wird in Tabelle \ref{hypergraphvanilsynth} dargestellt. \begin{table}[h!] \caption{Der für die Vanillinsynthese definierte Hypergraph in Gurobi mit den für jede Kante definierten Wahrscheinlichkeiten.} \begin{center} \scriptsize \begin{tabular}{ |c|c|c|c|c| } \label{hypergraphvanilsynth} Index & Edukt & Produkt & $n_{\text{p-Coumarsäure},e}$ & $n_{\text{4-Hydroxybenzaldehyd},e}$\\ \hline 1 & & Zimtsäure & 0.09 & 0.17\\ 2 &Zimtsäure & p-Coumarsäure & 0.38 & 0.44\\ 3 &Zimtsäure & m-Coumarsäure & 0.12 & 0.12\\ 4 &Zimtsäure & Benzaldehyd & 0.09 & 0.18\\ 5 &p-Coumarsäure & Kaffeesäure & 0.34 & 0.15\\ 6 &p-Coumarsäure & 4-Hydroxybenzaldehyd & 0.34 & 0.26\\ 7 &m-Coumarsäure & Kaffeesäure & 0.34 & 0.15\\ 8 &m-Coumarsäure & 3-Hydroxybenzaldehyd & 0.02 & 0.03\\ 9 &Benzaldehyd & 3-Hydroxybenzaldehyd & 0.02 & 0.03\\ 10 &Benzaldehyd & 4-Hydroxybenzaldehyd & 0.1 & 0.31\\ 11 &Kaffeesäure & 3,4-Dihydroxybenzaldehyd & 0.17 & 0.17\\ 12 &3-Hydroxybenzaldehyd & 3,4-Dihydroxybenzaldehyd & 0.17 & 0.17\\ 13 &4-Hydroxybenzaldehyd & 3,4-Dihydroxybenzaldehyd & 0.17 & 0.17\\ 14 &3,4-Dihydroxybenzaldehyd & & 0.0 & 0.0\\ \hline \end{tabular} \end{center} \end{table} Die beste gefundene Lösung ist für die Kanten 1, 2, 6, 13 und 14 jeweils ein flow $x_{e} = 1$ und alle anderen Kanten ein flow $x_{e} = 0$ (Tabelle \ref{ilpsolvan}). Dies entspricht der natürlich auftretenden Lösung. \begin{table}[h!] \caption{Lösung von Gurobi für die Vanillinsynthese.} \begin{center} \scriptsize \begin{tabular}{ |c|c|c|c| } \label{ilpsolvan} Index & Edukt & Produkt & flow $x_{e}$\\ \hline 1 & & Zimtsäure & 1\\ 2 &Zimtsäure & p-Coumarsäure & 1\\ 3 &Zimtsäure & m-Coumarsäure & 0\\ 4 &Zimtsäure & Benzaldehyd & 0\\ 5 &p-Coumarsäure & Kaffeesäure & 0\\ 6 &p-Coumarsäure & 4-Hydroxybenzaldehyd & 1\\ 7 &m-Coumarsäure & Kaffeesäure & 0\\ 8 &m-Coumarsäure & 3-Hydroxybenzaldehyd & 0\\ 9 &Benzaldehyd & 3-Hydroxybenzaldehyd & 0\\ 10 &Benzaldehyd & 4-Hydroxybenzaldehyd & 0\\ 11 &Kaffeesäure & 3,4-Dihydroxybenzaldehyd & 0\\ 12 &3-Hydroxybenzaldehyd & 3,4-Dihydroxybenzaldehyd & 0\\ 13 &4-Hydroxybenzaldehyd & 3,4-Dihydroxybenzaldehyd & 1\\ 14 &3,4-Dihydroxybenzaldehyd & & 1\\ \hline \end{tabular} \end{center} \end{table} Dabei ist jedoch anzumerken, dass die Lösung des ILP für das erste Moleküle 4-Hydroxybenzaldehyd und für das zweite Molekül p-Coumarsäure gewählt hat. Die Zuordnung erfolgte also genau anders herum. Zudem ist davon auszugehen, dass die zweitbeste Lösung mit erst Kaffeesäure und dann p-Coumarsäure genau den selben Zielwert aufweist. \FloatBarrier \section{Fazit: Ein-Molekül-Reaktionen} Bei den bisher betrachteten Ein-Molekül-Reaktionen wurde immer nur ein festgelegtes Startmolekül in das System gegeben. Dieses Startmolekül wurde durch Transformationsregeln verändert, wobei nach dem natürlichen Vorbild genau ein Zielmolekül synthetisiert werden sollte. In dem Koffein Beispiel gab es keine Enantiomere oder E/Z-Isomere. In dem Vanille Beispiel gab es keine Enantiomere. Es gab mit Zimtsäure, Coumarsäure und Kaffeesäure zwar E/Z-Isomere, welche aber durch die natürliche Vorgabe von E-Zimtsäure bereits vorgegeben waren. Außerdem entstand keine neue Doppelbindung mit einem E/Z-Isomer durch Graphtransformation. Der Schluss von dem Graphen auf eine Geometrie war deswegen immer eindeutig. \newline Die Berechnung der NMR-Spektren war trotzdem fehleranfällig. Dabei konnten bestimmte Molekülveränderungen besser in dem NMR-Berechnungstool abgebildet werden. Eine Methylierung erzeugte nur eine geringe Veränderung im Xanthingrundgerüst, weshalb die Zuordnung vor allem in dem Signal der Methylgruppe lag. Die Hydroxylierung und Oxidierung zum Aldehyd resultierten in deutlichen Veränderungen der Phenylgruppe, weshalb das Grundgerüst selber verändert wurde. \newline Für die Synthesen waren immer genau drei Reaktionen nötig, welche mittels des ILP nur bezüglich ihrer Reihenfolge variiert werden konnten. Das Problem ließe sich daher ebenfalls durch die einfache Addition zweier Moleküle entlang eines ausgewählten Pfades lösen. Der Raum möglicher Lösungen ist mit 6 pro Hypergraph stark begrenzt. Trotzdem eigneten sich beide Beispiele mit der natürlich auftretenden Selektivität der Reihenfolge, der größtenteils guten Differenzierbarkeit der NMR-Spektren und Kompaktheit der Synthese als Graphtransformationsregeln als Demonstration des Konzeptes. Es ist anzumerken, dass die erwartete Synthese über die Auswahl der Referenzspektren festgelegt wurde. Diese Synthese wurde über NMR-Spektren abstrahiert. Es lässt sich so sagen, dass die NMR-Spektren in dieser Analyse die Struktur hinreichend codieren, aber nicht eine unbekannte Synthese vollständig hergeleitet wurde. \newline Für ausgewählte Literaturspektren von nicht in der natürlichen Synthese erwarteten Molekülen wurden ebenfalls der $SI$ mit den berechneten Spektren ermittelt. Dabei konnte ebenfalls eine erhöhte Korrelation mit dem zugehörigen berechneten Spektrum festgestellt werden. Das ILP konnte dabei den neuen Syntheseplan ebenfalls zuverlässig zuordnen. Die Auswahl eines Syntheseplans basierend auf der Ähnlichkeit von NMR-Spektren ist daher gerechtfertigt und stellt eine Verbesserung gegenüber der Auswahl ohne selbige dar. \section{Multi-Molekül-Reaktionen}\label{chmultmore} Die bisherigen Ein-Molekül-Reaktionen lassen sich durch ein ILP lösen, erfordern die Verwendung jedoch nicht. Um die Grenzen des ILP weiter auszureizen, wird ein komplexerer Hypergraph entwickelt. Dieser Hypergraph wird mit wenigen, einfachen Transformationsregeln und mit einem Ausgangsmolekül erzeugt. Das Ziel ist ein hochvernetztes Reaktionsnetzwerk. \section{Olefinmetathese und Diels-Alder} Die zwei Reaktionen Olefinmetathese in Kombination mit der Diels-Alder-Reaktion erlaubt mit dem richtigen Ausgangsmolekülen die Entwicklung eines hochvernetzten Reaktionsnetzwerks und unterliegenden Hypergraphen. \subsection{Olefinmetathese} Die Olefinmetathese ist eine vergleichbar neuere Reaktion \cite{olefinmetathese}. Mit einem organischen Ruthenium-Katalysator werden dabei die Reste zweier Doppelbindungen rekombiniert. Der Ruthenium-Katalysator hat dabei eine offene Bindungsstelle an die eine Hälfte der C-C-Doppelbindung als Carben bindet (A und E). Ein Olefin (B und F) bindet an diesen Komplex (C und G). Aus diesem Komplex wird ein Olefin abgespalten, welches entweder das eingehende Olefin oder eine beliebige Kombination beider Olefine (D) ist (Abbildung \ref{olefinmetathese}). \newline In der Forschung wurde dieser katalytische Prozess zu immer gezielteren Ausbeuten durch spezialisierte Katalysatoren weiterentwickelt. Dabei wurde neben der generellen Selektivität der Produkte auch eine gezielte Selektivität hin zu E-Isomeren erzielt. Für die Anwendung auf ein Reaktionsnetzwerk ist eine Produktspezialisierung jedoch weniger interessant. Die E/Z-Selektivität ist dabei insbesondere im Hinblick auf das chemische Molekül als Graphen nicht umsetzbar, wenngleich sie bei der Geometrierealisierung berücksichtigt werden kann. Die Möglichkeiten der Olefinmetathese sind jedoch begrenzt, wenn es um die Queervernetzung geht. Bei einer Olefinmetathese bleibt die Gesamtzahl an Einfach und Doppelbindungen gleich. \begin{figure} [h!] \includegraphics[width=\textwidth]{olefinmetathese.png} \caption{Katalysezyklus der Olefinmetathese \cite{olefinmetathese}. } \label{olefinmetathese} \end{figure} \subsection{Diels-Alder} Die Diels-Alder-Reaktion ist eine ohne Katalysator natürlich auftretende Reaktion, bei der zwei Doppelbindungen enthaltenden Moleküle verbunden werden. Dabei wird eine Vinylgruppe und zwei mit einer Einfachbindung verbundenen Vinylgruppen zu einem Cyclohexen Ring verbunden (Abbildung \ref{dielsalder}). Wenn diese Bedingungen in einer Lösung vorhanden sind, kommt es zu dieser Kombination, mit nur geringer Rückreaktion \cite{dielsalder}. \begin{figure} [h!] \begin{center} \includegraphics[width=0.6\textwidth]{dielsalder.png} \caption{Die Diels-Alder-Reaktion \cite{dielsalder}. } \label{dielsalder} \end{center} \end{figure} Die Diels-Alder-Reaktion ist einfach und sorgt dafür, dass Doppelbindungen zu Einfachbindungen werden. In der Kombination mit der Olefinmetathese werden so viele verschiedene Produkte möglich. \subsection{Butadien} Das gewählte Startmolekül ist Butadien. Das Vorliegen von zwei Doppelbindungen erlaubt es der Olefinmetathese Kohlenstoffketten mit alternierenden Einfach und Doppelbindungen aufzubauen. Diese können wiederum mit der Diels-Alder-Reaktion untereinander kombiniert werden und so zu einer großen Menge an unterschiedlichen Molekülen führen, welche alle untereinander kombinierbar sind. \section{Graphtransformation: Multi-Molekül-Reaktionen} Bei der Definition der Transformationsregeln für Multi-Molekül-Reaktionen müssen auf der linken und rechten Seite, sowie im Kontext die gleichen Atome vorliegen. Bei den meisten mehrstufigen, chemischen Synthesen wird ein zentrales Molekül betrachtet. Es kommt zwar vor, dass für eine Reaktion zwei Moleküle zusammen geführt werden, aber Teile des Moleküls, die in Reaktionen abgespalten werden, werden in der restlichen Synthese nicht weiter erwähnt. Bei einer Multi-Molekül-Reaktion sollen diese abgespaltenen Moleküle potenziell wieder miteinander reagieren. Es ist also wichtig, dass eine Transformationsregel in Hinsicht auf Multi-Molekül-Reaktionen alle Atome vollständig weiterführt. \subsection{Olefinmetathese und Diels-Alder} In dem Graphtransformationstool wurde mit dem Molekül Butadien (Abbildung \ref{butadienmod}) gestartet. Auf dieses wurden wiederholt die Olefinmetathese (Abbildung \ref{olmetmod}) und die Diels-Alder Reaktion (Abbildung \ref{dialdmod}) angewendet. \begin{figure} [p] \begin{center} \includegraphics[width=0.5\textwidth]{butadienmod.png} \caption{Butadien in dem Graphtransformationstool.} \label{butadienmod} \end{center} \end{figure} \begin{figure} [p] \includegraphics[width=\textwidth]{olmetmod.png} \caption{Olefinmetathese in dem Graphtransformationstool.} \label{olmetmod} \end{figure} \begin{figure} [p] \includegraphics[width=\textwidth]{dialmod.png} \caption{Diels-Alder-Reaktion in dem Graphtransformationstool.} \label{dialdmod} \end{figure} Bei dieser Reaktion ist es anders als bei den bisherigen Reaktionen möglich, dass der Ergebnisraum unbegrenzt wachsen kann. Darum sind Beschränkungen des Ergebnisraums nötig. Die einfachste Beschränkung der Ergebnisse ist eine Obergrenze für die Anzahl an Atomen. Damit Butadien als Startmolekül diese Begrenzung erfüllt, müssen mindestens 4 Kohlenstoffe zugelassen werden. Damit die Diels-Alder-Reaktion stattfinden kann müssen 6 Kohlenstoffe erlaubt sein. Die Regeln zeigen einen diversen Raum an gefundenen Molekülen bei 10 Kohlenstoffen. Da dies die Computerleistung stark beansprucht, muss der Raum weiter eingeschränkt werden. Die Olefinmetathese ist bekannt dafür, Ringe mit großer Ringspannung aufzubrechen \cite{olefinmetathese}. Besonders strapaziert sind kleine Ringe mit der Größe 3 und 4, sowie große Ringe. Ein Kohlenstoff Ring der Größe 6 ist am stabilsten \cite{ringstrain}. Da die Olefinmetathese mit hoher Selektivität Ringe mit hoher Spannung aufbricht, ist davon auszugehen, dass zu große und zu kleine Ringe nicht Teil des Ergebnisraums sind. Insbesondere ist nur der 6 Ring als die stabilste Ringgröße zu erwarten. Deshalb wurden für jedes gefundene Molekül alle sehnenlosen Ringe gefunden und alle Ringe, die nicht die Größe 6 hatten aus dem Ergebnis Raum entfernt. Eine weitere Quelle der Ringspannung sind überlappende Ringe. Norbornan wird zum Beispiel aus drei überlagernden Ringen gebildet. Da auch in diesen Molekülen hohe Ringspannung vorliegt, wurden alle Moleküle mit Ringen, die mehr als zwei gemeinsame Atome haben aus dem Ergebnisraum entfernt. Die letzte Begrenzung für Ringspannung ist, dass kein Atom in vier Ringen sein kann. \newline Mit der Beschränkung wurden durch die Regelanwendung 19 Moleküle mit etwa 200 Transformationen aus Butadien entwickelt. Die Moleküle und ihre SMILES-Strings werden im Folgenden aufgeführt. \begin{itemize} \item Butadien: C=CC=C \item Ethylen: C=C \item Hexatrien: C=CC=CC=C oder C=C/C=C/C=C \item 4-Vinyl-Cyclohex-1-en: C1CCC(C=C)CC=1 \item 4-Vinyl-Octa-1,7-Dien: C(CCC(C=C)CC=C)=C \item Octatetraen: C=CC=CC=CC=C oder C=C/C=CC=C/C=C \item 4-Butadienyl-Cyclohex-1-en: C1CCC(C=CC=C)CC=1 \item 3-Vinyl-Cyclohex1en: C=CC1C=CCCC1 \item Cyclohexen: C1CCCCC=1 \item 4,5-Divinyl-Cyclohex-1-en: C=CC1CC=CCC1C=C \item 3,4-Divinyl-Cyclohex-1-en: C1CCC(C=C)C(C=C)C=1 \item 3,5-Divinyl-Cyclohex-1-en: C(C1CC(C=C)C=CC1)=C \item Benzol: C1C=CC=CC=1 \item 5-Butan-1-enyl-Cyclo-1,3-dien: C(CCC1C=CC=CC1)=C \item 3-Butadienyl-Cyclohex-1-en: C1C(C=CC=C)CCCC=1 \item 3-Vinyl-Octa-1,7-Dien: C=CC(C=C)CCCC=C \item Okta-1,7-dien: C=CCCCCC=C \item 3,6-Divinyl-Cyclohex1en: C=CC1C=CC(C=C)CC1 \item 1,2,3,4,4a,5,6,7,8,8a-Octahydronaphthalin: C1CC2CCCCC2CC=1 \item 1,4,4a,8a-Tetrahydronaphthalin: C1CC2C=CC=CC2CC=1 \end{itemize} Wenn Butadien mit einem generischen Olefinmetathese-Katalysator zum reagieren gebracht wird ist zu erwarten, dass alle Moleküle in verschiedenen Enantiomeren und E/Z-Isomeren synthetisiert werden. Daher ist es ohne das tatsächliche Durchführen dieses Experiments schwierig einen Reaktionspfad herauszusuchen, der von dem ILP bestätigt werden soll. Die Synthese von Benzol (Abbildung \ref{benzsynthmod}) stellt eine gute Option dar, da die hohe Symmetrie im Vergleich zu den anderen Spektren zu einem sehr markanten Spektrum führt. \begin{figure} [p] \includegraphics[width=\textwidth]{benzsynthmod.png} \caption{Benzolsynthese in dem Graphtransformationstool. Die Syntheseschritte lauten erst Olefinmetathese von zwei Butadien zu Hexatrien und Ethylen, dann Olefinmetathese von Hexatrien und Butdien zu Ethylen und Octatetraen und zuletzt eine Ringschluss Olefinmetathese von Octatetraen zu Benzol und Ethylen.} \label{benzsynthmod} \end{figure} Die direkteste Synthese stellen dabei die Olefinmetathese von zwei Butadien zu Hexatrien und Ethylen, die Olefinmetathese von Hexatrien mit Butadien zu Octatetraen und eine Ringschluss Olefinmetathese von Octatretraen zu Benzol und Ethylen dar. Es werden drei Butadien eingesetzt, um mit drei Olefinmetathesen ein Benzol und drei Ethylen zu erhalten. \section{NMR: Multi-Molekül-Reaktionen} \subsection{Olefinmetathese und Diels-Alder} Die generierten Graphen weisen anders als bei den vorherigen Beispielen verschiedene Enantiomere und E/Z-Isomere auf. Es ist daher notwendig diese zu unterscheiden. Für die Moleküle Hexatrien und Octatetraen wurden gezielt die Konfigurationen gewählt, bei der alle Doppelbindungen in der E-Konfiguration sind. Bei den Verbindungen mit Cyclohexanring und angehängten Vinylgruppen, wurde die Position so gewählt, dass die Reste äquatorial zum Ring stehen. Bei den beiden Hydronaphtalinen wurde die Konfiguration gewählt, bei der beide Ringe in der äquatorialen Konfiguration liegen. Verallgemeinert wurde für jedes Molekül die Konfiguration mit der größten Symmetrie und dadurch mit den meisten chemisch-äquivalenten Kohlenstoffen gewählt. \section{NMR $SI$: Multi-Molekül-Reaktionen} In der Literatur werden NMR-Spektren bevorzugt für einen vollständig aufgereinigten Stoff angegeben. Daher wurden die NMR-Spektren für jedes Molekül mit dem kombinierten NMR-Spektrum aller Produkte der ausgewählten Reaktionen verglichen. \subsection{Olefinmetathese und Diels-Alder} Die Referenzspektren der Benzolsynthese sind im Folgenden aufgeführt. \begin{itemize} \item Hexatrien und Ethylen \item Octatetraen und Ethylen \item Benzol und Ethylen \end{itemize} Für den Vergleich der NMR-Spektren wurden die gleichen Bingrößen $b$ wie für die bisherigen Ein-Molekül-Reaktionen (Gleichung \ref{binsize}) verwendet. Der Ähnlichkeitsindex $SI$ wurde als Mittelwert der $SI_{b}$ für alle Bingrößen berechnet. Die finale Wahrscheinlichkeit entspricht dem $SI$ normalisiert für jedes Referenzspektrum (Tabelle \ref{nmrsimbenz}). \begin{table} [h!] \caption{Resultierende Wahrscheinlichkeiten für den Vergleich aller aus Butadien Entwickelten Moleküle mit Hexatrien und Ethylen, Octatetraen und Ethylen und Benzol und Etylen.} \begin{center} \scriptsize \begin{tabular}{ |c|c|c|c| } \label{nmrsimbenz} Molekül & Hexatrien und Ethylen & Octatetraen und Ethylen & Benzol und Etylen \\ \hline Butadien & 0.11 & 0.11 & 0.0\\ Ethylen & 0.18 & 0.15 & 0.17\\ Hexatrien & 0.31 & 0.21 & 0.0\\ 4-Vinyl-Cyclohex-1-en & 0.0 & 0.0 & 0.0\\ 4-Vinyl-Octa-1,7-Dien & 0.04 & 0.04 & 0.0\\ Octatetraen & 0.21 & 0.31 & 0.0\\ 4-Butadienyl-Cyclohex-1-en & 0.02 & 0.02 & 0.03\\ 3-Vinyl-Cyclohex1en & 0.0 & 0.0 & 0.0\\ Cyclohexen & 0.0 & 0.0 & 0.0\\ 4,5-Divinyl-Cyclohex-1-en & 0.0 & 0.0 & 0.01\\ 3,4-Divinyl-Cyclohex-1-en & 0.0 & 0.0 & 0.0\\ 3,5-Divinyl-Cyclohex-1-en & 0.0 & 0.0 & 0.16\\ Benzol & 0.0 & 0.0 & 0.51\\ 5-Butan-1-enyl-Cyclo-1,3-dien & 0.0 & 0.0 & 0.09\\ 3-Butadienyl-Cyclohex-1-en & 0.03 & 0.03 & 0.0\\ 3-Vinyl-Octa-1,7-Dien & 0.0 & 0.0 & 0.0\\ Okta-1,7-dien & 0.0 & 0.0 & 0.0\\ 3,6-Divinyl-Cyclohex1en & 0.06 & 0.07 & 0.0\\ 1,2,3,4,4a,5,6,7,8,8a-Octahydronaphthalin & 0.0 & 0.0 & 0.0\\ 1,4,4a,8a-Tetrahydronaphthalin & 0.03 & 0.04 & 0.03\\ \hline \end{tabular} \end{center} \end{table} Das Referenzspekrtrum mit Hexatrien und Ethylen weist eine große Ähnlichkeit mit den beiden tatsächlichen Molekülen auf. Die Ähnlichkeit zum Referenzspektrum liegt für Ethylen bei 0.18 und für Hexatrien bei 0.31. Die Moleküle Butadien mit 0.11 und Octatetraen mit 0.21 weisen ebenfalls eine hohe Ähnlichkeit zum Referenzspektrum auf. Da es sich bei Octatetraen, Hexatrien und Butadien um Kohlenstoffketten mit alternierenden Einfach- und Doppelbindungen handelt, ist diese Ähnlichkeit zu erwarten. Alle drei haben ein Signal bei etwa 122 ppm und ein Signal bei etwa 145 ppm. Die längeren Ketten haben zusätzlich ein Signal bei etwa 140 ppm, weshalb der Überlapp bei Hexatrien und Octatetraen etwas größer ist. Das Referenzspektrum mit Octatetraen und Ethylen weist eine große Ähnlichkeit mit beiden tatsächlichen Molekülen auf. Die Ähnlichkeit zum Referenzspektrum liegt für Ethylen bei 0.15 und für Octatetraen bei 0.31. Die Moleküle Butadien mit 0.11 und Hexatrien mit 0.21 weisen ebenfalls eine hohe Ähnlichkeit zum Referenzspektrum auf. Die Gründe dafür sind die gleichen wie beim ersten Referenzsspektrum, da der Ähnlichkeitsindex $SI$ zweier Moleküle symmetrisch ist. Das Referenzspektrum mit Benzol und Etyhlen hat die größte Ähnlichkeit mit den beiden tatsächlich im Spektrum enthaltenden Molekülen. Ethylen hat einen $SI$ von 0.17 und Benzol von 0.51. Es gibt aber auch einige weitere Stoffe wie 5-Butan-1-enyl-Cyclo-1,3-dien und 3,5-Divinyl-Cyclohex-1-en, welche mit 0.09 und 0.16 ebenfalls eine deutliche Ähnlichkeit aufweisen. Diese Ähnlichkeit liegt daran, dass der eine Benzolpeak bei etwa 132.96 ppm mit manchen Peaks aus den anderen beiden Spektren überlagert. Beide Moleküle haben einen Peak bei 132.93 und 133.29 ppm respektive. Da der Peak von 5-Butan-1-enyl-Cyclo-1,3-dien näher liegt ist die Ähnlichkeit größer. Diese Vergleiche demonstrieren die Limitationen eines simpleren Ansatz wie dem Binning mit Crosskorrelation. Ähnliche Moleküle mit vergleichbaren chemischen Verschiebungen können einen hohen Ähnlichkeitsindex $SI$ aufweisen, auch wenn die Anzahl der Peaks nicht zum Molekül passt. Die tiefer gehenden chemischen Informationen können so verloren gehen. \section{ILP: Multi-Molekül-Reaktionen} Für die Definition des ILP Problems wurde der Hypergraph über alle mit dem Graphtransformationstool entwickelten Reaktionen definiert. Dabei wurde eine Kante, die den Ausgangsstoff in das System einführen kann und für jedes andere Molekül eine Kante, die es aus dem System entfernt, hinzugefügt. Die Kanten, bei der die Edukte und Produkte identisch zueinander sind, wurden aufgrund ihrer Redundanz entfernt. \subsection{Olefinmetathese und Diels-Alder} Für das ILP-Model wurde ein Ansatz entwickelt, welcher wie in den vorangegangenen Beispielen die Molekülwahrscheinlichkeiten den Kanten zuordnet. Zudem wurde ein Model entwickelt, welches zuerst die Wahrscheinlichsten Moleküle ermittelt und deren Produktion maximiert. \subsubsection{Kantenwahrscheinlichkeit} Bei diesem Model werden die Wahrscheinlichkeiten als konstante Werte $p_{i,e}$ den Kanten $e$ des Hypergraphen zugeordnet. In der Olefinmetathese fließt jede Kante entweder in keine, ein oder zwei Edukte. Ausgehende Kanten, die ins Nichts fließen, sollten in der Zielfunktion nicht berücksichtigt werden. Die chemischen Reaktionen sind also nur die Kanten mit ein oder zwei Produkten. Da den Produkten Wahrscheinlichkeitswerte zugewiesen sind liegt nahe, dass die Gesamtwahrscheinlichkeit dem Produkt beider Einzelwahrscheinlichkeiten entspricht. Da die Wahrscheinlichkeiten kleiner als 1 sind, wird dies bei Kanten mit zwei Edukten zu kleineren Wahrscheinlichkeiten führen als bei Kanten mit nur einem Edukt. Ein weiterer Ansatz ist es die Gesamtwahrscheinlichkeit als deb Mittelwert beider Produkte zu wählen. Dieser Ansatz sorgt dafür, dass eine Kante mit einem sehr unwahrscheinlichen und einem sehr wahrscheinlichen Produkt trotzdem eine mittlere Gesamtwahrscheinlichkeit hat. Zur Demonstration der beiden Ansätze zur Berechnung der Kantenwahrscheinlichkeit werden die Kantenwahrescheinlichkeiten von zwei Kanten im Hypergraphen für das dritte Referenzspektrum bestimmt. \begin{itemize} \item 42: Octatetraen $\rightarrow$ Benzol + Ethylen \item 44: 4-Butadienyl-Cyclohex-1-en $\rightarrow$ 5-Butan-1-enyl-Cyclo-1,3-dien \end{itemize} Die Wahrscheinlichkeiten für Kanten mit nur einem Produkt entsprechen der Wahrscheinlichkeit für dieses Produkt. Demnach hat die Kante 44 für das dritte Referenzspektrum die Wahrscheinlichkeit $p_{3,44}$ für Durchschnitt und Produkt. Kante 42 hingegen hat die Wahrscheinlichkeit 0.0876 für das Produkt (Gleichung \ref{pprod}) und 0.34 für den Durchschnitt (Gleichung \ref{pavg}). \begin{equation} \label{pprod} p_{3,42} = 0.51 \times 0.17 = 0.0876 \end{equation} \begin{equation} \label{pavg} p_{3,42} = \frac {0.51 + 0.17} {2} = 0.34 \end{equation} Für jede Kante werden fünf ganzzahlige Variablen $x_{e}$, $b_{e}$, $n_{1,e}$, $n_{2,e}$ und $n_{3,e}$ definiert. Die Variable $b_{e}$ ist dabei eine Wahrheitsvariable, die definiert ob eine Kante verwendet wird. Die Variable $x_{e}$ beschreibt wie oft die selbe Kante verwendet wird. Die Variablen $n_{1,e}$, $n_{2,e}$ und $n_{3,e}$ sind ebenfalls Wahrheitsvariablen und legen fest, welche drei Reaktionen den Referenzspektren entsprechen und deswegen zu der Zielfunktion dazuzählen. Dabei ist $n_{1,e}$ die Zuordnungsvariable zu dem ersten Referenzspektrum mit Hexatrien und Ethylen, $n_{2,e}$ die Zuordnungsvariable zu dem zweiten Referenzspektrum mit Octatetraen und Ethylen und $n_{3,e}$ die Zuordnungsvariable zum dritten Referenzspektrum mit Benzol und Ethylen. Da das ILP die Syntheseschritte durch den Vergleich der NMR-Spektren findet, werden alle synthetisierten Moleküle bestimmt. Um dahingehend die Lösung nicht einzuschränken ist kein Zielmolekül festgelegt. Die Variablen unterliegen den Beschränkungen, dass: \begin{itemize} \item $x_{e} \geq b_{e}$ \item $b_{e} \geq n_{1,e}$ \item $b_{e} \geq n_{2,e}$ \item $b_{e} \geq n_{3,e}$ \item s: ([], [Startmolekül]), $x_{s} = 3$ \item $\forall v \in V$ inflow = outflow \item $\forall i,j \in E$ i: head $\rightarrow$ tail und j: tail $\rightarrow$ $\sum b_{i} + b_{j} \leq 1$ \item $n_{1,e} + n_{2,e} + n_{3,e} \leq 1$ \item $\sum_{e} n_{1,e} = 1$ \item $\sum_{e} n_{2,e} = 1$ \item $\sum_{e} n_{3,e} = 1$ \item wenn die größte Wahrscheinlichkeit einer nicht ausgehenden Kante $p_{max,e} = 0$ ist, wird diese stattdessen auf $p_{max,e} = 0.0001$ gesetzt um einen Wertfehler zu verhindern \end{itemize} Unter diesen Beschränkungen soll ein Pfad durch den Hypergraphen bestimmt werden, welcher die drei wahrscheinlichsten Kanten für die drei Referenzspektren auswählt und dabei möglichst wenige unwahrscheinliche Kanten beinhaltet. Es ist nicht üblich, Kanten mit einer Wahrscheinlichkeit $p_{max,e} = 0.00$ auf $p_{max,e} = 0.0001$ zu setzen, da die Wahl des Wertes das Ergebnis beeinflussen kann. In dem ILP-Model ist es sehr restriktiv Kanten mit einer Wahrscheinlichkeit von $p_{max,e} = 0.00$ komplett zu verbieten, da die Differenz zu Kanten mit Wahrscheinlichkeit $p_{max,e} = 0.01$ nicht groß ist, jedoch $p_{max,e} = 0.00$ ausgeschlossen wird. Ein Ansatz der diese Entscheidung umgeht ist die Auswahl von Zielknoten (Kapitel \ref{nodechoice}). Die ausgehenden Kanten sind nicht in die Zielfunktion eingeflossen (Gleichung \ref{objectivebutadien}). \begin{equation} \label{objectivebutadien} max \: z = \sum_{e} 1000 \times (p_{1,e} \times n_{1,e} + p_{2,e} \times n_{2,e} + p_{3,e} \times n_{3,e}) \newline - \frac{x_{e}}{ max(p_{1,e}, p_{2,e}, p_{3,e})} \end{equation} Für beide Kantenwahrscheinlichkeitsdefinitionen wurden jeweils die 10 besten Lösungen gesucht. \subsubsection{Knoten Auswahl}\label{nodechoice} Ein Ansatz, der unabhängig von der Auswahl der Kantenwahrscheinlichkeitsberechnung ist, ist durch die NMR-Spektren Zielmoleküle zu wählen und einen Pfad zu finden, der diese synthetisiert. Dafür wurde jedem Molekül $v$ die Wahrscheinlichkeit $p_{i,v}$ direkt zugeordnet. Für jede Kante des Hypergraphen werden zwei Variablen $x_{e}$ und $b_{e}$ definiert. Die Variable $b_{e}$ ist dabei eine Wahrheitsvariable, die definiert ob eine Kante verwendet wird. Die Variable $x_{e}$ beschreibt wie oft die selbe Kante verwendet wird. Für jeden Knoten des Hypergraphen werden drei Variablen $n_{1,e}$, $n_{2,e}$ und $n_{3,e}$ definiert. Diese drei Variablen sind ebenfalls Wahrheitsvariablen und legen fest, welche drei Moleküle den drei Referenzspektren entsprechen. Dabei ist $n_{1,e}$ die Zuordnungsvariable zu dem ersten Referenz NMR-Spektrum mit Hexatrien und Ethylen, $n_{2,e}$ die Zuordnungsvariable zu dem zweiten Referenz NMR-Spektrum mit Octatetraen und Ethylen und $n_{3,e}$ die Zuordnungsvariable zum dritten Referenz NMR-Spektrum mit Benzol und Ethylen. Die Variablen unterliegen den Beschränkungen, dass: \begin{itemize} \item $x_{e} \geq b_{e}$ \item Hilfsvariable: $m_{v} = \sum_{e} b_{e}: v \in head(b_{e})$ \item $m_{v} \geq n_{1,v}$ \item $m_{v} \geq n_{2,v}$ \item $m_{v} \geq n_{3,v}$ \item s: ([], [Startmolekül]), $x_{s} = 3$ \item $\forall v \in V$ inflow = outflow \item $\forall i,j \in E$ i: head $\rightarrow$ tail und j: tail $\rightarrow$ $\sum b_{i} + b_{j} \leq 1$ \item $n_{1,v} + n_{2,v} + n_{3,v} \leq 1$ \item $\sum_{v} n_{1,v} = 1$ \item $\sum_{v} n_{2,v} = 1$ \item $\sum_{v} n_{3,v} = 1$ \end{itemize} Dabei soll die Wahrscheinlichkeit der ausgewählten Moleküle maximiert werden (Gleichung \ref{objspecmain}). Nachdem drei Zielmoleküle ausgewählt wurden wird der Fluss in diese Moleküle als Bedingung gesetzt. \begin{equation} \label{objspecmain} max \: z_{1} = \sum_{v}\sum_{i} n_{i,v} \times p_{i,v} \end{equation} Alle Lösungen, die positiven Fluss $b_{e} = 1$ in die Zielmoleküle haben, werden danach sortiert, welche Lösung die geringste Anzahl an Reaktionen hat (Gleichung \ref{objspecsec}). Wenn mehrere Lösungen die gleiche Anzahl an Reaktionen haben, wird die ausgewählt, die den geringsten flow hat (Gleichung \ref{objspectrd}). \begin{equation} \label{objspecsec} max \: z_{2} = \sum_{e} b_{e} \end{equation} \begin{equation} \label{objspectrd} max \: z_{3} = \sum_{e} x_{e} \end{equation} Es wurden die beiden besten Lösungen für das ILP bestimmt. \begin{table}[h!] \caption{Lösung 1 für die Knotenauswahl beim ILP vom Butadien Beispiel mit $z_{3} = 3$.} \begin{center} \scriptsize \begin{tabular}{ |c|c|c|c| } \label{ilpspec1} Flow & Kante & Edukt & Produkt \\ \hline 3 & 1 & Ethylen & \\ 1 & 4 & 2 Butadien & Hexatrien + Ethylen\\ 1 & 23 & Butadien + Hexatrien & Octatetraen + Ethylen\\ 1 & 41 & Benzol & \\ 1 & 42 & Octatetraen & Benzol + Ethylen\\ 3 & 213 & & Butadien\\ \hline \end{tabular} \end{center} \end{table} \begin{table}[h!] \caption{Lösung 2 für die Knotenauswahl beim ILP vom Butadien Beispiel mit $z_{3} = 4$.} \begin{center} \scriptsize \begin{tabular}{ |c|c|c|c| } \label{ilpspec2} Flow & Kante & Edukt & Produkt \\ \hline 3 & 1 & Ethylen & \\ 2 & 4 & 2 Butadien & Hexatrien + Ethylen\\ 1 & 18 & 2 Hexatrien & Octatetraen + Butadien\\ 1 & 41 & Benzol & \\ 1 & 42 & Octatetraen & Benzol + Ethylen\\ 3 & 213 & & Butadien\\ \hline \end{tabular} \end{center} \end{table} \newline Beide Lösungen die aus der Auswahl der besten Knoten resultieren ergeben chemisch Sinn. Die erste Lösung (Tabelle \ref{ilpspec1}) ist die erwartete Lösung. Die zweite Lösung (Tabelle \ref{ilpspec2}) verwendet viele der gleichen Kanten, aber produziert mit Kante 18, Octatetraen aus zwei Hexatrien. Diese Variation ist ebenfalls erwartbar. \section{Petrinetzzertifikat: Multi-Molekül-Reaktionen} Durch den einfachen Optimierungsprozess des ILP wird die Zielfunktion maximiert und eine Menge an Kanten zurückgegeben. Diese Kanten haben keine Reihenfolge oder Zusammenhang abseits der Definition des ILP. Darum kommt es oft vor, dass die Reaktionen in Sequenz ein Molekül für eine Reaktion benötigen, welches erst durch eine andere Reaktion erzeugt wird. Beide Reaktioen ließen sich so nur chemisch durchführen, wenn die jeweils andere zuvor stattgefunden hat. Um solche Lösungen des ILP-Models auszuschließen wurde für jede Lösung ein Petrinetzzertifikat erstellt. \subsection{Petrinetzzertifikat: Olefinmetathes und Diels-Alder} \subsubsection{Produkt} Die 10 Lösungen, welche durch die Bestimmung der Kantenwahrscheinlichkeit als Produkt beider Produkte erhalten wurden, wurden mit einem Petrinetz, basierend auf dem Hypergraphen, auf Erfüllbarkeit untersucht. Dabei wurden zu Beginn drei Token in Butadien angenommen. Nach der Untersuchung mit einem Petrinetz konnte nur eine Lösung zertifiziert werden (Tabelle \ref{ilpprodzert}). Dabei wird wie erwartet Hexatrien erzeugt. Aus Hexatrien erfolgt die Synthese von 3,5-Divinyl-Cyclohex-1-en. Es findet keine weitere Reaktion statt. \begin{table}[h!] \caption{Lösung 10 für die Kantenwahrscheinlichkeit als Produkt beim ILP vom Butadien Beispiel mit $z = 197.08$.} \begin{center} \scriptsize \begin{tabular}{ |c|c|c|c| } \label{ilpprodzert} Flow & Kante & Edukt & Produkt \\ \hline 1 & 1 & Ethylen & \\ 1 & 4 & 2 Butadien & Hexatrien + Ethylen\\ 1 & 36 & 3,5-Divinyl-Cyclohex-1-en & \\ 1 & 37 & Butadien + Ethylen & 3,5-Divinyl-Cyclohex-1-en\\ 3 & 213 & & Butadien\\ \hline \end{tabular} \end{center} \end{table} \newline Diese Lösung entspricht nicht der erwarteten Synthese von Benzol aus Butadien (Abbildung \ref{benzsynthmod}). \subsubsection{Durchschnitt} Die 10 Lösungen, welche durch die Bestimmung der Kantenwahrscheinlichkeit als Durchschnitt beider Produkte erhalten wurden, wurden mit einem Petrinetz, basierend auf dem Hypergraphen, auf Erfüllbarkeit untersucht. Dabei wurden zu Beginn drei Token in Butadien angenommen. Nach der Untersuchung mit einem Petrinetz konnte nur eine Lösung zertifiziert werden (Tabelle \ref{ilpavgzert}). \begin{table}[h!] \caption{Lösung 5 für die Kantenwahrscheinlichkeit als Durchschnitt beim ILP vom Butadien Beispiel mit $z = 798.54$.} \begin{center} \scriptsize \begin{tabular}{ |c|c|c|c| } \label{ilpavgzert} Flow & Kante & Edukt & Produkt \\ \hline 3 & 1 & Ethylen & \\ 1 & 4 & 2 Butadien & Hexatrien + Ethylen\\ 1 & 23 & Butadien + Hexatrien & Octatetraen + Ethylen\\ 1 & 41 & Benzol & \\ 1 & 42 & Octatetraen & Benzol + Ethylen\\ 3 & 213 & & Butadien\\ \hline \end{tabular} \end{center} \end{table} \newline Diese Lösung entspricht der erwarteten Synthese von Benzol aus Butadien (Abbildung \ref{benzsynthmod}). \subsubsection{Knoten Auswahl} Beide Lösungen der Knotenauswahlmethode erfüllen ein Petrinetzzertifikat. \section{Fazit: Multi-Molekül-Reaktionen} Im Vergleich zu einem kleinen Netzwerk mit festgelegten Start- und Zielmolekülen, ist die Formulierung der Zielfunktionen und Einschränkungen für ein Multi-Molekül-Reaktions-Netzwerk deutlich komplexer. Dazu kommt die Entscheidung nur die Ausgangsmoleküle zu definieren. Dabei kommt nur das ILP mit der Kantenwahrscheinlichkeit als Produkt gar nicht auf die erwartete Lösung. In beiden anderen Modellen ist die direkte Synthese von Benzol die erste Lösung mit Petrinetzzertifikat. Das ILP mit der Kantenwahrscheinlichkeit als Durchschnitt hatte trotzdem einige Lösungen, welche deutlich von der erwarteten Lösung abwichen und teilweise keine Kanten mit dieser teilten. Das ILP bei dem drei Knoten ausgewählt wurden die in der Synthese vorkommen müssen war dabei noch effizienter. Die zweitbeste Lösung bei diesem ILP ist ebenfalls eine realistische Möglichkeit unter den Reaktionsbedingungen, auch wenn sie in der Kohlenstoffeffizienz geringer ist. Die Lösung würde eigentlich vier Butadien zum Start benötigen. Da bei den Reaktionen aber ein neues Butadien entsteht, sind die ILP Bedingungen dennoch erfüllt. Da Butadien als Startmolekül theoretisch unbegrenzt vorliegt, ist die Synthese trotzdem chemisch möglich. Das überschüssige Butadien kann zu Beginn der nächsten Synthese erneut verbraucht werden. Zur Verifizierung des Butadien Beispiels wurden NMR-unabhängige Werte verwendet. Die betrachtete Synthese von Benzol aus Butadien stellt demnach auch ohne die Wahrscheinlichkeiten die optimale Lösung dar, wenn das Zielmolekül festgelegt wird. Ohne Wahrscheinlichkeiten oder mit zufälligen Wahrscheinlichkeiten und einem nicht definierten Zielmolekül, gibt es jedoch entweder eine Lösung ohne Reaktionen oder eine Lösung mit zufälligen Zielmolekül. Auch hier resultiert die Verwendung der NMR-basierten Wahrscheinlichkeiten in einer zielgerichteteren Lösungsauswahl. \newline Auch die Synthese von Benzol aus Butadien zeigt, dass die in den NMR-Spektren enthaltenen Informationen eine Rekonstruktion der vorgegeben Reaktionen ermöglicht. Mehr als in den natürlich auftretenden Synthesen muss angemerkt werden, dass die durch NMR-Spektren repräsentierten Reaktionen nicht auf in der Literatur vorgefundenen Synthesen oder im Labor gemessenen Spektren basiert. \chapter{Zusammenfassung und Ausblick} In dieser Masterarbeit wurde die Möglichkeit untersucht die Syntheseplanung und -analyse mittels NMR-Spektren zu automatisieren. Dabei wurden die charakteristischen Verschiebungen in $^{13}$C-NMR-Spektren verwendet, um die Wahrscheinlichkeit einzelner Moleküle und Reaktionen zu bestimmen. Diese Wahrscheinlichkeiten wurden über eine ILP-Berechnung optimiert. Noch öfter als $^{13}$C-NMR-Spektren in der organischen Chemie sind die $^{1}$H-NMR-Spektren. Diese haben zu den charakteristischen Verschiebungen noch zusätzliche Informationen durch Multiplettaufspaltungen. Da Wasserstoff im Gegensatz zu Kohlenstoffen als Protonen abgespalten werden können, kommt es vor allem bei Wasserstoff-Stickstoff- und Wasserstoff-Sauerstoff-Bindungen zu großen Schwankungen in der chemischen Verschiebung. Hinzu kommt, dass das quantenchemischen Rechnungsprogramm in der Berechnung der $^{1}$H-NMR-Spektren deutlich inkonstanter von den experimentellen Spektren abweicht. Zudem kommt es bei der Berechnung der Peakintensitäten und Multiplettaufspaltung zu größeren Fehlern, weshalb es für die verwendete Analyse nicht geeignet ist. Die Crosskorrelation bietet eine sehr einfache Methode zum Vergleich von NMR-Spektren. Da jedoch bei den $^{1}$H-NMR-Spektren die Multiplettaufspaltung oft mehr Informationen enthält als nur Intensität und chemische Verschiebung, geht bei der Analyse eine zentrale Komponente verloren. Beim Vergleich zwischen Spektren der gleichen Substanz, ist die Crosskorrelation auch bei $^{1}$H-NMR-Spektren geeignet, da nur leichte Verschiebungen der Peaks zu erwarten sind. Wenn jedoch zwei verschiedene Substanzen verglichen werden, kann die Verschiebung schon größer sein, weshalb die Ähnlichkeit gering ist, auch wenn die Spektren als Ganzes betrachtet ähnlich sind. Für eine solche Analyse könnte ein Machine-Learning-Model trainiert werden \cite{nmrsimml}, welches diese Details berücksichtigt. Zum Training dieses Models müssten jedoch deutlich mehr Beispiele verwendet werden und der NMR-Vergleich müsste eine deutlich zentralere Rolle in dieser Arbeit einnehmen. Mehr Beispiele können zudem verwendet werden, um den Schweratomeffekt in der NMR-Berechnung zu quantifizieren, sowie allgemeine Fehler in der Verschiebung zu ermitteln und so den idealen Verschiebungswert zu ermitteln. Die Kombination beider Ein-Molekül-Reaktionen legt jedoch nahe, dass dieser Fehler nicht durch eine Verschiebung mit einer Konstante zu beschreiben ist. In der NMR-Spektroskopie ist es auch möglich zwei Substanzen in einem NMR-Spektrum abzubilden. Dabei kann der Vergleich der Intensität der zugeordneten Peaks bestimmt werden, welches Stoffmengenverhältnis beide Substanzen zueinander haben \cite{nmrquant}. Es ist also möglich anstelle von mehreren NMR-Spektren für eine Synthese, nur ein gemeinsames NMR-Spektrum zu bestimmen und auf Basis dieses Spektrums neben den vorliegenden Substanzen, auch deren relative Stoffmenge zu bestimmen. Sowohl die natürlichen Synthesen im Organismus, als auch die Olefinmetathese und Diels-Alder-Reaktion für Butadien, eignen sich als Reaktion in einem Topf. Dadurch könnten die Startmoleküle mit den Enzymen und Katalysatoren in eine Lösung gebracht werden und nach einer festen Zeit NMR-Spektren aufgenommen werden. Basierend auf diesen Spektren könnten die relativen Stoffmengen und basierend auf diesen die Gleichgewichtskonstanten für diese Reaktionen bestimmt werden. Zudem erlaubt eine solche Analyse die Bestimmung von Nebenprodukten, welche in dem bisherigen Ansatz nicht definiert wurden. \newline Die Lösungen des ILP erfüllen die Erwartungen. Es ist aber denkbar, ein Model zu verwenden, dass in der Lage ist, die Tatsächlichen NMR-Spektren zu maximieren und so eine Abstraktion durch die Wahrscheinlichkeiten zu umgehen. Dabei könnte ein ILP-Model verwendet werden, dass zu einem Spektrum aus den Kandidatenmolekülen das am besten passende Molekül bestimmt. basierend auf diesem Molekül könnte ein weiteres ILP-Model diese Überlagerung für Reaktionen maximieren. Das verwendete ILP-Tool erlaubte jedoch nicht einzelnen Knoten größere Informationen zuzuordnen. Die bisherige Untersuchung eignet sich trotzdem als eine Demonstration der möglichen Kombination beider bisher getrennter Felder. \chapter*{Danksagung} Ich möchte die Gelegenheit nutzen, mich herzlichst bei dem Arbeitskreis Chemieinformatik zu bedanken. Besonderer Dank gilt Prof. Dr. Daniel Merkle, der mich sowohl tatkräftig bei inhaltlichen Fragen unterstützt hat als auch durch seine große Begeisterung eine konstante Quelle der Motivation war. Ebenso danke ich Dr. Yuri Kolcak, Addie Jordon und Cassandra König für die beratende Expertise zu theoretischen und inhaltlichen Fragen. Auch Dr. Marcel Friedrichs hat mich besonders bei technischen Problemen unterstützt. Zuletzt danke ich Verena Uhlig und Nicole Schellenberg, für die organisatorischen Unterstützung und das sie Prof. Merkle erinnert haben, wenn meine Fragen in den vielen Aufgaben etwas untergegangen sind. %%================================================================================== \addcontentsline{toc}{chapter}{Quellenverzeichnis} \renewcommand\bibname{Quellenverzeichnis} \fancyhead[R]{Quellenverzeichnis} \clearpage \printbibliography \clearpage \newpage %%================================================================================== \chapter*{Selbstständigkeitserklärung} \markboth{\textsc{Selbstständigkeitserklärung}}{\textsc{Selbstständigkeitserklärung}} \vspace{.5cm} \noindent \begin{tabular}{ll} {\bf Nachname:}& Seng \\ {\bf Vorname:} & Kilian \end{tabular} \vspace{1.2cm} \onehalfspacing \section*{Eigenständigkeitserklärung} Hiermit versichere ich, dass ich die vorliegende Arbeit selbstständig angefertigt und keine anderen als die angegebenen Hilfsmittel benutzt habe. Alle Stellen, die dem Wortlaut oder dem Sinn nach anderen Werken entnommen sind, habe ich in jedem einzelnen Fall unter genauer Angabe der Quelle (einschließlich des Internets sowie anderer digitaler Informationsquellen) als Entlehnung kenntlich gemacht. Dies gilt auch für eingefügte Zeichnungen, bildliche Darstellungen, Skizzen und Ähnliches. Ich nehme zur Kenntnis, dass die nachgewiesene Unterlassung der Herkunftsangabe als Plagiat bzw. Täuschungsversuch und auch die nachgewiesene Verwendung nicht angegebener KI- gestützter Werkzeuge als versuchte Täuschung gewertet und mit entsprechenden Maßnahmen geahndet wird. \newline \newline \begin{tabular}{lp{12em}l} Bielefeld, \today && \hspace{4cm} \\\cline{1-1}\cline{3-3} \\ \end{tabular} \vspace{-0.5 cm} \begin{tabular}{lp{20em}l} {\footnotesize Ort, Datum} && {\footnotesize Unterschrift } \end{tabular} \end{document}