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.

Perturbation theory for spectral subspaces

In der vorliegenden Arbeit wird die Variation abgeschlossener Unterräume eines Hilbertraumes untersucht, die mit isolierten Komponenten der Spektren von selbstadjungierten Operatoren unter beschränkten additiven Störungen assoziiert sind. Von besonderem Interesse ist hierbei die am wenigsten restriktive Bedingung an die Norm der Störung, die sicherstellt, dass die Differenz der zugehörigen orthogonalen Projektionen eine strikte Normkontraktion darstellt. Es wird ein Überblick über die bisher erzielten Resultate gegeben. Basierend auf einem Iterationsansatz wird eine allgemeine Schranke an die Variation der Unterräume für Störungen erzielt, die glatt von einem reellen Parameter abhängen. Durch Einführung eines Kopplungsparameters wird das Ergebnis auf den Fall additiver Störungen angewendet. Auf diese Weise werden zuvor bekannte Ergebnisse verbessert. Im Falle von additiven Störungen werden die Schranken an die Variation der Unterräume durch ein Optimierungsverfahren für die Stützstellen im Iterationsansatz weiter verschärft. Die zugehörigen Ergebnisse sind die besten, die bis zum jetzigen Zeitpunkt erzielt wurden.

Regularization methods for the solution of a nonlinear least-squares problem in tomography

In this work we study a polyenergetic and multimaterial model for the breast image reconstruction in Digital Tomosynthesis, taking into consideration the variety of the materials forming the object and the polyenergetic nature of the X-rays beam. The modelling of the problem leads to the resolution of a high-dimensional nonlinear least-squares problem that, due to its nature of inverse ill-posed problem, needs some kind of regularization. We test two main classes of methods: the Levenberg-Marquardt method (together with the Conjugate Gradient method for the computation of the descent direction) and two limited-memory BFGS-like methods (L-BFGS). We perform some experiments for different values of the regularization parameter (constant or varying at each iteration), tolerances and stop conditions. Finally, we analyse the performance of the several methods comparing relative errors, iterations number, times and the qualities of the reconstructed images.

Phase I/IIa study of cilengitide and temozolomide with concomitant radiotherapy followed by cilengitide and temozolomide maintenance therapy in patients with newly diagnosed glioblastoma

Invasion and migration are key processes of glioblastoma and are tightly linked to tumor recurrence. Integrin inhibition using cilengitide has shown synergy with chemotherapy and radiotherapy in vitro and promising activity in recurrent glioblastoma. This multicenter, phase I/IIa study investigated the efficacy and safety of cilengitide in combination with standard chemoradiotherapy in newly diagnosed glioblastoma.

The Politics of the "New North": Putting History and Geography at Stake in Arctic Futures

References to a “New North” have snowballed across popular media in the past 10 years. By invoking the phrase, scientists, policy analysts, journalists and others draw attention to the collision of global warming and global investment in the Arctic today and project a variety of futures for the region and the planet. While changes are apparent, the trope of a “New North” is not new. Discourses that appraised unfamiliar situations at the top of the world have recurred throughout the twentieth century. They have also accompanied attempts to cajole, conquer, civilize, consume, conserve and capitalize upon the far north. This article examines these politics of the “New North” by critically reading “New North” texts from the North American Arctic between 1910 and 2010. In each case, appeals to novelty drew from evaluations of the historical record and assessments of the Arctic’s shifting position in global affairs. “New North” authors pinpointed the ways science, state power, capital and technology transformed northern landscapes at different moments in time. They also licensed political and corporate influence in the region by delimiting the colonial legacies already apparent there. Given these tendencies, scholars need to approach the most recent iteration of the “New North” carefully without concealing or repeating the most troubling aspects of the Arctic’s past.

The quality of primary care in a country with universal health care coverage

Standard indicators of quality of care have been developed in the United States. Limited information exists about quality of care in countries with universal health care coverage.

Whither PQL?

Generalized linear mixed models (GLMM) are generalized linear models with normally distributed random effects in the linear predictor. Penalized quasi-likelihood (PQL), an approximate method of inference in GLMMs, involves repeated fitting of linear mixed models with “working” dependent variables and iterative weights that depend on parameter estimates from the previous cycle of iteration. The generality of PQL, and its implementation in commercially available software, has encouraged the application of GLMMs in many scientific fields. Caution is needed, however, since PQL may sometimes yield badly biased estimates of variance components, especially with binary outcomes. Recent developments in numerical integration, including adaptive Gaussian quadrature, higher order Laplace expansions, stochastic integration and Markov chain Monte Carlo (MCMC) algorithms, provide attractive alternatives to PQL for approximate likelihood inference in GLMMs. Analyses of some well known datasets, and simulations based on these analyses, suggest that PQL still performs remarkably well in comparison with more elaborate procedures in many practical situations. Adaptive Gaussian quadrature is a viable alternative for nested designs where the numerical integration is limited to a small number of dimensions. Higher order Laplace approximations hold the promise of accurate inference more generally. MCMC is likely the method of choice for the most complex problems that involve high dimensional integrals.

Origin of renal myofibroblasts in the model of unilateral ureter obstruction in the rat

Tubulo-interstitial fibrosis is a constant feature of chronic renal failure and it is suspected to contribute importantly to the deterioration of renal function. In the fibrotic kidney there exists, besides normal fibroblasts, a large population of myofibroblasts, which are supposedly responsible for the increased production of intercellular matrix. It has been proposed that myofibroblasts in chronic renal failure originate from the transformation of tubular cells via epithelial-mesenchymal transition (EMT) or from infiltration by bone marrow-derived precursors. Little attention has been paid to the possibility of a transformation of resident fibroblasts into myofibroblasts in renal fibrosis. Therefore we examined the fate of resident fibroblasts in the initial phase of renal fibrosis in the classical model of unilateral ureter obstruction (UUO) in the rat. Rats were perfusion-fixed on days 1, 2, 3 and 4 after ligature of the right ureter. Starting from 1 day of UUO an increasing expression of alpha-smooth muscle actin (alphaSMA) in resident fibroblasts was revealed by immunofluorescence and confirmed by the observation of bundles of microfilaments and webs of intermediate filaments in the electron microscope. Inversely, there was a decreased expression of 5'-nucleotidase (5'NT), a marker of renal cortical fibroblasts. The RER became more voluminous, suggesting an increased synthesis of matrix. Intercellular junctions, a characteristic feature of myofibroblasts, became more frequent. The mitotic activity in fibroblasts was strongly increased. Renal tubules underwent severe regressive changes but the cells retained their epithelial characteristics and there was no sign of EMT. In conclusion, after ureter ligature, resident peritubular fibroblasts proliferated and they showed progressive alterations, suggesting a transformation in myofibroblasts. Thus the resident fibroblasts likely play a central role in fibrosis in that model.

Automated optimization of polymer extrusion dies using finite element analysis

An extrusion die is used to continuously produce parts with a constant cross section; such as sheets, pipes, tire components and more complex shapes such as window seals. The die is fed by a screw extruder when polymers are used. The extruder melts, mixes and pressures the material by the rotation of either a single or double screw. The polymer can then be continuously forced through the die producing a long part in the shape of the die outlet. The extruded section is then cut to the desired length. Generally, the primary target of a well designed die is to produce a uniform outlet velocity without excessively raising the pressure required to extrude the polymer through the die. Other properties such as temperature uniformity and residence time are also important but are not directly considered in this work. Designing dies for optimal outlet velocity variation using simple analytical equations are feasible for basic die geometries or simple channels. Due to the complexity of die geometry and of polymer material properties design of complex dies by analytical methods is difficult. For complex dies iterative methods must be used to optimize dies. An automated iterative method is desired for die optimization. To automate the design and optimization of an extrusion die two issues must be dealt with. The first is how to generate a new mesh for each iteration. In this work, this is approached by modifying a Parasolid file that describes a CAD part. This file is then used in a commercial meshing software. Skewing the initial mesh to produce a new geometry was also employed as a second option. The second issue is an optimization problem with the presence of noise stemming from variations in the mesh and cumulative truncation errors. In this work a simplex method and a modified trust region method were employed for automated optimization of die geometries. For the trust region a discreet derivative and a BFGS Hessian approximation were used. To deal with the noise in the function the trust region method was modified to automatically adjust the discreet derivative step size and the trust region based on changes in noise and function contour. Generally uniformity of velocity at exit of the extrusion die can be improved by increasing resistance across the die but this is limited by the pressure capabilities of the extruder. In optimization, a penalty factor that increases exponentially from the pressure limit is applied. This penalty can be applied in two different ways; the first only to the designs which exceed the pressure limit, the second to both designs above and below the pressure limit. Both of these methods were tested and compared in this work.

Technical throughput of material flow systems required to achieve design flow – determined by simulation with aniteration algorithm

This article deals with complex material flow systems and series connections of conveyor and op-erating elements. These can be characterised by a specific availability. The thus resultant overall availabil-ity of necessary “technical throughput” of the individual elements for the achievement of a specified throughput. When the conveyor and operating elements are subjected to a stochastic distribution, the interposition of buffers is necessary but these can also lead to a reduction of the necessary throughput due faults. The system behaviour of complex installations can only be investigated by simulation. The parame-ter changes required in order to achieve specific target values can also be determined by simulation runs in iteration loops.

Approximating method of frames

In this paper, the well-known method of frames approach to the signal decomposition problem is reformulated as a certain bilevel goal-attainment linear least squares problem. As a consequence, a numerically robust variant of the method, named approximating method of frames, is proposed on the basis of a certain minimal Euclidean norm approximating splitting pseudo-iteration-wise method.

3D volumetric intensity reconsturction from 2D x-ray images using partial least squares regression

Reconstruction of shape and intensity from 2D x-ray images has drawn more and more attentions. Previously introduced work suffers from the long computing time due to its iterative optimization characteristics and the requirement of generating digitally reconstructed radiographs within each iteration. In this paper, we propose a novel method which uses a patient-specific 3D surface model reconstructed from 2D x-ray images as a surrogate to get a patient-specific volumetric intensity reconstruction via partial least squares regression. No DRR generation is needed. The method was validated on 20 cadaveric proximal femurs by performing a leave-one-out study. Qualitative and quantitative results demonstrated the efficacy of the present method. Compared to the existing work, the present method has the advantage of much shorter computing time and can be applied to both DXA images as well as conventional x-ray images, which may hold the potentials to be applied to clinical routine task such as total hip arthroplasty (THA).

Expectation Conditional Maximization-Based Deformable Shape Registration

This paper addresses the issue of matching statistical and non-rigid shapes, and introduces an Expectation Conditional Maximization-based deformable shape registration (ECM-DSR) algorithm. Similar to previous works, we cast the statistical and non-rigid shape registration problem into a missing data framework and handle the unknown correspondences with Gaussian Mixture Models (GMM). The registration problem is then solved by fitting the GMM centroids to the data. But unlike previous works where equal isotropic covariances are used, our new algorithm uses heteroscedastic covariances whose values are iteratively estimated from the data. A previously introduced virtual observation concept is adopted here to simplify the estimation of the registration parameters. Based on this concept, we derive closed-form solutions to estimate parameters for statistical or non-rigid shape registrations in each iteration. Our experiments conducted on synthesized and real data demonstrate that the ECM-DSR algorithm has various advantages over existing algorithms.

„Offene Grenzen“ oder „freiwillige Ghetttoisierung“? Aktuelle Kontroversen um die Sicherung der Kontinuität des Judentums

Die jüdische Gemeinschaft gilt generell als Musterbeispiel einer gut integrierten, religiösen Minderheit. Tatsächlich jedoch bewirken gerade die jüngsten gesellschaftlichen Entwicklungen − verstärkte Säkularisierung und Individualisierung verbunden mit steigenden Mischehenraten und einer Neudefinitionder Geschlechterrollen − eine Infragestellung der Kontinuität europäisch-jüdischer Existenz.Seit den 70iger Jahren des 20. Jahrhunderts bewegt sich die Mischehenrate fast überall in der Diaspora bei über 50%. Da die Weitergabe des Judentums religionsgesetzlich nur über die Mutter erfolgt,stellt der Umgang mit nichtjüdischen Familienmitgliedern einen hochsensiblen Bereich für die Gemeinschaft dar. Die soziale und religiöse Integration von nichtjüdischen Ehefrauen und vaterjüdischenKindern ist auf Grund einer nicht selten willkürlich erscheinenden Aufnahmepraxis ein häufig tabuisierter Aspekt des Gemeindelebens, der zu permanenten Spannungen führt. Konflikte bezüglich der Zugehörigkeitskriterien aber auch der religiösen Rolle der Frau führen zu Polarisierungs- und Pluralisierungstendenzen. Im Rahmen eines Projektes des NFP 58 wurden aktuelle innerjüdische Grenzziehungsdebatten im Kontext des Schweizer Judentums auch mit Methoden der Oral History festgehalten und analysiert. Die Auseinandersetzungen innerhalb der schweizerisch-jüdischen Gemeinschaft wurden zudem mit Entwicklungen in anderen Ländern der Diaspora und in Israel verglichen. Es ergab sich das Bild einer dynamischen und zugleich jedoch tief gespaltenen Religionsgemeinschaft, innerhalb der sich die verschiedenen Richtungen („liberal“ bis „ultra-orthodox“) die Verantwortung für eine zunehmende Schwächung und Spaltung des jüdischen Volkes zuweisen. Bibliographie Benbassa, Esther u. Jean-Christophe Attias. 2001. Les Juifs ont-t-ils un avenir? Paris. Lattés. Gerson, Daniel.2012. Ausbreitung und Bedeutung des Judentums in der Schweiz.in : Religionen in der Schweiz. Bulletin Schweizerische Akademie der Geistes- und Sozialwissenschaften, Nr 2.Bern. Gerson, Daniel.2011. Partizipation ohne Konversion? Grenzziehungsdebatten in neuen jüdischen Gemeinschaften der Schweiz,in: Chilufim. Zeitschrift für Jüdische Kulturgeschichte, Nr.11.Wien. Phoibos. Gerson, Daniel.2010. Gemeinschaftsbildung und «demokratischer» Antisemitismus: Das Entstehen eines Schweizer Judentums im Spannungsverhältnis von Akkulturation, Einwanderung und Ausgrenzung, in: Wyrwa, Ulrich (Hrsg.): Einspruch und Abwehr. Die Reaktion des europäischen Judentums und die Entstehung des Antisemitismus in Europa. Frankfurt am Main. Campus. Lambert, Nick.2008. Jews and Europe in the Twenty-First Century. London. Vallentine Mitchell. Picard, Jacques.2007. Judentum in der Schweiz: zwischen religiöser, kultureller und politischer Identität,in: Baumann, Martin u. Jörg Stolz (Hrsg.); Eine Schweiz - viele Religionen. Bielefeld. transcript. Wasserstein, Bernard.1996. Vanishing Diaspora. The Jews in Europa since 1945. New York.Harvard University Press.