Molecular dynamics simulations of sheared brush-like systems

Mit Hilfe von Molekulardynamik-Simulationen untersuchen wir bürstenartige Systeme unter guten Lösungsmittelbedingungen. Diese Systeme sind, dank ihren vielfältigen Beschaffenheiten, die von Molekularparametern und äußeren Bedingungen abhängig sind, wichtig für viele industrielle Anwendungen. Man vermutet, dass die Polymerbürsten eine entscheidende Rolle in der Natur wegen ihrer einzigartigen Gleiteigenschaften spielen. Ein vergröbertes Modell wird verwendet, um die strukturellen und dynamischen Eigenschaften zweier hochkomprimierter Polymerbürsten, die eine niedrige Reibung aufweisen, zu untersuchen. Allerdings sind die Lubrikationseigenschaften dieser Systeme, die in vielen biologischen Systemen vorhanden sind, beeinflußt. Wir untersuchen so-genannte "weiche Kolloide", die zwischen den beiden Polymerbürsten eingebettet sind, und wie diese Makroobjekte auf die Polymerbürsten wirken.rnrnNicht-Gleichgewichts-Molekulardynamik-Simulationen werden durchgeführt, in denen die hydrodynamischen Wechselwirkungen durch die Anwendung des DPD-Thermostaten mit expliziten Lösungsmittelmolekülen berücksichtigt werden. Wir zeigen, dass die Kenntnis der Gleichgewichtseigenschaften des Systems erlaubt, dynamische Nichtgleichgewichtsigenschaften der Doppelschicht vorherzusagen.rnrnWir untersuchen, wie die effektive Wechselwirkung zwischen kolloidalen Einschlüßen durch die Anwesenheit der Bürsten (in Abhängigkeit der Weichheit der Kolloide und der Pfropfdichte der Bürsten) beeinflußt wird. Als nächsten Schritt untersuchen wir die rheologische Antwort von solchen komplexen Doppelschichten auf Scherung. Wir entwickeln eine Skalen-Theorie, die die Abhängigkeit der makroskopischen Transporteigenschaften und der lateralen Ausdehnung der verankerten Ketten von der Weissenberg Zahl oberhalb des Bereichs, in dem die lineare Antwort-Theorie gilt, voraussagt. Die Vorhersagen der Theorie stimmen gut mit unseren und früheren numerischen Ergebnissen und neuen Experimenten überein. Unsere Theorie bietet die Möglichkeit, die Relaxationszeit der Doppelschicht zu berechnen. Wenn diese Zeit mit einer charakteristischen Längenskala kombiniert wird, kann auch das ''transiente'' (nicht-stationäre) Verhalten beschrieben werden.rnrnrnWir untersuchen die Antwort des Drucktensors und die Deformation der Bürsten während der Scherinvertierung für grosse Weissenberg Zahlen. Wir entwickeln eine Vorhersage für die charakteristische Zeit, nach der das System wieder den stationären Zustand erreicht.rnrnrnElektrostatik spielt eine bedeutende Rolle in vielen biologischen Prozessen. Die Lubrikationseigenschaften der Polymerbürsten werden durch die Anwesenheit langreichweitiger Wechselwirkungen stark beeinflusst. Für unterschiedliche Stärken der elektrostatischen Wechselwirkungen untersuchen wir rheologische Eigenschaften der Doppelschicht und vergleichen mit neutralen Systemen. Wir studieren den kontinuierlichen Übergang der Systemeigenschaften von neutralen zu stark geladenen Bürsten durch Variation der Bjerrumlänge und der Ladungsdichte.

Responsive Bürstenpolymere

'Responsive' Bürstenpolymere Bürstenpolymere sind definiert verzweigte Makromoleküle, die aus einer Hauptkette und vielen darauf (kovalent) gepfropften Seitenketten bestehen; ist der Pfropfungsgrad hoch und die Hauptkette wesentlich länger als die Seitenketten, dann haben sie die Form semiflexibler molekularer Zylinder. Lassen sich Form bzw. Ausdehnung eines solchen Zylinders gezielt ansteuern, dann könnten diese Moleküle entweder als (Nano-)Sensoren für die entsprechende Umgebungsbedingung oder als molekulare Motoren eingesetzt werden. Die Idee responsiver Bürstenpolymere beruht auf folgender Überlegung: Die gestreckte Konformation der Hauptkette ist entropisch gegenüber einem entsprechenden Knäuel benachteiligt, weshalb sie ,molekulare Federn‘ darstellen, die auf Änderung der repulsiven Wechselwirkung zwischen den Seitenketten reagieren. Dies wurde für den Wechsel zwischen gutem und schlechtem Lösungsmitteln untersucht. Ein zweites Konzept zur Änderung der Molekülform beruht auf der intramolekularen Phasentrennung (,Segmentbildung‘) miteinander unverträglicher Seitenketten in selektiven Lösungsmitteln, da die Hauptkette durch Ausbildung von Mikrophasen entlang des Moleküls ebenfalls aus ihrer gestreckten Form gebracht werden sollte. Die dritte Möglichkeit zur Änderung der Konformation ist die intramolekulare Vernetzung von Seitenketten, die ebenfalls zu verringerter Abstoßung und damit zur Verkürzung der Zylinder führen sollte. Eine weitere wichtige Untersuchung der Arbeit war der Übergang einer geknäuelten Hauptkette zu einer gestreckten Bürste als Funktion der Pfropfdichte. Zur Beantwortung dieser Fragestellungen wurden zylindrische Bürstenpolymere durch ,Grafting Trough‘ und ,Grafting Onto‘ synthetisiert (PS bzw. PI/PS und PnBMA/PMAA mit Kern/Schale- und ,Segment‘-Architektur) und systematisch Pfropfdichte, Vernetzungsgrad (Vernetzung durch gamma-Bestrahlung) und Lösungsbedingungen verändert. Die Möglichkeit gezielter Ansteuerung der Konformationsänderung durch Vernetzung konnte nach polymeranaloger Modifikation von PI/PS-Bürstenpolymeren durch Photovernetzung und vernetzende Komplexierung erfolgreich bestätigt werden. Zur Untersuchung der Probenreihen wurden AFM, Licht- und Neutronenstreuung herangezogen. Die Analysen bestätigten konsistent die Änderung von Steifigkeit, Zylinderquerschnitt und Streckung der Hauptkette durch Variation von Pfropfdichte, Vernetzung und Lösungsmittelqualität. Für die Änderung der Pfropfdichte gehorchen die Parameter dabei Potenzgesetzen.

Charge-transport simulations in organic semiconductors

In this thesis we have extended the methods for microscopic charge-transport simulations for organic semiconductors. In these materials the weak intermolecular interactions lead to spatially localized charge carriers, and the charge transport occurs as an activated hopping process between diabatic states. In addition to weak electronic couplings between these states, different electrostatic environments in the organic material lead to a broadening of the density of states for the charge energies which limits carrier mobilities.rnThe contributions to the method development includern(i) the derivation of a bimolecular charge-transfer rate,rn(ii) the efficient evaluation of intermolecular (outer-sphere) reorganization energies,rn(iii) the investigation of effects of conformational disorder on intramolecular reorganization energies or internal site energiesrnand (iv) the inclusion of self-consistent polarization interactions for calculation of charge energies.These methods were applied to study charge transport in amorphous phases of small molecules used in the emission layer of organic light emitting diodes (OLED).rnWhen bulky substituents are attached to an aromatic core in order to adjust energy levels or prevent crystallization, a small amount of delocalization of the frontier orbital to the substituents can increase electronic couplings between neighboring molecules. This leads to improved charge-transfer rates and, hence, larger charge-mobility. We therefore suggest using the mesomeric effect (as opposed to the inductive effect) when attaching substituents to aromatic cores, which is necessary for example in deep blue OLEDs, where the energy levels of a host molecule have to be adjusted to those of the emitter.rnFurthermore, the energy landscape for charges in an amorphous phase cannot be predicted by mesoscopic models because they approximate the realistic morphology by a lattice and represent molecular charge distributions in a multipole expansion. The microscopic approach shows that a polarization-induced stabilization of a molecule in its charged and neutral states can lead to large shifts, broadening, and traps in the distribution of charge energies. These results are especially important for multi-component systems (the emission layer of an OLED or the donor-acceptor interface of an organic solar cell), if the change in polarizability upon charging (or excitation in case of energy transport) is different for the components. Thus, the polarizability change upon charging or excitation should be added to the set of molecular parameters essential for understanding charge and energy transport in organic semiconductors.rnWe also studied charge transport in self-assembled systems, where intermolecular packing motives induced by side chains can increase electronic couplings between molecules. This leads to larger charge mobility, which is essential to improve devices such as organic field effect transistors, where low carrier mobilities limit the switching frequency.rnHowever, it is not sufficient to match the average local molecular order induced by the sidernchains (such as the pitch angle between consecutive molecules in a discotic mesophase) with maxima of the electronic couplings.rnIt is also important to make the corresponding distributions as narrow as possible compared to the window determined by the closest minima of thernelectronic couplings. This is especially important in one-dimensional systems, where charge transport is limited by the smallest electronic couplings.rnThe immediate implication for compound design is that the side chains should assist the self-assemblingrnprocess not only via soft entropic interactions, but also via stronger specific interactions, such as hydrogen bonding.rnrnrnrn

A molecular approach towards tethered bilayer lipid membranes: synthesis and characterization of novel anchor lipids

Tethered bilayer lipid membranes (tBLMs) are a promising model system for the natural cell membrane. They consist of a lipid bilayer that is covalently coupled to a solid support via a spacer group. In this study, we developed a suitable approach to increase the submembrane space in tBLMs. The challenge is to create a membrane with a lower lipid density in order to increase the membrane fluidity, but to avoid defects that might appear due to an increase in the lateral space within the tethered monolayers. Therefore, various synthetic strategies and different monolayer preparation techniques were examined. Synthetical attempts to achieve a large ion reservoir were made in two directions: increasing the spacer length of the tether lipids and increasing the lateral distribution of the lipids in the monolayer. The first resulted in the synthesis of a small library of tether lipids (DPTT, DPHT and DPOT) characterized by 1H and 13C NMR, FD-MS, ATR, DSC and TGA. The synthetic strategy for their preparation includes synthesis of precursor with a double bond anchor that can be easily modified for different substrates (e.g. metal and metaloxide). Here, the double bond was modified into a thiol group suitable for gold surface. Another approach towards the preparation of homogeneous monolayers with decreased two-dimensional packing density was the synthesis of two novel anchor lipids: DPHDL and DDPTT. DPHDL is “self-diluted” tether lipid containing two lipoic anchor moieties. DDPTT has an extended lipophylic part that should lead to the preparation of diluted, leakage free proximal layers that will facilitate the completion of the bilayer. Our tool-box of tether lipids was completed with two fluorescent labeled lipid precursors with respectively one and two phytanyl chains in the hydrophobic region and a dansyl group as a fluorophore. The use of such fluorescently marked lipids is supposed to give additional information for the lipid distribution on the air-water interface. The Langmuir film balance was used to investigate the monolayer properties of four of the synthesized thiolated anchor lipids. The packing density and mixing behaviour were examined. The results have shown that mixing anchor with free lipids can homogeneously dilute the anchor lipid monolayers. Moreover, an increase in the hydrophylicity (PEG chain length) of the anchor lipids leads to a higher packing density. A decrease in the temperature results in a similar trend. However, increasing the number of phytanyl chains per lipid molecule is shown to decrease the packing density. LB-monolayers based on pure and mixed lipids in different ratio and transfer pressure were tested to form tBLMs with diluted inner layers. A combination of the LB-monolayer transfer with the solvent exchange method accomplished successfully the formation of tBLMs based on pure DPOT. Some preliminary investigations of the electrical sealing properties and protein incorporation of self-assembled DPOT and DDPTT-based tBLMs were conducted. The bilayer formation performed by solvent exchange resulted in membranes with high resistances and low capacitances. The appearance of space beneath the membrane is clearly visible in the impedance spectra expressed by a second RC element. The latter brings the conclusion that the longer spacer in DPOT and the bigger lateral space between the DDPTT molecules in the investigated systems essentially influence the electrical parameters of the membrane. Finally, we could show the functional incorporation of the small ion carrier valinomycin in both types of membranes.

Spectroscopic properties from hybrid QM, MM molecular dynamics simulation

The central aim of this thesis work is the application and further development of a hybrid quantum mechanical/molecular mechanics (QM/MM) based approach to compute spectroscopic properties of molecules in complex chemical environments from electronic structure theory. In the framework of this thesis, an existing density functional theory implementation of the QM/MM approach is first used to calculate the nuclear magnetic resonance (NMR) solvent shifts of an adenine molecule in aqueous solution. The findings show that the aqueous solvation with its strongly fluctuating hydrogen bond network leads to specific changes in the NMR resonance lines. Besides the absolute values, also the ordering of the NMR lines changes under the influence of the solvating water molecules. Without the QM/MM scheme, a quantum chemical calculation could have led to an incorrect assignment of these lines. The second part of this thesis describes a methodological improvement of the QM/MM method that is designed for cases in which a covalent chemical bond crosses the QM/MM boundary. The development consists in an automatized protocol to optimize a so-called capping potential that saturates the electronic subsystem in the QM region. The optimization scheme is capable of tuning the parameters in such a way that the deviations of the electronic orbitals between the regular and the truncated (and "capped") molecule are minimized. This in turn results in a considerable improvement of the structural and spectroscopic parameters when computed with the new optimized capping potential within the QM/MM technique. This optimization scheme is applied and benchmarked on the example of truncated carbon-carbon bonds in a set of small test molecules. It turns out that the optimized capping potentials yield an excellent agreement of NMR chemical shifts and protonation energies with respect to the corresponding full molecules. These results are very promising, so that the application to larger biological complexes will significantly improve the reliability of the prediction of the related spectroscopic properties.

Systematic functionalization of molecules for molecular self-assembly

This thesis describes the investigation of systematically varied organic molecules for use in molecular self-assembly processes. All experiments were performed using high-resolution non-contact atomic force microscopy under UHV conditions and at room temperature. Using this technique, three different approaches for influencing intermolecular and molecule-surface interaction on the insulating calcite(10.4) surface were investigated by imaging the structure formation at the molecular scale. I first demonstrated the functionalization of shape-persistent oligo(p-benzamide)s that was engineered by introducing different functional groups and investigating their effect on the structural formation on the sample surface. The molecular core was designed to provide significant electrostatic anchoring towards the surface, while at the same time maintaining the flexibility to fine-tune the resulting structure by adjusting the intermolecular cohesion energy. The success of this strategy is based on a clear separation of the molecule-substrate interaction from the molecule-molecule interaction. My results show that sufficient molecule-surface anchoring can be achieved without restricting the structural flexibility that is needed for the design of complex molecular systems. Three derivatives of terephthalic acid (TPA) were investigated in chapter 7. Here, the focus was on changing the adhesion to the calcite surface by introducing different anchor functionalities to the TPA backbone. For all observed molecules, the strong substrate templating effect results in molecular structures that are strictly oriented along the calcite main crystal directions. This templating is especially pronounced in the case of 2-ATPA where chain formation on the calcite surface is observed in contrast to the formation of molecular layers in the bulk. At the same time, the amino group of 2-ATPA proved an efficient anchor functionality, successfully stabilizing the molecular chains on the sample surface. These findings emphasizes, once again, the importance of balancing and fine-tuning molecule-molecule and molecule-surface interactions in order to achieve stable, yet structurally flexible molecular arrangements on the sample surface. In the last chapter, I showed how the intrinsic property of molecular chirality decisively influences the structure formation in molecular self-assembly. This effect is especially pronounced in the case of the chiral heptahelicene-2-carboxylic acid. Deposition of the enantiopure molecules results in the formation of homochiral islands on the sample surface which is in sharp contrast to the formation of uni-directional double rows upon deposition of the racemate onto the same surface. While it remained uncertain from these previous experiments whether the double rows are composed of hetero- or homochiral molecules, I could clearly answer that question here and demonstrate that the rows are of heterochiral origin. Chirality, thus, proves to be another important parameter to steer the intermolecular interaction on surfaces. Altogether, the results of this thesis demonstrate that, in order to successfully control the structure formation in molecular self-assembly, the correct combination of molecule and surface properties is crucial. This is of special importance when working on substrates that exhibit a strong influence on the structure formation, such as the calcite(10.4) surface. Through the systematic variation of functional groups several important parameters that influence the balance between molecule-surface and molecule-molecule interaction were identified here, and the results of this thesis can, thus, act as a guideline for the rational design of molecules for use in molecular self-assembly.

Molecular plating of thin lanthanide layers with improved material properties for nuclear applications

In der vorliegenden Arbeit werden Experimente beschrieben, die zu einem vertieften Verständnis fundamentaler Prozesse bei der elektrochemischen Herstellung von Dünnschichten, sog. Targets, für kernphysikalische und -chemische Studien führten. Targets wurden mittels 'Molecular Plating' (MP) hergestellt, indem eine Elektrodeposition aus organischem Medium in der Regel bei konstantem Strom in Zwei-Elektroden-Zellen. Die Resultate erlaubten, optimierte Herstellungs-bedingungen zu ermitteln, welche die Produktion deutlich verbesserter Targets erlaubten. MP bei konstantem Strom ist ein massentransportkontrollierter Prozess. Der angelegte Strom wird durch einen konstanten Fluss elektroaktiver Spezies zur Kathode – auf der die Schicht wächst – und Anode aufrechterhalten. Die Untersuchungen zeigten, dass das Zellenpotential des Elektrodepositionsystems immer durch den Ohm'schen Spannungsabfall auf Grund des Widerstandes der verwendeten Lösung dominiert wurde. Dies erlaubte die Herleitung einer Beziehung zwischen dem Zellenpotential und der Konzentration der elektroaktiven Spezies. Die Beziehung erlaubt die Erklärung des gemessenen zeitlichen Verlaufs des Zellenpotentials während der Abscheidung als Funktion der Elektrolytkonzentration. Dies dient als Basis, auf der nun ein umfassenderes Bild der Prozesse, die für die charakteristischen Minima im Potentialverlauf einer Abscheidung verantwortlich sind, gewonnen werden kann. Es konnte gezeigt werden, dass die Minima mit der fast vollständigen Entfernung (durch Abscheidung) der aus einem gelösten Salz erzeugten Nd-Ionen korrespondieren. Die abgeschiedene Spezies wurde als Nd3+ identifiziert, vermutlich als Carboxylat, Oxid oder Hydroxid, was auf Grund der hohen negative Werte des Standardredoxpotentials der Lanthanide verständlich erscheint. Von den vorliegenden elektroaktiven Spezies tragen die Nd3+ Ionen nur zu knapp 20% zum Gesamtstrom bei. Durch Elektrolyse tragen auch die Lösungsmittelkomponenten zu diese Strom bei. Die Gegenwart von elektrolysiertem Lösungsmittel wurde in Analysen der Dünnschichten bestätigt. Diese waren immer mit chemi- und physisorbierten Lösungsmittelmolekülen bedeckt. Die Analyse der Dünnschichten zeigte, dass die Oberflächen von einem furchenartiges Netz durchzogen waren, und dass diese während des Trocknen der Schichten nach dem MP entstanden. Ob die Schichten an Luft oder in inerter Atmosphäre trockneten, hatte keinen Einfluss. Es wurden Experimente mit mehreren Lösungsmitteln durchgeführt, die sich deutlich in ihren physikalischen Eigenschaften, v.a. dem Siedepunkt, unterschieden. Furchenfreie Dünnschichten konnten insbesondere bei MP in N,N-dimethylformamide (DMF) erzeugt werden. Die Verwendung von DMF in Kombination mit einer Abscheidung auf sehr glatten Substraten erlaubte die Produktion von sehr homogenen, glatten und defektfreien Schichten. Diese waren vermutlich geringeren inneren Spannungen während des Trocknens ausgesetzt, als Schichten auf raueren Substraten oder solche, die aus flüchtigeren Lösungsmitteln hergestellt wurden. Die Oberflächenrauigkeit des Substrats und das gewählte Lösungsmittel wurden so als Schlüsselfaktoren für die Produktion hochqualitativer Schichten identifiziert. Es konnte gezeigt werden, dass mit MP eine sehr effiziente Methode zur Herstellung homogener Schichten mit exzellenter Ausbeute ist. In weiteren Experimenten mit dem primordialen Alpha-Emitter 147Sm als Modellisotop wurde die Eignung solcher Schichten als Alpha-Quelle untersucht. Sowohl die Energieauflösung als auch der Anteil der Alpha-Teilchen, die den Detektor erreichten, waren von den Quelleneigenschaften abhängig. Die Effekte wurden verschiedenen Variablen der Dünnschicht zugeordnet, welche die Alpha-Spektren beeinflussten. Dominant war die Wahl des Lösungsmittels und die Rauigkeit des Substrats. Dies beeinflusste Schichtdicke und -morphologie sowie die Art des Schichtwachstums und veränderte die Detektionseffizienz in Alpha-Messungen bis zu 15%. Nur homogene, ebene Schichten, die aus DMF auf glatten Substraten abgeschieden wurden, eignen sich optimal als Alpha-Quelle. Die gewonnenen Ergebnisse erlauben die optimierte Herstellung nuklearer Targets durch MP. Künftige Anwendungen beinhalten insbesondere die Herstellung von Targets für neutroneninduzierte Spaltexperimente und untergrundarmeAlpha-Messungen sehr kleiner Aktivitäten.

Novel algorithms for electronic structure based molecular dynamics with linear system-size scaling

Die vorliegende Arbeit behandelt die Entwicklung und Verbesserung von linear skalierenden Algorithmen für Elektronenstruktur basierte Molekulardynamik. Molekulardynamik ist eine Methode zur Computersimulation des komplexen Zusammenspiels zwischen Atomen und Molekülen bei endlicher Temperatur. Ein entscheidender Vorteil dieser Methode ist ihre hohe Genauigkeit und Vorhersagekraft. Allerdings verhindert der Rechenaufwand, welcher grundsätzlich kubisch mit der Anzahl der Atome skaliert, die Anwendung auf große Systeme und lange Zeitskalen. Ausgehend von einem neuen Formalismus, basierend auf dem großkanonischen Potential und einer Faktorisierung der Dichtematrix, wird die Diagonalisierung der entsprechenden Hamiltonmatrix vermieden. Dieser nutzt aus, dass die Hamilton- und die Dichtematrix aufgrund von Lokalisierung dünn besetzt sind. Das reduziert den Rechenaufwand so, dass er linear mit der Systemgröße skaliert. Um seine Effizienz zu demonstrieren, wird der daraus entstehende Algorithmus auf ein System mit flüssigem Methan angewandt, das extremem Druck (etwa 100 GPa) und extremer Temperatur (2000 - 8000 K) ausgesetzt ist. In der Simulation dissoziiert Methan bei Temperaturen oberhalb von 4000 K. Die Bildung von sp²-gebundenem polymerischen Kohlenstoff wird beobachtet. Die Simulationen liefern keinen Hinweis auf die Entstehung von Diamant und wirken sich daher auf die bisherigen Planetenmodelle von Neptun und Uranus aus. Da das Umgehen der Diagonalisierung der Hamiltonmatrix die Inversion von Matrizen mit sich bringt, wird zusätzlich das Problem behandelt, eine (inverse) p-te Wurzel einer gegebenen Matrix zu berechnen. Dies resultiert in einer neuen Formel für symmetrisch positiv definite Matrizen. Sie verallgemeinert die Newton-Schulz Iteration, Altmans Formel für beschränkte und nicht singuläre Operatoren und Newtons Methode zur Berechnung von Nullstellen von Funktionen. Der Nachweis wird erbracht, dass die Konvergenzordnung immer mindestens quadratisch ist und adaptives Anpassen eines Parameters q in allen Fällen zu besseren Ergebnissen führt.

Numerical analysis of the relation between interactions and structure in a molecular fluid

Coarse graining is a popular technique used in physics to speed up the computer simulation of molecular fluids. An essential part of this technique is a method that solves the inverse problem of determining the interaction potential or its parameters from the given structural data. Due to discrepancies between model and reality, the potential is not unique, such that stability of such method and its convergence to a meaningful solution are issues.rnrnIn this work, we investigate empirically whether coarse graining can be improved by applying the theory of inverse problems from applied mathematics. In particular, we use the singular value analysis to reveal the weak interaction parameters, that have a negligible influence on the structure of the fluid and which cause non-uniqueness of the solution. Further, we apply a regularizing Levenberg-Marquardt method, which is stable against the mentioned discrepancies. Then, we compare it to the existing physical methods - the Iterative Boltzmann Inversion and the Inverse Monte Carlo method, which are fast and well adapted to the problem, but sometimes have convergence problems.rnrnFrom analysis of the Iterative Boltzmann Inversion, we elaborate a meaningful approximation of the structure and use it to derive a modification of the Levenberg-Marquardt method. We engage the latter for reconstruction of the interaction parameters from experimental data for liquid argon and nitrogen. We show that the modified method is stable, convergent and fast. Further, the singular value analysis of the structure and its approximation allows to determine the crucial interaction parameters, that is, to simplify the modeling of interactions. Therefore, our results build a rigorous bridge between the inverse problem from physics and the powerful solution tools from mathematics. rn