Computer simulation methods to study interfacial tensions : from the Ising model to colloidal crystals

In condensed matter systems, the interfacial tension plays a central role for a multitude of phenomena. It is the driving force for nucleation processes, determines the shape and structure of crystalline structures and is important for industrial applications. Despite its importance, the interfacial tension is hard to determine in experiments and also in computer simulations. While for liquid-vapor interfacial tensions there exist sophisticated simulation methods to compute the interfacial tension, current methods for solid-liquid interfaces produce unsatisfactory results.rnrnAs a first approach to this topic, the influence of the interfacial tension on nuclei is studied within the three-dimensional Ising model. This model is well suited because despite its simplicity, one can learn much about nucleation of crystalline nuclei. Below the so-called roughening temperature, nuclei in the Ising model are not spherical anymore but become cubic because of the anisotropy of the interfacial tension. This is similar to crystalline nuclei, which are in general not spherical but more like a convex polyhedron with flat facets on the surface. In this context, the problem of distinguishing between the two bulk phases in the vicinity of the diffuse droplet surface is addressed. A new definition is found which correctly determines the volume of a droplet in a given configuration if compared to the volume predicted by simple macroscopic assumptions.rnrnTo compute the interfacial tension of solid-liquid interfaces, a new Monte Carlo method called ensemble switch method'' is presented which allows to compute the interfacial tension of liquid-vapor interfaces as well as solid-liquid interfaces with great accuracy. In the past, the dependence of the interfacial tension on the finite size and shape of the simulation box has often been neglected although there is a nontrivial dependence on the box dimensions. As a consequence, one needs to systematically increase the box size and extrapolate to infinite volume in order to accurately predict the interfacial tension. Therefore, a thorough finite-size scaling analysis is established in this thesis. Logarithmic corrections to the finite-size scaling are motivated and identified, which are of leading order and therefore must not be neglected. The astounding feature of these logarithmic corrections is that they do not depend at all on the model under consideration. Using the ensemble switch method, the validity of a finite-size scaling ansatz containing the aforementioned logarithmic corrections is carefully tested and confirmed. Combining the finite-size scaling theory with the ensemble switch method, the interfacial tension of several model systems, ranging from the Ising model to colloidal systems, is computed with great accuracy.

Comparison of three optical tracking systems in a complex navigation scenario

Three-dimensional rotational X-ray imaging with the SIREMOBIL Iso-C3D (Siemens AG, Medical Solutions, Erlangen, Germany) has become a well-established intra-operative imaging modality. In combination with a tracking system, the Iso-C3D provides inherently registered image volumes ready for direct navigation. This is achieved by means of a pre-calibration procedure. The aim of this study was to investigate the influence of the tracking system used on the overall navigation accuracy of direct Iso-C3D navigation. Three models of tracking system were used in the study: Two Optotrak 3020s, a Polaris P4 and a Polaris Spectra system, with both Polaris systems being in the passive operation mode. The evaluation was carried out at two different sites using two Iso-C3D devices. To measure the navigation accuracy, a number of phantom experiments were conducted using an acrylic phantom equipped with titanium spheres. After scanning, a special pointer was used to pinpoint these markers. The difference between the digitized and navigated positions served as the accuracy measure. Up to 20 phantom scans were performed for each tracking system. The average accuracy measured was 0.86 mm and 0.96 mm for the two Optotrak 3020 systems, 1.15 mm for the Polaris P4, and 1.04 mm for the Polaris Spectra system. For the Polaris systems a higher maximal error was found, but all three systems yielded similar minimal errors. On average, all tracking systems used in this study could deliver similar navigation accuracy. The passive Polaris system showed ? as expected ? higher maximal errors; however, depending on the application constraints, this might be negligible.

Accurate Predictions of Water Cluster Formation, (H2O)n=2-10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Structures of Cage, Prism, and Book Isomers of Water Hexamer from Broadband Rotational Spectroscopy

Theory predicts the water hexamer to be the smallest water cluster with a three-dimensional hydrogen-bonding network as its minimum energy structure. There are several possible low-energy isomers, and calculations with different methods and basis sets assign them different relative stabilities. Previous experimental work has provided evidence for the cage, book, and cyclic isomers, but no experiment has identified multiple coexisting structures. Here, we report that broadband rotational spectroscopy in a pulsed supersonic expansion unambiguously identifies all three isomers; we determined their oxygen framework structures by means of oxygen-18–substituted water (H218O). Relative isomer populations at different expansion conditions establish that the cage isomer is the minimum energy structure. Rotational spectra consistent with predicted heptamer and nonamer structures have also been identified.

Do women have impaired regional systolic function in hypertensive heart disease? A 3-dimensional echocardiography study

AIMS: In pressure overload left ventricular (LV) hypertrophy, gender-related differences in global LV systolic function have been previously reported. The goal of this study was to determine regional systolic function of the left ventricle in male and female patients with hypertensive heart disease. METHODS AND RESULTS: Regional LV function was analyzed from multiplane transesophageal echocardiography with three-dimensional (3D) reconstruction of the left ventricle. In 24 patients (13 males and 11 females), four parallel (2 basal and 2 apical) equidistant short axis cross-sections from base to apex were obtained from the reconstructed LV. In each short axis 24 wall-thickness measurements were carried out at 15 degrees intervals at end-diastole and end-systole. Thus, a total of 192 measurements were obtained in each patient. Wall thickening was calculated as difference of end-diastolic and end-systolic wall thickness, and fractional thickening as thickening divided by end-diastolic thickness. Fractional thickening and wall stress were inversely related to end-diastolic wall thickness in both, males and females. Females showed less LV systolic function when compared to males (p<0.001). However, when corrected for wall stress, which was higher in females, there was no gender difference in systolic function. CONCLUSION: There are regional differences in LV systolic function in females and males which are directly related to differences in wall stress. Thus, gender-related differences in LV regional function are load-dependent and not due to structural differences.

ANALYSIS AND VISUALIZATION OF FLOW FIELDS USING INFORMATION-THEORETIC TECHNIQUES AND GRAPH-BASED REPRESENTATIONS

Three-dimensional flow visualization plays an essential role in many areas of science and engineering, such as aero- and hydro-dynamical systems which dominate various physical and natural phenomena. For popular methods such as the streamline visualization to be effective, they should capture the underlying flow features while facilitating user observation and understanding of the flow field in a clear manner. My research mainly focuses on the analysis and visualization of flow fields using various techniques, e.g. information-theoretic techniques and graph-based representations. Since the streamline visualization is a popular technique in flow field visualization, how to select good streamlines to capture flow patterns and how to pick good viewpoints to observe flow fields become critical. We treat streamline selection and viewpoint selection as symmetric problems and solve them simultaneously using the dual information channel [81]. To the best of my knowledge, this is the first attempt in flow visualization to combine these two selection problems in a unified approach. This work selects streamline in a view-independent manner and the selected streamlines will not change for all viewpoints. My another work [56] uses an information-theoretic approach to evaluate the importance of each streamline under various sample viewpoints and presents a solution for view-dependent streamline selection that guarantees coherent streamline update when the view changes gradually. When projecting 3D streamlines to 2D images for viewing, occlusion and clutter become inevitable. To address this challenge, we design FlowGraph [57, 58], a novel compound graph representation that organizes field line clusters and spatiotemporal regions hierarchically for occlusion-free and controllable visual exploration. We enable observation and exploration of the relationships among field line clusters, spatiotemporal regions and their interconnection in the transformed space. Most viewpoint selection methods only consider the external viewpoints outside of the flow field. This will not convey a clear observation when the flow field is clutter on the boundary side. Therefore, we propose a new way to explore flow fields by selecting several internal viewpoints around the flow features inside of the flow field and then generating a B-Spline curve path traversing these viewpoints to provide users with closeup views of the flow field for detailed observation of hidden or occluded internal flow features [54]. This work is also extended to deal with unsteady flow fields. Besides flow field visualization, some other topics relevant to visualization also attract my attention. In iGraph [31], we leverage a distributed system along with a tiled display wall to provide users with high-resolution visual analytics of big image and text collections in real time. Developing pedagogical visualization tools forms my other research focus. Since most cryptography algorithms use sophisticated mathematics, it is difficult for beginners to understand both what the algorithm does and how the algorithm does that. Therefore, we develop a set of visualization tools to provide users with an intuitive way to learn and understand these algorithms.

Three members of the connective tissue growth factor family CCN are differentially regulated by mechanical stress

Expression of connective tissue growth factor (CTGF), a member of the CCN gene family, is known to be significantly induced by mechanical stress. We have therefore investigated whether other members of the CCN gene family, including Cyr61 and Nov, might reveal a similar stress-dependent regulation. Fibroblasts growing under stressed conditions within a three-dimensional collagen gel showed at least a 15 times higher level of Cyr61 mRNA than cells growing under relaxed conditions. Upon relaxation, the decline of the Cyr61 mRNA to a lower level occurred within 2 h, and was thus quicker than the response of CTGF. The regulation was fully reversible when stress was reapplied. Thus, Cyr61 represents another typical example of a stress-responsive gene. The level of the Nov mRNA was low in the stressed state, but increased in the relaxed state. This CCN gene therefore shows an inverted regulation relative to that of Cyr61 and CTGF. Inhibition of protein kinases by means of staurosporine suppressed the stress-induced expression of Cyr61 and CTGF. Elevated levels of cAMP induced by forskolin mimicked the effects of relaxation on the regulation of Cyr61, CTGF and Nov. Thus, adenylate cyclase as well as one or several protein kinases might be involved in the mechanoregulation of these CCN genes.

Stereological Modelling of Random Particles

Stochastic models for three-dimensional particles have many applications in applied sciences. Lévy–based particle models are a flexible approach to particle modelling. The structure of the random particles is given by a kernel smoothing of a Lévy basis. The models are easy to simulate but statistical inference procedures have not yet received much attention in the literature. The kernel is not always identifiable and we suggest one approach to remedy this problem. We propose a method to draw inference about the kernel from data often used in local stereology and study the performance of our approach in a simulation study.

The braneology of 3D dualities

In this paper, we study the reduction of four-dimensional Seiberg duality to three dimensions from a brane perspective. We reproduce the nonperturbative dynamics of three-dimensional field theory via a T–duality at a finite radius and the action of Euclidean D–strings. In this way, we also overcome certain issues regarding the brane description of Aharony duality. Moreover, we apply our strategy to more general dualities, such as toric duality for M2–branes and dualities with adjoint matter fields.

Two-Dimensional Crystallisation of Membrane Proteins and Structural Assessment

Two-dimensional (2D) crystallisation of Membrane proteins reconstitutes them into their native environment, the lipid bilayer. Electron crystallography allows the structural analysis of these regular protein–lipid arrays up to atomic resolution. The crystal quality depends on the protein purity, ist stability and on the crystallisation conditions. The basics of 2D crystallisation and different recent advances are reviewed and electron crystallography approaches summarised. Progress in 2D crystallisation, sample preparation, image detectors and automation of the data acquisition and processing pipeline makes 2D electron crystallography particularly attractive for the structural analysis of membrane proteins that are too small for single-particle analyses and too unstable to form three-dimensional (3D) crystals.

Newton-Cartan supergravity with torsion and Schrödinger supergravity

We derive a torsionfull version of three-dimensional N=2 Newton-Cartan supergravity using a non-relativistic notion of the superconformal tensor calculus. The “superconformal” theory that we start with is Schrödinger supergravity which we obtain by gauging the Schrödinger superalgebra. We present two non-relativistic N=2 matter multiplets that can be used as compensators in the superconformal calculus. They lead to two different off-shell formulations which, in analogy with the relativistic case, we call “old minimal” and “new minimal” Newton-Cartan supergravity. We find similarities but also point out some differences with respect to the relativistic case.

Planktonic foraminifera distribution during the mid-Piacenzian warm period in the North Atlantic

The mid-Piacenzian (MP) warm period (3.264-3.025 Ma) has been identified as the most recent time in geologic history during which mean global surface temperatures were considerably warmer than today for a sustained period. This interval has therefore been proposed as a potential (albeit imperfect) analog for future climate change and as such, has received much scientific attention over the past two decades. Central to this research effort is the Pliocene Research, Interpretation, and Synoptic Mapping (PRISM) project, an iterative paleoenvironmental reconstruction of the MP focused on increasing our understanding of warm-period climate forcings, dynamics, and feedbacks by providing three-dimensional data sets for general circulation models. A mainstay of the PRISM project has been the development of a global sea surface temperature (SST) data set based primarily upon quantitative analyses of planktic foraminifer assemblages, supplemented with geochemical SST estimates wherever possible. In order to improve spatial coverage of the PRISM faunal data set in the low and mid-latitude North Atlantic, this study provides a description of the MP planktic foraminifer assemblage from five Ocean Drilling Program sites (951, 958, 1006, 1062, and 1063) in the subtropical gyre, a region critical to Atlantic Ocean circulation and tropical heat advection. Assemblages from each core provide evidence for a temperature- and circulation-driven 5-10° northward displacement of MP faunal provinces, as well as regional shifts in planktic foraminifer populations linked to species ecology and interactions. General biogeographic trends also indicate that, relative to modern conditions, gyre circulation was stronger (particularly the Gulf Stream, North Atlantic Current, and North Equatorial Current) and meridionally broader. A comparison of mid-Piacenzian and modern North Atlantic planktic foraminifer assemblages suggests that low latitude western boundary currents were less than 1 °C warmer while eastern boundary currents were ~1-2 °C warmer, supporting the hypothesis of enhanced northward heat advection along western boundary currents and warming of high latitude Northeast Atlantic source regions for the Canary Current. These findings are consistent with a model of reduced meridional SST gradients, with little-to-no low latitude warming, and more vigorous ocean circulation. Results therefore support the theory that enhanced meridional overturn circulation and associated northward heat advection made an important contribution, in conjunction with elevated atmospheric CO2 concentrations, to the 2-3 °C global surface temperature increase (relative to today) and strong polar amplification of SST warmth during the MP warm period.

Modelos de elementos finitos con aplicación tridimensional de la analogía fluido-fibra, para predecir el comportamiento mecánico de piezas estructurales de madera bajo la presencia de los principales tipos de nudos

Los nudos son los defectos que más disminuyen la resistencia de piezas de madera en la escala estructural, al ocasionar no solo una discontinuidad material, sino también la desviación de las fibras que se encuentran a su alrededor. En la década de los 80 se introdujo la teoría de la analogía fluido-fibra, como un método que aproximaba adecuadamente todas estas desviaciones. Sin embargo en su aplicación tridimensional, nunca se consideraron las diferencias geométricas en el sentido perpendicular al eje longitudinal de las piezas estructurales, lo cual imposibilitaba la simulación numérica de algunos de los principales tipos de nudos, y disminuía la precisión obtenida en aquellos nudos en los que la modelización sí era viable. En este trabajo se propone un modelo programado en lenguaje paramétrico de un software de elementos finitos que, bajo una formulación en tres dimensiones más general, permitirá estudiar de forma automatizada el comportamiento estructural de la madera bajo la influencia de los principales tipos de nudos, a partir de la geometría visible de los mismos y la posición de la médula en la pieza, y el cual ha sido contrastado experimentalmente, simulando de forma muy precisa el comportamiento mecánico de vigas sometidas a ensayos de flexión a cuatro puntos. Knots are the defects that most reduce the strength of lumber at the structural level, by causing not only a material discontinuity but also the deviation of the fibers that surround them. In the 80's it was introduced the theory of the flow-grain analogy as a method to approximating adequately these deviations. However, in three-dimensional applications, geometrical differences in the direction perpendicular to the longitudinal axis of the structural specimens were never considered before, which prevented the numerical simulation of some of the main types of knots, and decreased the achieved precision in those kind of knots where modeling itself was possible. This paper purposes a parametric model programmed in a finite element software, in the way that with a more general three-dimensional formulation, an automated study of the structural behavior of timber under the influence of the main types of knots is allowed by only knowing the visible geometry of such defects, and the position of the pith. Furthermore that has been confirmed experimentally obtaining very accurately simulations of the mechanical behavior of beams under four points bending test.

Modeling of the Division Point of Different Propagation Mechanisms in the Near-Region Within Arched Tunnels

An accurate characterization of the near-region propagation of radio waves inside tunnels is of practical importance for the design and planning of advanced communication systems. However, there has been no consensus yet on the propagation mechanism in this region. Some authors claim that the propagation mechanism follows the free space model, others intend to interpret it by the multi-mode waveguide model. This paper clarifies the situation in the near-region of arched tunnels by analytical modeling of the division point between the two propagation mechanisms. The procedure is based on the combination of the propagation theory and the three-dimensional solid geometry. Three groups of measurements are employed to verify the model in different tunnels at different frequencies. Furthermore, simplified models for the division point in five specific application situations are derived to facilitate the use of the model. The results in this paper could help to deepen the insight into the propagation mechanism within tunnel environments.

Global Linear Instability at the Dawn of its 4th Decade: A List of Challenges (A Practical Guide on how to Contain the Euphoria and Avoid the Oversell)

Global linear instability theory is concerned with the temporal or spatial development of small-amplitude perturbations superposed upon laminar steady or time-periodic threedimensional flows, which are inhomogeneous in two (and periodic in one) or all three spatial directions.1 The theory addresses flows developing in complex geometries, in which the parallel or weakly nonparallel basic flow approximation invoked by classic linear stability theory does not hold. As such, global linear theory is called to fill the gap in research into stability and transition in flows over or through complex geometries. Historically, global linear instability has been (and still is) concerned with solution of multi-dimensional eigenvalue problems; the maturing of non-modal linear instability ideas in simple parallel flows during the last decade of last century2–4 has given rise to investigation of transient growth scenarios in an ever increasing variety of complex flows. After a brief exposition of the theory, connections are sought with established approaches for structure identification in flows, such as the proper orthogonal decomposition and topology theory in the laminar regime and the open areas for future research, mainly concerning turbulent and three-dimensional flows, are highlighted. Recent results obtained in our group are reported in both the time-stepping and the matrix-forming approaches to global linear theory. In the first context, progress has been made in implementing a Jacobian-Free Newton Krylov method into a standard finite-volume aerodynamic code, such that global linear instability results may now be obtained in compressible flows of aeronautical interest. In the second context a new stable very high-order finite difference method is implemented for the spatial discretization of the operators describing the spatial BiGlobal EVP, PSE-3D and the TriGlobal EVP; combined with sparse matrix treatment, all these problems may now be solved on standard desktop computers.