Coarse-graining and quantum-classical adaptive coupling in soft matter

Die vorliegende Arbeit untersucht den Zusammenhang zwischen Skalen in Systemen weicher Materie, der für Multiskalen-Simulationen eine wichtige Rolle spielt. Zu diesem Zweck wurde eine Methode entwickelt, die die Approximation der Separierbarkeit von Variablen für die Molekulardynamik und ähnliche Anwendungen bewertet. Der zweite und größere Teil dieser Arbeit beschäftigt sich mit der konzeptionellen und technischen Erweiterung des Adaptive Resolution Scheme'' (AdResS), einer Methode zur gleichzeitigen Simulation von Systemen mit mehreren Auflösungsebenen. Diese Methode wurde auf Systeme erweitert, in denen klassische und quantenmechanische Effekte eine Rolle spielen.rnrnDie oben genannte erste Methode benötigt nur die analytische Form der Potentiale, wie sie die meisten Molekulardynamik-Programme zur Verfügung stellen. Die Anwendung der Methode auf ein spezielles Problem gibt bei erfolgreichem Ausgang einen numerischen Hinweis auf die Gültigkeit der Variablenseparation. Bei nicht erfolgreichem Ausgang garantiert sie, dass keine Separation der Variablen möglich ist. Die Methode wird exemplarisch auf ein zweiatomiges Molekül auf einer Oberfläche und für die zweidimensionale Version des Rotational Isomer State (RIS) Modells einer Polymerkette angewandt.rnrnDer zweite Teil der Arbeit behandelt die Entwicklung eines Algorithmus zur adaptiven Simulation von Systemen, in denen Quanteneffekte berücksichtigt werden. Die Quantennatur von Atomen wird dabei in der Pfadintegral-Methode durch einen klassischen Polymerring repräsentiert. Die adaptive Pfadintegral-Methode wird zunächst für einatomige Flüssigkeiten und tetraedrische Moleküle unter normalen thermodynamischen Bedingungen getestet. Schließlich wird die Stabilität der Methode durch ihre Anwendung auf flüssigen para-Wasserstoff bei niedrigen Temperaturen geprüft.

Scale bridging concepts in molecular simulation : coarse-graining and thermodynamic coupling

In dieser Arbeit wurden Simulation von Flüssigkeiten auf molekularer Ebene durchgeführt, wobei unterschiedliche Multi-Skalen Techniken verwendet wurden. Diese erlauben eine effektive Beschreibung der Flüssigkeit, die weniger Rechenzeit im Computer benötigt und somit Phänomene auf längeren Zeit- und Längenskalen beschreiben kann.rnrnEin wesentlicher Aspekt ist dabei ein vereinfachtes (“coarse-grained”) Modell, welches in einem systematischen Verfahren aus Simulationen des detaillierten Modells gewonnen wird. Dabei werden ausgewählte Eigenschaften des detaillierten Modells (z.B. Paar-Korrelationsfunktion, Druck, etc) reproduziert.rnrnEs wurden Algorithmen untersucht, die eine gleichzeitige Kopplung von detaillierten und vereinfachten Modell erlauben (“Adaptive Resolution Scheme”, AdResS). Dabei wird das detaillierte Modell in einem vordefinierten Teilvolumen der Flüssigkeit (z.B. nahe einer Oberfläche) verwendet, während der Rest mithilfe des vereinfachten Modells beschrieben wird.rnrnHierzu wurde eine Methode (“Thermodynamische Kraft”) entwickelt um die Kopplung auch dann zu ermöglichen, wenn die Modelle in verschiedenen thermodynamischen Zuständen befinden. Zudem wurde ein neuartiger Algorithmus der Kopplung beschrieben (H-AdResS) der die Kopplung mittels einer Hamilton-Funktion beschreibt. In diesem Algorithmus ist eine zur Thermodynamischen Kraft analoge Korrektur mit weniger Rechenaufwand möglich.rnrnAls Anwendung dieser grundlegenden Techniken wurden Pfadintegral Molekulardynamik (MD) Simulationen von Wasser untersucht. Mithilfe dieser Methode ist es möglich, quantenmechanische Effekte der Kerne (Delokalisation, Nullpunktsenergie) in die Simulation einzubeziehen. Hierbei wurde zuerst eine Multi-Skalen Technik (“Force-matching”) verwendet um eine effektive Wechselwirkung aus einer detaillierten Simulation auf Basis der Dichtefunktionaltheorie zu extrahieren. Die Pfadintegral MD Simulation verbessert die Beschreibung der intra-molekularen Struktur im Vergleich mit experimentellen Daten. Das Modell eignet sich auch zur gleichzeitigen Kopplung in einer Simulation, wobei ein Wassermolekül (beschrieben durch 48 Punktteilchen im Pfadintegral-MD Modell) mit einem vereinfachten Modell (ein Punktteilchen) gekoppelt wird. Auf diese Weise konnte eine Wasser-Vakuum Grenzfläche simuliert werden, wobei nur die Oberfläche im Pfadintegral Modell und der Rest im vereinfachten Modell beschrieben wird.

Distamycin-NA: A DNA analog with an aromatic heterocyclic polyamide backbone. Part 2. Solid-phase synthesis of distamycin-NAs containing the nucleobase uracil: Unexpected solvent participation in the coupling step

10.1002/hlca.19980810512.abs The synthesis of the Fmoc-protected amino acid 2 is presented. First attempts of amide-bond formation to the homodimer 4 in solution showed only poor coupling yields indicative for the low reactivity of the amino and carboxy groups in the building blocks 1 and 2, respectively (Scheme 1). Best coupling yields were found using dicyclohexylcarbodiimide (DCC) without any additive. The oligomerization of building block 2 adopting the Fmoc ((9H-fluoren-9-ylmethoxy)carbonyl) solid-phase synthesis yielded a mixture of N-terminal-modified distamycin-NA derivatives. By combined HPLC and MALDI-TOF-MS analysis, the N-terminal functional groups could be identified as acetamide and N,N-dimethylformamidine functions, arising from coupling of the N-terminus of the growing chain with residual AcOH or DCC-activated solvent DMF. An improved preparation of building block 2 and coupling protocol led to the prevention of the N-terminal acetylation. However, ‘amidination’ could not be circumvented. A thus isolated tetramer of 2, containing a lysine unit at the C-terminus and a N,N-dimethylformamidine-modified N-terminus, not unexpectedly, showed no complementary base pairing to DNA and RNA, as determined by standard UV-melting-curve analysis.

Coupling light into graphene plasmons with the help of surface acoustic waves

We propose a scheme for coupling laser light into graphene plasmons with the help of electrically generated surface acoustic waves. The surface acoustic wave forms a diffraction grating which allows us to excite the long lived phononlike branch of the hybridized graphene plasmon-phonon dispersion with infrared laser light. Our approach avoids patterning the graphene sheet, does not rely on complicated optical near-field techniques, and allows us to electrically switch the coupling between far-field radiation and propagating graphene plasmons.

ISA S88 / Loose Coupling Based Solution for Production Control Software Reusability

A solution for the problem of reusability of software system for batch production systems is proposed. It is based on ISA S88 standard that prescribes the abstraction of elements in the manufacturing system that is equipment, processes and procedures abstraction, required to make a product batch. An easy to apply data scheme, compatible with the standard, is developed for management of production information. In addition to flexibility provided by the S88 standard, software system reusability requires a solution supporting manufacturing equipment reconfigurability. Toward this end a coupling mechanism is developed. A software tool, including these solutions, was developed and validated at laboratory level, using product manufacturing information of an actual plant.

Multiple solutions and numerical analysis to the dynamic and stationary models coupling a delayed energy balance model involving latent heat and discontinuous albedo with a deep ocean

We study a climatologically important interaction of two of the main components of the geophysical system by adding an energy balance model for the averaged atmospheric temperature as dynamic boundary condition to a diagnostic ocean model having an additional spatial dimension. In this work, we give deeper insight than previous papers in the literature, mainly with respect to the 1990 pioneering model by Watts and Morantine. We are taking into consideration the latent heat for the two phase ocean as well as a possible delayed term. Non-uniqueness for the initial boundary value problem, uniqueness under a non-degeneracy condition and the existence of multiple stationary solutions are proved here. These multiplicity results suggest that an S-shaped bifurcation diagram should be expected to occur in this class of models generalizing previous energy balance models. The numerical method applied to the model is based on a finite volume scheme with nonlinear weighted essentially non-oscillatory reconstruction and Runge–Kutta total variation diminishing for time integration.

A staggered approach for the coupling of Cahn–Hilliard type diffusion and finite strain elasticity

We develop an algorithm and computational implementation for simulation of problems that combine Cahn–Hilliard type diffusion with finite strain elasticity. We have in mind applications such as the electro-chemo- mechanics of lithium ion (Li-ion) batteries. We concentrate on basic computational aspects. A staggered algorithm is pro- posed for the coupled multi-field model. For the diffusion problem, the fourth order differential equation is replaced by a system of second order equations to deal with the issue of the regularity required for the approximation spaces. Low order finite elements are used for discretization in space of the involved fields (displacement, concentration, nonlocal concentration). Three (both 2D and 3D) extensively worked numerical examples show the capabilities of our approach for the representation of (i) phase separation, (ii) the effect of concentration in deformation and stress, (iii) the effect of Electronic supplementary material The online version of this article (doi:10.1007/s00466-015-1235-1) contains supplementary material, which is available to authorized users. B P. Areias pmaa@uevora.pt 1 Department of Physics, University of Évora, Colégio Luís António Verney, Rua Romão Ramalho, 59, 7002-554 Évora, Portugal 2 ICIST, Lisbon, Portugal 3 School of Engineering, Universidad de Cuenca, Av. 12 de Abril s/n. 01-01-168, Cuenca, Ecuador 4 Institute of Structural Mechanics, Bauhaus-University Weimar, Marienstraße 15, 99423 Weimar, Germany strain in concentration, and (iv) lithiation. We analyze con- vergence with respect to spatial and time discretization and found that very good results are achievable using both a stag- gered scheme and approximated strain interpolation.

Transition to complete synchronization in phase-coupled oscillators with nearest neighbor coupling

We investigate synchronization in a Kuramoto-like model with nearest neighbor coupling. Upon analyzing the behavior of individual oscillators at the onset of complete synchronization, we show that the time interval between bursts in the time dependence of the frequencies of the oscillators exhibits universal scaling and blows up at the critical coupling strength. We also bring out a key mechanism that leads to phase locking. Finally, we deduce forms for the phases and frequencies at the onset of complete synchronization.

Border figure detection using a phase oscillator network with dynamical coupling

Oscillator networks have been developed in order to perform specific tasks related to image processing. Here we analytically investigate the existence of synchronism in a pair of phase oscillators that are short-range dynamically coupled. Then, we use these analytical results to design a network able of detecting border of black-and-white figures. Each unit composing this network is a pair of such phase oscillators and is assigned to a pixel in the image. The couplings among the units forming the network are also dynamical. Border detection emerges from the network activity.

Bs0 meson and the Bs0BK coupling from QCD sum rules

We evaluate the mass of the B(s0) scalar meson and the coupling constant in the B(s0)BK vertex in the framework of QCD sum rules. We consider the B(s0) as a tetraquark state to evaluate its mass. We get m(Bs0) = (5.85 +/- 0.13) GeV, which is in agreement, considering the uncertainties, with predictions supposing it as a b (s) over bar state or a B (K) over bar bound state with J(P) = 0(+). To evaluate the g(Bs0BK) coupling, we use the three-point correlation functions of the vertex, considering B(s0) as a normal b (s) over bar state. The obtained coupling constant is: g(Bs0BK) = (16.3 +/- 3.2) GeV. This number is in agreement with light-cone QCD sum rules calculation. We have also compared the decay width of the B(s0) -> BK process considering the B(s0) to be a b (s) over bar state and a BK molecular state. The width obtained for the BK molecular state is twice as big as the width obtained for the b (s) over bar state. Therefore, we conclude that with the knowledge of the mass and the decay width of the B(s0) meson, one can discriminate between the different theoretical proposals for its structure.

Thermodynamic origin of cooperativity in actomyosin interactions: The coupling of short-range interactions with actin bending stiffness in an Ising-like model

We present Monte Carlo simulations for a molecular motor system found in virtually all eukaryotic cells, the acto-myosin motor system, composed of a group of organic macromolecules. Cell motors were mapped to an Ising-like model, where the interaction field is transmitted through a tropomyosin polymer chain. The presence of Ca(2+) induces tropomyosin to block or unblock binding sites of the myosin motor leading to its activation or deactivation. We used the Metropolis algorithm to find the transient and the equilibrium states of the acto-myosin system composed of solvent, actin, tropomyosin, troponin, Ca(2+), and myosin-S1 at a given temperature, including the spatial configuration of tropomyosin on the actin filament surface. Our model describes the short- and long-range cooperativity during actin-myosin binding which emerges from the bending stiffness of the tropomyosin complex. We found all transition rates between the states only using the interaction energy of the constituents. The agreement between our model and experimental data also supports the recent theory of flexible tropomyosin.

Continuum-continuum coupling and polarization potentials for weakly bound systems

We investigate the influence of couplings among continuum states in collisions of weakly bound nuclei. For this purpose, we compare cross sections for complete fusion, breakup, and elastic scattering evaluated by continuum discretized coupled channel (CDCC) calculations, including and not including these couplings. In our study, we discuss this influence in terms of the polarization potentials that reproduces the elastic wave function of the coupled channel method in single channel calculations. We find that the inclusion of couplings among continuum states renders the real part of the polarization potential more repulsive, whereas it leads to weaker absorption to the breakup channel. We show that the noninclusion of continuum-continuum couplings in CDCC calculations may lead to qualitative and quantitative wrong conclusions.

Spin-Orbit Coupling in n-Type PbTe/PbEuTe Quantum Wells

Magnetoresistance measurements were performed on an n-type PbTe/PbEuTe quantum well and weak antilocalization effects were observed. This indicates the presence of spin orbit coupling phenomena and we showed that the Rashba effect is the main mechanism responsible for this spin orbit coupling. Using the model developed by Iordanskii et al., we fitted the experimental curves and obtained the inelastic and spin orbit scattering times. Thus we could compare the zero field energy spin-splitting predicted by the Rashba theory with the energy spin-splitting obtained from the analysis of the experimental curves. The final result confirms the theoretical prediction of strong Rashba effect on IV-VI based quantum wells.

Ferromagnetic coupling in a Co-doped graphenelike ZnO sheet

Using first-principles calculations it is demonstrated that Co doped graphenelike ZnO sheet presents ferromagnetic coupling. The Co atoms are energetically barrierless absorbed in the Zn sites, suffering a Jahn-Teller distortion. The results reveal that the origin of the ferromagnetic coupling, different from the bulk 3D ZnO stacking, is mainly guided by a direct exchange interaction without any additional defect. This ferromagnetic coupling is due to the system topology, namely, it is a direct consequence of the two-dimensional character of the ZnO monolayer within graphenelike structure. Increasing the number of ZnO layers the ferromagnetic coupling vanishes.

Amorphous HfO(2) and Hf(1-x)Si(x)O via a melt-and-quench scheme using ab initio molecular dynamics

We have performed ab initio molecular dynamics simulations to generate an atomic structure model of amorphous hafnium oxide (a-HfO(2)) via a melt-and-quench scheme. This structure is analyzed via bond-angle and partial pair distribution functions. These results give a Hf-O average nearest-neighbor distance of 2.2 angstrom, which should be compared to the bulk value, which ranges from 1.96 to 2.54 angstrom. We have also investigated the neutral O vacancy and a substitutional Si impurity for various sites, as well as the amorphous phase of Hf(1-x)Si(x)O(2) for x=0.25, 0375, and 0.5.