Modulation of charge-density waves by superlattice structures

We discuss the interplay between electronic correlations and an underlying superlattice structure in determining the period of charge density waves (CDW's), by considering a one-dimensional Hubbard model with a repeated (nonrandom) pattern of repulsive (U > 0) and free (U=0) sites. Density matrix renormalization group diagonalization of finite systems (up to 120 sites) is used to calculate the charge-density correlation function and structure factor in the ground state. The modulation period can still be predicted through effective Fermi wave vectors k(F)(*) and densities, and we have found that it is much more sensitive to electron (or hole) doping, both because of the narrow range of densities needed to go from q(*)=0 to pi, but also due to sharp 2k(F)(*)-4k(F)(*) transitions; these features render CDW's more versatile for actual applications in heterostructures than in homogeneous systems.

Surface effects on moving vortices in superconducting stripes

Size and surface dynamical effects are investigated in thin superconducting stripes with variable width. We perform numerical simulations of the vortex dynamics, with the inclusion of the surface confining potential and a random distribution of pinning centers. To fully characterize the vortex flow, we calculate the differential resistance, the transverse diffusion coefficient, the structure factor and the intensity of the Bragg peaks, as functions of the transport force. We found that surface effects induce a premature ordering of the flux line lattice, and the system displays plastic and smectic behavior only in a very narrow range of forces. (c) 2007 Elsevier B.V. All rights reserved.

Influence of nearest-neighbor Coulomb interactions on the phase diagram of the ferromagnetic Kondo model

The influence of a nearest-neighbor Coulomb repulsion of strength V on the properties of the ferromagnetic Kondo model is analyzed using computational techniques. The Hamiltonian studied here is defined on a chain using localized S = 1/2 spins, and one orbital per site. Special emphasis is given to the influence of the Coulomb repulsion on the regions of phase separation recently discovered in this family of models, as well as on the double-exchange-induced ferromagnetic ground state. When phase separation dominates at V= 0, the Coulomb interaction breaks the large domains of the two competing phases into small islands of one phase embedded into the other. This is in agreement with several experimental results, as discussed in the text. Vestiges of the original phase separation regime are found in the spin structure factor as incommensurate peaks, even at large values of V. In the ferromagnetic regime close to density n = 0.5, the Coulomb interaction induces tendencies to charge ordering without altering the fully polarized character of the state. This regime of charge-ordered ferromagnetism may be related with experimental observations of a similar phase by Chen and Cheong [Phys. Rev. Lett. 76, 4042 (1996)]. Our results reinforce the recently introduced notion [see, e.g., S. Yunoki et al., Phys. Rev. Lett. 80, 845 (1998)] that in realistic models for manganites analyzed with unbiased many-body techniques, the ground state properties arise from a competition between ferromagnetism and phase-separation - charge-ordering tendencies. ©1999 The American Physical Society.

Multiperiodic magnetic structures in Hubbard superlattices

We consider fermions in one-dimensional superlattices (SL's), modeled by site-dependent Hubbard-U couplings arranged in a repeated pattern of repulsive (i.e., U>0) and free (U=0) sites. Density matrix renormalization group diagonalization of finite systems is used to calculate the local moment and the magnetic structure factor in the ground state. We have found four regimes for magnetic behavior: uniform local moments forming a spin-density wave (SDW), floppy local moments with short-ranged correlations, local moments on repulsive sites forming long-period SDW's superimposed with short-ranged correlations, and local moments on repulsive sites solely with long-period SDW's; the boundaries between these regimes depend on the range of electronic densities ρ and on the SL aspect ratio. Above a critical electronic density, ρ↑↓, the SDW period oscillates both with ρ and with the spacer thickness. The former oscillation allows one to reproduce all SDW wave vectors within a small range of electronic densities, unlike the homogeneous system. The latter oscillation is related to the exchange oscillation observed in magnetic multilayers. A crossover between regimes of thin to thick layers has also been observed.

Critical currents and melting temperature of a two-dimensional vortex lattice with periodic pinning

The critical current and melting temperature of a vortex system are analyzed. Calculations are made for a two-dimensional film at finite temperature with two kinds of periodic pinning: hexagonal and Kagomé. A transport current parallel and perpendicular to the main axis of the pinning arrays is applied and molecular dynamics simulations are used to calculate the vortex velocities to obtain the critical currents. The structure factor and displacements of vortices at zero transport current are used to obtain the melting temperature for both pinning arrays. The critical currents are higher for the hexagonal pinning lattice and anisotropic for both pinning arrays. This anisotropy is stronger with temperature for the hexagonal array. For the Kagomé pinning lattice, our analysis shows a multi stage phase melting; that is, as we increase the temperature, each different dynamic phase melts before reaching the melting temperature. Both the melting temperature and critical currents are larger for the hexagonal lattice, indicating the role for the interstitial vortices in decreasing the pinning strength. © 2012 Springer Science+Business Media New York.

Simulação computacional de materiais com elétrons fortemente interagentes: DMRG aplicado a super-redes Hubbard com modulação de condução entre camadas

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Charge dynamics in half-filled Hubbard chains with finite on-site interaction

We study the charge dynamic structure factor of the one-dimensional Hubbard model with finite on-site repulsion U at half-filling. Numerical results from the time-dependent density matrix renormalization group are analyzed by comparison with the exact spectrum of the model. The evolution of the line shape as a function of U is explained in terms of a relative transfer of spectral weight between the two-holon continuum that dominates in the limit U -> infinity and a subset of the two-holon-two-spinon continuum that reconstructs the electron-hole continuum in the limit U -> 0. Power-law singularities along boundary lines of the spectrum are described by effective impurity models that are explicitly invariant under spin and eta-spin SU(2) rotations. The Mott-Hubbard metal-insulator transition is reflected in a discontinuous change of the exponents of edge singularities at U = 0. The sharp feature observed in the spectrum for momenta near the zone boundary is attributed to a van Hove singularity that persists as a consequence of integrability.

Gaussian deconvolution: a useful method for a form-free modeling of scattering data from mono- and multilayered planar systems

A new method for analysis of scattering data from lamellar bilayer systems is presented. The method employs a form-free description of the cross-section structure of the bilayer and the fit is performed directly to the scattering data, introducing also a structure factor when required. The cross-section structure (electron density profile in the case of X-ray scattering) is described by a set of Gaussian functions and the technique is termed Gaussian deconvolution. The coefficients of the Gaussians are optimized using a constrained least-squares routine that induces smoothness of the electron density profile. The optimization is coupled with the point-of-inflection method for determining the optimal weight of the smoothness. With the new approach, it is possible to optimize simultaneously the form factor, structure factor and several other parameters in the model. The applicability of this method is demonstrated by using it in a study of a multilamellar system composed of lecithin bilayers, where the form factor and structure factor are obtained simultaneously, and the obtained results provided new insight into this very well known system.

Properties of a new poly-ether-glycol copolymer.

Triblock copolymers are made of monomer segments, being the central part usually hydrophobic and the outer parts hydrophilic. By varying sizes, molecular weights and monomer types of the segments one obtains different final molecules, with different physico-chemical properties, which are directly related to the performance of the final product. Looking for new products to be used, among other possibilities, in biological applications, a new polymer (Figure 1) was synthesized by the Dow Chemical and studied by Size Exclusion Chromatography, Fourier Transformed Infrared Spectrometry, Small-angle X-ray Scattering (SAXS) and its cloud point was determined by measuring light transmittance. The studies showed low molecular polydispersivety, but different polarities in the macromolecules fractions. Due to the low solubility of Diol in water, a mixture of water/butyl diglycol was used as solvent. An extensive analysis by SAXS was performed for concentrations from 50 wt% to 80 wt% of Diol in solution. Small concentrations showed very low signal to noise ratio, making it impossible to be analysed. The scattering intensity including the form factor of polydisperse non-homogeneous spheres, and the structure factor of interacting hard spheres was fitted to the curves. As the polymer concentration is high, the fitting of form factors of direct and reverse micelles were compared. The results for direct micelles were better up to 80 wt%, whereas at 90 wt% and 95 wt% the curves were better fitted by reverse micelles. It might seem odd that direct micelles are present up to such high concentrations, but it might have been caused by the presence of butyl diglycol, which increases the solubility of Diol in water. The inner and outer radius of the micelles, electron density distribution, and interaction radius of the micelles were obtained. The polydispersivety increases with Diol concentration. Besides, the interaction radius increases with solvent concentration, even when reversed micelles are present. In the last case, accompanied by an increase of inner radius (water content), as there are fewer Diol molecules to involve the water nuclei, which become larger, further apart, and in less number.

Simulation studies of correlation functions and relaxation in polymeric systems

We investigate the statics and dynamics of a glassy,non-entangled, short bead-spring polymer melt with moleculardynamics simulations. Temperature ranges from slightlyabove the mode-coupling critical temperature to the liquidregime where features of a glassy liquid are absent. Ouraim is to work out the polymer specific effects on therelaxation and particle correlation. We find the intra-chain static structure unaffected bytemperature, it depends only on the distance of monomersalong the backbone. In contrast, the distinct inter-chainstructure shows pronounced site-dependence effects at thelength-scales of the chain and the nearest neighbordistance. There, we also find the strongest temperaturedependence which drives the glass transition. Both the siteaveraged coupling of the monomer and center of mass (CM) andthe CM-CM coupling are weak and presumably not responsiblefor a peak in the coherent relaxation time at the chain'slength scale. Chains rather emerge as soft, easilyinterpenetrating objects. Three particle correlations arewell reproduced by the convolution approximation with theexception of model dependent deviations. In the spatially heterogeneous dynamics of our system weidentify highly mobile monomers which tend to follow eachother in one-dimensional paths forming ``strings''. Thesestrings have an exponential length distribution and aregenerally short compared to the chain length. Thus, arelaxation mechanism in which neighboring mobile monomersmove along the backbone of the chain seems unlikely.However, the correlation of bonded neighbors is enhanced. When liquids are confined between two surfaces in relativesliding motion kinetic friction is observed. We study ageneric model setup by molecular dynamics simulations for awide range of sliding speeds, temperatures, loads, andlubricant coverings for simple and molecular fluids. Instabilities in the particle trajectories are identified asthe origin of kinetic friction. They lead to high particlevelocities of fluid atoms which are gradually dissipatedresulting in a friction force. In commensurate systemsfluid atoms follow continuous trajectories for sub-monolayercoverings and consequently, friction vanishes at low slidingspeeds. For incommensurate systems the velocity probabilitydistribution exhibits approximately exponential tails. Weconnect this velocity distribution to the kinetic frictionforce which reaches a constant value at low sliding speeds. This approach agrees well with the friction obtaineddirectly from simulations and explains Amontons' law on themicroscopic level. Molecular bonds in commensurate systemslead to incommensurate behavior, but do not change thequalitative behavior of incommensurate systems. However,crossed chains form stable load bearing asperities whichstrongly increase friction.

Molekulardynamiksimulationen zur Untersuchung des Mischalkali-Effektes in silikatischen Gläsern

In der vorliegenden Arbeit wird mittels Molekulardynamik(MD)-Computersimulationen die Dynamik von verschiedenen Alkalisilikaten in der Schmelze und im Glas untersucht. Es ist bekannt, daß diese Systeme ionenleitend sind, was auf eine hohe Mobilität der Alkaliionen im Vergleich zu den glasbildenden Komponenten Si und O zurückzuführen ist. Im Mittelpunkt des Interesses steht der sog. Mischalkalieffekt (MAE), der in ternären Mischungen aus Siliziumdioxid mit zwei Alkalioxiden auftritt. Gegenüber Mischungen mit nur einer Alkaliionensorte weisen letztere Systeme eine signifikante Verlangsamung der Alkaliionendiffusion auf. Zunächst werden zwei binäre Alkalisilikate simuliert, nämlich Lithiumdisilikat (LS2) und Kaliumdisilikat (KS2). Die Simulationen zeigen, daß der Ursprung der hohen Mobilität der Alkaliionen in der Struktur begründet ist. KS2 und LS2 weisen auf intermediären Längenskalen Ordnung auf, die in partiellen statischen Strukturfaktoren durch Prepeaks reflektiert ist. Die den Prepeaks zugrundeliegende Struktur erklärt sich durch perkolierende Netzwerke aus alkalioxidreichen Kanälen, die als Diffusionskanäle für die mobilen Alkaliionen fungieren. In diesen Kanälen bewegen sich die Ionen mittels Sprüngen (Hopping) zwischen ausgezeichneten Plätzen. In der Simulation beobachtet man für die hohen Temperaturen (4000K>=1500K) eine ähnliche Aktivierungsenergie wie im Experiment. Im Experiment findet allerdings unterhalb von ca.1200K ein Crossover in ein Arrheniusverhalten mit höherer Aktivierungsenergie statt, welches von der Simulation nicht nachvollzogen wird. Das kann mit der in der Simulation nicht im Gleichgewicht befindlichen Si-O-Matrix erklärt werden, bei der Alterungseffekte beobachtet werden. Am stärksten ist der MAE für eine Alkalikomponente, wenn deren Konzentrationsanteil in einem ternären Mischalkalisystem gegen 0 geht. Daher wird ein LS2-System untersucht, in dem ein Li-Ion gegen ein K-Ion getauscht wird. Der Einfluß des K-Ions ist sowohl lokal in den charakteristischen Abständen zu den ersten nächsten Nachbarn (NN) zu sehen, als auch in der ortsaufgelösten Koordinationszahlverteilung bis zu Längenskalen von ca. 8,5 Angstrom. Die Untersuchung der Dynamik des eingesetzten K-Ions zeigt, daß die Sprungwahrscheinlichkeit nicht mit der Lokalisierung, einem Maß für die Bewegung eines Teilchens um seine Ruheposition, korreliert ist, aber daß eine chemische Umgebung mit wenig Li- und vielen O-NN oder vielen Li- und wenig O-NN ein Sprungereignis begünstigt. Zuletzt wird ein ternäres Alkalisilikat (LKS2) untersucht, dessen Struktur alle charakteristischen Längenskalen von LS2 und KS2 aufweist. Es stellt sich also eine komplexe Struktur mit zwei perkolierenden Subnetzwerken für Alkaliionen ein. Die Untersuchung der Dynamik zeigt eine geringe Wahrscheinlichkeit dafür auf, daß Ionen in ein Subnetzwerk andersnamiger Ionen springen. Auch kann gezeigt werden, daß das Modellpotential den MAE reproduzieren kann, daß also die Diffusionskonstanten in LKS2 bei bis zu einer Größenordnung langsamer sind als in KS2 bzw. LS2. Der beobachtete Effekt stellt sich zudem vom funktionalen Verlauf her so dar, wie er beim MAE erwartet wird. Es wurde auch festgestellt, daß trotz der zeitlichen Verzögerung in den dynamischen Größen die Anzahl der Sprünge pro Zeit nicht geringer ist und daß für niedrige Temperaturen (d.h.im Glas) Sprünge auf den Nachbarplatz mit anschließendem Rücksprung auf die vorherige Position deutlich wahrscheinlicher sind als bei hohen Temperaturen (also in der Schmelze). Die vorliegenden Resultate geben Aufschluß über die Details der Mechanismen mikroskopischer Ionenleitung in binären und ternären Alkalisilikaten sowie dem MAE.

Dynamics of the halo gas in disc galaxies

This Thesis studies the dynamics of hot and cold gas outside the plane in galaxies like the Milky-Way (extra-planar gas) and focuses on the interaction between disc and halo material. Stationary models for the cold phase of the extra-planar gas are presented. They show that the kinematics of this phase must be influenced by the interaction with an ambient medium that we identify as the hot cosmological corona that surrounds disc galaxies. To study this interaction a novel hydrodynamical code has been implemented and a series of hydrodynamical simulations has been run to investigate the mass and momentum exchange between the cold extra-planar gas clouds and the hot corona. These simulations show that the coronal gas can condense efficiently in the turbulent wakes that form behind the cold clouds and it can be accreted by the disc to sustain star formation. They also predict that the corona cannot be a static structure but it must rotate and lag by approximately 80-120 km/s with respect to the disc. Implications of the results of this Thesis for the evolution of star-forming galaxies and for the large-scale dynamics of galactic coronae are also briefly discussed.

Strong Dynamical Heterogeneity and Universal Scaling in Driven Granular Fluids

Large-scale simulations of two-dimensional bidisperse granular fluids allow us to determine spatial correlations of slow particles via the four-point structure factor S-4 (q, t). Both cases, elastic (epsilon = 1) and inelastic (epsilon < 1) collisions, are studied. As the fluid approaches structural arrest, i.e., for packing fractions in the range 0.6 <= phi <= 0.805, scaling is shown to hold: S-4 (q, t)/chi(4)(t) = s(q xi(t)). Both the dynamic susceptibility chi(4)(tau(alpha)) and the dynamic correlation length xi(tau(alpha)) evaluated at the alpha relaxation time tau(alpha) can be fitted to a power law divergence at a critical packing fraction. The measured xi(tau(alpha)) widely exceeds the largest one previously observed for three-dimensional (3d) hard sphere fluids. The number of particles in a slow cluster and the correlation length are related by a robust power law, chi(4)(tau(alpha)) approximate to xi(d-p) (tau(alpha)), with an exponent d - p approximate to 1.6. This scaling is remarkably independent of epsilon, even though the strength of the dynamical heterogeneity at constant volume fraction depends strongly on epsilon.

Supporting Task-oriented Navigation in IDEs with Configurable HeatMaps

Mainstream IDEs generally rely on the static structure of a software project to support browsing and navigation. We propose HeatMaps, a simple but highly configurable technique to enrich the way an IDE displays the static structure of a software system with additional kinds of information. A heatmap highlights software artifacts according to various metric values, such as bright red or pale blue, to indicate their potential degree of interest. We present a prototype system that implements heatmaps, and we describe an initial study that assesses the degree to which different heatmaps effectively guide developers in navigating software.

Querying Runtime Information in the IDE

Code queries focus mainly on the static structure of a system. To comprehend the dynamic behavior of a system however, a software engineer needs to be able to reason about the dynamics of this system, for instance by querying a database of dynamic information. Such a querying mechanism should be directly available in the IDE where the developers implements, navigates and reasons about the software system. We propose (i) concepts to gather dynamic information, (ii) the means to query this information, and (iii) tools and techniques to integrate querying of dynamic information in the IDE, including the presentation of results generated by queries.